版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

shidoudou

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 94.1
帖子: 21
在線: 11小時(shí)
蟲號(hào): 1410410
注冊(cè): 2011-09-21
專業(yè): 生物化工與食品化工

[求助] 請(qǐng)教高手一個(gè)小問題，多謝指點(diǎn)哈

我要用到matlab處理一部分?jǐn)?shù)據(jù)，只是我不會(huì)從擬合的曲線上取幾個(gè)散點(diǎn)。我的主程序是：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];
我只是想把對(duì)于t1 的各點(diǎn)在擬合的曲線上的數(shù)值顯示出來(lái)，用disp顯示數(shù)組的形式。謝謝各位高手的點(diǎn)撥，解救我于困境，謝謝！

回復(fù)此樓

» 猜你喜歡

290求調(diào)劑已經(jīng)有3人回復(fù)
求調(diào)劑已經(jīng)有3人回復(fù)
327求調(diào)劑已經(jīng)有3人回復(fù)
085600 286分材料求調(diào)劑已經(jīng)有11人回復(fù)
一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑已經(jīng)有15人回復(fù)
085601材料工程找調(diào)劑已經(jīng)有15人回復(fù)
材料與化工272求調(diào)劑已經(jīng)有15人回復(fù)
南京大學(xué)化學(xué)調(diào)劑已經(jīng)有14人回復(fù)
一志愿北京工業(yè)大學(xué)，324分求調(diào)劑已經(jīng)有6人回復(fù)
求調(diào)劑已經(jīng)有5人回復(fù)

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

急求一個(gè)過柱子分離的問題，求高手助陣，萬(wàn)分感謝！已經(jīng)有27人回復(fù)
請(qǐng)教涂膜高手幫小弟解決一個(gè)問題，困擾我很久了。謝謝！已經(jīng)有12人回復(fù)
c語(yǔ)言高手請(qǐng)進(jìn)，請(qǐng)勿灌水，謝謝！（關(guān)于讀寫文件的問題）已經(jīng)有6人回復(fù)
【求助】關(guān)于PMMA的幾個(gè)小問題，請(qǐng)教各位高手已經(jīng)有25人回復(fù)
【求助】請(qǐng)教高手解決一道行測(cè)問題謝謝已經(jīng)有15人回復(fù)
投AC文章回來(lái)，遇棘手問題，請(qǐng)教高手指點(diǎn)！已經(jīng)有4人回復(fù)
【求助】請(qǐng)教高手關(guān)于核磁譜中積分的一個(gè)問題已經(jīng)有10人回復(fù)

點(diǎn)擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 shidoudou 的主題更新