1stopt公式參數(shù)擬合
各位蟲友,我嘗試用1stopt對下面的公式進(jìn)行擬合,想求解出里面的6個參數(shù),但是運行一直出錯,本人是新手,望不吝賜教。謝謝。
Parameter E;A;N;B;C;D;
Variable x,y;
y=(1+E*x)exp[(A+N*E)*x+B*x^2+C*x^3+D*x^4]
data
0 0
5 0.35
10 0.7
20 1.85
50 6
100 12
200 54
返回小木蟲查看更多
今日熱帖
京公網(wǎng)安備 11010802022153號
我知道m(xù)atlab擬合法
可以賜教下嗎
%
x1=temperature(:,1);
x2=temperature(:,4);
x3=temperature(:,5);
y=temperature(:,3); % 目標(biāo) 獲得 f=g(x1,x2,x3)的函數(shù)
X = [ones(size(x1)) x1 x2 x1.*x3 x2.*x3 ]; % 選取的函數(shù)具有形式: y=b0+b1*x1+b2*x2+b3*x1*x3+b4*x2*x3
[b,bint,r,rint,stats] = regress(y,X);
b' % 上述展開系數(shù) b0-b4
stats(1) % 相關(guān)系數(shù) 反映擬合函數(shù)的精度 接近1為優(yōu)
x1fit = min(x1)max(x1)-min(x1))/20:max(x1);
x2fit = min(x2)max(x2)-min(x2))/20:max(x2);
x3fit = min(x3)max(x3)-min(x3))/20:max(x3);
YFIT = b(1)+b(2)*X1FIT+b(3)*(X1FIT.*X3FIT)+ b(4)*(X2FIT.*X3FIT); %目標(biāo)函數(shù)形式,
具體百度regress函數(shù)。 剛才發(fā)的是我的算例函數(shù)
小毛病不少:
模型公式有問題:“(A+N*E)”,A和N兩個參數(shù)完全可以用一個表示即可,否則過擬合(A和N不唯一)。
函數(shù)前要加Function ,要不然它識別不了啊,函數(shù)中需要加乘號的地方也不能吝嗇,比如說exp前面。。。
Parameter E;A;N;B;C;D;
Variable x,y;
Function y=(1+E*x)*exp[(A+N*E)*x+B*x^2+C*x^3+D*x^4];
data;
0 0
5 0.35
10 0.7
20 1.85
50 6
100 12
200 54
還請賜教