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

xxz903

應(yīng)助: 33 (小學(xué)生)
金幣: 1045.6
散金: 190
紅花: 4
帖子: 253
在線: 66.7小時(shí)
蟲(chóng)號(hào): 1360336
注冊(cè): 2011-08-03
性別: GG
專業(yè): 化工系統(tǒng)工程

[求助] matlab數(shù)據(jù)擬合求助�。∠Ｍ呤謳兔懸幌耺atlab的程序謝謝了��！

自變量 m,因變量y1,y2;
需要擬合的參數(shù)t12,t21,q1,q2,r1,r2,h
其中z=10;
x1=m/(m+h);
x2=1-x1;
o1=x1*r1/(x1*r1+x2*r2);
o2=x2*r2/(x1*r1+x2*r2);
b1=x1*q1/(x1*q1+x2*q2);
b2=x2*q2/(x1*q1+x2*q2);
l1=0.5*z*(r1-q1)-(r1-1);
l2=0.5*z*(r2-q2)-(r2-1);
需要擬合的方程：
y1=ln(o1/x1)+0.5*z*q1*ln(o1/b1)+o2*(l1-r1*l2/r2)-q1*ln(b1+b2*t21)+b2*q1*(t21/(b1+b2*t21)-t12/(b2+b1*t12))-ln(r1/r2)-05*z*q1*ln(q1*r2/(q2*r1))-(l1-r1*l2/r2)+q1*ln(t21)-q1*(1-t12);

y2=ln(o2/x2)+0.5*z*q2*ln(o2/b2)+o1*(l2-r2*l1/r1)-q2*ln(b2+b1*t12)+b1*q2*(t12/(b2+b1*t12)-t21/(b1+b2*t21))-ln(r2/r1)-05*z*q2*ln(q2*r1/(q1*r2))-(l2-r2*l1/r1)+q2*ln(t12)-q2*(1-t21);
數(shù)據(jù)：
0 0 0
0.1 -0.673344553 -0.003214304
0.2 -0.783071888 -0.005447496
0.3 -0.841647189 -0.008713882
0.4 -0.869884359 -0.010990339
0.5 -0.884307686 -0.014309599
0.6 -0.889162064 -0.018687569
0.7 -0.889162064 -0.02207116
0.8 -0.88673193 -0.025484701
0.9 -0.879476759 -0.029975034
1 -0.869884359 -0.038731848
1.2 -0.850971266 -0.042636996
1.4 -0.830113036 -0.051995295
1.6 -0.807436327 -0.061533872
1.8 -0.783071888 -0.071257081
2 -0.759286983 -0.081169452
2.5 -0.701179352 -0.107972916
3 -0.650087691 -0.136153773
3.5 -0.599656837 -0.165810242
4 -0.553385238 -0.197052419
4.5 -0.512493681 -0.227286842
5 -0.47320876 -0.260553522

回復(fù)此樓

» 猜你喜歡

光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有32人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有102人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有79人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職已經(jīng)有0人回復(fù)
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)

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

1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？已經(jīng)有13人回復(fù)
matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進(jìn)行數(shù)值比較已經(jīng)有4人回復(fù)
需要用matlab擬合數(shù)據(jù)的外行人急需幫助已經(jīng)有10人回復(fù)
求助大俠幫忙，有關(guān)matlab模擬數(shù)據(jù) 已經(jīng)有29人回復(fù)
【求助】使用MATLAB擬合出按照公式的系數(shù) 已經(jīng)有19人回復(fù)
【求助】用matlab編程，需要將數(shù)據(jù)自檢分類（高手進(jìn)來(lái)指點(diǎn)下）已經(jīng)有5人回復(fù)
【求助】請(qǐng)教一個(gè)matlab程序，謝謝幫忙！！已經(jīng)有15人回復(fù)
【求助】matlab曲面擬合表達(dá)式已經(jīng)有10人回復(fù)
【求助】如何使得matlab擬合的曲線強(qiáng)制經(jīng)過(guò)一個(gè)點(diǎn) 已經(jīng)有8人回復(fù)
【求助】使用MATLAB怎么實(shí)現(xiàn)擬合力場(chǎng)參數(shù)的程序？【已解決】已經(jīng)有15人回復(fù)

1樓 2012-08-21 11:08:27

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

xxz903

金蟲(chóng) (小有名氣)

應(yīng)助: 33 (小學(xué)生)
金幣: 1045.6
散金: 190
紅花: 4
帖子: 253
在線: 66.7小時(shí)
蟲(chóng)號(hào): 1360336
注冊(cè): 2011-08-03
性別: GG
專業(yè): 化工系統(tǒng)工程

引用回帖:

2樓: Originally posted by csgt0 at 2012-08-21 11:39:11
程序大概這樣，初值得你自己調(diào)了

function xxz
data=;
xdata=data(:,1); %Im
ydata=data(:,2:3); %
x0=ones(1,7); %t12,t21,q1,q2,r1,r2,h

options= optimset('MaxFunEvals',5000,'Maxiter',1000) ...

如果我想讓這幾個(gè)參數(shù)隨機(jī)在-10到10取100次初值然后進(jìn)行計(jì)算，然后在結(jié)果中找到最優(yōu)解再進(jìn)行畫圖我改怎么修改這個(gè)程序？？謝謝你了

贊一下

回復(fù)此樓

4樓2012-08-21 14:30:14

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

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

專家經(jīng)驗(yàn): +2
計(jì)算強(qiáng)帖: 2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時(shí)
蟲(chóng)號(hào): 449046
注冊(cè): 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計(jì)算模擬

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
xxz903: 金幣+10, ★★★★★最佳答案 2012-08-21 14:05:54

程序大概這樣，初值得你自己調(diào)了

CODE:

function xxz

data=[

 0.1        -0.673344553        -0.003214304

 0.2        -0.783071888        -0.005447496

 0.3        -0.841647189        -0.008713882

 0.4        -0.869884359        -0.010990339

 0.5        -0.884307686        -0.014309599

 0.6        -0.889162064        -0.018687569

 0.7        -0.889162064        -0.02207116

 0.8        -0.88673193        -0.025484701

 0.9        -0.879476759        -0.029975034

 1        -0.869884359        -0.038731848

 1.2        -0.850971266        -0.042636996

 1.4        -0.830113036        -0.051995295

 1.6        -0.807436327        -0.061533872

 1.8        -0.783071888        -0.071257081

 2        -0.759286983        -0.081169452

 2.5        -0.701179352        -0.107972916

 3        -0.650087691        -0.136153773

 3.5        -0.599656837        -0.165810242

 4        -0.553385238        -0.197052419

 4.5        -0.512493681        -0.227286842

 5        -0.47320876        -0.260553522];

xdata=data(:,1);  %Im

ydata=data(:,2:3);  % 

x0=ones(1,7);    %t12,t21,q1,q2,r1,r2,h

options= optimset('MaxFunEvals',5000,'Maxiter',1000);

%lb=[10,20,0,0];

%ub=[10000,10000,100,1000];

lb=[];

ub=[];

[x,resnorm]=lsqcurvefit(@myfun,x0,xdata,ydata,lb,ub,options)

y=myfun(x,xdata);

plot(xdata,ydata,'g',xdata,y,'r')

end

function F=myfun(x,xdata)  %x=[]

t12=x(1);

t21=x(2);

q1=x(3);

q2=x(4);

r1=x(5);

r2=x(6);

h=x(7);

m=xdata;

z=10;

 x1=m./(m+h);

 x2=1-x1;

 o1=x1*r1./(x1*r1+x2*r2);

 o2=x2*r2./(x1*r1+x2*r2);

 b1=x1*q1./(x1*q1+x2*q2);

 b2=x2*q2./(x1*q1+x2*q2);

 l1=0.5*z*(r1-q1)-(r1-1);

 l2=0.5*z*(r2-q2)-(r2-1);

 y1=log(o1./x1)+0.5*z*q1*log(o1./b1)+o2*(l1-r1*l2/r2)-q1*log(b1+b2*t21)...

     +b2*q1.*(t21./(b1+b2*t21)-t12./(b2+b1*t12))-log(r1/r2)...

     -05*z*q1*log(q1*r2/(q2*r1))-(l1-r1*l2/r2)+q1*log(t21)-q1*(1-t12);

 y2=log(o2./x2)+0.5*z*q2*log(o2./b2)+o1*(l2-r2*l1/r1)-q2*log(b2+b1*t12)...

     +b1*q2.*(t12./(b2+b1*t12)-t21./(b1+b2*t21))-log(r2/r1)...

     -05*z*q2*log(q2*r1/(q1*r2))-(l2-r2*l1/r1)+q2*log(t12)-q2*(1-t21);

F=[y1,y2];

end

贊一下

回復(fù)此樓

showmethemoney

2樓2012-08-21 11:39:11

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

xxz903

金蟲(chóng) (小有名氣)

應(yīng)助: 33 (小學(xué)生)
金幣: 1045.6
散金: 190
紅花: 4
帖子: 253
在線: 66.7小時(shí)
蟲(chóng)號(hào): 1360336
注冊(cè): 2011-08-03
性別: GG
專業(yè): 化工系統(tǒng)工程

引用回帖:

牛��！終于能算了啊我試試啊~有問(wèn)題再聯(lián)系你啊謝謝了

贊一下

回復(fù)此樓

3樓2012-08-21 14:05:43

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

專家經(jīng)驗(yàn): +2
計(jì)算強(qiáng)帖: 2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時(shí)
蟲(chóng)號(hào): 449046
注冊(cè): 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計(jì)算模擬

引用回帖:

4樓: Originally posted by xxz903 at 2012-08-21 14:30:14
如果我想讓這幾個(gè)參數(shù)隨機(jī)在-10到10取100次初值然后進(jìn)行計(jì)算，然后在結(jié)果中找到最優(yōu)解再進(jìn)行畫圖我改怎么修改這個(gè)程序？？謝謝你了...

兩種方法
1、你可以使用matlab里的globalsearch或者multistart來(lái)進(jìn)行，但是需要先構(gòu)建問(wèn)題結(jié)構(gòu)，你可以看幫助里的例子。
2、第2種方法你可以做100次循環(huán)，每次算完后保留x和殘差，最后取殘差較小的值。
如果你的參數(shù)有范圍，可以在lb里設(shè)置下限，在ub里設(shè)置上限

贊一下

回復(fù)此樓

showmethemoney

5樓2012-08-21 15:30:13

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

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

[考研] 289求調(diào)劑 +9	懷瑾握瑜l 2026-03-20	9/450	2026-03-25 11:02 by userper
[考研] 生物技術(shù)與工程 +3	1294608413 2026-03-25	3/150	2026-03-25 10:45 by hxsm
[考研] 318求調(diào)劑 +3	plum李子 2026-03-23	3/150	2026-03-25 09:42 by 霧散后相遇lc
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +7	天天好運(yùn)來(lái)上岸?/a> 2026-03-24	7/350	2026-03-24 20:26 by peike
[考研] 0854AI CV方向招收調(diào)劑 +3	章小魚567 2026-03-23	3/150	2026-03-24 20:25 by 汪！？！
[有機(jī)交流] 有機(jī)合成求助 20+3	FENGSHUJEI 2026-03-23	5/250	2026-03-24 19:31 by 88817753
[考研] 070300化學(xué)求調(diào)劑 +9	苑豆豆 2026-03-20	9/450	2026-03-24 17:15 by licg0208
[考研] 321求調(diào)劑 +4	Ymlll 2026-03-24	4/200	2026-03-24 14:44 by sprinining
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑 +7	vv迷 2026-03-22	7/350	2026-03-23 23:44 by Txy@872106
[考研] 一志愿陜師大生物學(xué)071000，298分，求調(diào)劑 +3	SYA！ 2026-03-23	3/150	2026-03-23 19:09 by macy2011
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過(guò) +6	幸運(yùn)哩哩 2026-03-22	10/500	2026-03-22 20:10 by edmund7
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +5	困困困困坤坤 2026-03-20	6/300	2026-03-22 17:41 by hxsm
[考研] 生物學(xué)調(diào)劑 +5	Surekei 2026-03-21	5/250	2026-03-22 14:39 by tcx007
[考研] 296求調(diào)劑 +4	www_q 2026-03-20	4/200	2026-03-21 17:26 by 學(xué)員8dgXkO
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7	墨墨漠 2026-03-20	7/350	2026-03-21 16:36 by barlinike
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9	Ncdx123456 2026-03-19	9/450	2026-03-20 23:41 by lovewei0727
[考研] 材料與化工 322求調(diào)劑 +4	然11 2026-03-19	4/200	2026-03-20 22:12 by luoyongfeng
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] 材料學(xué)碩318求調(diào)劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166