wgy318

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1161.7
散金: 1024
紅花: 1
帖子: 84
在線: 26.5小時(shí)
蟲號(hào): 587490
注冊(cè): 2008-08-22
專業(yè): 金屬材料的制備科學(xué)與跨學(xué)

[交流] 【求助】數(shù)據(jù)擬合

我想擬合一個(gè)兩個(gè)自變量的模型的六個(gè)參數(shù)，編輯了如下的程序，可是系統(tǒng)總是報(bào)錯(cuò)，請(qǐng)大家看看程序編寫的是否正確,謝謝！
我是新手，請(qǐng)大家多給點(diǎn)意見。

M文件如下：

CODE:

function F=myfun(a,data)

x(1)=data(1,:);

x(2)=data(2,:);

F=a(3)*exp((-a(5))*(x(2)-a(4))/(a(6)+(x(2)-a(4))))/(1+(a(3)*exp((-a(5))*(x(2)-a(4))/(a(6)+(x(2)-a(4))))*x(1)/a(1))^(1-a(2))) 

命令窗口程序如下：

>>  x1=[40 60 80 100 150 200 300 400 500 600 800 1000 1200 40 60 80 100 150 200 300 400 500 600 800 1000 1200 40 60 80 100 150 200 300 400 500 600 800 1000 1200 40 60 80 100 150 200 300 400 500 600 800 1000 1200 40 60 80 100 150 200 300 400 500 600 800 1000 1200];

x2=[423.15,423.15,423.15,425.15,425.15,425.15,425.15,425.15,425.15,425.15,425.15,425.15,425.15,431.15,431.15,431.15,431.15,431.15,431.15,431.15,431.15,431.15,431.15,431.15,431.15,431.15,437.15,437.15,437.15,437.15,437.15,437.15,437.15,437.15,437.15,437.15,437.15,437.15,437.15,443.15,443.15,443.15,443.15,443.15,443.15,443.15,443.15,443.15,443.15,443.15,443.15,443.15,449.15,449.15,449.15,449.15,449.15,449.15,449.15,449.15,449.15,449.15,449.15,449.15,449.15];

 ydata=[500 367.4 316.4 351 272 249 207 180 128 121 86.2 80.8 69.1 172.9 142.4 122.1 114 95 84 69.2 60 53 48 41.2 36.2 31.6 152.9 122.4 107.1 98 84 71 58.2 51 45.6 41.8 36.4 32.7 30.0 132.5 115.4 101.9 90 76 69 57.1 51 44.9 40.8 35.2 31.0 28.2 112.1 101.9 86.7 82 68 61 50.3 45 41 38 32.6 29 26.2];

data=[x1;x2];

a0=[5000 0.3 3.5*10e13 300 20 1.5];

[a,resnorm]=lsqcurvefit(@myfun,a0,data,ydata)

[ Last edited by nono2009 on 2010-9-14 at 06:02 ]

回復(fù)此樓

» 猜你喜歡

337求調(diào)劑已經(jīng)有3人回復(fù)
329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601）已經(jīng)有10人回復(fù)
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑已經(jīng)有16人回復(fù)
291求調(diào)劑已經(jīng)有6人回復(fù)
299求調(diào)劑已經(jīng)有10人回復(fù)
面上5B能上會(huì)嗎？已經(jīng)有4人回復(fù)
322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過六級(jí)。已經(jīng)有4人回復(fù)
085600 材料與化工 329分求調(diào)劑已經(jīng)有14人回復(fù)
材料科學(xué)與工程求調(diào)劑已經(jīng)有3人回復(fù)
289求調(diào)劑已經(jīng)有5人回復(fù)

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

【求助】線性擬合已經(jīng)有2人回復(fù)
【求助】阻抗擬合已經(jīng)有9人回復(fù)
【求助】Tafel 擬合數(shù)據(jù)輸出辦法已經(jīng)有1人回復(fù)
【求助】matlab 擬合問題已經(jīng)有1人回復(fù)
【求助】關(guān)于數(shù)據(jù)擬合，請(qǐng)教數(shù)學(xué)達(dá)人已經(jīng)有4人回復(fù)
【求助】XPS擬合軟件及方法已經(jīng)有12人回復(fù)
【求助】數(shù)據(jù)擬合，函數(shù)系數(shù)對(duì)函數(shù)的影響已經(jīng)有0人回復(fù)
【求助】幫忙擬合一下已經(jīng)有7人回復(fù)
【求助】數(shù)據(jù)擬合軟件及教程已經(jīng)有4人回復(fù)
【求助】關(guān)于數(shù)據(jù)擬合轉(zhuǎn)換成origin圖已經(jīng)有5人回復(fù)
【求助】SiTF數(shù)據(jù)擬合已經(jīng)有0人回復(fù)
【求助】MRTD數(shù)據(jù)擬合已經(jīng)有0人回復(fù)
【求助】XPS數(shù)據(jù)擬合已經(jīng)有2人回復(fù)
【求助】ORIGIN 擬合已經(jīng)有5人回復(fù)
【求助】關(guān)于擬合已經(jīng)有4人回復(fù)
【求助】數(shù)據(jù)擬合已經(jīng)有1人回復(fù)
【求助】數(shù)據(jù)擬合已經(jīng)有11人回復(fù)

1樓 2010-09-07 10:53:22

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

lidaxue

木蟲 (正式寫手)

之乎者也

應(yīng)助: 10 (幼兒園)
金幣: 1952.1
散金: 612
紅花: 2
帖子: 527
在線: 856.1小時(shí)
蟲號(hào): 259962
注冊(cè): 2006-06-17
性別: GG
專業(yè): 幾何學(xué)

wgy318(金幣+3):問題已經(jīng)解決了，不過還是要謝謝你 2010-09-09 08:53:57

請(qǐng)樓主給出報(bào)錯(cuò)信息！

贊一下

回復(fù)此樓

Comeon！

2樓2010-09-08 10:27:17

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

孤鴻飄逸

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 5861.4
散金: 499
紅花: 17
帖子: 617
在線: 432.5小時(shí)
蟲號(hào): 864807
注冊(cè): 2009-10-07
性別: GG
專業(yè): 機(jī)器人學(xué)及機(jī)器人技術(shù)

★
nono2009:編輯內(nèi)容 2010-09-13 06:18
nono2009(金幣+1):用[code][/code]貼程序的效果比較好。 2010-09-13 06:19:15
wgy318(金幣+6):問題已經(jīng)解決，非常感謝 2010-09-13 08:43:58

函數(shù)寫為：（注意是點(diǎn)乘，點(diǎn)除）

CODE:

function F=myfun(a,data)

xx=data(1,:);

xxx=data(2,:);

F=a(3).*exp((-a(5))*(xxx-a(4))./(a(6)+(xxx-a(4))))./(1+(a(3).*exp((-a(5)).*(xxx-a(4))./(a(6)+(xxx-a(4)))).*xx/a(1)).^(1-a(2)))

輸入?yún)?shù)值運(yùn)行結(jié)果為：

Maximum number of function evaluations exceeded;
increase options.MaxFunEvals

a =
1.0e+014 *

0.0000 0.0000 3.5000 0.0000 0.0000 0.0000

resnorm =

3.9906e+004
我這臺(tái)機(jī)子是實(shí)驗(yàn)室的matlab6.1，估計(jì)版本太差，你用高版本的運(yùn)行看看

[ Last edited by nono2009 on 2010-9-13 at 06:18 ]

贊一下(1人)

回復(fù)此樓

3樓2010-09-12 18:11:56

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

lijinfeng042

木蟲 (小有名氣)

Matlab

仿真EPI: 2
應(yīng)助: 1 (幼兒園)
金幣: 2156.1
散金: 115
帖子: 291
在線: 31.5小時(shí)
蟲號(hào): 1019062
注冊(cè): 2010-05-15
性別: GG
專業(yè): 通信理論與系統(tǒng)

★ ★
robert2020(金幣+2):辛苦了！鼓勵(lì)繼續(xù)應(yīng)助！ 2010-09-13 18:49:53
wgy318(金幣+2): 2010-09-13 21:07:35

這個(gè)函數(shù)F擬合的效果不太理想啊 resnorm = 3.9906e+004

贊一下(1人)

回復(fù)此樓

工作了,偶爾會(huì)上來~可以關(guān)注新浪微博 @云是風(fēng)的夢(mèng)_Matlab

4樓2010-09-13 09:18:16

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

wgy318

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1161.7
散金: 1024
紅花: 1
帖子: 84
在線: 26.5小時(shí)
蟲號(hào): 587490
注冊(cè): 2008-08-22
專業(yè): 金屬材料的制備科學(xué)與跨學(xué)

resnorm要多少才算理想

贊一下

回復(fù)此樓

5樓2010-09-13 21:09:18

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

lijinfeng042

木蟲 (小有名氣)

Matlab

仿真EPI: 2
應(yīng)助: 1 (幼兒園)
金幣: 2156.1
散金: 115
帖子: 291
在線: 31.5小時(shí)
蟲號(hào): 1019062
注冊(cè): 2010-05-15
性別: GG
專業(yè): 通信理論與系統(tǒng)

wgy318(金幣+2): 2010-09-14 09:06:49

引用回帖:

Originally posted by wgy318 at 2010-09-13 21:09:18:
resnorm要多少才算理想

resnorm越小是表示擬合精度越高。你試試用別的優(yōu)化函數(shù)進(jìn)行擬合，或者就是你的模型沒建立好,殘差太大啊

贊一下

回復(fù)此樓

工作了,偶爾會(huì)上來~可以關(guān)注新浪微博 @云是風(fēng)的夢(mèng)_Matlab

6樓2010-09-13 23:10:04

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

wgy318

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1161.7
散金: 1024
紅花: 1
帖子: 84
在線: 26.5小時(shí)
蟲號(hào): 587490
注冊(cè): 2008-08-22
專業(yè): 金屬材料的制備科學(xué)與跨學(xué)

能不能給你一個(gè)例子或者一個(gè)大體的程序框架啊？matlab里面應(yīng)該有優(yōu)化函數(shù)的程序可以調(diào)用吧。我現(xiàn)在函數(shù)已經(jīng)有了，只是想優(yōu)化下擬合的程度。但是我matlab又不怎么懂，不知道該怎么進(jìn)一步優(yōu)化。請(qǐng)多多指點(diǎn)

贊一下

回復(fù)此樓

7樓2010-09-14 09:09:29

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

lijinfeng042

木蟲 (小有名氣)

Matlab

仿真EPI: 2
應(yīng)助: 1 (幼兒園)
金幣: 2156.1
散金: 115
帖子: 291
在線: 31.5小時(shí)
蟲號(hào): 1019062
注冊(cè): 2010-05-15
性別: GG
專業(yè): 通信理論與系統(tǒng)

wgy318(金幣+1): 2010-09-15 15:52:47

引用回帖:

Originally posted by wgy318 at 2010-09-14 09:09:29:
能不能給你一個(gè)例子或者一個(gè)大體的程序框架��？matlab里面應(yīng)該有優(yōu)化函數(shù)的程序可以調(diào)用吧。我現(xiàn)在函數(shù)已經(jīng)有了，只是想優(yōu)化下擬合的程度。但是我matlab又不怎么懂，不知道該怎么進(jìn)一步優(yōu)化。請(qǐng)多多指點(diǎn)

就是這個(gè)函數(shù)?
F=a(3).*exp((-a(5))*(xxx-a(4))./(a(6)+(xxx-a(4))))./(1+(a(3).*exp((-a(5)).*(xxx-a(4))./(a(6)+(xxx-a(4)))).*xx/a(1)).^(1-a(2)))

贊一下

回復(fù)此樓

工作了,偶爾會(huì)上來~可以關(guān)注新浪微博 @云是風(fēng)的夢(mèng)_Matlab

8樓2010-09-14 19:07:58

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

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

[考研] 一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑 +16	稻妻小編 2026-03-26	16/800	2026-03-29 18:15 by Reek buddy
[基金申請(qǐng)] 面上5B能上會(huì)嗎？ +4	redcom 2026-03-29	4/200	2026-03-29 17:39 by 最愛編纂辭典
[考研] 南京大學(xué)化學(xué)調(diào)劑 +6	景隨風(fēng) 2026-03-29	9/450	2026-03-29 15:27 by 唐沐兒
[考研] 一志愿鄭州大學(xué)，080500學(xué)碩，總分317分求調(diào)劑 +8	舉個(gè)栗子oi 2026-03-24	9/450	2026-03-29 13:08 by peike
[考研] 數(shù)一英一271專碩（085401）求調(diào)劑，可跨 +7	前行必有光 2026-03-28	8/400	2026-03-28 23:22 by 小木蟲tim
[考研] 295求調(diào)劑 +4	wei-5 2026-03-26	4/200	2026-03-28 23:20 by 小木蟲tim
[考研] 317求調(diào)劑 +6	十閑wx 2026-03-24	6/300	2026-03-28 13:27 by Iveryant
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +9	丹青奶蓋 2026-03-26	10/500	2026-03-28 07:45 by barnett0632
[考研] 085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富 +4	來財(cái)aa 2026-03-25	4/200	2026-03-28 07:41 by 棒棒球手
[考研] 265求調(diào)劑 +8	小木蟲085600 2026-03-27	8/400	2026-03-27 22:16 by 無際的草原
[考研] 化學(xué)308分求調(diào)劑 +8	你好明天你好 2026-03-23	9/450	2026-03-27 14:01 by 楊光于青云
[考研] 298調(diào)劑 +3	jiyingjie123 2026-03-27	3/150	2026-03-27 11:57 by wxiongid
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4	崔wj 2026-03-26	4/200	2026-03-27 08:04 by chemisry
[考研] 342求調(diào)劑 +3	加油a李zs 2026-03-26	3/150	2026-03-27 00:29 by wxiongid
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5	星沉uu 2026-03-26	6/300	2026-03-26 19:02 by macy2011
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3	江往賣魚 2026-03-26	3/150	2026-03-26 09:42 by 王小欠i
[考研] 293求調(diào)劑 +7	加一一九 2026-03-24	7/350	2026-03-25 12:02 by userper
[考研] 化工專碩求調(diào)劑 +3	question挽風(fēng) 2026-03-24	3/150	2026-03-24 18:48 by jhhcooi
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5	番茄頭--- 2026-03-22	6/300	2026-03-23 23:45 by Txy@872106
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長(zhǎng)期招生公告 +3	李子目 2026-03-22	3/150	2026-03-22 21:01 by 怎么釋懷