| 6 | 1/1 | 返回列表 |
| 查看: 1335 | 回復: 5 | ||
Jounry-Jo銀蟲 (小有名氣)
|
[求助]
急~~~~~~~求利用MATLAB或1STOPT非線性擬合兩組數(shù)據(jù) 已有1人參與
|
| 1stopt5.0是試用版不能擬合三個參數(shù)以上的非線性擬合,matlab試了好久都沒能出正確數(shù)據(jù),不知是否編程有問題。著急得數(shù)據(jù)分析,希望有軟件的大俠幫幫忙。不勝感激 |
鐵桿木蟲 (職業(yè)作家)
|
用1.5老版本的都行吧: 第一組: 均方差(RMSE):1.07205511441145 殘差平方和(SSE):48.2706910701012 相關系數(shù)(R): 0.939291118069207 相關系數(shù)之平方(R^2): 0.8822678044837 確定系數(shù)(DC): 0.882147234885844 卡方系數(shù)(Chi-Square): 1.7441087621366 F統(tǒng)計(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 相關系數(shù)(R): 0.97164462582794 相關系數(shù)之平方(R^2): 0.944093278900318 確定系數(shù)(DC): 0.944088417392726 卡方系數(shù)(Chi-Square): 1.25185088936219 F統(tǒng)計(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 |
銀蟲 (小有名氣)
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 資源與環(huán)境 調劑申請(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿西安交通大學材料工程專業(yè) 282分求調劑 +10 | 楓橋ZL 2026-03-18 | 12/600 |
|
|
[考研] 考研化學學碩調劑,一志愿985 +5 | 張vvvv 2026-03-15 | 7/350 |
|
|
[考研] 材料學碩333求調劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
[考研] 材料 271求調劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[考研] 085601調劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 324分 085600材料化工求調劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
[考研] 藥學383 求調劑 +3 | 藥學chy 2026-03-15 | 5/250 |
|
|
[考研] 中南大學化學學碩337求調劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
[考研] 295材料求調劑,一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
[考研] 319求調劑 +3 | 小力氣珂珂 2026-03-20 | 3/150 |
|
|
[考研] 求調劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[考研] 材料考研調劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 312求調劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 334求調劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
[考研] 308求調劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 321求調劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 070303 總分349求調劑 +3 | LJY9966 2026-03-15 | 5/250 |
|