| 5 | 1/1 | 返回列表 |
| 查看: 1337 | 回復(fù): 5 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
Jounry-Jo銀蟲 (小有名氣)
|
[求助]
急~~~~~~~求利用MATLAB或1STOPT非線性擬合兩組數(shù)據(jù) 已有1人參與
|
||
| 1stopt5.0是試用版不能擬合三個(gè)參數(shù)以上的非線性擬合,matlab試了好久都沒(méi)能出正確數(shù)據(jù),不知是否編程有問(wèn)題。著急得數(shù)據(jù)分析,希望有軟件的大俠幫幫忙。不勝感激 |
鐵桿木蟲 (職業(yè)作家)
|
用1.5老版本的都行吧: 第一組: 均方差(RMSE):1.07205511441145 殘差平方和(SSE):48.2706910701012 相關(guān)系數(shù)(R): 0.939291118069207 相關(guān)系數(shù)之平方(R^2): 0.8822678044837 確定系數(shù)(DC): 0.882147234885844 卡方系數(shù)(Chi-Square): 1.7441087621366 F統(tǒng)計(jì)(F-Statistic): 61.0931783656144 參數(shù) 最佳估算 ---------- ------------- a21 11.9497208389958 a31 17.4598494080349 a12 -6.19190121843902 a32 0.00659902374601305 a13 4.12126008280832 a23 -0.296670034331323 第二組: 均方差(RMSE):0.973759323354068 殘差平方和(SSE):34.135459913483 相關(guān)系數(shù)(R): 0.97164462582794 相關(guān)系數(shù)之平方(R^2): 0.944093278900318 確定系數(shù)(DC): 0.944088417392726 卡方系數(shù)(Chi-Square): 1.25185088936219 F統(tǒng)計(jì)(F-Statistic): 107.312290588236 參數(shù) 最佳估算 ---------- ------------- a21 7.4558436152277 a31 13.1441376631054 a12 6.26895600341922 a32 3.27183196434597 a13 -3.59060566306702 a23 3.09901050411686 |
銀蟲 (小有名氣)
木蟲 (小有名氣)
|
clc clear % http://www.gaoyang168.com/bbs/viewthread.php?tid=9102854&fpage=4&target=blank data = xlsread('data.xlsx'); x1 = data(:,1); x2 = data(:,2); x3 = data(:,3); y = data(:,4); beta0 = ones(6,1)*1e-5; options = optimset('TolFun' , 1e-15, ... 'TolX' , 1e-15, ... 'MaxFunEvals', 1e5, ... 'MaxIter' , 1e4 ); [X,RESNORM,RESIDUAL,EXITFLAG,OUTPUT,LAMBDA,JACOBIAN]= ... lsqnonlin(@OptObjFunc,beta0,[],[],options,x1,x2,x3,y); y_err = OptObjFunc(X,x1,x2,x3,y); y_cal = y-y_err; plot(y_cal,y,'r*') axis([0 20 0 20]) % =================================== function ans = OptObjFunc(param,x1,x2,x3,y) % [a12,a21,a13,a31,a23,a32] = beta; m = size(x1,1); y_cal = zeros(m,1); for i = 1:m y_cal(i) = nrtl3comp(x1(i),x2(i),x3(i),param); end ans = y-y_cal; % ============================================ function y = nrtl3comp(x1,x2,x3,param) a12 = param(1); a21 = param(2); a13 = param(3); a31 = param(4); a23 = param(5); a32 = param(6); y=(a21*exp(-0.3*a21)*x2+a31*exp(-0.3*a31)*x3)/(x1+exp(-0.3*a21)*x2+exp(-0.3*a31)*x3) ... -(a21*exp(-0.3*a21)*x2+a31*exp(-0.3*a31)*x3)/(x1+exp(-0.3*a21)*x2+exp(-0.3*a31)*x3) ... *x1/(x1+exp(-0.3*a21)*x2+exp(-0.3*a31)*x3)+ ... a12-(a12*exp(-0.3*a12)*x1+a32*exp(-0.3*a32)*x3)/ ... (x2+exp(-0.3*a12)*x1+exp(-0.3*a32)*x3)*x2* ... exp(-0.3*a12)/(x2+exp(-0.3*a12)*x1+exp(-0.3*a32)*x3)+a13 ... -(a13*exp(-0.3*a13)*x1+a23*exp(-0.3*a23)*x2)/ ... (x3+exp(-0.3*a13)*x1+exp(-0.3*a23)*x2)*exp(-0.3*a13)*x3/ ... (x3+exp(-0.3*a13)*x1+exp(-0.3*a23)*x2) ; 5.4478;6.184;-3.5576;13.645;-14.644;2.5835 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 化學(xué)調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
[考研] 311求調(diào)劑 +13 | 冬十三 2026-03-15 | 14/700 |
|
|
[考研] 326求調(diào)劑 +5 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 8/400 |
|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
[考研] 268求調(diào)劑 +9 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 9/450 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂(lè)29 2026-03-18 | 5/250 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +5 | tcxiaoxx 2026-03-20 | 5/250 |
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機(jī)用戶 2026-03-17 | 9/450 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專利,省、國(guó)級(jí)證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|