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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

咔咔小男人

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 650
散金: 3
紅花: 2
帖子: 88
在線: 53.6小時
蟲號: 2725088
注冊: 2013-10-15
性別: GG
專業(yè): 環(huán)境化工

[求助] matlab求模型參數(shù)有錯誤，求指導(dǎo) 已有1人參與

以下的程序是我照葫蘆畫瓢改的，不知道哪里錯了。。。。。。本人剛開始接觸matlab，求大神指導(dǎo)��！
想用方程dy/dx=c*(A-y)^a*(B-y)進行擬合，求出c和a，其中A=0.9，B=8.5。

CODE:

function KineticsEst1_Int_modified_by_Yuezhilan

clear all;clc

format long

tspan=[0,1,2,3,4,5,7,10,15,20,30,45,90,120]; 

yexp=[0,0.3425,0.61125,0.6125,0.67,0.8025,0.85,0.8375,0.93,0.955,0.95125,0.95375,0.95625,0.96125]';

k0=[1 1];   %%%請注意這里，初值的選取

y0=0.2;

lb=[0 0];

ub=[1000 10];   

yy=[y0 yexp'];

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

    lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,y0,yexp);

fprintf('\n\n使用函數(shù)lsqnonlin()估計得到的參數(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:max(tspan);

[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('計算值','實驗值','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-');

%---------------------------------------------------------

function f = ObjFunc(k,tspan,y0,yexp)           

[t Xsim] = ode45(@KineticsEqs,tspan,y0,[],k)

ysim = Xsim(2:end);

size(ysim);

size(yexp);

f=ysim-yexp;

%----------------------------------------------------------

function dydt = KineticsEqs(t,y,k)

beta(1)=k(1);

beta(2)=k(2);

dydt =beta(1)*(0.9-y)^beta(2)*(7.5-y);

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有218人回復(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ù)
散金幣，求好運，祝面上順利！已經(jīng)有72人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

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

Matlab程序一直有錯誤，崩潰了已經(jīng)有9人回復(fù)
求matlab大神幫忙看下一個計算公式，怎么一直運行錯誤啊已經(jīng)有12人回復(fù)
comsol導(dǎo)出模型到matlab 已經(jīng)有4人回復(fù)
求高手相助！matlab求反應(yīng)動力學(xué)模型參數(shù) 已經(jīng)有6人回復(fù)
使用MATLAB擬合動力學(xué)參數(shù)問題已經(jīng)有18人回復(fù)
matlab求取循環(huán)中最大值對應(yīng)的參數(shù)值已經(jīng)有12人回復(fù)
動力學(xué)方程組參數(shù)求解問題已經(jīng)有9人回復(fù)
logit模型求解求助已經(jīng)有4人回復(fù)
matlab求模型參數(shù)有錯誤，求指導(dǎo) 已經(jīng)有5人回復(fù)
想用遺傳算法求解動力學(xué)參數(shù) 已經(jīng)有13人回復(fù)
matlab 運行ode45出錯，不知道什么原因已經(jīng)有4人回復(fù)
使用matlab最優(yōu)化方法擬合獲得多個動力學(xué)參數(shù)中的問題已經(jīng)有4人回復(fù)
求助利用補償最小二乘法辨識半?yún)?shù)模型的matlab程序已經(jīng)有3人回復(fù)
matlab動力學(xué)模型編程求助已經(jīng)有6人回復(fù)
偏微分方程Matlab求解已經(jīng)有7人回復(fù)
實在是看不出這段matlab程序里哪里出錯了，請大家?guī)兔纯矗?/a> 已經(jīng)有10人回復(fù)

MATLAB關(guān)于lsqcurvefit程序幾個錯誤，求大神指導(dǎo)一下已經(jīng)有6人回復(fù)
matlab求解非線性方程組，錯誤提示怎么解決已經(jīng)有5人回復(fù)
鋰離子電池matlab仿真的問題已經(jīng)有6人回復(fù)
matlab 擬合反應(yīng)動力學(xué)參數(shù)結(jié)果很差。大家?guī)兔匆幌?/a> 已經(jīng)有14人回復(fù)

請教matlab反應(yīng)動力學(xué)參數(shù)估計遇到的問題，謝謝已經(jīng)有15人回復(fù)
求助matlab-simulink運行時報錯問題已經(jīng)有3人回復(fù)

一個人知識的界限就是一個人所能達到的最大范圍的界限。對一種事物的理解、認識深度決定了其利用這種事物上所能實現(xiàn)的價值的有限性。

1樓 2015-06-01 14:44:43

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

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
計算強帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
咔咔小男人: 金幣+20, ★★★★★最佳答案 2015-06-02 18:30:50

CODE:

function KineticsEst1_Int_modified_by_Yuezhilan

clear all;clc

format long

tspan=[0,1,2,3,4,5,7,10,15,20,30,45,90,120]; 

yexp=[0.3425,0.61125,0.6125,0.67,0.8025,0.85,0.8375,0.93,0.955,0.95125,0.95375,0.95625,0.96125]';

k0=[1 1];   %%%請注意這里，初值的選取

y0=0.2;

lb=[-inf -inf];

ub=[+inf +inf];   

yy=[y0 yexp'];

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

    lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,y0,yexp);

fprintf('\n\n使用函數(shù)lsqnonlin()估計得到的參數(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:max(tspan);

[ts ys]=ode45(@KineticsEqs,ts,y0,[],k);

[ttt XXsim] = ode45(@KineticsEqs,tspan,y0,[],k); 

y=XXsim(2:end);

xexp=yexp;

[ttt(2:end) y xexp]

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('計算值','實驗值','Location','best');

%---------------------------------------------------------

function f = ObjFunc(k,tspan,y0,yexp)           

[t Xsim] = ode45(@KineticsEqs,tspan,y0,[],k);

ysim = Xsim(2:end);

size(ysim);

size(yexp);

f=ysim-yexp;

%----------------------------------------------------------

function dydt = KineticsEqs(t,y,k)

beta(1)=k(1);

beta(2)=k(2);

dydt =beta(1)*(0.9-y).^beta(2)*(7.5-y);

初值k0的取值很重要，自己再改改吧。

附圖1.png

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

2樓2015-06-01 18:52:09

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

咔咔小男人

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 650
散金: 3
紅花: 2
帖子: 88
在線: 53.6小時
蟲號: 2725088
注冊: 2013-10-15
性別: GG
專業(yè): 環(huán)境化工

引用回帖:

2樓: Originally posted by 月只藍 at 2015-06-01 18:52:09
function KineticsEst1_Int_modified_by_Yuezhilan
clear all;clc
format long
tspan=;
yexp=';

k0=; %%%請注意這里，初值的選取
y0=0.2;
lb=;
ub=;
yy=;

= ...
lsqnonlin(@ObjFunc,k0,l ...

沒想到作者來了

，感謝您

贊一下

回復(fù)此樓

一個人知識的界限就是一個人所能達到的最大范圍的界限。對一種事物的理解、認識深度決定了其利用這種事物上所能實現(xiàn)的價值的有限性。

3樓2015-06-02 18:31:17

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

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

返回列表

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

[考研] 289材料與化工（085600）B區(qū)求調(diào)劑 +3	這么名字咋樣 2026-03-22	4/200	2026-03-22 17:56 by 云民大李老師
[考研] 317求調(diào)劑 +10	申子申申 2026-03-19	16/800	2026-03-22 17:03 by ColorlessPI
[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +9	niko- 2026-03-19	10/500	2026-03-22 16:08 by ColorlessPI
[考研] 324求調(diào)劑 +6	lucky呀呀呀鴨 2026-03-20	6/300	2026-03-22 16:01 by ColorlessPI
[考研] 319求調(diào)劑 +4	小力氣珂珂 2026-03-20	4/200	2026-03-22 15:53 by ColorlessPI
[考研] 298求調(diào)劑一志愿211 +3	上岸6666@ 2026-03-20	3/150	2026-03-22 15:50 by ColorlessPI
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike
[考研] 269專碩求調(diào)劑 +6	金恩貝 2026-03-21	6/300	2026-03-22 14:31 by ColorlessPI
[考研] 材料求調(diào)劑 +5	@taotao 2026-03-21	5/250	2026-03-21 20:55 by lbsjt
[考研] 一志愿南大，0703化學(xué)，分數(shù)336，求調(diào)劑 +3	收到VS 2026-03-21	3/150	2026-03-21 18:42 by 學(xué)員8dgXkO
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3	邱gl 2026-03-21	3/150	2026-03-21 13:47 by lature00
[考研] 330求調(diào)劑0854 +3	assdll 2026-03-21	3/150	2026-03-21 13:01 by 搏擊518
[考研] 一志愿南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5	@taotao 2026-03-20	5/250	2026-03-20 20:16 by JourneyLucky
[考研] 環(huán)境工程調(diào)劑 +9	大可digkids 2026-03-16	9/450	2026-03-20 17:38 by 醉在風(fēng)里
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4	xbxudjdn 2026-03-18	4/200	2026-03-20 09:07 by 不168
[考研] 求調(diào)劑 +3	暗涌afhb 2026-03-16	3/150	2026-03-20 00:28 by 河南大學(xué)校友
[考研] 材料學(xué)碩318求調(diào)劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

咔咔小男人

[求助] matlab求模型參數(shù)有錯誤，求指導(dǎo) 已有1人參與

» 猜你喜歡

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

月只藍

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

咔咔小男人

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频