testfit.m
796ebcd0f67bTekRSA-Freq-vs-Time

testfit.m
View Options

			%Testing tool for MyFit
			clear
			close(figure(1))
			x_vec=linspace(0,200,1000);

			% testFit=MyFit('fit_name','Linear','enable_gui',1);
			testFit=MyGorodetsky2000Fit();
			params=cell(1,testFit.n_params);
			switch testFit.fit_name
			case 'Lorentzian'
			for i=1:testFit.n_params
			params{i}=5*rand;
			params{3}=200*rand;
			end
			case 'Exponential'
			params{1}=+5*rand;
			params{2}=-0.1*rand;
			params{3}=0.01*rand;
			case 'Linear'
			params{1}=10*rand;
			params{2}=-50*rand+25;
			case 'DoubleLorentzian'
			for i=1:testFit.n_params
			params{i}=5*rand;
			end
			params{3}=50+10*rand;
			params{2}=params{5};
			params{6}=50+100*rand;
			case 'Gorodetsky2000'
			for i=1:testFit.n_params
			params{i}=5*rand;
			end
			params{3}=0.001*rand;
			params{2}=80+20*rand;
			params{4}=5+10*rand;
			otherwise
			for i=1:testFit.n_params
			params{i}=5*rand;
			end
			end

			params
			y_vec=testFit.anon_fit_fun(x_vec,params{:}).*normrnd(1,0.04,size(x_vec));
			figure(1)
			ax=gca;
			plot(x_vec,y_vec,'x');
			axis([min(x_vec),max(x_vec),0.5min(y_vec),1.5max(y_vec)]);
			hold on
			testFit.plot_handle=ax;
			testFit.enable_plot=1;
			testFit.Data.x=x_vec;
			testFit.Data.y=y_vec;
			% testFit.genInitParams;
			% testFit.init_params
			% testFit.fitTrace;
			% testFit.init_params