lanlan0822

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 21.5
帖子: 19
在線: 3.3小時
蟲號: 856446

[交流] 【求助】急求關(guān)于matlab做非線性擬合的函數(shù)

我這里有三十組X1、X2與Y值，想擬合出它們之間的二元非線性回歸方程，不知道如何下手，望知道的童鞋指點。

X1=[367 357 348 338 333360 351 343 335 329 363 355 344 337 330 358 351 343 334
327 358 352 340 332 326 362 356 344 337 332]

X2=[600 630 650 680 710 625 645 678 712 720 640 635 688 679 727 642 664 693 67
5 740 630 648 689 710 742 620 643 667 701 710]

Y=[10.60 17.01 25.40 38.52 46.99 14.47 21.77 31.30 41.89 52.05 12.65 19.49 30.
28 40.32 51.33 16.11 22.22 32.02 44.25 57.64 15.68 21.09 34.85 47.01 59.45 13.
37 17.39 30.03 39.29 47.97]

如果用matlaB回歸的話，希望可以轉(zhuǎn)告程序！
不勝感激，多謝！我一共就32枚金幣，但是系統(tǒng)提示送不了30枚，所以就送25枚了
急求！

回復(fù)此樓

» 猜你喜歡

【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】已經(jīng)有7人回復(fù)
08工學調(diào)劑已經(jīng)有12人回復(fù)
材料與化工328分調(diào)劑已經(jīng)有4人回復(fù)
279分求調(diào)劑一志愿211 已經(jīng)有20人回復(fù)
293求調(diào)劑已經(jīng)有3人回復(fù)
一志愿華東理工大學081700，初試分數(shù)271 已經(jīng)有6人回復(fù)
085404電子信息284分求調(diào)劑已經(jīng)有3人回復(fù)
306求0703調(diào)劑一志愿華中師范已經(jīng)有9人回復(fù)
一志愿吉大化學322求調(diào)劑已經(jīng)有6人回復(fù)
求調(diào)劑已經(jīng)有5人回復(fù)

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

在matlab中擬合曲面的表達式已經(jīng)有29人回復(fù)
matlab擬合擬均相動力學參數(shù) 已經(jīng)有6人回復(fù)
需要用matlab擬合數(shù)據(jù)的外行人急需幫助已經(jīng)有10人回復(fù)
matlab非線性參數(shù)擬合問題已經(jīng)有7人回復(fù)
求助Matlab lsqcurve函數(shù)曲線擬合已經(jīng)有3人回復(fù)
【求助】使用MATLAB擬合出按照公式的系數(shù) 已經(jīng)有19人回復(fù)
如何用MATLAB擬合指數(shù)函數(shù) 已經(jīng)有3人回復(fù)
【原創(chuàng)】多項式和非多項式曲線擬合的MATLAB實現(xiàn)視頻教程已經(jīng)有10人回復(fù)
【求助】如何使得matlab擬合的曲線強制經(jīng)過一個點已經(jīng)有8人回復(fù)

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

查看全部散金貼

1樓 2011-01-04 09:18:05

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

lijinfeng042

木蟲 (小有名氣)

lanlan0822(金幣+1): 2011-01-05 09:01:46

方程的形式呢

回復(fù)此樓

2樓2011-01-04 13:03:28

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

rewind

至尊木蟲 (正式寫手)

lanlan0822(金幣+1): 2011-01-05 09:01:52

回歸的話您的方程形式要您設(shè)定不知道您的方程是什么

贊一下

回復(fù)此樓

3樓2011-01-04 13:12:18

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

藍云天

木蟲 (著名寫手)

應(yīng)助: 2 (幼兒園)
金幣: 4414.6
帖子: 2832
在線: 290.2小時
蟲號: 775959

★ ★
lanlan0822(金幣+5): 2011-01-05 09:01:31
xiegangmai(金幣+2):辛苦了，謝謝應(yīng)助！ 2011-01-07 22:07:27

樓主給你個例子，你自己琢磨琢磨吧。我覺得最好的學習方式就是按照一個例子自己摸索試驗。
例程：方程為F=K*X+C*Y+Z*（1-exp(-B*y^2)）-F0
F,X,Y,F0,都是已知的列向量，擬合參數(shù)K，C，Z，B
解：在Matlab下輸入：edit，然后將下面兩行百分號之間的內(nèi)容，復(fù)制進去，保存
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function F=zhidao_31(a,x)
%其中K,C,Z,B分別用a(1),a(2),a(3),a(4)代替
%X,Y,F0分別用x(:,1)，x(:,2)，x(:,3)代替
F=a(1)*x(:,1)+a(2)*x(:,2)+a(3)*(1-exp(-a(4)*x(:,2).^2))-x(:,3);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

以下是測試：在Matlab下面輸入：
%先由已經(jīng)的x,y,F0和已知的k,c,z,b來產(chǎn)生F，然后通過x,y,F0,F
%驗證一下：能否得到我們剛才所選用的參數(shù)k,c,z,b
%下面以隨機數(shù)為例
x=rand(100,1);
y=rand(100,1);
F0=rand(100,1);

k=1;
c=1.5;
z=2.5;
b=3.5;

X=[x,y,F0];
F=zhidao_31([k,c,z,b],X);
[A,res]=lsqcurvefit('zhidao_31',ones(1,4),X,F);
A

得到的結(jié)果：
A =
1.0033 1.4446 2.5520 3.4853
跟k,c,z,b的值大致一樣。
------------------------------------------
例題中的問題：
設(shè)好x,y,F0,F后，在Matlab下輸入：
[A,res]=lsqcurvefit('zhidao_31',ones(1,4),[x,y,F0],F); %注意x,y大小寫。
A

就可以了。

贊一下(1人)

回復(fù)此樓

5樓2011-01-04 13:46:08

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

lanlan0822

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 21.5
帖子: 19
在線: 3.3小時
蟲號: 856446

關(guān)于matlab做非線性擬合的函數(shù)

我想得出一個二元二次的回歸方程，但是找到的都是只能進行線性回歸，或者是二元二次回歸的。

贊一下

回復(fù)此樓

6樓2011-01-04 15:29:37

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

lanlan0822

銅蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 21.5
帖子: 19
在線: 3.3小時
蟲號: 856446

[關(guān)于matlab做非線性擬合的函數(shù)

我想得出一個二元二次的回歸方程，但是找到的都是只能進行線性回歸，或者是二元二次回歸的。

贊一下

回復(fù)此樓

7樓2011-01-04 15:30:35

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

feixiaolin

榮譽版主 (文壇精英)

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

y=p0+p1*x1+p2*x2;p1
p0 = 352.718330203856
p1 = -1.03742355830517
p2 = 0.0544705855433249
相關(guān)系數(shù)(R): 0.988109884761729

y=p0+p1*x1+p2*x2+p3*x1^2+p4*x2^2+p5*x1*x2;
p0 = 5409.56927884948
p1 = -22.0861948342262
p2 = -4.06582793544753
p3 = 0.024111003125518
p4 = 0.00140539387129161
p5 = 0.0063262492366287
相關(guān)系數(shù)(R): 0.999477545777094

贊一下

回復(fù)此樓

8樓2013-11-01 20:38:21

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

cobrasq

金蟲 (小有名氣)

應(yīng)助: 9 (幼兒園)
金幣: 570.6
帖子: 83
在線: 33.9小時
蟲號: 144382

★ ★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
1592203609: 金幣+1, 謝謝回帖 2013-11-02 12:18:48

假設(shè)你預(yù)期的二元二次方程為
y = b1 + b2*X1 + b3*X2 + b4*X1*X2 + b5*X1^2 + b6*X2^2

程序如下
X1=[367 357 348 338 333360 351 343 335 329 363 355 344 337 330 358 351 343 334 327 358 352 340 332 326 362 356 344 337 332 357];

X2=[600 630 650 680 710 625 645 678 712 720 640 635 688 679 727 642 664 693 675 740 630 648 689 710 742 620 643 667 701 710];

Y=[10.60 17.01 25.40 38.52 46.99 14.47 21.77 31.30 41.89 52.05 12.65 19.49 30.28 40.32 51.33 16.11 22.22 32.02 44.25 57.64 15.68 21.09 34.85 47.01 59.45 13.37 17.39 30.03 39.29 47.97];

modelfun = @(b,x)b(1)+b(2)*x(1,

+b(3)*x(2,

+b(4)*x(1,

.*x(2,

+b(5)*x(1,

.^2+b(6)*x(2,

.^2;
beta0 = [1 2 3 4 5 6];
mdl = nlinfit([X1;X2],Y,modelfun,beta0)

如果效果不好，可以考慮加權(quán)函數(shù)或啟用robust.

贊一下(1人)

回復(fù)此樓

9樓2013-11-02 05:09:24

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

簡單回復(fù)

zyxme4樓

2011-01-04 13:20 回復(fù)

robert2020(金幣-2):為了他人的方便，請勿在求助帖中純表無意義回復(fù)！ 2011-01-11 10:34:40

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

[考研] 求調(diào)劑 +4	研研，接電話 2026-03-24	5/250	2026-03-24 10:17 by 研研，接電話
[考研] 環(huán)境學碩288求調(diào)劑 +8	皮皮皮123456 2026-03-22	8/400	2026-03-23 23:47 by 熱情沙漠
[考研] 一志愿武理材料工程348求調(diào)劑 +6	￣^￣゜汗 2026-03-19	9/450	2026-03-23 19:53 by pswait
[考研] 招08考數(shù)學 +6	laoshidan 2026-03-20	14/700	2026-03-23 14:37 by 15614359529
[考研] 298求調(diào)劑 +8	上岸6666@ 2026-03-20	8/400	2026-03-23 11:02 by laoshidan
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調(diào)劑 +4	幸運的醬醬 2026-03-22	5/250	2026-03-22 20:15 by edmund7
[考研] 280求調(diào)劑 +11	咕嚕曉曉 2026-03-18	12/600	2026-03-21 22:40 by ACS Nano——
[考研] 330求調(diào)劑0854 +3	assdll 2026-03-21	3/150	2026-03-21 13:01 by 搏擊518
[考研] 332求調(diào)劑 +3	鳳凰院丁真 2026-03-20	3/150	2026-03-21 10:27 by luoyongfeng
[考研] 華東師范大學-071000生物學-293分-求調(diào)劑 +3	研究生何瑤明 2026-03-18	3/150	2026-03-21 01:30 by JourneyLucky
[考研] 一志愿西南交大，求調(diào)劑 +5	材化逐夢人 2026-03-18	5/250	2026-03-21 00:26 by JourneyLucky
[考研] 一志愿南昌大學，327分，材料與化工085600 +9	Ncdx123456 2026-03-19	9/450	2026-03-20 23:41 by lovewei0727
[考研] 一志愿西北大學，070300化學學碩，總分287，雙非一本，求調(diào)劑。 +4	晨昏線與星海 2026-03-19	4/200	2026-03-20 22:15 by JourneyLucky
[考研] 求調(diào)劑，一志愿:南京航空航天大學大學，080500材料科學與工程學碩，總分289分 +4	@taotao 2026-03-19	4/200	2026-03-20 22:14 by JourneyLucky
[考研] 材料與化工 322求調(diào)劑 +4	然11 2026-03-19	4/200	2026-03-20 22:12 by luoyongfeng
[考研] 295復(fù)試調(diào)劑 +8	簡木ChuFront 2026-03-19	8/400	2026-03-20 20:44 by zhukairuo
[論文投稿] 申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[考研] 一志愿中國海洋大學，生物學，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	6/300	2026-03-19 23:46 by zcl123
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 一志愿福大288有機化學，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu