 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 848 | 回復: 3 | ||
[求助]
最小二乘 平方和初值小于函數(shù)公差的默認值（動力學方程擬合） 已有2人參與
|
|
本人剛開始接觸matlab在做發(fā)酵動力學，運行結果都可以出來，可是殘差總是出不來，每次都報錯。 系統(tǒng)提示：lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the default value of the function tolerance. (相對平方和的最終變化的初始值小于該函數(shù)公差的默認值。) 我看幫助，好像是說局部最優(yōu)，但不一定是整體最優(yōu)，不太明白，還請知道的大神幫幫我啊 請問是我程序哪里有問題嗎，以下是我的程序 M文件 function f=ObjFunc4LNL8(k,x0,yexp) tspan=[0 2 4 6 8 10 12 14 16 18 20 22 24 26 28]'; [t1,x]=ode45(@kineticseqs8,tspan,x0,[],k); y(:,1)=x(:,1);y(:,2)=x(:,2);y(:,3)=x(:,3); f1=y(:,1)-yexp(:,1);f2=y(:,2)-yexp(:,2);f3=y(:,3)-yexp(:,3); f=[f1 f2 f3]; M文件 function dxdt=kineticseqs8(t,x,k) %模型方程 dxdt=[k(1)*x(1)*(1.0-x(1)/25.022) -k(1)/k(2)*x(1)*(1.0-x(1)/25.022) k(3)*x(1)]; 程序： k0=[0.3654 0.4 0.321]; x0=[1.964 95.668 0.059]; t1=[0 2 4 6 8 10 12 14 16 18 20 22 24 26 28]'; tspan=[0 2 4 6 8 10 12 14 16 18 20 22 24 26 28]'; yexp=[[1.964 3.281 6.250 10.714 15.625 20.758 23.660 24.553 24.843 25.022 24.970 24.703 24.3013 24.078 24.078]; [95.668 90.514 76.770 70.658 63.541 60.228 59.442 55.442 54.215 51.883 49.306 46.116 44.398 44.275 42.925]; [0.059 0.933 16.794 25.060 42.533 59.105 71.308 84.649 104.210 119.331 134.731 150.211 164.022 185.013 196.408]]'; lb=[0.2 0.3 0.3];ub=[0.5 2 0.6]; [k,resnorm,resid,exitflag,output,lambda,jacobian]=lsqnonlin(@ObjFunc4LNL8,k0,lb,ub,[],x0,yexp);%非線性擬合 ci=nlparci(k,resid,jacobian);%計算非線性模型中參數(shù)估計值的95%置信區(qū)間 y1=[yexp(:,1)]'; y2=[yexp(:,2)]'; y3=[yexp(:,3)]'; [t4plot,x4plot]=ode45(@kineticseqs8,[tspan(1) tspan(end)],x0,[],k); plot(t1,y1,'bo',t1,y2,'g*',t1,y3,'r*',t4plot,x4plot,'k-'); fprintf('Estimated Parameter\n'), fprintf('\tk=%.4f±%.4f\n',k(1),ci(1,2)-k(1)) fprintf('\tk=%.4f±%.4f\n',k(2),ci(2,2)-k(2)) fprintf('\tk=%.4f±%.4f\n',k(3),ci(3,2)-k(3)) [t,x]=ode45(@kineticseqs8,tspan,x0,[],k); |
回復此樓» 猜你喜歡
求調(diào)劑
已經(jīng)有3人回復
-
265求調(diào)劑
已經(jīng)有4人回復
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有6人回復
-
一志愿吉林大學材料學碩321求調(diào)劑
已經(jīng)有12人回復
-
286分人工智能專業(yè)請求調(diào)劑愿意跨考！
已經(jīng)有3人回復
-
329求調(diào)劑
已經(jīng)有5人回復
-
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復
已經(jīng)有4人回復
-
材料學碩318求調(diào)劑
已經(jīng)有5人回復
-
一志愿中國海洋大學，生物學，301分，求調(diào)劑
已經(jīng)有6人回復
-
081700化工學碩調(diào)劑
已經(jīng)有3人回復
張曉騰58
新蟲 (正式寫手)
- 應助: 1 (幼兒園)
- 金幣: 5481.5
- 紅花: 5
- 帖子: 560
- 在線: 91.6小時
- 蟲號: 4006790
- 注冊: 2015-08-05
- 性別: GG
- 專業(yè): 微生物生理與生物化學
2樓2016-01-09 19:55:14
dingd
鐵桿木蟲 (職業(yè)作家)
- 應助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
3樓2016-01-09 21:51:10
張曉騰58
新蟲 (正式寫手)
- 應助: 1 (幼兒園)
- 金幣: 5481.5
- 紅花: 5
- 帖子: 560
- 在線: 91.6小時
- 蟲號: 4006790
- 注冊: 2015-08-05
- 性別: GG
- 專業(yè): 微生物生理與生物化學
4樓2016-01-10 16:30:13
相關版塊跳轉
我要訂閱樓主 微冷 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 081700化工學碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 085600材料與化工 求調(diào)劑 +14 | enenenhui 2026-03-13 | 15/750 |
|
|
|
[考研] 復試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
|
[考研] 【考研調(diào)劑】化學專業(yè) 281分，一志愿四川大學，誠心求調(diào)劑 +5 | 吃吃吃才有意義 2026-03-19 | 5/250 |
|
|
|
[考研]
 271材料工程求調(diào)劑
+6
|
.6lL 2026-03-18 | 6/300 |
|
|
|
[考研] 一志愿 西北大學 ，070300化學學碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
|
[考研] 085600材料與化工 +5 | 安全上岸！ 2026-03-16 | 5/250 |
|
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑 +7 | 手機用戶 2026-03-17 | 7/350 |
|
|
|
[考研] 331求調(diào)劑（0703有機化學 +7 | ZY-05 2026-03-13 | 8/400 |
|
|
|
[考研] 070300化學319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
|
[考研] 0703化學336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 藥學383 求調(diào)劑 +3 | 藥學chy 2026-03-15 | 4/200 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|
|
|
[考研] 080500，材料學碩302分求調(diào)劑學校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|