wgy318

金蟲 (小有名氣)

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

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

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

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ù)此樓

» 猜你喜歡

調(diào)劑310 已經(jīng)有12人回復(fù)
化學(xué)0703 調(diào)劑 306分一志愿211 已經(jīng)有6人回復(fù)
26申博自薦已經(jīng)有4人回復(fù)
化學(xué)工程085602 305分求調(diào)劑已經(jīng)有20人回復(fù)
一志愿南昌大學(xué)324求調(diào)劑已經(jīng)有5人回復(fù)
289求調(diào)劑已經(jīng)有3人回復(fù)
085601材料工程找調(diào)劑已經(jīng)有14人回復(fù)
一志愿211，335分，0856，求調(diào)劑院校和導(dǎo)師已經(jīng)有6人回復(fù)
334分一志愿武理材料求調(diào)劑已經(jīng)有7人回復(fù)
275求調(diào)劑已經(jīng)有18人回復(fù)

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

【求助】線性擬合已經(jīng)有2人回復(fù)
【求助】阻抗擬合已經(jīng)有9人回復(fù)
【求助】Tafel 擬合數(shù)據(jù)輸出辦法已經(jīng)有1人回復(fù)
【求助】matlab 擬合問題已經(jīng)有1人回復(fù)
【求助】關(guān)于數(shù)據(jù)擬合，請教數(shù)學(xué)達人已經(jīng)有4人回復(fù)
【求助】XPS擬合軟件及方法已經(jīng)有12人回復(fù)
【求助】數(shù)據(jù)擬合，函數(shù)系數(shù)對函數(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小時
蟲號: 259962
注冊: 2006-06-17
性別: GG
專業(yè): 幾何學(xué)

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

請樓主給出報錯信息！

贊一下

回復(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小時
蟲號: 864807
注冊: 2009-10-07
性別: GG
專業(yè): 機器人學(xué)及機器人技術(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ù)寫為：（注意是點乘，點除）

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ù)值運行結(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
我這臺機子是實驗室的matlab6.1，估計版本太差，你用高版本的運行看看

[ 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小時
蟲號: 1019062
注冊: 2010-05-15
性別: GG
專業(yè): 通信理論與系統(tǒng)

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

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

贊一下(1人)

回復(fù)此樓

工作了,偶爾會上來~可以關(guān)注新浪微博 @云是風(fē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小時
蟲號: 587490
注冊: 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小時
蟲號: 1019062
注冊: 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ù)進行擬合，或者就是你的模型沒建立好,殘差太大啊

贊一下

回復(fù)此樓

工作了,偶爾會上來~可以關(guān)注新浪微博 @云是風(fē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小時
蟲號: 587490
注冊: 2008-08-22
專業(yè): 金屬材料的制備科學(xué)與跨學(xué)

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

贊一下

回復(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小時
蟲號: 1019062
注冊: 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:
能不能給你一個例子或者一個大體的程序框架��？matlab里面應(yīng)該有優(yōu)化函數(shù)的程序可以調(diào)用吧。我現(xiàn)在函數(shù)已經(jīng)有了，只是想優(yōu)化下擬合的程度。但是我matlab又不怎么懂，不知道該怎么進一步優(yōu)化。請多多指點

就是這個函數(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ù)此樓

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

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

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

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

[考研] 085601材料工程找調(diào)劑 +13	oatmealR 2026-03-29	14/700	2026-03-29 16:41 by 學(xué)員8dgXkO
[考研] 085600 286分材料求調(diào)劑 +8	麻辣魷魚 2026-03-27	9/450	2026-03-29 15:36 by Serene1974
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5	大饞小子 2026-03-28	5/250	2026-03-29 08:56 by qingfeng258
[考研] 0856求調(diào)劑 +7	楒桉 2026-03-28	7/350	2026-03-29 08:28 by fmesaito
[考研] 一志愿華北電力大學(xué)能動專碩，293，求調(diào)劑 +3	15537177284 2026-03-23	5/250	2026-03-28 16:11 by xxxsssccc
[考研] 291求調(diào)劑 +15	hhhhxn.. 2026-03-23	21/1050	2026-03-28 11:26 by self2008
[考研] 材料與化工考研調(diào)劑 +17	孅華 2026-03-22	17/850	2026-03-28 08:35 by WYUMater
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8	曼111 2026-03-24	9/450	2026-03-28 07:58 by YYYYX1234
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級，六級521，求調(diào)劑 +4	美味可樂雞翅 2026-03-26	4/200	2026-03-27 15:27 by 星空星月
[考研] 085600材料與化工調(diào)劑 +10	A-哆啦Z夢 2026-03-23	16/800	2026-03-27 15:13 by caszguilin
[考研] 322求調(diào)劑 +4	我真的很想學(xué)習(xí) 2026-03-23	4/200	2026-03-27 13:51 by 楊楊楊紫
[考研] 材料求調(diào)劑 +5	.m.. 2026-03-25	5/250	2026-03-27 11:08 by 不吃魚的貓
[考研] 325求調(diào)劑 +3	Aoyijiang 2026-03-23	3/150	2026-03-26 20:46 by 不吃魚的貓
[考研] 081700 調(diào)劑 267分 +11	迷人的哈哈 2026-03-23	11/550	2026-03-26 15:41 by zzll406
[考研] 263求調(diào)劑 +6	yqdszhdap－ 2026-03-22	10/500	2026-03-26 13:11 by 公瑾逍遙
[考研] 296求調(diào)劑 +4	汪��？！ 2026-03-25	7/350	2026-03-25 16:41 by 汪��？！
[考研] 282求調(diào)劑 +3	wcq131415 2026-03-24	3/150	2026-03-25 12:16 by userper
[考研] 277分求調(diào)劑，跨調(diào)材料 +3	考研調(diào)劑lxh 2026-03-24	3/150	2026-03-24 13:52 by JourneyLucky
[考研] 一志愿國科過程所081700，274求調(diào)劑 +3	三水研0水立方 2026-03-23	3/150	2026-03-23 23:11 by MajorWen
[論文投稿] 急發(fā)核心期刊論文 +3	賢達問津 2026-03-23	5/250	2026-03-23 17:13 by 妹子不好惹