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

xxz903

金蟲 (小有名氣)

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

[求助] matlab數(shù)據(jù)擬合求助��！希望高手幫忙寫一下matlab的程序謝謝了�。�

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

» 猜你喜歡

計算與實(shí)驗(yàn)科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有4人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有110人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有28人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有26人回復(fù)
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
UJN物理學(xué)專業(yè)調(diào)劑已經(jīng)有3人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有0人回復(fù)

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

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)來指點(diǎn)下）已經(jīng)有5人回復(fù)
【求助】請教一個matlab程序，謝謝幫忙�。� 已經(jīng)有15人回復(fù)
【求助】matlab曲面擬合表達(dá)式已經(jīng)有10人回復(fù)
【求助】如何使得matlab擬合的曲線強(qiáng)制經(jīng)過一個點(diǎn) 已經(jīng)有8人回復(fù)
【求助】使用MATLAB怎么實(shí)現(xiàn)擬合力場參數(shù)的程序？【已解決】已經(jīng)有15人回復(fù)

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

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

csgt0

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

彩色掛圖

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

引用回帖:

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

兩種方法
1、你可以使用matlab里的globalsearch或者multistart來進(jìn)行，但是需要先構(gòu)建問題結(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的回帖

查看全部 16 個回答

csgt0

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

彩色掛圖

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

【答案】應(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

金蟲 (小有名氣)

應(yīng)助: 33 (小學(xué)生)
金幣: 1045.6
散金: 190
紅花: 4
帖子: 253
在線: 66.7小時
蟲號: 1360336
注冊: 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) ...

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

贊一下

回復(fù)此樓

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

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

xxz903

金蟲 (小有名氣)

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

引用回帖:

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

贊一下

回復(fù)此樓

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

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

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

[考研] 環(huán)境工程297分求調(diào)劑一志愿杭高院 +10	GENJIOW 2026-03-31	11/550	2026-04-01 21:38 by GENJIOW
[考研] 292求調(diào)劑 +17	木蟲er12138 2026-04-01	17/850	2026-04-01 21:37 by 七度不信任
[考研] 339求調(diào)劑，想調(diào)回江蘇 +7	烤麥芽 2026-03-27	10/500	2026-04-01 21:35 by 495374996
[考研] 材料科學(xué)與工程調(diào)劑 +16	深V宿舍吧 2026-03-30	17/850	2026-04-01 21:00 by 賴春艷
[考研] 362求調(diào)劑 +13	西南交材料專碩3 2026-03-31	13/650	2026-04-01 17:38 by JYD2011
[考研] 314求調(diào)劑 +4	溪云珂 2026-03-26	4/200	2026-04-01 17:00 by oooqiao
[考研] 318求調(diào)劑 +8	七憶77 2026-04-01	8/400	2026-04-01 10:37 by Jaylen.
[考研] 358求調(diào)劑 +3	王向陽花 2026-03-31	3/150	2026-04-01 09:56 by zzchen2000
[考研] 材料工程085601數(shù)二英一335求調(diào)劑 +5	雙馬尾痞老板2 2026-03-31	5/250	2026-03-31 19:07 by Wang200018
[考研] 物理學(xué)調(diào)劑 +4	小羊36 2026-03-30	4/200	2026-03-31 16:16 by lishahe
[考研] 0703化學(xué) +20	妮妮ninicgb 2026-03-27	20/1000	2026-03-31 13:33 by 無際的草原
[考研] 一志愿西電085401數(shù)一英一299求調(diào)劑六級521 +4	愛吃大鴨梨 2026-03-31	4/200	2026-03-31 11:51 by 搏擊518
[考博] 材料專業(yè)申博 +5	杜雨婷dyt 2026-03-29	5/250	2026-03-31 11:19 by oooqiao
[考研] 085701求調(diào)劑初試286分 +5	secret0328 2026-03-28	5/250	2026-03-30 12:54 by fangnagu
[考研] 本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑 +3	AZMK 2026-03-27	5/250	2026-03-28 16:19 by xxxsssccc
[考研] 材料求調(diào)劑一志愿哈工大324 +7	閆旭東 2026-03-28	9/450	2026-03-28 08:51 by Xu de nuo
[考研] 352分化工與材料 +5	海納百川Ly 2026-03-27	5/250	2026-03-28 03:39 by fmesaito
[考研] 一志愿上海理工能源動力（085800）310分求調(diào)劑 +3	zhangmingc 2026-03-27	4/200	2026-03-27 19:01 by 給你你注意休息
[考研] 266求調(diào)劑 +11	陽陽哇塞 2026-03-27	12/600	2026-03-27 17:56 by yu221
[考研] 中國科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5	YangTyu1 2026-03-26	5/250	2026-03-26 18:27 by 貓咪貓咪呀