| 5 | 1/1 | 返回列表 |
| 查看: 2100 | 回復(fù): 12 | ||||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | ||||
[求助]
用MATLAB擬合如下數(shù)據(jù)(求程序代碼和擬合結(jié)果) 已有1人參與
|
||||
|
現(xiàn)有如下數(shù)據(jù) x=[2.07931 ,1.80644,1.61,1.48184,1.4116]; y=[1.74024,1.63,1.49,1.4,1.36]; 根據(jù)如下關(guān)系式擬合 y=atan((3.14^2*459*1.449e-29*x*A/6.626E-34*3753^2)^0.5)*x/(3.14^2*459*1.449e-29*x*A/6.626E-34*3753^2)^0.5) 辛苦了 |
主管區(qū)長(zhǎng) (職業(yè)作家)
![]() |
專(zhuān)家經(jīng)驗(yàn): +1059 |

主管區(qū)長(zhǎng) (職業(yè)作家)
![]() |
專(zhuān)家經(jīng)驗(yàn): +1059 |
|
感覺(jué)你給出的公式有點(diǎn)問(wèn)題,一下結(jié)果僅供參考: 均方差(RMSE): 0.0590588265039104 殘差平方和(SSE): 0.0174397249400949 相關(guān)系數(shù)(R): 0.993602844977662 相關(guān)系數(shù)之平方(R^2): 0.987246613547703 決定系數(shù)(DC): 0.828021511100056 卡方系數(shù)(Chi-Square): 0.00560692501155114 F統(tǒng)計(jì)(F-Statistic): 232.231639158851 參數(shù) 最佳估算 ---------- ------------- a 1.52138772416843E-16 ====== 結(jié)果輸出 ===== No 實(shí)測(cè)值y 計(jì)算值y 1 1.7402 1.8354974 2 1.63 1.6175395 3 1.49 1.4570489 4 1.4 1.3505577 5 1.36 1.2916440 |

|
function fit_nonl clear all;clc format long data=[ 2.07931 1.74024 1.80644 1.63 1.61 1.49 1.48184 1.4 1.4116 1.36 ]; xspan=data(:,2); %x的數(shù)據(jù),在此輸入 Texp=data(:,1); %T的數(shù)據(jù),在此輸入 k0=[0.1 0.1 1 100]; lb=-[1 1 1 1]*1e9; ub=[1 1 1 1]*1e9; %------------------------------------------------------------------------- % 使用函數(shù)lsqnonlin()進(jìn)行參數(shù)估計(jì) OPTIONS=optimset('MaxFunEvals',1000); [k,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@ObjFunc,k0,lb,ub,OPTIONS,xspan,Texp); ci = nlparci(k,residual,jacobian); %residual; fprintf('\n\n擬合結(jié)果:\n') fprintf('\n\t參數(shù) E1 = %.16f',k(1)) y=KineticsEqs(xspan,k); R2=1-sum((Texp-y).^2)./sum((Texp-mean(y)).^2); fprintf('\n\tR^2 = %.16f',R2); figure plot(xspan,KineticsEqs(xspan,k),'b',xspan,Texp,'or'),legend('計(jì)算值','實(shí)驗(yàn)值','Location','Best') %------------------------------------------------------------------------- function f = ObjFunc(k,xspan,Texp) f=KineticsEqs(xspan,k)-Texp; %------------------------------------------------------------------------ function xt = KineticsEqs(x,k) E1=k(1); xt=atan((3.14^2*459*1.449e-29*x*E1/6.626E-34*3753^2)^0.5)*x/(3.14^2*459*1.449e-29*x*E1/6.626E-34*3753^2)^0.5); 我這個(gè)為啥擬合不出結(jié)果,系統(tǒng)老提示有問(wèn)題,謝謝你的幫忙 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +3 | 閆旭東 2026-03-28 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 340求調(diào)劑 +4 | jhx777 2026-03-27 | 4/200 |
|
|
[考研] 328求調(diào)劑 +6 | 嗯滴的基本都 2026-03-27 | 6/300 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +7 | 邱gl 2026-03-27 | 7/350 |
|
|
[考研] 一志愿北京化工大學(xué)材料與化工 264分各科過(guò)A區(qū)國(guó)家線(xiàn) +10 | 哈哈157349 2026-03-21 | 10/500 |
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
[考研]
|
平樂(lè)樂(lè)樂(lè) 2026-03-26 | 4/200 |
|
|
[考研] 總分293求調(diào)劑 +6 | 加一一九 2026-03-25 | 8/400 |
|
|
[考研] 一志愿哈工大,085400,320,求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
[考研] 332求調(diào)劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑。。! +4 | Acaciad 2026-03-25 | 4/200 |
|
|
[考研] 284求調(diào)劑 +15 | Zhao anqi 2026-03-22 | 15/750 |
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5 | 番茄頭--- 2026-03-22 | 6/300 |
|
|
[考研] 材料/農(nóng)業(yè)專(zhuān)業(yè),07/08開(kāi)頭均可,過(guò)線(xiàn)就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|
|
[考研] 070300,一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|