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

小木蟲(chóng)zb

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 11 (小學(xué)生)
金幣: 1560.9
紅花: 5
帖子: 345
在線: 142.1小時(shí)
蟲(chóng)號(hào): 1227944
注冊(cè): 2011-03-10
性別: GG
專業(yè): 分離過(guò)程

[求助] 求助matlab擬合方程組

遇到個(gè)方程組，不知道怎么用matlab擬合其中的參數(shù)，求助大神。k1=0.45967,n1=1/0.33824,k2=0.40308,n2=1/0.38477，擬合求解n11,n12,k12,n22,n21,k21
方程組如下
q1=(k1*c1^(1/n11 +1/n1 ))/(c1^(1/n11 )+k12*c2^(1/n12 ) )
q2=(k2*c2^(1/n22 +1/n2 ))/(c2^(1/n22 )+k21*c1^(1/n21 ) )
源數(shù)據(jù)如下
c1 q1 c2 q2
0.04556 0.04352 0.099 0.03592
0.10253 0.08542 0.12602 0.08207
0.2056 0.17075 0.27154 0.16137
0.6889 0.30201 0.87702 0.27525
1.51941 0.38387 1.81968 0.34115
2.40038 0.45855 2.72298 0.41266
3.37753 0.51955 3.82469 0.45595
6.32386 0.60045 6.74548 0.54047
9.51436 0.64661 9.89625 0.59228
12.75747 0.68528 13.15723 0.62842
16.07689 0.7131 16.59229 0.63979

回復(fù)此樓

» 猜你喜歡

調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有162人回復(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ù)
北工大宋曉艷教授團(tuán)隊(duì)招收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ù)

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

matlab數(shù)據(jù)擬合問(wèn)題求助已經(jīng)有7人回復(fù)
matlab數(shù)據(jù)擬合求助��！希望高手幫忙寫(xiě)一下matlab的程序謝謝了！！已經(jīng)有15人回復(fù)
matlab擬合求助已經(jīng)有11人回復(fù)
matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進(jìn)行數(shù)值比較已經(jīng)有4人回復(fù)
高金求助matlab解微分方程組已經(jīng)有12人回復(fù)
非線性方程組matlab擬合已經(jīng)有9人回復(fù)
急切求助用matlab曲線擬合色散方程已經(jīng)有10人回復(fù)
求助matlab擬合參數(shù)的問(wèn)題已經(jīng)有2人回復(fù)
【求助】matlab擬合二元線性方程已經(jīng)有3人回復(fù)
【求助】使用MATLAB擬合出按照公式的系數(shù) 已經(jīng)有19人回復(fù)
【求助】求助matlab數(shù)值擬合已經(jīng)有8人回復(fù)
【求助】求教多元指數(shù)方程如何解有軟件否？已經(jīng)有6人回復(fù)
【求助】matlab函數(shù)擬合，擬合結(jié)果并保存，方便調(diào)用已經(jīng)有4人回復(fù)
【求助】matlab曲面擬合表達(dá)式已經(jīng)有10人回復(fù)
【求助】如何使得matlab擬合的曲線強(qiáng)制經(jīng)過(guò)一個(gè)點(diǎn) 已經(jīng)有8人回復(fù)
【求助】指數(shù)函數(shù)方程求解已經(jīng)有5人回復(fù)
【求助】求MAtlab求解程序！��！已經(jīng)有20人回復(fù)
【求助】求教matlab解非線性方程組已經(jīng)有9人回復(fù)
【求助】使用Matlab擬合反應(yīng)動(dòng)力學(xué)方程問(wèn)題已經(jīng)有7人回復(fù)

1樓 2013-01-09 11:07:22

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

小木蟲(chóng)zb

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 11 (小學(xué)生)
金幣: 1560.9
紅花: 5
帖子: 345
在線: 142.1小時(shí)
蟲(chóng)號(hào): 1227944
注冊(cè): 2011-03-10
性別: GG
專業(yè): 分離過(guò)程

引用回帖:

2樓: Originally posted by csgt0 at 2013-01-09 11:45:59
function zbzbzbzb
data=;
xdata=;
ydata=;
x0=; %n11,n12,k12,n22,n21,k21
= lsqcurvefit(@myfun,x0,xdata,ydata)
y=myfun(x,xdata)
end

function y=myfun(x,xdata)
c1=xdata(:,1);
c2=xdata(:,2); ...

運(yùn)行不出來(lái)，出現(xiàn)錯(cuò)誤

360截圖20130110092051698.jpg

贊一下

回復(fù)此樓

4樓2013-01-10 09:21:45

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

查看全部 13 個(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
小木蟲(chóng)zb: 金幣+15, ★★★★★最佳答案 2013-01-10 09:52:54

CODE:

function zbzbzbzb

data=[ 0.04556        0.04352        0.099        0.03592   %  c1,q1,c2,q2

 0.10253        0.08542        0.12602        0.08207

 0.2056        0.17075        0.27154        0.16137

 0.6889        0.30201        0.87702        0.27525

 1.51941        0.38387        1.81968        0.34115

 2.40038        0.45855        2.72298        0.41266

 3.37753        0.51955        3.82469        0.45595

 6.32386        0.60045        6.74548        0.54047

 9.51436        0.64661        9.89625        0.59228

 12.75747        0.68528        13.15723        0.62842

 16.07689        0.7131        16.59229        0.63979];

xdata=[data(:,1),data(:,3)];

ydata=[data(:,2),data(:,4)];

x0=[0.2,0.2,0.2,9314,106,0.6]; %n11,n12,k12,n22,n21,k21

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

y=myfun(x,xdata)

end

function y=myfun(x,xdata)

c1=xdata(:,1);

c2=xdata(:,2);

n11=x(1);

n12=x(2);

k12=x(3);

n22=x(4);

n21=x(5);

k21=x(6);

k1=0.45967;

n1=1/0.33824;

k2=0.40308;

n2=1/0.38477;

q1=(k1*c1.^(1/n11 +1/n1 ))./(c1.^(1/n11 )+k12*c2.^(1/n12 ) );

q2=(k2*c2.^(1/n22 +1/n2 ))./(c2.^(1/n22 )+k21*c1.^(1/n21 ) );

y=[q1,q2];

end

贊一下(1人)

回復(fù)此樓

showmethemoney

2樓2013-01-09 11:45:59

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

dingd

鐵桿木蟲(chóng) (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲(chóng)號(hào): 291104
注冊(cè): 2006-10-28

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

★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
csgt0: 金幣+2, 謝謝 2013-01-09 15:58:11
小木蟲(chóng)zb: 金幣+5, ★★★很有幫助 2013-01-10 08:53:05

用1stOpt求解試試：

CODE:

Constant  k1=0.45967,n1=1/0.33824,k2=0.40308,n2=1/0.38477;

Variable c1,q1,c2,q2;

SharedModel;

Function q1=(k1*c1^(1/n11 +1/n1 ))/(c1^(1/n11 )+k12*c2^(1/n12 ));

         q2=(k2*c2^(1/n22 +1/n2 ))/(c2^(1/n22 )+k21*c1^(1/n21 ));

Data;

  c1            q1          c2            q2

0.04556        0.04352        0.099        0.03592

0.10253        0.08542        0.12602        0.08207

0.2056        0.17075        0.27154        0.16137

0.6889        0.30201        0.87702        0.27525

1.51941        0.38387        1.81968        0.34115

2.40038        0.45855        2.72298        0.41266

3.37753        0.51955        3.82469        0.45595

6.32386        0.60045        6.74548        0.54047

9.51436        0.64661        9.89625        0.59228

12.75747        0.68528        13.15723        0.62842

16.07689        0.7131        16.59229        0.63979

均方差(RMSE): 0.0335802018362698
殘差平方和(SSE): 0.0248078590180215
相關(guān)系數(shù)(R): 0.997283667098149
相關(guān)系數(shù)之平方(R^2): 0.994574712660733
決定系數(shù)(DC): 0.994566132239987
F統(tǒng)計(jì)(F-Statistic): 36.2326353309996

參數(shù) 最佳估算
-------------------- -------------
n11 0.241571501600743
k12 0.158571836379447
n12 0.218728458244284
n22 -0.168403298883327
k21 0.136718576047343
n21 -0.186531568585624

贊一下(2人)

回復(fù)此樓

3樓2013-01-09 13:56:50

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

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)助回帖

引用回帖:

4樓: Originally posted by 小木蟲(chóng)zb at 2013-01-10 09:21:45
運(yùn)行不出來(lái)，出現(xiàn)錯(cuò)誤

360截圖20130110092051698.jpg
...

你改程序了吧？記得用點(diǎn)乘和點(diǎn)除

贊一下

回復(fù)此樓

showmethemoney

5樓2013-01-10 09:26:59

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

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

[考研] 一志愿武漢理工，總分321，英一數(shù)二，求老師收留。 +6	nnnnnnn5 2026-03-25	6/300	2026-03-28 13:33 by 果果媽咪
[考研] 0856求調(diào)劑 +11	zhn03 2026-03-25	12/600	2026-03-28 13:32 by 唐沐兒
[考研] 299求調(diào)劑 +7	嗯嗯嗯嗯2 2026-03-27	7/350	2026-03-28 13:09 by 唐沐兒
[考研] 328求調(diào)劑 +7	嗯滴的基本都 2026-03-27	7/350	2026-03-28 04:19 by fmesaito
[考研] 265求調(diào)劑 +8	小木蟲(chóng)085600 2026-03-27	8/400	2026-03-27 22:16 by 無(wú)際的草原
[考研] 266求調(diào)劑 +11	陽(yáng)陽(yáng)哇塞 2026-03-27	12/600	2026-03-27 17:56 by yu221
[考研] 305求調(diào)劑 +5	哇盧卡庫(kù) 2026-03-26	5/250	2026-03-27 14:01 by laoshidan
[考研] 322求調(diào)劑 +4	我真的很想學(xué)習(xí) 2026-03-23	4/200	2026-03-27 13:51 by 楊楊楊紫
[考研] 333求調(diào)劑 +3	question挽風(fēng) 2026-03-23	3/150	2026-03-27 11:29 by 不吃魚(yú)的貓
[考研] 286求調(diào)劑 +4	lim0922 2026-03-26	4/200	2026-03-27 10:28 by guoweigw
[考研] 材料學(xué)碩333求調(diào)劑 +8	北道巷 2026-03-24	8/400	2026-03-27 10:18 by 我是小康
[考研] 276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷 +10	材料學(xué)257求調(diào)劑 2026-03-23	11/550	2026-03-27 10:13 by YCIT- LHL
[考研] 333求調(diào)劑 +6	wfh030413@ 2026-03-23	6/300	2026-03-26 22:45 by 學(xué)員8dgXkO
[考研] 085602 289分求調(diào)劑 +8	WWW西西弗斯 2026-03-24	8/400	2026-03-26 16:33 by 不吃魚(yú)的貓
[考研] 考研一志愿蘇州大學(xué)初始315（英一）求調(diào)劑 +3	sbdksD 2026-03-24	4/200	2026-03-25 18:16 by xcjcqu
[考研] 293求調(diào)劑 +7	加一一九 2026-03-24	7/350	2026-03-25 12:02 by userper
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +7	天天好運(yùn)來(lái)上岸?/a> 2026-03-24	7/350	2026-03-24 20:26 by peike
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4	17501029541 2026-03-23	6/300	2026-03-24 10:21 by 戴圍脖的小蚊子
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7	jhybd 2026-03-23	12/600	2026-03-24 09:03 by jhybd
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑 +7	墨墨漠 2026-03-23	8/400	2026-03-23 20:36 by Creta