 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 2046 | 回復(fù): 0 | ||
shidoudou新蟲 (初入文壇)
|
[求助]
matlab 擬合優(yōu)度的判斷語句？
|
|
用matlab軟件中的lsqnonlin 擬合了幾條曲線，想判斷擬合的好壞程度，可是lsqnonlin輸出的是誤差平方和及殘差等，加上什么樣的語句可以顯示類似擬合優(yōu)度R2，及平均相對誤差這類的數(shù)據(jù)呢？謝謝～ 主程序： format short global Umax Ks Kp Y a Uo Ko D Ux Kx Ku Ki; Ki=0.0005;Ku=0.0001; k0=[0.2 0.5 0.1 10 0.01 0.01 0.01 0.1 0.01 0.01]; %參數(shù)的初始值 x0=[0.0514 0.250 0]; %菌體濃度、底物濃度和產(chǎn)物濃度的初始值 t1=[0 6 12 18 24 30 36 42 48 54 60 66]';%發(fā)酵時間的實驗數(shù)據(jù) tspan=[0 6 12 18 24 30 36 42 48 54 60 66 ]';%時間間隔 %yexp實驗數(shù)據(jù)[菌體濃度x1、底物濃度x2和產(chǎn)物濃度x3] yexp=[[0.0514 0.0711 0.0997 0.1887 0.3478 0.6634 0.8711 1.1023 1.2336 1.4239 1.5849 1.6739];[0.25 0.2447 0.2366 0.2258 0.2196 0.203 0.1872 0.1698 0.1453 0.1124 0.0855 0.0608];[0 0.0012 0.0021 0.0039 0.0088 0.0137 0.0183 0.0268 0.0361 0.0438 0.0557 0.0643]]'; lb=[0 0 0 0 0.0001 0.0116 0.001 1 0.001 0.001];ub=[1 10 13 50 50 2 1 1.001 1 1];%參數(shù)的下、上限 [k,resnorm,residual,exitflag,output,lambda,jacobian]=lsqnonlin(@ObjFunc4LNL83,k0,lb,ub,[],x0,yexp);%非線性擬合 y1=[yexp(:,1)];y2=[yexp(:,2)];y3=[yexp(:,3)]'; [t4plot,x4plot]=ode45(@kineticseqs83,[0 100],x0,[],k); plot(t1,y1,'b*',t4plot,x4plot(:,1),'k-'),xlabel('T(h)'),ylabel('Cx/（g/L）'); figure plot(t1,y2,'g*',t4plot,x4plot(:,2),'k-'),xlabel('T(h)'),ylabel('Cs/(mmol/L)'); figure plot(t1,y3,'r*',t4plot,x4plot(:,3),'k-'),xlabel('T(h)'),ylabel('Cp/(mmol/L)'); disp(k) ======================================== function dxdt=kineticseqs83(t,x,k) %模型方程 global Umax Ks Kp Y a Uo Ko D Ux Kx Ku Ki Umax=k(1); Ks=k(2); Kp=k(3); Y=k(4); a=k(5); Uo=k(6); Ko=k(7); D=k(8); Ux=k(9); Kx=k(10); dxdt1=Umax*x(2)*x(1)*(1-Kp*x(3))/(Ks+x(2)); if dxdt1<0 dxdt1=0; end dxdt3=a*dxdt1+Uo*x(2)*x(1)/(Ko+x(2))*(Ki*x(1)/(Ki+x(3))); if dxdt3<0 dxdt3=0; end dxdt2=-1/Y*dxdt1-Ux*x(2)*x(1)/(Kx+x(2))*(Ku*x(1)/(Ku+x(3)))-D*dxdt3; dxdt=[dxdt1;dxdt2;dxdt3]; ========================================== function f=ObjFunc4LNL83(k,x0,yexp) tspan=[0 6 12 18 24 30 36 42 48 54 60 66 ]; [t1,x]=ode45(@kineticseqs83,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*5 f3*10]; |
回復(fù)此樓» 猜你喜歡
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有137人回復(fù)
-
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團隊招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
找到一些相關(guān)的精華帖子，希望有用哦~
遞推最小二乘法在什么書本中有詳細(xì)的講解和matlab仿真介紹？
已經(jīng)有5人回復(fù)
- Matlab數(shù)學(xué)手冊
已經(jīng)有88人回復(fù)
- 關(guān)于Matlab作圖的問題
已經(jīng)有14人回復(fù)
- 聯(lián)苯類手性的判斷
已經(jīng)有22人回復(fù)
- 在網(wǎng)絡(luò)層如何判斷HTTP協(xié)議的請求
已經(jīng)有7人回復(fù)
- matlab 求指點 動力學(xué)方程擬合過程中導(dǎo)數(shù)的獲取
已經(jīng)有10人回復(fù)
- matlab作業(yè)哪位大神麻煩給做一下,十萬火急,周五之間要交
已經(jīng)有3人回復(fù)
- 熱力學(xué)第二定律的疑問——自發(fā)過程的判斷？
已經(jīng)有7人回復(fù)
- MATLAB 數(shù)據(jù)導(dǎo)入
已經(jīng)有13人回復(fù)
- 【求助】matlab 中，幾個m文件調(diào)用
已經(jīng)有14人回復(fù)
- MATLAB安裝后出現(xiàn)一堆java的問題
已經(jīng)有12人回復(fù)
- 【求助】幫我編個程啊~C語言，C++，Matlab，VB...都行~
已經(jīng)有14人回復(fù)
點擊這里搜索更多相關(guān)資源科研從小木蟲開始，人人為我，我為人人
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 shidoudou 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級，六級521，求調(diào)劑 +4 | 美味可樂雞翅 2026-03-26 | 4/200 |
|
|
|
[考研] 299求調(diào)劑 +6 | 嗯嗯嗯嗯2 2026-03-27 | 6/300 |
|
|
|
[考研] 一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師 +4 | 傾____蕭 2026-03-27 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +7 | 丹青奶蓋 2026-03-26 | 8/400 |
|
|
|
[考研] 312求調(diào)劑 +9 | 上岸吧ZJY 2026-03-22 | 13/650 |
|
|
|
[考研] 284求調(diào)劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
|
[考研] 281求調(diào)劑 +3 | 亞克西good 2026-03-26 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | QiMing7 2026-03-25 | 3/150 |
|
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
|
[考研] 0854AI CV方向招收調(diào)劑 +4 | 章小魚567 2026-03-23 | 4/200 |
|
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
|
[考研] 化工專碩求調(diào)劑 +3 | question挽風(fēng) 2026-03-24 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 6/300 |
|
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|