版塊導航: 正在加載中...

登錄注冊

應《網絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

>論壇更新日志 (6849)
>考研 (1943)
>導師招生 (1600)
>文獻求助 (518)
>考博 (316)
>蟲友互識 (304)
>碩博家園 (165)
>休閑灌水 (140)
>博后之家 (127)
>基金申請 (115)
>招聘信息布告欄 (113)
>論文投稿 (99)
>第一性原理 (79)
>教師之家 (74)
>公派出國 (50)
>外文書籍求助 (32)

返回列表

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

小木蟲zb

木蟲 (正式寫手)

應助: 11 (小學生)
金幣: 1560.9
紅花: 5
帖子: 345
在線: 142.1小時
蟲號: 1227944
注冊: 2011-03-10
性別: GG
專業(yè): 分離過程

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

遇到個方程組，不知道怎么用matlab擬合其中的參數，求助大神。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 ) )
源數據如下
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

回復此樓

» 猜你喜歡

光學工程學碩調劑信息已經有32人回復
歡迎加入課題組已經有0人回復
物理學I論文潤色/翻譯怎么收費? 已經有221人回復
散金幣，求好運，祝面上順利！已經有79人回復
華東師范大學芯片設計徐瓏真實水平咋樣已經有1人回復
調劑已經有0人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經有0人回復
中德博士后2026項目5月30日截止- 國內博士后入站、德國兩年全職已經有0人回復
面向全球申請人的中德博士后2026項目5月30日截止- 國內博后入站德國兩年全職已經有0人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經有2人回復
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結構與器件課題組。已經有0人回復

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

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

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

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

小木蟲zb

木蟲 (正式寫手)

應助: 11 (小學生)
金幣: 1560.9
紅花: 5
帖子: 345
在線: 142.1小時
蟲號: 1227944
注冊: 2011-03-10
性別: GG
專業(yè): 分離過程

引用回帖:

7樓: Originally posted by csgt0 at 2013-01-10 09:42:08
公式里的矩陣之間的除法也要用點的...

matlab做出的結果貌似和三樓用1stOpt做的結果相差很大

贊一下

回復此樓

9樓2013-01-10 09:57:24

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

查看全部 13 個回答

csgt0

榮譽版主 (著名寫手)

彩色掛圖

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

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數 +1
小木蟲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人)

回復此樓

showmethemoney

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

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

dingd

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

計算強帖: 4
應助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時
蟲號: 291104
注冊: 2006-10-28

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數 +1
csgt0: 金幣+2, 謝謝 2013-01-09 15:58:11
小木蟲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
相關系數(R): 0.997283667098149
相關系數之平方(R^2): 0.994574712660733
決定系數(DC): 0.994566132239987
F統(tǒng)計(F-Statistic): 36.2326353309996

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

贊一下(2人)

回復此樓

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

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

小木蟲zb

木蟲 (正式寫手)

應助: 11 (小學生)
金幣: 1560.9
紅花: 5
帖子: 345
在線: 142.1小時
蟲號: 1227944
注冊: 2011-03-10
性別: GG
專業(yè): 分離過程

引用回帖:

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

運行不出來，出現(xiàn)錯誤

360截圖20130110092051698.jpg

贊一下

回復此樓

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

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

查看全部 13 個回答

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

[考研] 求調劑，一志愿:南京航空航天大學大學，080500材料科學與工程學碩，總分289分 +5	@taotao 2026-03-19	5/250	2026-03-24 21:07 by greychen00
[考研] 求調劑323材料與化工 +3	1124361 2026-03-24	3/150	2026-03-24 20:24 by peike
[考研] 0854人工智能方向招收調劑 +3	章小魚567 2026-03-24	3/150	2026-03-24 19:46 by zhouxuan..
[考研] 招08考數學 +7	laoshidan 2026-03-20	16/800	2026-03-24 17:52 by 烏拉兒山脈
[考研] 資源與環(huán)境調劑申請(333分) +7	holy J 2026-03-21	7/350	2026-03-24 17:24 by xiaohai104
[考研] 【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調劑】 +4	Higraduate 2026-03-22	7/350	2026-03-24 11:23 by 種大樹
[考研] 材料專碩英一數二306 +8	z1z2z3879 2026-03-18	8/400	2026-03-23 20:49 by baobaoye
[考研] 一志愿重慶大學085700資源與環(huán)境，總分308求調劑 +7	墨墨漠 2026-03-23	8/400	2026-03-23 20:36 by Creta
[考研] 考研化學308分求調劑 +7	你好明天你好 2026-03-23	8/400	2026-03-23 18:39 by macy2011
[論文投稿] 急發(fā)核心期刊論文 +3	賢達問津 2026-03-23	5/250	2026-03-23 17:13 by 妹子不好惹
[考研] 一志愿南京理工大學085701資源與環(huán)境302分求調劑 +5	葵梓衛(wèi)隊 2026-03-18	7/350	2026-03-23 16:26 by lingjue
[考研] 北科281學碩材料求調劑 +8	tcxiaoxx 2026-03-20	9/450	2026-03-23 12:16 by tcxiaoxx
[考研] 298求調劑 +8	上岸6666@ 2026-03-20	8/400	2026-03-23 11:02 by laoshidan
[考研] 求調劑材料學碩080500，總分289分 5+3	@taotao 2026-03-19	21/1050	2026-03-23 10:17 by 冠c哥
[考研] 276求調劑 +3	YNRYG 2026-03-21	4/200	2026-03-23 08:31 by 醉在風里
[考研] 一志愿東華大學化學070300，求調劑 +7	2117205181 2026-03-21	8/400	2026-03-22 22:55 by chixmc
[考研] 308求調劑 +3	墨墨漠 2026-03-21	3/150	2026-03-22 16:54 by i_cooler
[考研] 一志愿華中科技大學071000，求調劑 +4	沿岸有貝殼6 2026-03-21	4/200	2026-03-22 07:21 by ilovexiaobin
[考研] 299求調劑 +5	shxchem 2026-03-20	7/350	2026-03-21 17:09 by ColorlessPI
[考研] 一志愿蘇州大學材料求調劑，總分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

小木蟲zb

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

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

小木蟲zb

csgt0

【答案】應助回帖

dingd

【答案】應助回帖

小木蟲zb

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

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