| 5 | 1/1 | 返回列表 |
| 查看: 2106 | 回復(fù): 12 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
[求助]
用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) 辛苦了 |
|
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)題,謝謝你的幫忙 |
主管區(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 |

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

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856求調(diào)劑 +5 | 楒桉 2026-03-28 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +7 | 爭(zhēng)取九點(diǎn)睡 2026-03-28 | 8/400 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動(dòng)285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
[考研] 311(085601)求調(diào)劑 +4 | liziyeyeye 2026-03-28 | 4/200 |
|
|
[考研] 339求調(diào)劑,想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研] 311求調(diào)劑 +3 | 希望上岸阿小楊 2026-03-23 | 3/150 |
|
|
[考研] 291求調(diào)劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
[考研] 085600材料與化工調(diào)劑 +10 | A-哆啦Z夢(mèng) 2026-03-23 | 16/800 |
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[考研] 0703化學(xué)338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
[考研] 一志愿武理085500機(jī)械專(zhuān)業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
[考研] 300分,材料,求調(diào)劑,英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
[考研]
|
黃粱一夢(mèng)千年 2026-03-24 | 3/150 |
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程(085601) +5 | WW.' 2026-03-23 | 7/350 |
|
|
[基金申請(qǐng)] 請(qǐng)教下大家 2026年國(guó)家基金申請(qǐng)是雙盲審嗎? +3 | lishucheng1 2026-03-22 | 5/250 |
|
|
[考研] 一志愿國(guó)科過(guò)程所081700,274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|