版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

zfym

銅蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 884.4
帖子: 556
在線: 283.1小時(shí)
蟲(chóng)號(hào): 854647

[交流] 三個(gè)變量之間的非線性回歸問(wèn)題求助

三個(gè)變量之間的非線性回歸問(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é)工具非線性回歸原函數(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代碼
謝謝！

回復(fù)此樓

» 猜你喜歡

【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有102人回復(fù)
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！已經(jīng)有1人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)
計(jì)算與實(shí)驗(yàn)科學(xué)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院已經(jīng)有10人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【公式擬合問(wèn)題】如何擬合下面這樣的一個(gè)公式，多謝大牛。已經(jīng)有6人回復(fù)
非線性回歸參數(shù)..急求幫助��！謝謝已經(jīng)有8人回復(fù)
【求助】注冊(cè)化工工程師的報(bào)考條件，報(bào)考時(shí)間以及考試科目都是什么？已經(jīng)有26人回復(fù)
【求助】關(guān)于SPSS多元非線性回歸的問(wèn)題已經(jīng)有3人回復(fù)
【求助】求MAtlab求解程序�。。� 已經(jīng)有20人回復(fù)
【求助】大家談?wù)剬?duì)注冊(cè)環(huán)保工程師的了解已經(jīng)有15人回復(fù)
【求助】什么樣的數(shù)據(jù)適合主成分分析已經(jīng)有8人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-03-16 20:57:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

change0618

鐵桿木蟲(chóng) (著名寫(xiě)手)

計(jì)算強(qiáng)帖: 1
應(yīng)助: 44 (小學(xué)生)
金幣: 17724.5
帖子: 2413
在線: 546.7小時(shí)
蟲(chóng)號(hào): 496517

你的數(shù)組太少，而且你總要給出個(gè)參數(shù)的范圍吧，有物理意義吧

贊一下

回復(fù)此樓

5樓2011-03-17 11:35:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 6 個(gè)回答

change0618

鐵桿木蟲(chóng) (著名寫(xiě)手)

計(jì)算強(qiáng)帖: 1
應(yīng)助: 44 (小學(xué)生)
金幣: 17724.5
帖子: 2413
在線: 546.7小時(shí)
蟲(chóng)號(hào): 496517

★
zfym(金幣+1):謝謝參與
zfym(金幣+9): 謝謝了，太好了，麻煩你了！ 2011-03-17 10:13:44

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ì)角線
figure(2)
plot(z,residual./z*100,'ro')       % 實(shí)驗(yàn)值與誤差值繪圖
ylim([-10 10])
xlabel('z')
ylabel('相對(duì)誤差 %')
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

贊一下(1人)

回復(fù)此樓

2樓2011-03-17 08:50:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zfym

銅蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 884.4
帖子: 556
在線: 283.1小時(shí)
蟲(chóng)號(hào): 854647

引用回帖:

Originally posted by change0618 at 2011-03-17 08:50:49:
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= ...

??? function ccccc
|
Error: Function definitions are not permitted in this context.

該怎么辦

贊一下

回復(fù)此樓

3樓2011-03-17 10:25:06

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zfym

銅蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 884.4
帖子: 556
在線: 283.1小時(shí)
蟲(chóng)號(hào): 854647

引用回帖:

你算的結(jié)果和1stOpt差別較大，是什么原因呢，并且結(jié)果和你選的初始值比較接近。

贊一下

回復(fù)此樓

4樓2011-03-17 10:39:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 6 個(gè)回答

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 08工科求調(diào)劑286 +4	tgs_001 2026-03-28	4/200	2026-03-30 21:03 by dophin1985
[考研] 317求調(diào)劑 +8	十閑wx 2026-03-24	8/400	2026-03-30 21:00 by peike
[考研] 085600 295分求調(diào)劑 +7	W55j 2026-03-30	9/450	2026-03-30 20:36 by dick_runner
[考研] 0703化學(xué)321分求調(diào)劑 +10	三dd. 2026-03-30	11/550	2026-03-30 19:24 by markhwc
[考研] 289求調(diào)劑 +16	新時(shí)代材料 2026-03-27	16/800	2026-03-30 19:04 by Wang200018
[考研] 一志愿南開(kāi)大學(xué)0710生物學(xué)359求調(diào)劑 +5	兔兔兔111223314 2026-03-29	7/350	2026-03-30 18:29 by 兔兔兔111223314
[考研] 293求調(diào)劑 +3	末未mm 2026-03-30	5/250	2026-03-30 17:23 by 王保杰33
[考研] 284求調(diào)劑 +14	junqihahaha 2026-03-26	15/750	2026-03-30 14:12 by 探123
[考研] 329求調(diào)劑 +8	星野? 2026-03-26	8/400	2026-03-30 13:41 by chemdavid
[考研] 南京大學(xué)化學(xué)調(diào)劑 +10	景隨風(fēng) 2026-03-29	15/750	2026-03-30 11:21 by limeifeng
[考研] 348求調(diào)劑 +6	小懶蟲(chóng)不懶了 2026-03-28	6/300	2026-03-30 10:29 by Evan_Liu
[考研] 085600，材料與化工321分求調(diào)劑 +10	大饞小子 2026-03-28	10/500	2026-03-29 23:35 by 飛行日記西
[考研] 343求調(diào)劑 +6	愛(ài)羈絆 2026-03-29	6/300	2026-03-29 12:00 by 無(wú)際的草原
[考研] 298求調(diào)劑 +4	種圣賜 2026-03-28	4/200	2026-03-29 08:42 by q1092522407
[考研] 316求調(diào)劑 +7	江辭666 2026-03-26	7/350	2026-03-28 21:28 by sanrepian
[考研] 394求調(diào)劑 +3	好事多磨靜候佳?/a> 2026-03-26	5/250	2026-03-28 14:24 by 唐沐兒
[考研] 266求調(diào)劑 +11	陽(yáng)陽(yáng)哇塞 2026-03-27	12/600	2026-03-27 17:56 by yu221
[考研] 085602 289分求調(diào)劑 +8	WWW西西弗斯 2026-03-24	8/400	2026-03-26 16:33 by 不吃魚(yú)的貓
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4	鋼鐵大炮 2026-03-24	4/200	2026-03-26 16:28 by dick_runner
[考研] 材料專碩 335 分求調(diào)劑 +4	拒絕冷暴力 2026-03-25	4/200	2026-03-25 18:45 by haxia