| 5 | 1/1 | 返回列表 |
| 查看: 1192 | 回復(fù): 5 | |||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | |||
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/186
+1/81
+1/73
+1/43
+2/40
+1/39
+1/39
+1/18
+1/11
+1/11
+1/10
+1/7
+1/6
+1/6
+1/4
+1/4
+1/3
+1/3
+1/2
+1/1
銅蟲(chóng) (正式寫(xiě)手)
鐵桿木蟲(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ě)手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 266分,求材料相關(guān)專(zhuān)業(yè)調(diào)劑 +6 | 哇呼哼呼哼 2026-03-30 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 各位老師好,我的一志愿為北京科技大學(xué)085601材料專(zhuān)碩 +10 | Koxui 2026-03-28 | 10/500 |
|
|
[考研] 287求調(diào)劑 +14 | land xuxu 2026-03-26 | 14/700 |
|
|
[考研] 08工科,295,接受跨專(zhuān)業(yè)調(diào)劑 +3 | lmnlzy 2026-03-30 | 3/150 |
|
|
[考研] 材料科學(xué)與工程調(diào)劑 +8 | 深V宿舍吧 2026-03-30 | 8/400 |
|
|
[考研] 求調(diào)劑 +7 | 青春裁為三截 2026-03-29 | 7/350 |
|
|
[考研] 一志愿北京工業(yè)大學(xué),324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +8 | 柳溪邊 2026-03-26 | 8/400 |
|
|
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書(shū)。 2026-03-28 | 3/150 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研]
|
18419759900 2026-03-25 | 8/400 |
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書(shū)讀?/a> 2026-03-26 | 6/300 |
|
|
[考研] 08開(kāi)頭275求調(diào)劑 +4 | 拉誰(shuí)不重要 2026-03-26 | 4/200 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 一志愿鄭大085600,310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|