我有一组数据,认为是正态分布的。
我用 fit函数,希望得到其中的参数。
参数为 y0,a,w,xc。
正态函数我定义为fun_futing_normal.m
数据以及fit 的程序是如下定义的。
clear all
clc
clf
x=[-402.50;-397.50;-392.50;-387.50;-382.50;-377.50;-372.50;-367.50;-362.50;-357.50;-352.50;-347.50;-342.50;-337.50;-332.50;-327.50;-322.50;-317.50;-312.50;-307.50;-302.50;-297.50;-292.50;-287.50;-282.50;-277.50;-272.50;-267.50;-262.50;-257.50;-252.50;-247.50;-242.50;-237.50;-232.50;-227.50;-222.50;-217.50;-212.50;-207.50;-202.50;-197.50;-192.50;-187.50;-182.50;-177.50;-172.50;-167.50;-162.50;];
y=[0;0;0;0;1;2;2;9;5;14;25;32;49;69;98;127;112;162;160;153;149;131;112;73;80;63;65;51;49;47;36;39;29;25;13;10;4;4;0;0;0;0;0;0;0;0;0;0;0;];
% define fit
ft = fittype( 'fun_futing_normal(x,y0,a,w,xc)' )
f = fit( x, y, ft, 'StartPoint', [-0.4, 10000, 70, -300] );
% function generate points