| 5 | 1/1 | 返回列表 |
| 查看: 2101 | 回復(fù): 12 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
[求助]
用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ū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗: +1059 |
|
感覺你給出的公式有點問題,一下結(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)計(F-Statistic): 232.231639158851 參數(shù) 最佳估算 ---------- ------------- a 1.52138772416843E-16 ====== 結(jié)果輸出 ===== No 實測值y 計算值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 |

主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗: +1059 |

|
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ù)估計 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('計算值','實驗值','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); 我這個為啥擬合不出結(jié)果,系統(tǒng)老提示有問題,謝謝你的幫忙 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +10 | Micky11223 2026-03-25 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 一志愿上海理工能源動力(085800)310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
[考研] 289求調(diào)劑 +7 | 新時代材料 2026-03-27 | 7/350 |
|
|
[考研] 292求調(diào)劑 +13 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 14/700 |
|
|
[考研] 085600,材料與化工321分,求調(diào)劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實習(xí)經(jīng)歷 +10 | 材料學(xué)257求調(diào)劑 2026-03-23 | 11/550 |
|
|
[考研] 材料調(diào)劑 +8 | 匹克i 2026-03-23 | 8/400 |
|
|
[考研] 284求調(diào)劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
[考研] 329求調(diào)劑 +7 | 鈕恩雪 2026-03-25 | 7/350 |
|
|
[考研] 271求調(diào)劑 +6 | 生如夏花… 2026-03-22 | 6/300 |
|
|
[考研] 297求調(diào)劑 +6 | 田洪有 2026-03-26 | 6/300 |
|
|
[考研] 求調(diào)劑 +3 | QiMing7 2026-03-25 | 3/150 |
|
|
[考研] 求b區(qū)院校調(diào)劑 +4 | 周56 2026-03-24 | 5/250 |
|
|
[考研] 0854AI CV方向招收調(diào)劑 +4 | 章小魚567 2026-03-23 | 4/200 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +10 | 紙魚ly 2026-03-21 | 11/550 |
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|
|
[考研] 材料/農(nóng)業(yè)專業(yè),07/08開頭均可,過線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|