云淡风轻小刚 发表于 2012-12-21 12:35

微分方程系数拟合

参考了一些文献，拟合三个微分方程组的7个系数，拟合的效果不是很好，麻烦高人帮我看看！
function Dy=ggg(t,y,k)
  Dy=zeros(3,1);
  Dy(1)=k(1)*y(1)-k(2)*y(1).*y(1);
  Dy(2)=k(3)*Dy(1)+k(4)*y(1);
  Dy(3)=-k(5)*Dy(1)-k(6)*Dy(2)-k(7)*y(1);
end

function f=ww(k,y0,yexp)
y0=';
tspan=';
=ode45(@ggg,tspan,y0,[],k);
f1=y(:,1)-yexp(:,1);
f2=y(:,2)-yexp(:,2);
f3=y(:,3)-yexp(:,3);
f=;
end

k0=;
t=;
y0=;
yexp=xlsread('sample.xls');
lb=;ub=[];
=lsqnonlin(@ww,k0,lb,ub,[],y0,yexp);
y1=;
y2=;
y3=;
=ode45(@ggg,,y0,[],k);
plot(t,y1,'bo',t,y2,'g+',t,y3,'r*',t4,y4,'k-')
legend('x','p','s','Model'),
xlabel('time(h)'),ylabel('x(g/l),p(g/l)s(g/l)')


数据
t        X        P        S
0        3.50         0.000         0.0000
12        4.55         1.435         0.7310
24        8.75         1.793         0.8310
36        7.00         2.260         0.8321
48        6.30         2.475         0.8372
60        5.60         2.626         0.8882
72        4.90         2.601         0.8905
84        4.55         2.601         0.8913
96        4.20         2.618         0.9067

squirrelcs 发表于 2012-12-21 18:19

学习学习~~~~~~~~~~~~~~

frankteller 发表于 2012-12-22 10:18

好东西啊！谢谢分享！！

霁雨晴蓝 发表于 2015-2-6 21:30

{:3_59:}{:3_59:}

小舒同学 发表于 2018-5-28 19:17

请问你最后有修改好程序吗？
页: [1]
查看完整版本: 微分方程系数拟合