| 查看: 258 | 回复: 0 | |||
[求助]
请教高手一个小问题，多谢指点哈
|
|
我要用到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显示数组的形式。 谢谢各位高手的点拨，解救我于困境，谢谢！ |
回复此楼» 猜你喜欢
274求调剂
已经有14人回复
-
北京林业大学硕导招生广告
已经有3人回复
-
309求调剂
已经有5人回复
-
292求调剂
已经有8人回复
-
求调剂
已经有4人回复
-
一志愿 西北大学 总分282 英语一62 求调剂
已经有3人回复
-
调剂310
已经有4人回复
-
北科281学硕材料求调剂
已经有16人回复
-
一志愿上海交大生物与医药专硕324分，求调剂
已经有6人回复
-
一志愿哈工大，085400，320，求调剂
已经有4人回复
找到一些相关的精华帖子，希望有用哦~
急求一个过柱子分离的问题，求高手助阵，万分感谢！
已经有27人回复
- 请教涂膜高手帮小弟解决一个问题，困扰我很久了。谢谢！
已经有12人回复
- c语言高手请进，请勿灌水，谢谢！（关于读写文件的问题）
已经有6人回复
- 【求助】关于PMMA的几个小问题，请教各位高手
已经有25人回复
- 【求助】请教高手解决一道行测问题 谢谢
已经有15人回复
- 投AC文章回来，遇棘手问题，请教高手指点！
已经有4人回复
- 【求助】请教高手关于核磁谱中积分的一个问题
已经有10人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
20