chj787

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 933
帖子: 104
在線: 42.2小時(shí)
蟲號(hào): 798564
注冊(cè): 2009-06-24
專業(yè): 原子和分子物理

[求助] matlab擬合求助

用matlab擬合一個(gè)方程：
y=(a+bx+cx^2+dx^3)exp(-bx)
x=[0.8,0.85,0.9,0.94,1,1.1,1.14,1.144,1.146,1.148,1.149,1.15,1.151,1.152,1.153,1.157,1.25,1.3,1.35,1.4,1.45,1.5,1.55,1.6,1.65,1.7,1.75,1.8,1.85,1.9,1.95,2,2.05,2.1,2.15,2.2,2.25,2.3,2.35,2.4];
y=[129698.6607,81517.14907,48462.38869,30160.54842,12625.39888,830.7904906,0,2.343989057,8.695584894,18.42928479,24.54823752,31.4946096,39.26181681,47.84108015,57.22581538,102.7141275,687.9234793,3872.493796,7589.20011,11578.99871,15654.88623,19681.73873,23751.25051,27388.0944,30756.13673,33816.1486,36542.19252,38921.37873,40955.22819,42660.53072,44067.23581,45213.10408,46138.07991,46880.36067,47474.18002,47948.97387,48329.25318,48634.92208,48881.88591,49082.65691];

變量：x,y
擬合常數(shù)： a,b,c,d
求高手幫助！給出matlab擬合程序。

[ Last edited by chj787 on 2012-5-20 at 16:29 ]

回復(fù)此樓

» 猜你喜歡

303求調(diào)劑已經(jīng)有3人回復(fù)
有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～已經(jīng)有7人回復(fù)
梁成偉老師課題組歡迎你的加入已經(jīng)有10人回復(fù)
一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) 已經(jīng)有8人回復(fù)
070300化學(xué)319求調(diào)劑已經(jīng)有3人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有4人回復(fù)
268求調(diào)劑已經(jīng)有7人回復(fù)
290求調(diào)劑已經(jīng)有11人回復(fù)
材料與化工求調(diào)劑已經(jīng)有4人回復(fù)
293求調(diào)劑已經(jīng)有11人回復(fù)

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

matlab擬合求值已經(jīng)有17人回復(fù)
在matlab中擬合曲面的表達(dá)式已經(jīng)有29人回復(fù)
matlab 求指點(diǎn) 動(dòng)力學(xué)方程擬合過程中導(dǎo)數(shù)的獲取已經(jīng)有10人回復(fù)
求助大俠幫忙，有關(guān)matlab模擬數(shù)據(jù) 已經(jīng)有29人回復(fù)
【求助】用matlab最優(yōu)化方法進(jìn)行參數(shù)擬合已經(jīng)有17人回復(fù)
【求助】使用MATLAB擬合出按照公式的系數(shù) 已經(jīng)有19人回復(fù)
【求助】請(qǐng)教一個(gè)matlab程序，謝謝幫忙��！已經(jīng)有15人回復(fù)
【求助】matlab曲面擬合表達(dá)式已經(jīng)有10人回復(fù)
【求助】如何使得matlab擬合的曲線強(qiáng)制經(jīng)過一個(gè)點(diǎn) 已經(jīng)有8人回復(fù)
【求助】使用MATLAB怎么實(shí)現(xiàn)擬合力場參數(shù)的程序？【已解決】已經(jīng)有15人回復(fù)
【求助】求MAtlab求解程序�。。� 已經(jīng)有20人回復(fù)

1樓 2012-05-20 16:13:46

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

回帖支持 ( 顯示支持度最高的前 50 名 )

cactus586

金蟲 (著名寫手)

流浪過的人總是想再次踏上旅途

應(yīng)助: 95 (初中生)
貴賓: 0.036
金幣: 143.8
散金: 5509
紅花: 35
帖子: 2929
在線: 534.9小時(shí)
蟲號(hào): 484017
注冊(cè): 2007-12-26
性別: GG
專業(yè): 自然地理學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
xzhdty: 金幣+1, 歡迎常來程序語言看看 2012-05-21 05:57:43
chj787: 金幣+30, ★有幫助, 1 2012-05-24 21:55:56

昨天看到一個(gè)用MatLab演奏卡儂的可以站短問問他
http://www.gaoyang168.com/bbs/viewthread.php?tid=4519096

贊一下(1人)

回復(fù)此樓

靜水流深

2樓2012-05-20 16:40:48

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

普通回帖

chj787

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 933
帖子: 104
在線: 42.2小時(shí)
蟲號(hào): 798564
注冊(cè): 2009-06-24
專業(yè): 原子和分子物理

引用回帖:

2樓: Originally posted by cactus586 at 2012-05-20 16:40:48:
昨天看到一個(gè)用MatLab演奏卡儂的可以站短問問他
http://www.gaoyang168.com/bbs/viewthread.php?tid=4519096

謝謝！

回復(fù)此樓

3樓2012-05-20 17:59:36

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

lijie169

銅蟲 (著名寫手)

應(yīng)助: 159 (高中生)
金幣: 161.5
散金: 2018
紅花: 7
帖子: 1338
在線: 491.5小時(shí)
蟲號(hào): 1180408
注冊(cè): 2010-12-30
專業(yè): 分離過程

【答案】應(yīng)助回帖

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
xzhdty: 金幣+2, 歡迎常來程序語言 2012-05-21 05:58:00
chj787: 金幣+3, ★★★很有幫助 2012-05-29 15:30:03

CODE:

function f=fit(x,y)

f=lsqnonlin(@(a)fun(a,x,y),[1,1,1,1]');

function f1=fun(a,x,y)

 f1=y-(a(1)+a(2)*x+a(3)*x.^2+a(4)*x.^3).*exp(-a(2)*x);

end

end

x,y就是樓主輸入的,a是一個(gè)矩陣，對(duì)應(yīng)著abcd

贊一下

回復(fù)此樓

4樓2012-05-20 18:08:21

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

lijie169

銅蟲 (著名寫手)

應(yīng)助: 159 (高中生)
金幣: 161.5
散金: 2018
紅花: 7
帖子: 1338
在線: 491.5小時(shí)
蟲號(hào): 1180408
注冊(cè): 2010-12-30
專業(yè): 分離過程

【答案】應(yīng)助回帖

額，輸入fit(x,y)就可以了

贊一下

回復(fù)此樓

5樓2012-05-20 18:09:06

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

coolxh

木蟲 (小有名氣)

應(yīng)助: 6 (幼兒園)
金幣: 3261.8
散金: 969
紅花: 4
帖子: 184
在線: 208小時(shí)
蟲號(hào): 318919
注冊(cè): 2007-03-07
專業(yè): 交通工程

【答案】應(yīng)助回帖

★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
xzhdty: 金幣+1, 歡迎常來程序語言看看 2012-05-21 05:58:22
chj787: 金幣+2, ★★★很有幫助 2012-05-29 15:29:52

若不想自己寫代碼，亦可用曲線擬合工具，
cftool(x,y)

然后采用自定義公式，按你想要的公式定義表達(dá)式，進(jìn)行擬合

然后再自動(dòng)生成擬合代碼

贊一下

回復(fù)此樓

6樓2012-05-20 21:43:41

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

lgycjpcqu

金蟲 (正式寫手)

應(yīng)助: 35 (小學(xué)生)
金幣: 566.8
散金: 411
帖子: 371
在線: 176.9小時(shí)
蟲號(hào): 1327328
注冊(cè): 2011-06-20
性別: GG
專業(yè): 電化學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
xzhdty: 金幣+1, 謝謝了 2012-05-22 06:01:48
chj787: 金幣+15, ★★★★★最佳答案, 非常感謝！評(píng)分晚了，萬分抱歉~！ 2012-05-29 15:29:20

x=[0.8,0.85,0.9,0.94,1,1.1,1.14,1.144,1.146,1.148,1.149,1.15,1.151,1.152,1.153,1.157,1.25,1.3,1.35,1.4,1.45,1.5,1.55,1.6,1.65,1.7,1.75,1.8,1.85,1.9,1.95,2,2.05,2.1,2.15,2.2,2.25,2.3,2.35,2.4];
y=[129698.6607,81517.14907,48462.38869,30160.54842,12625.39888,830.7904906,0,2.343989057,8.695584894,18.42928479,24.54823752,31.4946096,39.26181681,47.84108015,57.22581538,102.7141275,687.9234793,3872.493796,7589.20011,11578.99871,15654.88623,19681.73873,23751.25051,27388.0944,30756.13673,33816.1486,36542.19252,38921.37873,40955.22819,42660.53072,44067.23581,45213.10408,46138.07991,46880.36067,47474.18002,47948.97387,48329.25318,48634.92208,48881.88591,49082.65691];
ybeta=inline('(beta(1)+beta(2)*x+beta(3)*x.^2+beta(4)*x.^3).*exp(-beta(2)*x)','beta','x');
beta=nlinfit(x,y,ybeta,[1e+004,0,-1e+004,1e+004]);
a=beta(1);
b=beta(2);
c=beta(3);
d=beta(4);
yy=(a+b*x+c*x.^2+d*x.^3).*exp(-b*x);
plot(x,yy,x,y,'*')
但是會(huì)顯示
Warning: The Jacobian at the solution is ill-conditioned, and some
model parameters may not be estimated well (they are not identifiable).
Use caution in making predictions.所以可能還需要進(jìn)一步手動(dòng)優(yōu)化，我做你何時(shí)用nlinfit也經(jīng)常出現(xiàn)這個(gè)警告，期待有高手解答。
擬合結(jié)果是
3866678.20479525
2.55862906421524
-8281766.51676900
4642226.52606192
希望能幫到你

贊一下

回復(fù)此樓

7樓2012-05-21 10:31:55

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

lijie169

銅蟲 (著名寫手)

應(yīng)助: 159 (高中生)
金幣: 161.5
散金: 2018
紅花: 7
帖子: 1338
在線: 491.5小時(shí)
蟲號(hào): 1180408
注冊(cè): 2010-12-30
專業(yè): 分離過程

引用回帖:

4樓: Originally posted by lijie169 at 2012-05-20 18:08:21
function f=fit(x,y)
f=lsqnonlin(@(a)fun(a,x,y),');
function f1=fun(a,x,y)
f1=y-(a(1)+a(2)*x+a(3)*x.^2+a(4)*x.^3).*exp(-a(2)*x);
end
end

x,y就是樓主輸入的,a是一個(gè)矩陣，對(duì)應(yīng)著abcd...

很有幫助才給3個(gè)金幣

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

chj787

8樓2012-05-29 17:14:23

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

chj787

金蟲 (小有名氣)

送鮮花一朵

引用回帖:

8樓: Originally posted by lijie169 at 2012-05-29 17:14:23
很有幫助才給3個(gè)金幣

...

9樓2012-05-30 22:10:20

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

lijie169

銅蟲 (著名寫手)

引用回帖:

9樓: Originally posted by chj787 at 2012-05-30 22:10:20

...

10樓2012-05-31 07:04:05

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

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

[考研] 070300化學(xué)319求調(diào)劑 +3	錦鯉0909 2026-03-17	3/150	2026-03-17 15:01 by 我的船我的海
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +4	yangfz 2026-03-17	4/200	2026-03-17 15:00 by 我的船我的海
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 341求調(diào)劑 +5	搗蛋豬豬 2026-03-11	5/250	2026-03-17 06:56 by Muziwl
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進(jìn)步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[考研] 279求調(diào)劑 +3	Dizzy123@ 2026-03-10	3/150	2026-03-13 23:02 by JourneyLucky
[考研] 304求調(diào)劑 +6	Mochaaaa 2026-03-12	7/350	2026-03-13 22:18 by 星空星月
[考研] 304求調(diào)劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[考研] 329求調(diào)劑 +3	miaodesi 2026-03-12	4/200	2026-03-13 20:53 by 18595523086
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 材料工程調(diào)劑 +4	咪咪空空 2026-03-11	4/200	2026-03-13 19:57 by JourneyLucky
[考研] 0703化學(xué)求調(diào)劑 +7	綠豆芹菜湯 2026-03-12	7/350	2026-03-13 17:25 by njzyff
[考研] 085600材料與化工 309分請(qǐng)求調(diào)劑 +7	dtdxzxx 2026-03-12	8/400	2026-03-13 14:43 by jxchenghu
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-12	3/150	2026-03-13 14:30 by 求調(diào)劑zz
[考研] 一志愿山大07化學(xué) 332分四六級(jí)已過本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-12	3/150	2026-03-13 14:18 by JourneyLucky
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike