 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 866 | 回復: 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); |
回復此樓» 猜你喜歡
求b區(qū)院校調劑
已經有4人回復
-
材料與化工考研調劑
已經有9人回復
-
求調劑
已經有3人回復
-
347求調劑
已經有3人回復
-
總分293求調劑
已經有4人回復
-
284求調劑
已經有15人回復
-
359求調劑
已經有3人回復
-
各位老師您好：本人初試372分
已經有5人回復
-
285求調劑
已經有3人回復
-
282求調劑
已經有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ā)表 | |
|---|---|---|---|---|
|
[考研] 各位老師您好：本人初試372分 +4 | jj涌77 2026-03-25 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 材料277求調劑 +3 | min3 2026-03-24 | 3/150 |
|
|
[考研] 274求調劑 +6 | 顧九笙要謙虛 2026-03-24 | 6/300 |
|
|
|
[考研] 289求調劑 +9 | 懷瑾握瑜l 2026-03-20 | 9/450 |
|
|
|
[考研] 085601求調劑總分293英一數(shù)二 +3 | 鋼鐵大炮 2026-03-24 | 3/150 |
|
|
|
[考研] 0703化學調劑，求導師收 +7 | 天天好運來上岸?/a> 2026-03-24 | 7/350 |
|
|
|
[考研] 0854人工智能方向招收調劑 +3 | 章小魚567 2026-03-24 | 3/150 |
|
|
|
[考研] 070300化學求調劑 +9 | 苑豆豆 2026-03-20 | 9/450 |
|
|
|
[考研] 環(huán)境學碩288求調劑 +8 | 皮皮皮123456 2026-03-22 | 8/400 |
|
|
|
[考研] 335分 | 材料與化工專碩 | GPA 4.07 | 有科研經歷 +4 | cccchenso 2026-03-23 | 4/200 |
|
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調劑 +5 | jiajunX 2026-03-22 | 5/250 |
|
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達問津 2026-03-23 | 5/250 |
|
|
|
[考研] 316求調劑 +7 | 梁茜雯 2026-03-19 | 7/350 |
|
|
|
[考研] 070300，一志愿北航320求調劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
|
[考研] 一志愿華中農業(yè)071010，總分320求調劑 +5 | 困困困困坤坤 2026-03-20 | 6/300 |
|
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
|
[考研] 求調劑 +3 | 13341 2026-03-20 | 3/150 |
|
|
|
[考研] 南京大學化學376求調劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調劑院校 +3 | 葵梓衛(wèi)隊 2026-03-20 | 3/150 |
|