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

返回列表

下一頁(yè)

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

258190169

銅蟲(chóng) (小有名氣)

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

[求助] 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 以及相對(duì)偏差
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論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有123人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？已經(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à)值貼推薦，對(duì)您同樣有幫助:

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

高手進(jìn)來(lái)用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求解方程出問(wèn)題，請(qǐng)幫忙看看已經(jīng)有3人回復(fù)
求回歸方程參數(shù)估計(jì)的最優(yōu)化算法matlab代碼已經(jīng)有9人回復(fù)
[求助]matlab畫(huà)一個(gè)非參數(shù)方程的三維圖形已經(jīng)有3人回復(fù)
matlab中dydt的編寫(xiě)要求已經(jīng)有12人回復(fù)
【求助】用matlab最優(yōu)化方法進(jìn)行參數(shù)擬合已經(jīng)有17人回復(fù)
【求助】求助一個(gè)一階常微分方程的初值問(wèn)題已經(jīng)有16人回復(fù)
【求助】用MATLAB寫(xiě)程序求解簡(jiǎn)單的正弦函數(shù) 已經(jīng)有4人回復(fù)
【求助】使用Matlab預(yù)估動(dòng)力學(xué)方程問(wèn)題已經(jīng)有13人回復(fù)

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

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

回帖置頂 ( 共有1個(gè) )

dbb627

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

專家經(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í)
蟲(chóng)號(hào): 149791
注冊(cè): 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:
本人編寫(xiě)的程序如下但是無(wú)法運(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的回帖

普通回帖

258190169

銅蟲(chóng) (小有名氣)

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

本人編寫(xiě)的程序如下但是無(wú)法運(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的回帖

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

★
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

木蟲(chóng) (職業(yè)作家)

應(yīng)助: 5 (幼兒園)
金幣: 6188.4
散金: 8
紅花: 3
帖子: 3600
在線: 809小時(shí)
蟲(chóng)號(hào): 382693
注冊(cè): 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

銅蟲(chóng) (小有名氣)

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

引用回帖:

樓: 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
但是沒(méi)有計(jì)算出結(jié)果，我的選擇的優(yōu)化算法是第一種，版本是1.5，請(qǐng)問(wèn)為什么

贊一下

回復(fù)此樓

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

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

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

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

贊一下

回復(fù)此樓

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

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

258190169

銅蟲(chóng) (小有名氣)

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

內(nèi)容已刪除

回復(fù)此樓

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

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

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

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

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

贊一下

回復(fù)此樓

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

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

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

上面“起始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 的主題更新

返回列表

下一頁(yè)

普通表情龍兔虎貓

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

[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +5	生物工程調(diào)劑 2026-03-17	9/450	2026-03-21 23:32 by zhujy1982
[考研] 資源與環(huán)境調(diào)劑申請(qǐng)(333分) +5	holy J 2026-03-21	5/250	2026-03-21 22:42 by Catalysis25
[考研] 311求調(diào)劑 +13	冬十三 2026-03-15	14/700	2026-03-21 22:10 by peike
[考研] 0805 316求調(diào)劑 +3	大雪深藏 2026-03-18	3/150	2026-03-21 18:55 by 學(xué)員8dgXkO
[考研] 求調(diào)劑 +5	十三加油 2026-03-21	5/250	2026-03-21 18:48 by 學(xué)員8dgXkO
[考研] 材料工程專碩 348分求調(diào)劑 +3	冬辭. 2026-03-17	5/250	2026-03-21 18:47 by 學(xué)員8dgXkO
[考研] 279分求調(diào)劑一志愿211 +14	chaojifeixia 2026-03-19	15/750	2026-03-21 13:24 by zhukairuo
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4	順順順mr 2026-03-18	5/250	2026-03-21 10:22 by luoyongfeng
[考研] 346求調(diào)劑[0856] +4	WayneLim327 2026-03-16	7/350	2026-03-21 04:02 by JourneyLucky
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3	hisfailed 2026-03-19	6/300	2026-03-20 23:43 by hisfailed
[考研] 考研調(diào)劑求學(xué)校推薦 +3	伯樂(lè)29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] A區(qū)線材料學(xué)調(diào)劑 +5	周周無(wú)極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan
[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	3/150	2026-03-20 11:58 by 呱呱呱呱叫
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 334求調(diào)劑 +3	志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 304求調(diào)劑 +5	素年祭語(yǔ) 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82