 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2102 | 回復(fù): 12 | |||||
[求助]
用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) 辛苦了  |
回復(fù)此樓» 猜你喜歡
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有292人回復(fù)
-
中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團(tuán)隊(duì)招收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)價(jià)值貼推薦，對(duì)您同樣有幫助:
- matlab 多元函數(shù)擬合
已經(jīng)有9人回復(fù)
- 使用matlab最優(yōu)化方法擬合獲得多個(gè)動(dòng)力學(xué)參數(shù)中的問題
已經(jīng)有4人回復(fù)
- 求助：用matlab或其它軟件進(jìn)行如下計(jì)算模擬，多謝
已經(jīng)有3人回復(fù)
- 拜求！matlab擬合動(dòng)力學(xué)參數(shù)
已經(jīng)有14人回復(fù)
- ]緊急求助復(fù)介電常數(shù)的模型擬合-Matlab or 1stOpt
已經(jīng)有9人回復(fù)
- 求大神用matlab給我處理一組數(shù)據(jù)，有人沒有��？
已經(jīng)有15人回復(fù)
- 求Origin或Matlab高手非線性擬合公式
已經(jīng)有5人回復(fù)
- matlab 擬合求助
已經(jīng)有8人回復(fù)
- 請(qǐng)幫忙擬合一組數(shù)據(jù)，求回歸系數(shù)
已經(jīng)有3人回復(fù)
- 三維圖譜中，用matlab如何實(shí)現(xiàn)去除一部分?jǐn)?shù)據(jù)并實(shí)現(xiàn)曲面擬合
已經(jīng)有10人回復(fù)
- 用matlab里的nlinfit如何擬合分段函數(shù)
已經(jīng)有8人回復(fù)
- matlab對(duì)數(shù)據(jù)進(jìn)行最小二乘擬合，若果數(shù)據(jù)為負(fù)數(shù)怎么擬合？急急急！
已經(jīng)有3人回復(fù)
- matlab數(shù)據(jù)擬合問題求助
已經(jīng)有7人回復(fù)
- 用MATLAB或者ORIGIN怎么擬合復(fù)數(shù)數(shù)據(jù)
已經(jīng)有10人回復(fù)
- 坐等----MAtlab中如何有三維數(shù)據(jù)擬合曲面，并得到曲面方程？？？？？？？
已經(jīng)有10人回復(fù)
- matlab擬合求值
已經(jīng)有17人回復(fù)
- 需要用matlab擬合數(shù)據(jù)的外行人急需幫助
已經(jīng)有10人回復(fù)
- 【求助】matlab曲面擬合表達(dá)式
已經(jīng)有10人回復(fù)
月只藍(lán)
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68128.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
感覺你給出的公式有點(diǎ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í)測值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 |
2樓2014-10-28 17:00:36
3樓2014-10-28 17:02:51
月只藍(lán)
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68128.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
qj2006: 金幣+20 2014-10-28 21:23:55
qj2006: 金幣+20 2014-10-28 21:23:55
4樓2014-10-28 17:07:52
|
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)老提示有問題，謝謝你的幫忙 |
5樓2014-10-28 17:27:27
6樓2014-10-28 18:24:58
月只藍(lán)
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68128.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
7樓2014-10-28 19:15:53
8樓2014-10-28 19:53:04
月只藍(lán)
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗(yàn): +1059 - 計(jì)算強(qiáng)帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68128.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時(shí)
- 蟲號(hào): 1122189
- 注冊(cè): 2010-10-14
- 專業(yè): 宇宙學(xué)
- 管轄: 計(jì)算模擬區(qū)
9樓2014-10-29 08:28:13
|
你好幫我把剩下幾組擬合一下，我實(shí)在不會(huì)，謝謝你 第一組 x=[2.07931 ,1.80644,1.61,1.48184,1.4116]; y1=[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)； 第二組 x=[2.07931 ,1.80644,1.61,1.48184,1.4116]; y2=[1.73395,1.58,1.45,1.37,1.34]; 根據(jù)如下關(guān)系式擬合 y=atan((3.14^2*460*1.444e-29*x*A/6.626E-34*3756^2)^0.5)*x/(3.14^2*460*1.444e-29*x*A/6.626E-34*3756^2)^0.5)； 第三組 x=[2.07931 ,1.80644,1.61,1.48184,1.4116]; y3=[1.55275,1.51,1.41,1.34,1.3]; 根據(jù)如下關(guān)系式擬合 y=atan((3.14^2*460*1.431-29*x*A/6.626E-34*3750^2)^0.5)*x/(3.14^2*460*1.431e-29*x*A/6.626E-34*3750^2)^0.5)； 第四組 x=[2.07931 ,1.80644,1.61,1.48184,1.4116]; y4=[1.55959,1.52,1.42,1.35,1.33]; 根據(jù)如下關(guān)系式擬合 y=atan((3.14^2*461*1.418-29*x*A/6.626E-34*3742^2)^0.5)*x/(3.14^2*461*1.418e-29*x*A/6.626E-34*3742^2)^0.5)； |
10樓2014-10-29 09:30:25
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 qj2006 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過 2026-03-27 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 求調(diào)劑 +6 | 白QF 2026-03-21 | 7/350 |
|
|
|
[有機(jī)交流]
高溫高壓反應(yīng)求助
 10+4
|
chibby 2026-03-25 | 4/200 |
|
|
[考研] 283求調(diào)劑（080500） +4 | A child 2026-03-27 | 4/200 |
|
|
|
[考研] 08開頭275求調(diào)劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
|
[考研] 一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271 +6 | kotoko_ik 2026-03-23 | 7/350 |
|
|
|
[考研] 071000生物學(xué)求調(diào)劑，初試成績343 +6 | 小小甜面團(tuán) 2026-03-25 | 6/300 |
|
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
|
[考研] 調(diào)劑310 +3 | 溫柔的晚安 2026-03-25 | 4/200 |
|
|
|
[考研] 生物技術(shù)與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
|
[考研] 0854電子信息求調(diào)劑 +7 | α____ 2026-03-22 | 9/450 |
|
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
|
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境） +4 | 我愛學(xué)電池 2026-03-23 | 4/200 |
|
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
|
[考研]
 材料考研調(diào)劑生
+3
|
黃粱一夢千年 2026-03-24 | 3/150 |
|
|
|
[考研]
 一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
+7
|
2117205181 2026-03-21 | 8/400 |
|
|
|
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
|
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|