 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 847 | 回復(fù): 3 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
微冷新蟲 (初入文壇)
|
[求助]
最小二乘 平方和初值小于函數(shù)公差的默認(rèn)值（動(dòng)力學(xué)方程擬合） 已有2人參與
|
|
|
本人剛開始接觸matlab在做發(fā)酵動(dòng)力學(xué)，運(yùn)行結(jié)果都可以出來，可是殘差總是出不來，每次都報(bào)錯(cuò)。 系統(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ù)公差的默認(rèn)值。) 我看幫助，好像是說局部最優(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);%計(jì)算非線性模型中參數(shù)估計(jì)值的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); |
回復(fù)此樓» 猜你喜歡
271材料工程求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料專碩274一志愿陜西師范大學(xué)求調(diào)劑
已經(jīng)有8人回復(fù)
-
288求調(diào)劑
已經(jīng)有9人回復(fù)
-
一志愿武漢理工材料工程專碩調(diào)劑
已經(jīng)有4人回復(fù)
-
085601材料工程專碩求調(diào)劑
已經(jīng)有10人回復(fù)
-
321求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑
已經(jīng)有3人回復(fù)
-
307求調(diào)劑
已經(jīng)有6人回復(fù)
-
266求調(diào)劑
已經(jīng)有10人回復(fù)
-
本人考085602 化學(xué)工程 專碩
已經(jīng)有19人回復(fù)
張曉騰58
新蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 5481.5
- 紅花: 5
- 帖子: 560
- 在線: 91.6小時(shí)
- 蟲號(hào): 4006790
- 注冊: 2015-08-05
- 性別: GG
- 專業(yè): 微生物生理與生物化學(xué)
4樓2016-01-10 16:30:13
張曉騰58
新蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 5481.5
- 紅花: 5
- 帖子: 560
- 在線: 91.6小時(shí)
- 蟲號(hào): 4006790
- 注冊: 2015-08-05
- 性別: GG
- 專業(yè): 微生物生理與生物化學(xué)
2樓2016-01-09 19:55:14
dingd
鐵桿木蟲 (職業(yè)作家)
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號(hào): 291104
- 注冊: 2006-10-28
3樓2016-01-09 21:51:10
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 321求調(diào)劑 +8 | 何潤采123 2026-03-18 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑 +3 | 一瓶苯甲酸 2026-03-14 | 3/150 |
|
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +17 | 不知道叫什么！ 2026-03-15 | 19/950 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-18 | 3/150 |
|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +11 | 陌の森林 2026-03-18 | 11/550 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +8 | 妮妮ninicgb 2026-03-15 | 12/600 |
|
|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
|
[考研] 材料專業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 268求調(diào)劑 +7 | 好運(yùn)連綿不絕 2026-03-12 | 8/400 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|
|
|
[考研] 321求調(diào)劑（食品/專碩） +3 | mxcz321 2026-03-12 | 6/300 |
|