當(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ù)擬合求助！！希望高手幫忙寫(xiě)一下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ù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有187人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有68人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(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)工程

引用回帖:

8樓: Originally posted by csgt0 at 2012-08-21 17:34:26
中間幾句改成這樣可以試試。
沒(méi)有給你在-10到10上隨機(jī)，只在0-1上隨機(jī)初值了。每個(gè)初值都要保證你的ln里面的恒正，不然就錯(cuò)了。解決這個(gè)問(wèn)題得轉(zhuǎn)化為其他的函數(shù)用不等式約束的方法，那將十分麻煩。

xdata=dat ...

版主為什么算出來(lái)的是復(fù)數(shù)？？我自己修改的隨機(jī)取初值的算出的也是有復(fù)數(shù)的這個(gè)問(wèn)題應(yīng)該怎么解決？

贊一下

回復(fù)此樓

9樓2012-08-21 19:33:48

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

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

csgt0

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

彩色掛圖

專家經(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)工程

引用回帖:

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) ...

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

贊一下

回復(fù)此樓

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

已閱回復(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)工程

引用回帖:

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

贊一下

回復(fù)此樓

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

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

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

[考研] 333求調(diào)劑 +5	87639 2026-03-21	5/250	2026-03-21 17:54 by JourneyLucky
[考研] 280求調(diào)劑 +10	咕嚕曉曉 2026-03-18	11/550	2026-03-21 17:40 by 學(xué)員8dgXkO
[考研] 313求調(diào)劑 +4	肆叁貳壹22 2026-03-19	4/200	2026-03-21 17:33 by ColorlessPI
[考研] 材料 271求調(diào)劑 +5	展信悅_ 2026-03-21	5/250	2026-03-21 17:29 by 學(xué)員8dgXkO
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +8	楓橋ZL 2026-03-18	10/500	2026-03-21 15:29 by Shawn0911
[考研] 08工科 320總分求調(diào)劑 +6	梨花珞晚風(fēng) 2026-03-17	6/300	2026-03-21 03:40 by JourneyLucky
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +12	yangfz 2026-03-17	12/600	2026-03-21 03:30 by JourneyLucky
[考研] 274求調(diào)劑 +10	S.H1 2026-03-18	10/500	2026-03-20 23:51 by JourneyLucky
[考研] 304求調(diào)劑 +7	司空. 2026-03-18	7/350	2026-03-20 23:08 by JourneyLucky
[考研] 考研調(diào)劑求學(xué)校推薦 +3	伯樂(lè)29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4	晨昏線與星海 2026-03-19	4/200	2026-03-20 22:15 by JourneyLucky
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4	@taotao 2026-03-19	4/200	2026-03-20 22:14 by JourneyLucky
[考研] 工科材料085601 279求調(diào)劑 +7	困于星晨 2026-03-17	9/450	2026-03-20 17:38 by 無(wú)懈可擊111
[考研] 招收調(diào)劑碩士 +4	lidianxing 2026-03-19	12/600	2026-03-20 12:25 by lidianxing
[考研] 材料與化工求調(diào)劑 +7	為學(xué)666 2026-03-16	7/350	2026-03-19 14:48 by 盡舜堯1
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國(guó)級(jí)證書(shū)40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 283求調(diào)劑 +3	聽(tīng)風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客