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

登錄注冊

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

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

返回列表

本帖產(chǎn)生 1 個(gè) 計(jì)算強(qiáng)帖，點(diǎn)擊這里進(jìn)行查看

258190169

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 72
紅花: 1
帖子: 61
在線: 40.8小時(shí)
蟲號: 581262
注冊: 2008-07-18
專業(yè): 分離過程

[求助] matlab-常微分方程參數(shù)估計(jì)

初始數(shù)據(jù)濃度和時(shí)間
t=[0,10,30,50,70,90,110,130,150,160];
c=[0,0.23211,0.45906,0.68601,0.92328,1.21213,1.32561,1.34624,1.39782,1.398];
微分方程，dc/dt=[4.41/96485-(4.41*k+L*4.41/96485)c]/(1+4.41*k*t)

要求： 1.得到擬合參數(shù)：k 和L 以及相對偏差
2.得到擬合曲線和數(shù)據(jù)點(diǎn)的圖
3.最好附上院程序

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有231人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有68人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

參數(shù)估計(jì)問題，本人Matlab初學(xué)者，求教大神！已經(jīng)有3人回復(fù)
緊急求助，利用Matlab對實(shí)驗(yàn)數(shù)據(jù)進(jìn)行擬合求解參數(shù)。已經(jīng)有27人回復(fù)
微分方程組求參數(shù)問題，求高人指教，感謝已經(jīng)有11人回復(fù)
關(guān)于matlab的參數(shù)估計(jì) 已經(jīng)有15人回復(fù)
微分方程組參數(shù)擬合的問題（多參數(shù)動力系統(tǒng)參數(shù)估計(jì)）已經(jīng)有11人回復(fù)
matlab處理常微分方程作圖問題已經(jīng)有7人回復(fù)
多變量齊次方程解已經(jīng)有4人回復(fù)
關(guān)于matlab線性規(guī)劃背包問題的解法求救呀~~方程知道但不知道如何解答···· 已經(jīng)有4人回復(fù)
matlab 擬合反應(yīng)動力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a> 已經(jīng)有14人回復(fù)

高手進(jìn)來用matlab求解常微分方程，十分感謝！已經(jīng)有17人回復(fù)
求助：幫忙做個(gè)matlab優(yōu)化計(jì)算已經(jīng)有11人回復(fù)
MATLAB用于三元回歸參數(shù) 已經(jīng)有3人回復(fù)
matlab擬合方程參數(shù)時(shí)初值的選擇已經(jīng)有15人回復(fù)
用matlab求解方程出問題，請幫忙看看已經(jīng)有3人回復(fù)
求回歸方程參數(shù)估計(jì)的最優(yōu)化算法matlab代碼已經(jīng)有9人回復(fù)
[求助]matlab畫一個(gè)非參數(shù)方程的三維圖形已經(jīng)有3人回復(fù)
matlab中dydt的編寫要求已經(jīng)有12人回復(fù)
【求助】用matlab最優(yōu)化方法進(jìn)行參數(shù)擬合已經(jīng)有17人回復(fù)
【求助】求助一個(gè)一階常微分方程的初值問題已經(jīng)有16人回復(fù)
【求助】用MATLAB寫程序求解簡單的正弦函數(shù) 已經(jīng)有4人回復(fù)
【求助】使用Matlab預(yù)估動力學(xué)方程問題已經(jīng)有13人回復(fù)

1樓 2011-12-16 15:34:33

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

258190169

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 72
紅花: 1
帖子: 61
在線: 40.8小時(shí)
蟲號: 581262
注冊: 2008-07-18
專業(yè): 分離過程

本人編寫的程序如下但是無法運(yùn)行，由于是新手，希望高手幫忙調(diào)試一下：
function PenicilliumEst
clear all;
t=[0,10,30,50,70,90,110,130,150,160];
y=[0,0.23211,0.45906,0.68601,0.92328,1.21213,1.32561,1.34624,1.39782,1.398];
y0=0;

% Nonlinear least square estimate using lsqnonlin()
beta0=[0.005 0.001];
lb=[0 0];ub=[inf inf];
[beta,resnorm,residual,exitflag,output,lambda,jacobian] = ...
lsqnonlin(@Func,beta0,lb,ub,[],t,y);
ci = nlparci(beta,residual,jacobian);

% =======================================
function f = Func(beta,t,y,y0)       % Define objective function
tspan = [0  max(x)];
[tt yy] = ode45(@ModelEqs,tspan,y0,[],beta);
yc= spline(tt,yy,x);
f1=y-yc
% ==================================
function dydt = ModelEqs(t,y,beta)       % Model equations
dydt = [4.41/96485-(4.41*beta(1)+beta(2)*4.41/96485)*y]/(1+4.41*beta(1)*t)

% result
fprintf('\n Estimated Parameters by Lsqnonlin():\n')
fprintf('\t k1 = %.4f ± %.4f\n',beta(1),ci(1,2)-beta(1))
fprintf('\t k2 = %.4f ± %.4f\n',beta(2),ci(2,2)-beta(2))
fprintf('  The sum of the residual squares is: %.1e\n\n',sum(residual.^2))

% plot of fit results
tspan = [0  max(t)];
[tt yc] = ode45(@modeleqs,tspan,c0,[],beta);
tc=linspace(0,max(t),200);
yc = spline(tt,yc,tc);
plot(t,c,'ro',tc,yca,'r-');
hold on
xlabel('Time');
ylabel('Concentration');
hold off

贊一下

回復(fù)此樓

2樓2011-12-16 16:49:37

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

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗(yàn): +4
計(jì)算強(qiáng)帖: 12
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時(shí)
蟲號: 149791
注冊: 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

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

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
258190169(金幣+20): 多謝大俠你的QQ是多少可以和你聯(lián)系一下嘛 2011-12-16 22:12:03
cenwanglai(金幣+5, 計(jì)算強(qiáng)帖+1): 謝謝給予幫助~ 2011-12-20 09:07:46

引用回帖:

2樓: Originally posted by 258190169 at 2011-12-16 16:49:37:
本人編寫的程序如下但是無法運(yùn)行，由于是新手，希望高手幫忙調(diào)試一下：
function PenicilliumEst
clear all;
t=[0,10,30,50,70,90,110,130,150,160];
y=[0,0.23211,0.45906,0.68601,0.92328,1.21213,1.32561, ...

CODE:

function PenicilliumEst

clear all;

t=[0,10,30,50,70,90,110,130,150,160];

y=[0,0.23211,0.45906,0.68601,0.92328,1.21213,1.32561,1.34624,1.39782,1.398];

y0=0;

% Nonlinear least square estimate using lsqnonlin()

beta0=[0.001 0.001];

lb=[0 0];ub=[inf inf];

[beta,resnorm,residual,exitflag,output,lambda,jacobian] = ...

    lsqnonlin(@Func,beta0,lb,ub,[],t,y,y0);         

ci = nlparci(beta,residual,jacobian);

beta

% result

fprintf('\n Estimated Parameters by Lsqnonlin():\n')

fprintf('\t k1 = %.4f ± %.4f\n',beta(1),ci(1,2)-beta(1))

fprintf('\t k2 = %.4f ± %.4f\n',beta(2),ci(2,2)-beta(2))

fprintf('  The sum of the residual squares is: %.1e\n\n',sum(residual.^2))

% plot of fit results

tspan = [0  max(t)];

[tt yc] = ode45(@ModelEqs,tspan,y0,[],beta);

tc=linspace(0,max(t),200);

yca = spline(tt,yc,tc);

plot(t,y,'ro',tc,yca,'r-');

hold on

xlabel('Time');

ylabel('Concentration');

hold off 

% =======================================

function f1 = Func(beta,t,y,y0)        % Define objective function

tspan =t;

[tt yy] = ode45(@ModelEqs,tspan,y0,[],beta);

yc= spline(tt,yy,t);

f1=y-yc;

% ==================================

function dydt = ModelEqs(t,y,beta)          % Model equations

dydt = (4.41/96485-(4.41*beta(1)+beta(2)*4.41/96485)*y)/(1+4.41*beta(1)*t);

改的可以運(yùn)行了，但是初值不合適

贊一下(1人)

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

3樓2011-12-16 17:31:51

已閱回復(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í)
蟲號: 291104
注冊: 2006-10-28

★
dbb627(金幣+1): 感謝參與 2011-12-17 11:40:34

用1stOpt試試：

CODE:

Variable t,c;

ODEFunction c'=(4.41/96485-(4.41*k+L*4.41/96485)*c)/(1+4.41*k*t);

Data;

0        0

10        0.23211

30        0.45906

50        0.68601

70        0.92328

90        1.21213

110        1.32561

130        1.34624

150        1.39782

160        1.398

均方差(RMSE): 0.272257391318038
殘差平方和(SSE): 0.66711678414573
相關(guān)系數(shù)(R): 0.932733204157737
相關(guān)系數(shù)之平方(R^2): 0.869991230138359
決定系數(shù)(DC): 0.576237758605965

參數(shù) 最佳估算
-------------------- -------------
k -2.32798002359073
l 514537.340812352

贊一下(1人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

258190169

4樓2011-12-17 09:52:26

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

sui2066

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

應(yīng)助: 5 (幼兒園)
金幣: 6188.4
散金: 8
紅花: 3
帖子: 3600
在線: 809小時(shí)
蟲號: 382693
注冊: 2007-05-25
性別: GG
專業(yè): 聚合物共混與復(fù)合材料

謝謝老大！

回復(fù)此樓

氟硅(富貴)http://www.dowpont.com/bbs/index.php^_^

5樓2011-12-17 19:37:49

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

258190169

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 72
紅花: 1
帖子: 61
在線: 40.8小時(shí)
蟲號: 581262
注冊: 2008-07-18
專業(yè): 分離過程

引用回帖:

樓: Originally posted by dingd at 2011-12-17 09:52:26:
用1stOpt試試：
[code]
Variable t,c;
ODEFunction c'=(4.41/96485-(4.41*k+L*4.41/96485)*c)/(1+4.41*k*t);
Data;
0 0
10 0.23211
30 0.45906
50 0.68601
70 0.92328
90 1.21213
110 1.32561
130 1. ...

Title "model";
Parameters k, L;
Variable t,c;
ODEFunction c'=(4.41/96485-(4.41*k+L*4.41/96485)*c)/(1+4.41*k*t);
Data;
0 0.02063
600 0.23211
1800 0.45906
3000 0.68601
4200 0.92328
5400 1.21213
6600 1.32561
7800 1.34624
9000 1.39782
9600 1.398
但是沒有計(jì)算出結(jié)果，我的選擇的優(yōu)化算法是第一種，版本是1.5，請問為什么

贊一下

回復(fù)此樓

6樓2011-12-19 21:45:05

已閱回復(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í)
蟲號: 291104
注冊: 2006-10-28

1.5版好像還沒有微分方程計(jì)算功能，換個(gè)新的吧，5.0都快出來了。

贊一下

回復(fù)此樓

7樓2011-12-19 22:12:17

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

258190169

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 72
紅花: 1
帖子: 61
在線: 40.8小時(shí)
蟲號: 581262
注冊: 2008-07-18
專業(yè): 分離過程

內(nèi)容已刪除

回復(fù)此樓

8樓2011-12-20 09:27:31

已閱回復(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í)
蟲號: 291104
注冊: 2006-10-28

現(xiàn)在有點(diǎn)明白國產(chǎn)軟件發(fā)展艱難的原因了

，起始Origin的創(chuàng)始人好像也是中國人，但在美國發(fā)展，大家有用正版Origin的嗎？

贊一下

回復(fù)此樓

9樓2011-12-20 14:10:09

已閱回復(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í)
蟲號: 291104
注冊: 2006-10-28

上面“起始Origin”應(yīng)該為“其實(shí)Origin”

贊一下

回復(fù)此樓

10樓2011-12-20 14:18:57

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

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

返回列表

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

[考研] 材料學(xué)碩301分求調(diào)劑 +7	Liyouyumairs 2026-03-21	7/350	2026-03-21 22:31 by peike
[考研] 初試 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 化學(xué)調(diào)劑 +5	yzysaa 2026-03-21	5/250	2026-03-21 22:12 by peike
[考研] 材料求調(diào)劑 +5	@taotao 2026-03-21	5/250	2026-03-21 20:55 by lbsjt
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3	收到VS 2026-03-21	3/150	2026-03-21 18:42 by 學(xué)員8dgXkO
[考研] 0703化學(xué)297求調(diào)劑 +3	Daisy☆ 2026-03-20	3/150	2026-03-21 17:45 by ColorlessPI
[考研] 299求調(diào)劑 +4	某某某某位 2026-03-21	4/200	2026-03-21 16:30 by barlinike
[考研] 310求調(diào)劑 +3	baibai1314 2026-03-16	3/150	2026-03-21 03:56 by JourneyLucky
[考研] 265求調(diào)劑 +9	梁梁校校 2026-03-17	9/450	2026-03-21 02:17 by JourneyLucky
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5	守候夕陽CF 2026-03-18	5/250	2026-03-21 01:06 by JourneyLucky
[考研] 330求調(diào)劑 +4	小材化本科 2026-03-18	4/200	2026-03-20 23:13 by JourneyLucky
[考研] 324求調(diào)劑 +5	lucky呀呀呀鴨 2026-03-20	5/250	2026-03-20 22:30 by 促天成
[考研] A區(qū)線材料學(xué)調(diào)劑 +5	周周無極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3	我想要讀研究生 2026-03-20	3/150	2026-03-20 20:13 by JourneyLucky
[考研] 材料與化工專碩調(diào)劑 +7	heming3743 2026-03-16	7/350	2026-03-20 19:31 by zhukairuo
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10	Liwangman 2026-03-15	10/500	2026-03-19 10:25 by 無際的草原
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy

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

[求助] matlab-常微分方程參數(shù)估計(jì)

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

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

» 本帖已獲得的紅花（最新10朵）

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助: