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

xxz903

金蟲 (小有名氣)

應(yīng)助: 33 (小學(xué)生)
金幣: 1045.6
散金: 190
紅花: 4
帖子: 253
在線: 66.7小時(shí)
蟲號(hào): 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ù)此樓

» 猜你喜歡

光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有32人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有96人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有78人回復(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日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請人的中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(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)來指點(diǎn)下）已經(jīng)有5人回復(fù)
【求助】請教一個(gè)matlab程序，謝謝幫忙��！已經(jīng)有15人回復(fù)
【求助】matlab曲面擬合表達(dá)式已經(jīng)有10人回復(fù)
【求助】如何使得matlab擬合的曲線強(qiáng)制經(jīng)過一個(gè)點(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的回帖

xxz903

金蟲 (小有名氣)

應(yīng)助: 33 (小學(xué)生)
金幣: 1045.6
散金: 190
紅花: 4
帖子: 253
在線: 66.7小時(shí)
蟲號(hào): 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的回帖

查看全部 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í)
蟲號(hào): 449046
注冊: 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

金蟲 (小有名氣)

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

引用回帖:

如果我想讓這幾個(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的回帖

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í)
蟲號(hào): 449046
注冊: 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來進(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的回帖

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

[考研] 材料與化工328分調(diào)劑 +3	。，。，。，。i 2026-03-23	3/150	2026-03-24 10:34 by Pegasus0429
[考研] 一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271 +5	kotoko_ik 2026-03-23	6/300	2026-03-24 10:29 by 學(xué)術(shù)搬磚er
[考研] 【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】 +3	Higraduate 2026-03-22	6/300	2026-03-24 09:55 by AZMK
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑 +7	vv迷 2026-03-22	7/350	2026-03-23 23:44 by Txy@872106
[考研] 284求調(diào)劑 +3	yanzhixue111 2026-03-23	6/300	2026-03-23 22:58 by pswait
[考研] 材料專碩英一數(shù)二306 +8	z1z2z3879 2026-03-18	8/400	2026-03-23 20:49 by baobaoye
[考研] 303求調(diào)劑 +4	元夕元 2026-03-20	4/200	2026-03-23 19:00 by macy2011
[考研] 336化工調(diào)劑 +4	王大坦1 2026-03-23	5/250	2026-03-23 18:32 by allen-yin
[考研] 291 求調(diào)劑 +4	化工2026屆畢業(yè)?/a> 2026-03-21	5/250	2026-03-23 16:46 by 化工2026屆畢業(yè)?/a>
[考研] 08工學(xué)調(diào)劑 +7	用戶573181 2026-03-20	11/550	2026-03-23 15:47 by 我愛學(xué)習(xí)學(xué)習(xí)使?/a>
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4	allen-yin 2026-03-23	6/300	2026-03-23 15:04 by 汪�。�！
[考研] 293求調(diào)劑 +3	濤濤Wjt 2026-03-22	5/250	2026-03-22 22:21 by jiangpengfei
[考研] 289求調(diào)劑 +7	懷瑾握瑜l 2026-03-20	7/350	2026-03-22 15:57 by ColorlessPI
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5	vv迷 2026-03-21	8/400	2026-03-22 14:20 by ColorlessPI
[考研] 初試 317 +7	半拉月丙 2026-03-20	7/350	2026-03-21 22:26 by peike
[考研] 297求調(diào)劑 +3	喜歡還是不甘心 2026-03-20	3/150	2026-03-21 18:33 by 學(xué)員8dgXkO
[考研] 303求調(diào)劑 +5	睿08 2026-03-17	7/350	2026-03-21 03:11 by JourneyLucky
[考研] 330求調(diào)劑 +4	小材化本科 2026-03-18	4/200	2026-03-20 23:13 by JourneyLucky
[考研] 320求調(diào)劑0856 +3	不想起名字112 2026-03-19	3/150	2026-03-19 22:53 by 學(xué)員8dgXkO
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438