 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2223 | 回復: 10 | ||
| 【懸賞金幣】回答本帖問題，作者Evoly_z將贈送您 88 個金幣 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
Evoly_z銀蟲 (小有名氣)
|
[求助]
求助代謝動力學系數模擬代碼 1stOpt或者MATLAB 已有3人參與
|
|
 求代謝動力學系數模擬一級代謝物的動力學方程：dC2/dt=k1*C1-k2*C2 初始條件：t=0，C2=0 且C1=exp（-A*t）。A=0.2779 我嘗試了用1stOpt(破解版)和MATLAB ODE方法，都沒成功，想請教一下大神。 另外t不是嚴格的等差數列，取值如：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使用函數lsqnonlin()估計得到的參數值為:\n') fprintf('\t待擬合參數 k1 = %.6f\n',k(1)) fprintf('\t待擬合參數 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決定系數R-Square = %.6f',R2); figure(1) plot(ts,ys,'b',tspan,yy,'or'),legend('計算值','實驗值','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的，沒有做修改) %--------------------------------------------------------- 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; 求求啦，被這個問題卡了兩個多月了，不知道怎么解出k1 k2 |
回復此樓» 猜你喜歡
考研交流
已經有0人回復
-
【新加坡】納米電子器件項目組有“聯合培養(yǎng)博士生”名額
已經有0人回復
-
物理學I論文潤色/翻譯怎么收費?
已經有119人回復
-
【新加坡】納米電子器件項目組有“聯合培養(yǎng)博士生”名額
已經有0人回復
-
幫我的英語口語老師找學生
已經有3人回復
-
計算三階力常數用thirdorder擴胞時參數如何選擇？臨近原子數怎么取值
已經有2人回復
-
2026年循環(huán)經濟功能材料國際會議（ICFMCE 2026）
已經有0人回復
-
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經有0人回復
Evoly_z
銀蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 109.5
- 帖子: 51
- 在線: 30.9小時
- 蟲號: 2116918
- 注冊: 2012-11-09
- 性別: GG
- 專業(yè): 環(huán)境污染化學
3樓2021-03-24 11:13:13
獨孤神宇
版主 (知名作家)
- 應助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時
- 蟲號: 3522474
- 注冊: 2014-11-06
- 性別: GG
- 專業(yè): 機械動力學
- 管轄: 計算模擬
2樓2021-03-24 07:20:06
dingd
鐵桿木蟲 (職業(yè)作家)
- 計算強帖: 4
- 應助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
4樓2021-03-24 16:12:00
Evoly_z
銀蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 109.5
- 帖子: 51
- 在線: 30.9小時
- 蟲號: 2116918
- 注冊: 2012-11-09
- 性別: GG
- 專業(yè): 環(huán)境污染化學
5樓2021-03-24 20:15:22
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] A區(qū)一本交叉課題組，低分調劑，招收機械電子信息通信等交叉方向 +28 | lisimayy 2026-03-04 | 40/2000 |
|
|---|---|---|---|---|
|
|
[考研] 中科大材料299求調劑 +7 | DAIjiayo 2026-03-05 | 10/500 |
|
|
|
[考博] 2026申博自薦 六級440電催化方向 +4 | 櫻落成影花成雙 2026-03-05 | 4/200 |
|
|
[考研] 334求調劑 +6 | Trying] 2026-03-05 | 8/400 |
|
|
|
[考研] 考研282分求調劑，接受跨專業(yè) +4 | 劉淄博 2026-03-04 | 7/350 |
|
|
|
[考研] 085701環(huán)境工程 求調劑 +7 | xiiiia 2026-03-04 | 7/350 |
|
|
|
[考研] 材料328求調劑 +10 | 一個蘿卜02 2026-03-03 | 10/500 |
|
|
|
[考研] 22408-273求調劑-擔任3個項目的負責人-1篇國際期刊論文（一作）1篇核心期刊論文在投。 +3 | 沒想好取什么名 2026-03-03 | 3/150 |
|
|
|
[考研] 307求調劑 +4 | xuyuan111 2026-03-04 | 5/250 |
|
|
|
[考研] 一志愿中科大080500總分324求調劑 +3 | jorna 2026-03-03 | 6/300 |
|
|
|
[考研] 0703化學求調劑 +3 | 將未有期 2026-03-04 | 3/150 |
|
|
|
[考研] 0857求調劑 +4 | CDDD1013 2026-03-04 | 4/200 |
|
|
|
[考研] 本科太原理工采礦工程，求調劑 +3 | onlx 2026-03-01 | 3/150 |
|
|
|
[考研] 306求調劑 +7 | 張張張張oo 2026-03-03 | 7/350 |
|
|
|
[考研] 289求調劑 +7 | BrightLL 2026-03-02 | 9/450 |
|
|
|
[考研] 268求調劑 +10 | 簡單點0 2026-03-02 | 14/700 |
|
|
|
[考研] 291求調劑 +3 | MuoLuo1312 2026-03-02 | 6/300 |
|
|
|
[考研] 298求調劑 +10 | 人間唯你是清歡 2026-02-28 | 14/700 |
|
|
|
[考研] 求調劑 +11 | yunziaaaaa 2026-03-01 | 13/650 |
|
|
|
[考研] 272求調劑 +6 | 田智友 2026-02-28 | 6/300 |
|