 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學術科研互動平臺 » 計算模擬區(qū) » 程序語言 » 其它 » 請教高手一個小問題，多謝指點哈
| 1 | 1/1 | 返回列表 |
| 查看: 260 | 回復: 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顯示數組的形式。 謝謝各位高手的點撥，解救我于困境，謝謝！ |
回復此樓» 猜你喜歡
085600 材料與化工 329分求調劑
已經有13人回復
-
一志愿南昌大學324求調劑
已經有4人回復
-
085601材料工程找調劑
已經有12人回復
-
求調劑
已經有4人回復
-
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經歷
已經有7人回復
-
求調劑
已經有10人回復
-
化學工程085602 305分求調劑
已經有18人回復
-
一志愿雙一流機械285分求調劑
已經有5人回復
-
一志愿：西北大學，英一數一408-284分求調劑
已經有4人回復
-
291求調劑
已經有5人回復
找到一些相關的精華帖子，希望有用哦~
急求一個過柱子分離的問題，求高手助陣，萬分感謝！
已經有27人回復
- 請教涂膜高手幫小弟解決一個問題，困擾我很久了。謝謝！
已經有12人回復
- c語言高手請進，請勿灌水，謝謝�。P于讀寫文件的問題）
已經有6人回復
- 【求助】關于PMMA的幾個小問題，請教各位高手
已經有25人回復
- 【求助】請教高手解決一道行測問題 謝謝
已經有15人回復
- 投AC文章回來，遇棘手問題，請教高手指點！
已經有4人回復
- 【求助】請教高手關于核磁譜中積分的一個問題
已經有10人回復
點擊這里搜索更多相關資源科研從小木蟲開始，人人為我，我為人人
相關版塊跳轉
我要訂閱樓主 shidoudou 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 300求調劑，材料科學英一數二 +9 | leaflight 2026-03-24 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 298求調劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
[考研] 318一志愿吉林大學生物與醫(yī)藥 求調劑 +4 | 篤行致遠. 2026-03-28 | 4/200 |
|
|
|
[考研] 352分-085602-一志愿985 +5 | 海納百川Ly 2026-03-29 | 5/250 |
|
|
|
[考研] 調劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
|
[考研] 295材料工程專碩求調劑 +7 | 1428151015 2026-03-27 | 7/350 |
|
|
|
[考研] 320分，材料與化工專業(yè)，求調劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
|
[考研] 347求調劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
|
[考研] 求調劑 +6 | 蘆lty 2026-03-25 | 7/350 |
|
|
|
[考研] 286求調劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
|
[考研] 一志愿上海理工能源動力（085800）310分求調劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
|
[考研] 08開頭275求調劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
|
[考研] 314求調劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
|
[考研] 一志愿鄭大085600，310分求調劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
|
[考研]
材料學碩，求調劑
 6+5
|
糖葫蘆888ll 2026-03-22 | 10/500 |
|
|
|
[考研] 281求調劑 +3 | 亞克西good 2026-03-26 | 5/250 |
|
|
|
[考研] 347求調劑 +4 | L when 2026-03-25 | 4/200 |
|
|
|
[考研] 求調劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
|
[考研] 341求調劑(一志愿湖南大學070300) +5 | 番茄頭--- 2026-03-22 | 6/300 |
|
|
|
[考研] 333求調劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|