| 6 | 1/1 | 返回列表 |
| 查看: 1191 | 回復(fù): 5 | |||
[交流]
三個(gè)變量之間的非線性回歸問題求助
|
|
三個(gè)變量之間的非線性回歸問題求助 已知三個(gè)變量之間的函數(shù)關(guān)系為z=(p1*p2*p3*x*y^0.5)/[(1+p2*x)(1+p3*y^0.5)],且知道n組x,y,z的值,能否利用matlab等數(shù)學(xué)工具非線性回歸原函數(shù),求出其中的參數(shù)p1,p2,p3呢? 數(shù)據(jù)如下 x y z 9.9621 89.9837 0.00585 9.0558 90.8959 0.00574 8.3014 91.6560 0.00552 7.6633 92.2989 0.00531 7.1158 92.8498 0.00520 請告訴我matlab代碼 謝謝! |
» 搶金幣啦!回帖就可以得到:
+1/89
+1/81
+1/71
+2/68
+1/57
+1/40
+1/35
+1/34
+2/32
+1/31
+1/30
+3/28
+2/20
+1/17
+1/11
+1/7
+1/6
+1/3
+1/2
+1/1
鐵桿木蟲 (著名寫手)
|
function ccccc clear clc format long data=[9.9621 89.9837 0.00585 9.0558 90.8959 0.00574 8.3014 91.6560 0.00552 7.6633 92.2989 0.00531 7.1158 92.8498 0.00520]; x=data(:,1:2); z=data(:,3); beta0 = [0.0088 0.2012 45.7356]; % 參數(shù)初值 options=optimset('TolFun',1e-8,'display','off'); [beta,resnorm,residual,EXITFLAG] = lsqnonlin(@Objfun,beta0,[],[],options,x,z); zcal=myfun(x,beta); % 模型計(jì)算值 % 結(jié)果輸出 fprintf('\tp1=%10.6f\n\tp2=%10.6f\n\tp3=%10.6f\n',beta) fprintf('\tresidual=\n') fprintf('\t%10.6f\n',residual) fprintf('\tresnorm=%10.6f\n',resnorm) fprintf('\n\tEXITFLAG=%d\n',EXITFLAG) figure(1) plot(z,zcal,'b.') % 實(shí)驗(yàn)值與模型計(jì)算值比較 xlabel('z 實(shí)驗(yàn)值') ylabel('z 計(jì)算值') hold on refline([1 0]) % 對角線 figure(2) plot(z,residual./z*100,'ro') % 實(shí)驗(yàn)值與誤差值繪圖 ylim([-10 10]) xlabel('z') ylabel('相對誤差 %') hold on refline([0 0]) % 誤差為0的線 function f=Objfun(beta,lamda,z) % 目標(biāo)函數(shù),sum((T(i)-Tcal(i))^2)最小化 zcal=myfun(lamda,beta); f=z-zcal; function z=myfun(x,beta) % 模型 p1=beta(1); p2=beta(2); p3=beta(3); z=(p1.*p2.*p3.*x(:,1).*x(:,2).^0.5)./((1+p2.*x(:,1)).*(1+p3*x(:,2).^0.5)); p1= 0.008833 p2= 0.201290 p3= 45.735600 residual= -0.000030 0.000049 0.000007 -0.000037 0.000010 resnorm= 0.000000 EXITFLAG=1 |
鐵桿木蟲 (著名寫手)
|
本帖內(nèi)容被屏蔽 |
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 295材料工程專碩求調(diào)劑 +9 | 1428151015 2026-03-27 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +8 | cccchenso 2026-03-23 | 8/400 |
|
|
[考研] 296求調(diào)劑 +10 | 彼岸t 2026-03-29 | 10/500 |
|
|
[考研] 375求調(diào)劑 +6 | 雨夏整夜 2026-03-29 | 6/300 |
|
|
[考研] 一志愿211,335分,0856,求調(diào)劑院校和導(dǎo)師 +7 | 傾____蕭 2026-03-27 | 8/400 |
|
|
[考研] 327求調(diào)劑 +4 | 小卡不卡. 2026-03-29 | 4/200 |
|
|
[考博] 26申博自薦 +6 | whh869393 2026-03-24 | 6/300 |
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
[考研] 085602 307分 求調(diào)劑 +7 | 不知道叫什么! 2026-03-26 | 7/350 |
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
[考研] 359求調(diào)劑 +4 | 王了個(gè)楠 2026-03-25 | 4/200 |
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 中國科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|