 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 845 | 回復(fù): 3 | ||
微冷新蟲 (初入文壇)
|
[求助]
最小二乘 平方和初值小于函數(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. (相對(duì)平方和的最終變化的初始值小于該函數(shù)公差的默認(rèn)值。) 我看幫助，好像是說局部最優(yōu)，但不一定是整體最優(yōu)，不太明白，還請(qǐng)知道的大神幫幫我啊 請(qǐng)問是我程序哪里有問題嗎，以下是我的程序 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ù)此樓» 猜你喜歡
294求調(diào)劑材料與化工專碩
已經(jīng)有10人回復(fù)
-
288求調(diào)劑
已經(jīng)有6人回復(fù)
-
287求調(diào)劑
已經(jīng)有4人回復(fù)
-
招收調(diào)劑碩士
已經(jīng)有7人回復(fù)
-
085601材料工程專碩求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿武漢理工材料工程專碩調(diào)劑
已經(jīng)有3人回復(fù)
-
【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑
已經(jīng)有3人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有19人回復(fù)
-
材料考研調(diào)劑
已經(jīng)有3人回復(fù)
-
346求調(diào)劑[0856]
已經(jīng)有6人回復(fù)
張曉騰58
新蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 5481.5
- 紅花: 5
- 帖子: 560
- 在線: 91.6小時(shí)
- 蟲號(hào): 4006790
- 注冊(cè): 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
- 注冊(cè): 2006-10-28
3樓2016-01-09 21:51:10
張曉騰58
新蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 5481.5
- 紅花: 5
- 帖子: 560
- 在線: 91.6小時(shí)
- 蟲號(hào): 4006790
- 注冊(cè): 2015-08-05
- 性別: GG
- 專業(yè): 微生物生理與生物化學(xué)
4樓2016-01-10 16:30:13
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 微冷 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085601材料工程專碩求調(diào)劑 +8 | 慕寒mio 2026-03-16 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 5/250 |
|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
|
[考研] 一志愿天大材料與化工（085600）總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +9 | rare12345 2026-03-18 | 9/450 |
|
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
|
[考研]
 0854，計(jì)算機(jī)類招收調(diào)劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[基金申請(qǐng)] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-12 | 3/150 |
|
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|
|
|
[考研] 321求調(diào)劑（食品/專碩） +3 | xc321 2026-03-12 | 6/300 |
|