 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
- >論壇更新日志 (897)
- >考研 (114)
- >導(dǎo)師招生 (111)
- >蟲(chóng)友互識(shí) (38)
- >考博 (20)
- >基金申請(qǐng) (18)
- >教師之家 (10)
- >公派出國(guó) (9)
- >休閑灌水 (9)
- >論文投稿 (8)
- >找工作 (7)
- >碩博家園 (5)
- >論文道賀祈福 (4)
- >博后之家 (3)
- >文獻(xiàn)求助 (3)
- >藥品生產(chǎn) (2)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2214 | 回復(fù): 10 | |||
| 【懸賞金幣】回答本帖問(wèn)題，作者Evoly_z將贈(zèng)送您 88 個(gè)金幣 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
Evoly_z銀蟲(chóng) (小有名氣)
|
[求助]
求助代謝動(dòng)力學(xué)系數(shù)模擬代碼 1stOpt或者M(jìn)ATLAB 已有3人參與
|
||
 求代謝動(dòng)力學(xué)系數(shù)模擬一級(jí)代謝物的動(dòng)力學(xué)方程：dC2/dt=k1*C1-k2*C2 初始條件：t=0，C2=0 且C1=exp（-A*t）。A=0.2779 我嘗試了用1stOpt(破解版)和MATLAB ODE方法，都沒(méi)成功，想請(qǐng)教一下大神。 另外t不是嚴(yán)格的等差數(shù)列，取值如：t=0，1，2，4，6，10，15，24 1stOpt代碼： Title Kinetic_ave Parameters k1[0,100], k2[0,100]; Variable t, C; StartProgram var i:integer; begin for i:=0 to DataLength -1 do begin if i ==0 C=0; else C:=C[i-1]+k1*(t-t[i-1])*exp(-0.2779 *t) - k2*C*(t-t[i-1]);end; EndProgram; Data; //t C 0 xxx 1 xxx 2 xxx 4 xxx 6 xxx 10 xxx 15 xxx 24 xxx Matlab代碼： function ODE_ave clear all;clc format long aveall; t=T_h(  ;yexp=OLEave( ;k0=[1 1]; y0=0; lb=[0 0]; ub=[+inf +inf]; yy=[y0 yexp']; tspan=0:1:24; [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,y0,yexp); fprintf('\n\n使用函數(shù)lsqnonlin()估計(jì)得到的參數(shù)值為:\n') fprintf('\t待擬合參數(shù) k1 = %.6f\n',k(1)) fprintf('\t待擬合參數(shù) k2 = %.6f\n',k(2)) fprintf(' \t殘差平方和= %.6f\n\n',resnorm) ts=0:1:24; [ts ys]=ode45(@KineticsEqs,ts,y0,[],k); [ttt XXsim] = ode45(@KineticsEqs,tspan,y0,[],k); y=XXsim(2:end); xexp=yexp; R2=1-sum((xexp-y).^2)./sum((xexp-mean(y)).^2); fprintf('\n\t決定系數(shù)R-Square = %.6f',R2); figure(1) plot(ts,ys,'b',tspan,yy,'or'),legend('計(jì)算值','實(shí)驗(yàn)值','Location','best'); yr=y-yexp; figure(2) plot(tspan(2:end),yr,'r*',[-1 15],[0 0]),axis([-1 15 -0.5 0.5]); figure(3) plot(yexp,y,'ro',[21 29],[21 29],'b-'); (作圖這塊兒是copy的，沒(méi)有做修改) %--------------------------------------------------------- function f = ObjFunc(k,tspan,y0,yexp) [t Xsim] = ode45(@KineticsEqs,tspan,y0,[],k) ; ysim = Xsim(2:end); size(ysim); size(yexp); f=ysim(1,1)+ysim(2,1)+ysim(4,1)+ysim(6,1)+ysim(10,1)+ysim(15,1)+ysim(24,1) - sum(yexp(:,1)); %---------------------------------------------------------- function dydt = KineticsEqs(t,y,k) beta(1)=k(1); beta(2)=k(2); dydt = beta(1)*exp(-0.2779*t)-beta(2)*y; 求求啦，被這個(gè)問(wèn)題卡了兩個(gè)多月了，不知道怎么解出k1 k2 |
回復(fù)此樓» 猜你喜歡
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有2人回復(fù)
-
一志愿鄭大材料學(xué)碩298分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有170人回復(fù)
-
尋合作：應(yīng)力腐蝕多尺度模擬
已經(jīng)有6人回復(fù)
-
考研交流
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
幫我的英語(yǔ)口語(yǔ)老師找學(xué)生
已經(jīng)有3人回復(fù)
Evoly_z
銀蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 109.5
- 帖子: 51
- 在線: 30.9小時(shí)
- 蟲(chóng)號(hào): 2116918
- 注冊(cè): 2012-11-09
- 性別: GG
- 專業(yè): 環(huán)境污染化學(xué)
|
我想買一個(gè)5.0以上版本的，結(jié)果入手了一個(gè)5.0的破解版，這個(gè)程序好像也運(yùn)行不了。 哥你能幫我算一次嗎？ Variable x, y; ODEFunction y'= k1*(0.701*exp(-3.211*x)+0.299*exp(-0.067*x))-k2*y; Data; //x, y 0 0 0.33 0.061043 1 0.03675 2 0.05932 4 0.095993 6 0.072057 10 0.05085 15 0.04678 24 0.047673 30 0.034973 48 0.030375 72 0 我按教程寫成這樣了，應(yīng)該沒(méi)問(wèn)題吧？就是算不出結(jié)果。 |
9樓2021-03-30 01:33:28
獨(dú)孤神宇
版主 (知名作家)
- 應(yīng)助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時(shí)
- 蟲(chóng)號(hào): 3522474
- 注冊(cè): 2014-11-06
- 性別: GG
- 專業(yè): 機(jī)械動(dòng)力學(xué)
- 管轄: 計(jì)算模擬
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
參數(shù)擬合C2缺少數(shù)據(jù) 1stOpt容易實(shí)現(xiàn)，1.5不支持微分方程擬合，需要下載5.0版本的 發(fā)自小木蟲(chóng)Android客戶端 |
2樓2021-03-24 07:20:06
Evoly_z
銀蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 109.5
- 帖子: 51
- 在線: 30.9小時(shí)
- 蟲(chóng)號(hào): 2116918
- 注冊(cè): 2012-11-09
- 性別: GG
- 專業(yè): 環(huán)境污染化學(xué)
3樓2021-03-24 11:13:13
dingd
鐵桿木蟲(chóng) (職業(yè)作家)
- 計(jì)算強(qiáng)帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲(chóng)號(hào): 291104
- 注冊(cè): 2006-10-28
4樓2021-03-24 16:12:00
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 334求調(diào)劑 +6 | Trying] 2026-03-05 | 8/400 |
|
|---|---|---|---|---|
|
[考博] 2026年博士名額撿漏 +4 | 科研ya 2026-03-04 | 5/250 |
|
|
|
[考研] 274求調(diào)劑 +6 | cgyzqwn 2026-03-01 | 11/550 |
|
|
|
[考研] 一志愿蘇大材料工程專碩293求調(diào)劑 +3 | 鋼鐵大炮 2026-03-04 | 3/150 |
|
|
|
[考研] 274求調(diào)劑 +9 | 一個(gè)學(xué)習(xí)者 2026-03-04 | 9/450 |
|
|
|
[考研] 求調(diào)劑，學(xué)校研究所都可以，材料與化工267分 +6 | wmx1 2026-03-05 | 6/300 |
|
|
|
[考研] 求材料調(diào)劑 +4 | berdmond 2026-03-05 | 4/200 |
|
|
|
[考研] 接收調(diào)劑 +19 | 津萌津萌 2026-03-02 | 29/1450 |
|
|
|
[考研] 沒(méi)上岸的看過(guò)來(lái) +3 | tangxiaotian 2026-03-01 | 5/250 |
|
|
|
[考研] 0856材料求調(diào)劑 +17 | hyf hyf hyf 2026-02-28 | 18/900 |
|
|
|
[考研] 295求調(diào)劑 +4 | 小賽不吃香菜 2026-03-04 | 4/200 |
|
|
|
[考研] 266材料化工求調(diào)劑 +3 | 哇塞王帥 2026-03-03 | 3/150 |
|
|
|
[考研] 化工專碩348，一志愿985求調(diào)劑 +8 | 弗格個(gè) 2026-02-28 | 11/550 |
|
|
|
[考研]
 計(jì)算機(jī)學(xué)碩 分?jǐn)?shù)285求調(diào)劑
+4
|
glwshine 2026-03-02 | 5/250 |
|
|
|
[考研] 278求調(diào)劑 +3 | 滿天星11_22 2026-03-02 | 3/150 |
|
|
|
[考研] 材料工程求調(diào)劑 +3 | 1431251 2026-03-03 | 3/150 |
|
|
|
[考研] 271求調(diào)劑 +4 | Ricardo1113 2026-03-02 | 4/200 |
|
|
|
[考研] 一志愿華南理工大學(xué)材料與化工326分，求調(diào)劑 +3 | wujinrui1 2026-02-28 | 3/150 |
|
|
|
[考研] 化工專碩342，一志愿大連理工大學(xué)，求調(diào)劑 +6 | kyf化工 2026-02-28 | 7/350 |
|
|
|
[考研] 調(diào)劑 +3 | 13853210211 2026-03-02 | 4/200 |
|