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

» 猜你喜歡

中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有222人回復(fù)
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團隊招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！已經(jīng)有1人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人已經(jīng)有0人回復(fù)
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人已經(jīng)有0人回復(fù)
計算與實驗科學(xué)工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學(xué)換能器分會S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有1人回復(fù)

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

1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？已經(jīng)有13人回復(fù)
matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進行數(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īng)有5人回復(fù)
【求助】請教一個matlab程序，謝謝幫忙！！已經(jīng)有15人回復(fù)
【求助】matlab曲面擬合表達式已經(jīng)有10人回復(fù)
【求助】如何使得matlab擬合的曲線強制經(jīng)過一個點已經(jīng)有8人回復(fù)
【求助】使用MATLAB怎么實現(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小時
蟲號: 1360336
注冊: 2011-08-03
性別: GG
專業(yè): 化工系統(tǒng)工程

引用回帖:

5樓: Originally posted by csgt0 at 2012-08-21 15:30:13
兩種方法
1、你可以使用matlab里的globalsearch或者multistart來進行，但是需要先構(gòu)建問題結(jié)構(gòu)，你可以看幫助里的例子。
2、第2種方法你可以做100次循環(huán)，每次算完后保留x和殘差，最后取殘差較小的值。
如果你的 ...

版主~~再幫我改一下吧我改了好久了也沒畫出最優(yōu)解的圖形~~感謝你了啊

贊一下

回復(fù)此樓

7樓2012-08-21 16:34:50

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

查看全部 16 個回答

csgt0

榮譽版主 (著名寫手)

彩色掛圖

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

【答案】應(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ù)隨機在-10到10取100次初值然后進行計算，然后在結(jié)果中找到最優(yōu)解再進行畫圖我改怎么修改這個程序？？謝謝你了

贊一下

回復(fù)此樓

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

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

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

[考研] 085601材料工程找調(diào)劑 +14	oatmealR 2026-03-29	15/750	2026-03-29 21:48 by plmuchong
[考研] 一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑 +17	稻妻小編 2026-03-26	17/850	2026-03-29 19:07 by 學(xué)員8dgXkO
[考研] 085600 材料與化工 329分求調(diào)劑 +14	Mr. Z 2026-03-25	14/700	2026-03-29 17:27 by 學(xué)員8dgXkO
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5	hanamiko 2026-03-29	5/250	2026-03-29 16:46 by 學(xué)員8dgXkO
[考研] 22408 359分調(diào)劑 +4	Qshers 2026-03-27	6/300	2026-03-29 14:16 by Qshers
[考研] 327求調(diào)劑 +6	汲亦昊 2026-03-29	6/300	2026-03-29 13:40 by peike
[考研] 343求調(diào)劑 +6	愛羈絆 2026-03-29	6/300	2026-03-29 12:00 by 無際的草原
[考研] 332求調(diào)劑 +8	蕉蕉123 2026-03-28	8/400	2026-03-29 10:46 by 周梓丹
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +9	天天好運來上岸?/a> 2026-03-24	10/500	2026-03-28 22:17 by chemzp
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +5	0857求調(diào)劑. 2026-03-24	6/300	2026-03-28 22:04 by zhq0425
[考研] 316求調(diào)劑 +7	江辭666 2026-03-26	7/350	2026-03-28 21:28 by sanrepian
[考研] 320分，材料與化工專業(yè)，求調(diào)劑 +9	一定上岸aaa 2026-03-27	13/650	2026-03-28 15:00 by 神馬都不懂
[考研] 材料277求調(diào)劑 +7	min3 2026-03-24	7/350	2026-03-28 11:39 by xuxiang
[考研] 311求調(diào)劑 +3	希望上岸阿小楊 2026-03-23	3/150	2026-03-28 07:57 by 熱情沙漠
[考研] 330一志愿中國海洋大學(xué) 化學(xué)工程 085602 有讀博意愿求調(diào)劑 +3	wywy.. 2026-03-27	4/200	2026-03-28 03:32 by fmesaito
[考研] 考研調(diào)劑 +4	Sanmu-124 2026-03-26	4/200	2026-03-27 17:49 by kiokin
[考研] 287求調(diào)劑 +10	land xuxu 2026-03-26	10/500	2026-03-27 15:33 by 帕爾馬拉特
[考研] 290分調(diào)劑求助 +3	吉祥止止陳 2026-03-25	3/150	2026-03-25 19:58 by barlinike
[考研] 0854AI CV方向招收調(diào)劑 +4	章小魚567 2026-03-23	4/200	2026-03-25 17:04 by CoderLoser
[考研] 293求調(diào)劑 +7	加一一九 2026-03-24	7/350	2026-03-25 12:02 by userper