Greatfit.m (Matlab program fot fitting) - Revision history

Oborn at 04:53, 30 October 2008

2008-10-30T04:53:27Z

Oborn at 07:39, 29 October 2008

2008-10-29T07:39:09Z

Oborn at 07:38, 29 October 2008

2008-10-29T07:38:39Z

Oborn at 07:37, 29 October 2008

2008-10-29T07:37:33Z

Oborn at 07:37, 29 October 2008

2008-10-29T07:37:09Z

Oborn at 07:33, 29 October 2008

2008-10-29T07:33:56Z

New page

fid = fopen('datafile.m','r');
z = fscanf(fid, '%g');
fclose(fid);
x=[0:4095];
time=[979:1999];
counts=([z(979:1999)])';

plot(x,z,'r');
xlabel('Channel'); ylabel('Number of Counts');
hold on
modelFun = @(p,x) p(1).*exp((-(x-p(2)).^2)/(2*p(3).^2))+p(4).*exp((-(x-p(5)).^2)/(2*p(6).^2))+p(7).*exp((-(x-p(8)).^2)/(2*p(9).^2));
startingVals = [800 1150 100 500 1250 100 100 1450 300];
coefEsts = nlinfit(time, counts, modelFun, startingVals);
xgrid = linspace(900,2000,1000);
line(xgrid, modelFun(coefEsts, xgrid), 'Color','m');
plot(x, coefEsts(1).*exp((-(x-coefEsts(2)).^2)/(2*coefEsts(3).^2))+coefEsts(4).*exp((-(x-coefEsts(5)).^2)/(2*coefEsts(6).^2))+coefEsts(7).*exp((-(x-coefEsts(8)).^2)/(2*coefEsts(9).^2)),'g','LineWidth',1.5)
axis([700 2500 0 1000])
xsmallarea=0:976;
xlargearea=976:4095;
ysmallarea=coefEsts(1).*exp((-(xsmallarea-coefEsts(2)).^2)/(2*coefEsts(3).^2))+coefEsts(4).*exp((-(xsmallarea-coefEsts(5)).^2)/(2*coefEsts(6).^2))+coefEsts(7).*exp((-(xsmallarea-coefEsts(8)).^2)/(2*coefEsts(9).^2));
ylargearea=coefEsts(1).*exp((-(xlargearea-coefEsts(2)).^2)/(2*coefEsts(3).^2))+coefEsts(4).*exp((-(xlargearea-coefEsts(5)).^2)/(2*coefEsts(6).^2))+coefEsts(7).*exp((-(xlargearea-coefEsts(8)).^2)/(2*coefEsts(9).^2));
SmallArea=trapz(xsmallarea,ysmallarea)
LargeArea=trapz(xlargearea,ylargearea)
Ratio_percent=100*SmallArea/(SmallArea+LargeArea)

← Older revision		Revision as of 04:53, 30 October 2008
Line 44:		Line 44:

	Ratio_percent=100*SmallArea/(SmallArea+LargeArea)		Ratio_percent=100*SmallArea/(SmallArea+LargeArea)
		+
		+	modelFun1 = @(s,x) s(1).exp((-(x-s(2)).^2)/(2s(3).^2));
		+
		+	startingVals = [800 1150 100 ];
		+
		+	coefEsts = nlinfit(timegamma, countsgamma, modelFun1, startingVals);
		+
		+	xgrid = linspace(0,900,1000);
		+
		+	line(xgrid, modelFun1(coefEsts, xgrid), 'Color','m');
		+
		+	plot(x, coefEsts(1).exp((-(x-coefEsts(2)).^2)/(2coefEsts(3).^2)),'b','LineWidth',1.1)
		+

	[http://www.iac.isu.edu/mediawiki/index.php/Peaks_fitting Go Back]		[http://www.iac.isu.edu/mediawiki/index.php/Peaks_fitting Go Back]

← Older revision		Revision as of 07:39, 29 October 2008
Line 45:		Line 45:
	Ratio_percent=100*SmallArea/(SmallArea+LargeArea)		Ratio_percent=100*SmallArea/(SmallArea+LargeArea)

−	[http://www.iac.isu.edu/mediawiki/index.php/~~HRRL~~ Go Back]	+	[http://www.iac.isu.edu/mediawiki/index.php/Peaks_fitting Go Back]

← Older revision		Revision as of 07:38, 29 October 2008
Line 44:		Line 44:

	Ratio_percent=100*SmallArea/(SmallArea+LargeArea)		Ratio_percent=100*SmallArea/(SmallArea+LargeArea)
		+
		+	[http://www.iac.isu.edu/mediawiki/index.php/HRRL Go Back]

← Older revision		Revision as of 07:37, 29 October 2008
Line 18:		Line 18:

	modelFun = @(p,x) p(1).exp((-(x-p(2)).^2)/(2p(3).^2))+p(4).exp((-(x-p(5)).^2)/(2p(6).^2))+p(7).exp((-(x-p(8)).^2)/(2p(9).^2));		modelFun = @(p,x) p(1).exp((-(x-p(2)).^2)/(2p(3).^2))+p(4).exp((-(x-p(5)).^2)/(2p(6).^2))+p(7).exp((-(x-p(8)).^2)/(2p(9).^2));
		+
	startingVals = [800 1150 100 500 1250 100 100 1450 300];		startingVals = [800 1150 100 500 1250 100 100 1450 300];

@@ Line 1: / Line 1: @@
 z = fscanf(fid, '%g');
 fclose(fid);
 x=[0:4095];
 time=[979:1999];
 counts=([z(979:1999)])';
 plot(x,z,'r');
 xlabel('Channel'); ylabel('Number of Counts');
 hold on
 modelFun =  @(p,x) p(1).*exp((-(x-p(2)).^2)/(2*p(3).^2))+p(4).*exp((-(x-p(5)).^2)/(2*p(6).^2))+p(7).*exp((-(x-p(8)).^2)/(2*p(9).^2));
 startingVals = [800 1150 100 500 1250 100 100 1450 300];
 coefEsts = nlinfit(time, counts, modelFun, startingVals);
 xgrid = linspace(900,2000,1000);
 line(xgrid, modelFun(coefEsts, xgrid), 'Color','m');
 plot(x, coefEsts(1).*exp((-(x-coefEsts(2)).^2)/(2*coefEsts(3).^2))+coefEsts(4).*exp((-(x-coefEsts(5)).^2)/(2*coefEsts(6).^2))+coefEsts(7).*exp((-(x-coefEsts(8)).^2)/(2*coefEsts(9).^2)),'g','LineWidth',1.5)
 axis([700 2500 0 1000])
 xsmallarea=0:976;
 xlargearea=976:4095;
 ysmallarea=coefEsts(1).*exp((-(xsmallarea-coefEsts(2)).^2)/(2*coefEsts(3).^2))+coefEsts(4).*exp((-(xsmallarea-coefEsts(5)).^2)/(2*coefEsts(6).^2))+coefEsts(7).*exp((-(xsmallarea-coefEsts(8)).^2)/(2*coefEsts(9).^2));
 ylargearea=coefEsts(1).*exp((-(xlargearea-coefEsts(2)).^2)/(2*coefEsts(3).^2))+coefEsts(4).*exp((-(xlargearea-coefEsts(5)).^2)/(2*coefEsts(6).^2))+coefEsts(7).*exp((-(xlargearea-coefEsts(8)).^2)/(2*coefEsts(9).^2));
 SmallArea=trapz(xsmallarea,ysmallarea)
 LargeArea=trapz(xlargearea,ylargearea)
 Ratio_percent=100*SmallArea/(SmallArea+LargeArea)