版塊導(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í)熱門版塊排行榜

返回列表

【懸賞金幣】回答本帖問題，作者Evoly_z將贈(zèng)送您 88 個(gè)金幣

Evoly_z

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 109.5
帖子: 51
在線: 30.9小時(shí)
蟲號(hào): 2116918
注冊(cè): 2012-11-09
性別: GG
專業(yè): 環(huán)境污染化學(xué)

[求助] 求助代謝動(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方法，都沒成功，想請(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的，沒有做修改)
%---------------------------------------------------------
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è)問題卡了兩個(gè)多月了，不知道怎么解出k1 k2

回復(fù)此樓

» 猜你喜歡

德國亥姆霍茲Hereon中心招收兩位醫(yī)用鎂合金腐蝕與LPSO相變方向2026公派博士生已經(jīng)有0人回復(fù)
推薦一款可以AI輔助寫作的Latex編輯器SmartLatexEditor，超級(jí)好用，推薦試試已經(jīng)有11人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有92人回復(fù)
2026-CJ開始申報(bào)了已經(jīng)有1人回復(fù)
西安電子科學(xué)大學(xué)杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲(chǔ)能器件方向博士已經(jīng)有13人回復(fù)
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息已經(jīng)有1人回復(fù)
一志愿鄭大材料學(xué)碩298分，求調(diào)劑已經(jīng)有6人回復(fù)
尋合作：應(yīng)力腐蝕多尺度模擬已經(jīng)有3人回復(fù)
考研交流已經(jīng)有0人回復(fù)

RELEASE/ME/FROM/THE/PRISON!

1樓 2021-03-24 02:33:08

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

獨(dú)孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時(shí)
蟲號(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ā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

數(shù)值計(jì)算

2樓2021-03-24 07:20:06

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

Evoly_z

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 109.5
帖子: 51
在線: 30.9小時(shí)
蟲號(hào): 2116918
注冊(cè): 2012-11-09
性別: GG
專業(yè): 環(huán)境污染化學(xué)

引用回帖:

2樓: Originally posted by 獨(dú)孤神宇 at 2021-03-24 07:20:06
參數(shù)擬合C2缺少數(shù)據(jù)
1stOpt容易實(shí)現(xiàn)，1.5不支持微分方程擬合，需要下載5.0版本的

謝謝回復(fù)。

有C2的數(shù)據(jù)，因?yàn)镃1服從指數(shù)方程，直接把方程寫在程序里了。C2的數(shù)據(jù)就寫在了下面的數(shù)據(jù)表里。

能分享5.0版的1stOpt下載么？我在網(wǎng)上實(shí)在找不到了，我這里也沒有校內(nèi)bbs之類的東西。。

謝謝了

贊一下

回復(fù)此樓

RELEASE/ME/FROM/THE/PRISON!

3樓2021-03-24 11:13:13

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

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

C2數(shù)據(jù)給出來看看

贊一下

回復(fù)此樓

4樓2021-03-24 16:12:00

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

Evoly_z

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 109.5
帖子: 51
在線: 30.9小時(shí)
蟲號(hào): 2116918
注冊(cè): 2012-11-09
性別: GG
專業(yè): 環(huán)境污染化學(xué)

引用回帖:

4樓: Originally posted by dingd at 2021-03-24 16:12:00
C2數(shù)據(jù)給出來看看

不是數(shù)據(jù)的問題，是程序不能運(yùn)行的問題。點(diǎn)運(yùn)行之后就轉(zhuǎn)圈，輸出那里也沒反應(yīng)。

贊一下

回復(fù)此樓

RELEASE/ME/FROM/THE/PRISON!

5樓2021-03-24 20:15:22

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

獨(dú)孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時(shí)
蟲號(hào): 3522474
注冊(cè): 2014-11-06
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)
管轄: 計(jì)算模擬

【答案】應(yīng)助回帖

將C 的數(shù)據(jù)補(bǔ)上就可以運(yùn)行了。

Parameters k1[0,100], k2[0,100];
Variable t, C;
InitialODEValue t=0,c=0;
ODEFunction c'=k1*exp(-0.2779*t)-k2*c;
Data;
//t C
0 xxx
1 xxx
2 xxx
4 xxx
6 xxx
10  xxx
15  xxx
24  xxx

贊一下

回復(fù)此樓

數(shù)值計(jì)算

6樓2021-03-25 09:23:53

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

Evoly_z

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 109.5
帖子: 51
在線: 30.9小時(shí)
蟲號(hào): 2116918
注冊(cè): 2012-11-09
性別: GG
專業(yè): 環(huán)境污染化學(xué)

引用回帖:

6樓: Originally posted by 獨(dú)孤神宇 at 2021-03-25 09:23:53
將C 的數(shù)據(jù)補(bǔ)上就可以運(yùn)行了。

Parameters k1, k2;
Variable t, C;
InitialODEValue t=0,c=0;
ODEFunction c'=k1*exp(-0.2779*t)-k2*c;
Data;
//t C
0 xxx
1 xxx
2 xxx
4 x ...

好像不對(duì)頭。

大哥加個(gè)qq啊，我好請(qǐng)教你啊

贊一下

回復(fù)此樓

RELEASE/ME/FROM/THE/PRISON!

7樓2021-03-25 13:31:45

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

獨(dú)孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時(shí)
蟲號(hào): 3522474
注冊(cè): 2014-11-06
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)
管轄: 計(jì)算模擬

【答案】應(yīng)助回帖

引用回帖:

7樓: Originally posted by Evoly_z at 2021-03-25 13:31:45
好像不對(duì)頭。

大哥加個(gè)qq啊，我好請(qǐng)教你啊...

1stOpt需要4.0以上版本，你將數(shù)據(jù)附上就行了

發(fā)自小木蟲Android客戶端

贊一下

回復(fù)此樓

數(shù)值計(jì)算

8樓2021-03-25 13:37:38

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

Evoly_z

銀蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 109.5
帖子: 51
在線: 30.9小時(shí)
蟲號(hào): 2116918
注冊(cè): 2012-11-09
性別: GG
專業(yè): 環(huán)境污染化學(xué)

引用回帖:

8樓: Originally posted by 獨(dú)孤神宇 at 2021-03-25 13:37:38
1stOpt需要4.0以上版本，你將數(shù)據(jù)附上就行了
...

我想買一個(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)該沒問題吧？就是算不出結(jié)果。

贊一下

回復(fù)此樓

RELEASE/ME/FROM/THE/PRISON!

9樓2021-03-30 01:33:28

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

獨(dú)孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時(shí)
蟲號(hào): 3522474
注冊(cè): 2014-11-06
性別: GG
專業(yè): 機(jī)械動(dòng)力學(xué)
管轄: 計(jì)算模擬

【答案】應(yīng)助回帖

引用回帖:

9樓: Originally posted by Evoly_z at 2021-03-30 01:33:28
我想買一個(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;
...

這是免費(fèi)的試用版，你這被人忽悠了啊，最新版已經(jīng)到9.0了，而且只能去官網(wǎng)購買。

均方差(RMSE): 0.0171691634594409
殘差平方和(SSR): 0.00324258191286701
相關(guān)系數(shù)(R): 0.802484090026993
相關(guān)系數(shù)之平方(R^2): 0.64398071474645
修正R平方(Adj. R^2): 0.554975893433063
確定系數(shù)(DC): 0.468028249120648
F統(tǒng)計(jì)(F-Statistic): 23.5213130305514

參數(shù) 最佳估算
-------------------- -------------
k1 0.13256115104059
k2 0.358504426638858

====== 結(jié)果輸出 ======

文件: 數(shù)據(jù)文件-1
No x 目標(biāo) y 計(jì)算 y
1 0.33 0.061043 0.0298491702527491
2 1 0.03675 0.053602315170092
3 2 0.05932 0.0683875960100214
4 4 0.095993 0.079360928514414
5 6 0.072057 0.0789300799888936
6 10 0.05085 0.0667093831255373
7 15 0.04678 0.0492934543907035
8 24 0.047673 0.02721397228202
9 30 0.034973 0.0182160615899268
10 48 0.030375 0.00545427669955797
11 72 0 0.00109238646097193

贊一下

回復(fù)此樓

數(shù)值計(jì)算

10樓2021-03-30 08:47:04

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

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

返回列表

不應(yīng)助 確定回帖應(yīng)助 (注意：應(yīng)助才可能被獎(jiǎng)勵(lì)，但不允許灌水，必須填寫15個(gè)字符以上)

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 085600材料與化工（高分子）290分求調(diào)劑 +3	wengyujian 2026-03-04	3/150	2026-03-04 13:08 by zhukairuo
[考研] 一志愿西交大材料學(xué)碩 346 求調(diào)劑 +3	zju51 2026-03-04	3/150	2026-03-04 13:01 by zhukairuo
[考研] 一志愿中南大學(xué)理學(xué)化學(xué) +5	15779376950 2026-03-01	7/350	2026-03-04 11:03 by 一志愿京區(qū)211
[考研] 0703化學(xué)學(xué)碩 288分 +4	旅行中的紫葡萄 2026-03-03	4/200	2026-03-03 23:41 by ch1979
[考研] 中國林科院林化所（南京）2026年招收化學(xué)/材料/環(huán)境工程等背景碩士研究生3名 +3	realstar2006 2026-02-27	3/150	2026-03-03 21:39 by AC.
[考研] 306求調(diào)劑 +7	張張張張oo 2026-03-03	7/350	2026-03-03 21:35 by L135790
[考研] 298求調(diào)劑 +3	人間唯你是清歡 2026-03-03	4/200	2026-03-03 18:09 by ms629
[考研] 江蘇省農(nóng)科院招調(diào)劑1名 +5	Qwertyuop 2026-03-01	5/250	2026-03-03 17:50 by 梅子鴨鴨
[考研] 環(huán)境調(diào)劑 +5	柒槿levana 2026-03-01	5/250	2026-03-03 15:47 by stewie_jz
[考研] 材料學(xué)碩318求調(diào)劑 +15	February_Feb 2026-03-01	17/850	2026-03-03 14:47 by xiaomc_gzh
[考研] 267求調(diào)劑 +6	釣魚佬as 2026-03-02	6/300	2026-03-03 13:59 by 13589
[考研] 清華大學(xué) 材料與化工 353分求調(diào)劑 +5	awaystay 2026-03-02	6/300	2026-03-03 09:03 by houyaoxu
[考研] 298求調(diào)劑 +7	axyz3 2026-02-28	8/400	2026-03-03 06:11 by tgxtgxtgx9
[考研] 高分子化學(xué)與物理調(diào)劑 +6	好好好1233 2026-02-28	16/800	2026-03-03 00:36 by 好好好1233
[考研] 0856材料求調(diào)劑 +12	hyf hyf hyf 2026-02-28	13/650	2026-03-02 20:19 by hypershenger
[考研] 化學(xué)，材料，環(huán)境類求調(diào)劑 +7	考研版棒棒 2026-03-02	7/350	2026-03-02 19:56 by hypershenger
[考研] 一志愿山東大學(xué)材料與化工325求調(diào)劑 +5	半截的詩0927 2026-03-02	5/250	2026-03-02 18:37 by 明亮9527
[考研] 295求調(diào)劑 +8	19171856320 2026-02-28	8/400	2026-03-02 11:19 by yuchj
[考研] 材料調(diào)劑 +6	愛擦汗的可樂冰 2026-02-28	7/350	2026-03-02 10:42 by Jy?
[考研] 317一志愿華南理工電氣工程求調(diào)劑 +6	Soliloquy_Q 2026-02-28	11/550	2026-03-01 11:14 by 歌liekkas