在利用origin等商业软件对某些复杂的公式进行拟合时,常常会出现未知原因的拟合失败。这并不是说明我们的公式本身有问题,很有可能是因为商业软件本身的算法有问题。这时候就可以应用matlab中的曲线拟合功能进行拟合。下面附上案例代码,为方便后文叙述,不同部分之间用空行隔开。
syms t;
x=[50; 200; 400; 600; 800];
y=[0.00000209; 0.00000267; 0.00000825; 0.000748; 0.0216];
f=fittype('C*(1-(0.0005)^(0.00013-0.000000372*t+0.0000000036*t^2))^(n+1)','independent','t','coefficients',{'C','n'});
[cfun,rsquare]=fit(x,y,f,'Lower',[580,1.4],'Upper',[3000000,3],'StartPoint',[600, 1.5]);
xi=0:1:850;
yi=cfun(xi);
plot(x,y,'r*',xi,yi,'b-');
第一部分进行的操作是定义变量及赋已知点
需要注意的点:
1、自变量t要单独定义
2、已知点必须保存为列向量。
3、列向量要一 一对应(数据量大的时候一定要检查一下)。
第二部分进行的操作是定义拟合函数f
这里调用了fittype函数来进行函数的自定义。其语法可简述为:
f=fittype('公式具体表达','ind