| 6 | 1/1 | 返回列表 |
| 查看: 1190 | 回復(fù): 5 | |||
zfym銅蟲(chóng) (正式寫(xiě)手)
|
[交流]
三個(gè)變量之間的非線(xiàn)性回歸問(wèn)題求助
|
|
三個(gè)變量之間的非線(xiàn)性回歸問(wèn)題求助 已知三個(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é)工具非線(xiàn)性回歸原函數(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 請(qǐng)告訴我matlab代碼 謝謝! |
» 搶金幣啦!回帖就可以得到:
+1/185
+1/89
+2/86
+1/85
+1/78
+1/73
+2/52
+3/48
+1/33
+1/13
+1/11
+1/10
+1/10
+1/7
+1/6
+1/6
+1/5
+1/4
+1/3
+1/1
鐵桿木蟲(chóng) (著名寫(xiě)手)
|
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]) % 對(duì)角線(xiàn) figure(2) plot(z,residual./z*100,'ro') % 實(shí)驗(yàn)值與誤差值繪圖 ylim([-10 10]) xlabel('z') ylabel('相對(duì)誤差 %') hold on refline([0 0]) % 誤差為0的線(xiàn) 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 |
銅蟲(chóng) (正式寫(xiě)手)
銅蟲(chóng) (正式寫(xiě)手)
鐵桿木蟲(chóng) (著名寫(xiě)手)
|
本帖內(nèi)容被屏蔽 |
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +14 | 新時(shí)代材料 2026-03-27 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +11 | 芋泥小鈴鐺 2026-03-28 | 11/550 |
|
|
[考研] 材料工程專(zhuān)碩求調(diào)劑 +6 | hyl3153942 2026-03-29 | 6/300 |
|
|
[考研] 285求調(diào)劑 +5 | AZMK 2026-03-29 | 8/400 |
|
|
[考研] 282求調(diào)劑 +4 | wcq131415 2026-03-24 | 4/200 |
|
|
[考研] 375求調(diào)劑 +6 | 雨夏整夜 2026-03-29 | 6/300 |
|
|
[考研] 329求調(diào)劑,一志愿西北工業(yè)大學(xué),材料工程(085601) +4 | 小小機(jī)靈蟲(chóng) 2026-03-29 | 10/500 |
|
|
[考研] 291求調(diào)劑 +5 | Y-cap 2026-03-29 | 6/300 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +7 | @taotao 2026-03-29 | 7/350 |
|
|
[考研] 329求調(diào)劑 +7 | 星野? 2026-03-26 | 7/350 |
|
|
[考研] 083000學(xué)碩274求調(diào)劑 +8 | Li李魚(yú) 2026-03-26 | 8/400 |
|
|
[考研] 295材料工程專(zhuān)碩求調(diào)劑 +7 | 1428151015 2026-03-27 | 7/350 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動(dòng)285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
[考研] 347求調(diào)劑 +3 | 山頂見(jiàn)α 2026-03-25 | 3/150 |
|
|
[考研] 331環(huán)境科學(xué)與工程求調(diào)劑 +3 | 熠然好運(yùn)氣 2026-03-27 | 3/150 |
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
[考研] 330一志愿中國(guó)海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
[考研] 333求調(diào)劑 +3 | question挽風(fēng) 2026-03-23 | 3/150 |
|
|
[考研] 調(diào)劑推薦 +5 | 清酒714 2026-03-26 | 6/300 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|