24小时热门版块排行榜

返回列表

shidoudou

新虫 (初入文坛)

应助: 1 (幼儿园)
金币: 94.1
帖子: 21
在线: 11小时
虫号: 1410410
注册: 2011-09-21
专业: 生物化工与食品化工

[求助] 请教高手一个小问题，多谢指点哈

我要用到matlab处理一部分数据，只是我不会从拟合的曲线上取几个散点。我的主程序是：format long
global Umax Ks Kp Y m a Ko Uo Ki He Hs Hp D
Umax=1.000;Ks=2.7071; Kp=7.5535; Y=50; m=0.0000; a=0.0329; Ki=0.00005; D=2.5;Ko=0.9315;Uo=0.0532;He=1.480;Hs=20317;Hp=500000;
x0=[0.0514 0.250 0];
t1=[0 6 12 18 24 30 36 42 48 54 60 66]';
yexp=[[0.0514 0.0643 0.0722 0.1081 0.1441 0.2395 0.3688 0.4807 0.5726 0.7192 0.7854 0.8632];[0.25 0.2413 0.2324 0.223 0.216 0.1921 0.1777 0.1574 0.0993 0.0625 0.0447 0.0338];[0 0.0021 0.0034 0.0059 0.0072 0.0106 0.0183 0.0244 0.0375 0.0643 0.0693 0.0721]]';
y1=[yexp(:,1)];y2=[yexp(:,2)];y3=[yexp(:,3)]'
tspan=[0 6 12 18 24 30 36 42 48 54 60 66]';
[t4plot,x4plot]=ode45(@kineticseqs56,tspan,x0);
plot(t1,y1,'b*',t4plot,x4plot(:,1),'k-'),xlabel('T(h)'),ylabel('DOW（g/L）');
figure
plot(t1,y2,'g*',t4plot,x4plot(:,2),'k-'),xlabel('T(h)'),ylabel('Remainder of DBT/(mmol/L)');
figure
plot(t1,y3,'r*',t4plot,x4plot(:,3),'k-'),xlabel('T(h)'),ylabel('Production of 2-HBP/(mmol/L)');

子程序是：
function dxdt=kineticseqs56(t,x)  %模型方程
global Umax Ks Kp Y m a Ko Uo Ki He Hs Hp D
  dxdt1=Umax*x(2)*x(1)*(1-Kp*x(3)/Hp)/(He*Ks+x(2));
  dxdt3=Hp/(1+Hp)*(a*dxdt1+Uo*x(2)*x(1)/(He*Ko+x(2))*(Ki*x(1)/(Ki+x(3)/Hp)));
  dxdt2=Hs/(1+Hs)*(-1/Y*dxdt1-m*x(1)-D*dxdt3*(1+Hp)/Hp);
  dxdt=[dxdt1;dxdt2;dxdt3];
我只是想把对于t1 的各点在拟合的曲线上的数值显示出来，用disp显示数组的形式。谢谢各位高手的点拨，解救我于困境，谢谢！

回复此楼

» 猜你喜欢

1楼 2011-09-21 09:29:00

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主 shidoudou 的主题更新

返回列表