| 查看: 308 | 回复: 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显示数组的形式。 谢谢各位高手的点拨，解救我于困境，谢谢！ |
回复此楼» 猜你喜欢
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
-
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有69人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
找到一些相关的精华帖子，希望有用哦~
急求一个过柱子分离的问题，求高手助阵，万分感谢！
已经有27人回复
- 请教涂膜高手帮小弟解决一个问题，困扰我很久了。谢谢！
已经有12人回复
- c语言高手请进，请勿灌水，谢谢！（关于读写文件的问题）
已经有6人回复
- 【求助】关于PMMA的几个小问题，请教各位高手
已经有25人回复
- 【求助】请教高手解决一道行测问题 谢谢
已经有15人回复
- 投AC文章回来，遇棘手问题，请教高手指点！
已经有4人回复
- 【求助】请教高手关于核磁谱中积分的一个问题
已经有10人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
