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

Oborn at 04:53, 30 October 2008

2008-10-30T04:53:27Z

← 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]

Oborn at 07:39, 29 October 2008

2008-10-29T07:39:09Z

← 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]

Oborn at 07:38, 29 October 2008

2008-10-29T07:38:39Z

← 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]

Oborn at 07:37, 29 October 2008

2008-10-29T07:37:33Z

← 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];

Oborn at 07:37, 29 October 2008

2008-10-29T07:37:09Z

← Older revision		Revision as of 07:37, 29 October 2008
Line 1:		Line 1:
			fid = fopen('datafile.m','r');

	~~fid = fopen('datafile.m','r');~~
	z = fscanf(fid, '%g');		z = fscanf(fid, '%g');

	fclose(fid);		fclose(fid);

	x=[0:4095];		x=[0:4095];

	time=[979:1999];		time=[979:1999];

	counts=([z(979:1999)])';		counts=([z(979:1999)])';

	plot(x,z,'r');		plot(x,z,'r');

	xlabel('Channel'); ylabel('Number of Counts');		xlabel('Channel'); ylabel('Number of Counts');

	hold on		hold on

	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];

	coefEsts = nlinfit(time, counts, modelFun, startingVals);		coefEsts = nlinfit(time, counts, modelFun, startingVals);

	xgrid = linspace(900,2000,1000);		xgrid = linspace(900,2000,1000);

	line(xgrid, modelFun(coefEsts, xgrid), 'Color','m');		line(xgrid, modelFun(coefEsts, xgrid), 'Color','m');

	plot(x, coefEsts(1).exp((-(x-coefEsts(2)).^2)/(2coefEsts(3).^2))+coefEsts(4).exp((-(x-coefEsts(5)).^2)/(2coefEsts(6).^2))+coefEsts(7).exp((-(x-coefEsts(8)).^2)/(2coefEsts(9).^2)),'g','LineWidth',1.5)		plot(x, coefEsts(1).exp((-(x-coefEsts(2)).^2)/(2coefEsts(3).^2))+coefEsts(4).exp((-(x-coefEsts(5)).^2)/(2coefEsts(6).^2))+coefEsts(7).exp((-(x-coefEsts(8)).^2)/(2coefEsts(9).^2)),'g','LineWidth',1.5)

	axis([700 2500 0 1000])		axis([700 2500 0 1000])

	xsmallarea=0:976;		xsmallarea=0:976;

	xlargearea=976:4095;		xlargearea=976:4095;

	ysmallarea=coefEsts(1).exp((-(xsmallarea-coefEsts(2)).^2)/(2coefEsts(3).^2))+coefEsts(4).exp((-(xsmallarea-coefEsts(5)).^2)/(2coefEsts(6).^2))+coefEsts(7).exp((-(xsmallarea-coefEsts(8)).^2)/(2coefEsts(9).^2));		ysmallarea=coefEsts(1).exp((-(xsmallarea-coefEsts(2)).^2)/(2coefEsts(3).^2))+coefEsts(4).exp((-(xsmallarea-coefEsts(5)).^2)/(2coefEsts(6).^2))+coefEsts(7).exp((-(xsmallarea-coefEsts(8)).^2)/(2coefEsts(9).^2));

	ylargearea=coefEsts(1).exp((-(xlargearea-coefEsts(2)).^2)/(2coefEsts(3).^2))+coefEsts(4).exp((-(xlargearea-coefEsts(5)).^2)/(2coefEsts(6).^2))+coefEsts(7).exp((-(xlargearea-coefEsts(8)).^2)/(2coefEsts(9).^2));		ylargearea=coefEsts(1).exp((-(xlargearea-coefEsts(2)).^2)/(2coefEsts(3).^2))+coefEsts(4).exp((-(xlargearea-coefEsts(5)).^2)/(2coefEsts(6).^2))+coefEsts(7).exp((-(xlargearea-coefEsts(8)).^2)/(2coefEsts(9).^2));

	SmallArea=trapz(xsmallarea,ysmallarea)		SmallArea=trapz(xsmallarea,ysmallarea)

	LargeArea=trapz(xlargearea,ylargearea)		LargeArea=trapz(xlargearea,ylargearea)

	Ratio_percent=100*SmallArea/(SmallArea+LargeArea)		Ratio_percent=100*SmallArea/(SmallArea+LargeArea)

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)