版塊導(dǎo)航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

3209的煙屁

金蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 789.3
紅花: 1
帖子: 131
在線: 57.4小時
蟲號: 1003622
注冊: 2010-04-23
專業(yè): 鋼鐵冶金

[求助] 解帶參數(shù)方程，參數(shù)在變

function f=fafun(x)
syms M01 M02 M03 M04 M05 M06 T N0 R
M01=30;
M02=30;
M03=5;
M04=5;
M05=0;
M06=30;
M11=M01/60;
M12=M02/56;
M13=M03/102;
M14=M04/40.3;
M15=M05/62;
M16=M06/78;
M21=M11/(M11+M12+M13+M14+M15+M16);
M22=M12/(M11+M12+M13+M14+M15+M16);
M23=M13/(M11+M12+M13+M14+M15+M16);
M24=M14/(M11+M12+M13+M14+M15+M16);
M25=M15/(M11+M12+M13+M14+M15+M16);
M26=M16/(M11+M12+M13+M14+M15+M16);
N31=0.5*M21/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N32=0.6875*M22/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N33=0.3542*M23/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N34=0.4583*M24/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N35=0.6736*M25/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N36=0.7444*M26/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
T=1623.15;
N11=(243.2+0.031*T)/1000;
N12=(791-0.0935*T)/1000;
N13=(1024-0.177*T)/1000;
N14=(1770-0.636*T)/1000;
N15=(438-0.116*T)/1000;
N16=(1604.6-0.72*T)/1000;
N0=6.02*10^23;
R=8.314;
N21=R*T/( N0^(1/3)*(27.516*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N22=R*T/( N0^(1/3)*(20.7*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N23=R*T/(N0^(1/3)*(28.3*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N24=R*T/(N0^(1/3)*(16.1*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N25=R*T/(N0^(1/3)*(33.0*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N26=R*T/(N0^(1/3)*(31.3*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
f=[N31*exp((x(7)-N11)/N21)-x(1);
N32*exp((x(7)-N12)/N22)-x(2);
N33*exp((x(7)-N13)/N23)-x(3);
N34*exp((x(7)-N14)/N24)-x(4);
N35*exp((x(7)-N15)/N25)-x(5);
N36*exp((x(7)-N16)/N26)-x(6);
N31*exp((x(7)-N11)/N21)+N32*exp((x(7)-N12)/N22)+N33*exp((x(7)-N13)/N23)+N34*exp((x(7)-N14)/N24)+N35*exp((x(7)-N15)/N25)+N36*exp((x(7)-N16)/N26)-1];
M文件如上所示，我想讓M05+M06=30，所以就是M05=30,29，28....不想一個一個去改，有沒有類似for M05=30 if M05>1 M05=M05-1這樣的辦法直接能輸出30個結(jié)果？還有可能變?nèi)齻€量，所以想請教怎么辦？

回復(fù)此樓

» 猜你喜歡

【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有201人回復(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)有4人回復(fù)
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院已經(jīng)有10人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑本科期間發(fā)表一篇Sci 已經(jīng)有4人回復(fù)

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

微分方程組參數(shù)擬合的問題（多參數(shù)動力系統(tǒng)參數(shù)估計）已經(jīng)有11人回復(fù)
如何用得到的吸附數(shù)據(jù)進行Langmuir-Freundlich方程擬合，而得出相關(guān)參數(shù)的具體值已經(jīng)有5人回復(fù)
請問這個4階的非線性偏微分方程組在PDE中怎么解？等待高手解答。已經(jīng)有17人回復(fù)
求高版本1stopt，解決一個問題，微分方程組參數(shù)擬合已經(jīng)有3人回復(fù)
動力學(xué)方程參數(shù)估計方法已經(jīng)有14人回復(fù)
用matlab解帶參數(shù)的一元三次方程--求解已經(jīng)有3人回復(fù)
拜求，這種帶參數(shù)的方程怎么解比較好？已經(jīng)有4人回復(fù)
origin8.0擬合方程中參數(shù)怎么求已經(jīng)有5人回復(fù)
一階微分方程求解，并優(yōu)化參數(shù)值已經(jīng)有14人回復(fù)
matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進行數(shù)值比較已經(jīng)有4人回復(fù)
matlab擬合方程參數(shù)時初值的選擇已經(jīng)有15人回復(fù)
matlab解微分方程已經(jīng)有10人回復(fù)
matlab-常微分方程參數(shù)估計已經(jīng)有12人回復(fù)
求回歸方程參數(shù)估計的最優(yōu)化算法matlab代碼已經(jīng)有9人回復(fù)
求各位前輩幫忙！matlab解復(fù)數(shù)方程，出現(xiàn)圖中的警告，怎來修改呢？謝謝了啊已經(jīng)有5人回復(fù)
請教Matlab高手解多元非線性方程組中的參數(shù) 已經(jīng)有6人回復(fù)
【求助】matlab 求解微分方程中的未知參數(shù) 已經(jīng)有20人回復(fù)
求非線性方程中的參數(shù)問題已經(jīng)有13人回復(fù)
【求助】如何從Murnaghan方程的結(jié)果中得到優(yōu)化的晶格參數(shù)？已經(jīng)有11人回復(fù)
【求助】不僅控制方程里含有未知參數(shù)，邊界條件里也有已經(jīng)有5人回復(fù)

1樓 2013-04-27 10:44:04

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

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
3209的煙屁: 金幣+10, ★★★很有幫助 2013-04-28 11:28:23

1、用矩陣M06=0:30,M01=30-M06;后面都加上點
2、用循環(huán)，用數(shù)組存儲結(jié)果

贊一下

回復(fù)此樓

showmethemoney

2樓2013-04-27 17:15:11

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

3209的煙屁

金蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 789.3
紅花: 1
帖子: 131
在線: 57.4小時
蟲號: 1003622
注冊: 2010-04-23
專業(yè): 鋼鐵冶金

引用回帖:

2樓: Originally posted by csgt0 at 2013-04-27 17:15:11
1、用矩陣M06=0:30,M01=30-M06;后面都加上點
2、用循環(huán)，用數(shù)組存儲結(jié)果

大哥不要生氣，我想再問一下具體應(yīng)該怎么寫，自己改了幾下也不行，希望把改的地方寫一下好嗎，跪謝！

贊一下

回復(fù)此樓

3樓2013-04-28 11:30:48

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

3209的煙屁

金蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 789.3
紅花: 1
帖子: 131
在線: 57.4小時
蟲號: 1003622
注冊: 2010-04-23
專業(yè): 鋼鐵冶金

引用回帖:

2樓: Originally posted by csgt0 at 2013-04-27 17:15:11
1、用矩陣M06=0:30,M01=30-M06;后面都加上點
2、用循環(huán)，用數(shù)組存儲結(jié)果

M01=30;
M02=30;
M03=5;
M04=5;
for M05=0:30
if M05<=30
M05=M05+1;
end
M06=30-M05;
這是我改的，運行不了啊，問題是這些代碼是M文件中的，解這個方程還需要在命令窗口中輸入
x0=[0.4;0.2;0.007;0.005;0.16;0.14;0.4];
options=optimset('Display','off');
[x,fval,exitflag,output,Jacobian]=fsolve(@fafun,x0,options)
我需要改變M01-M05得到很多數(shù)據(jù)點，這樣一個一個來太慢了

贊一下

回復(fù)此樓

4樓2013-04-28 16:13:46

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

tracuer2011

木蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 1745.2
紅花: 1
帖子: 223
在線: 93.2小時
蟲號: 1260856
注冊: 2011-04-10
性別: GG
專業(yè): 光學(xué)

★
csgt0: 金幣+1, 謝謝 2013-05-02 09:10:58

M01=30;
M02=30;
M03=5;
M04=5;
x0=[0.4;0.2;0.007;0.005;0.16;0.14;0.4];
options=optimset('Display','off');
for M05=0:30
   M06=30-M05;
   options=optimset('Display','off');
   [x,fval,exitflag,output,Jacobian]=fsolve(@fafun,x0,options)
   Fval(M05+1)=fval; %Fval 即為所求數(shù)組
end

%%類似于這樣行不行

贊一下

回復(fù)此樓

5樓2013-04-30 00:09:22

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

3209的煙屁

金蟲 (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 789.3
紅花: 1
帖子: 131
在線: 57.4小時
蟲號: 1003622
注冊: 2010-04-23
專業(yè): 鋼鐵冶金

引用回帖:

5樓: Originally posted by tracuer2011 at 2013-04-30 00:09:22
M01=30;
M02=30;
M03=5;
M04=5;
x0=;
options=optimset('Display','off');
for M05=0:30
   M06=30-M05;
   options=optimset('Display','off');
   =fsolve(@fafun,x0,options)
   Fval(M05+1)= ...

試了一下，還是不行。M01-m06是定義在M文件中的，這樣在命令窗口給行嗎，是不是這個循環(huán)要放在M文件的代碼中�。�

贊一下

回復(fù)此樓

6樓2013-04-30 17:08:43

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

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)、準晶和納米晶
管轄: 計算模擬

function f3209
global M05 M06
options=optimset('Display','off');
x0=[0.4;0.2;0.007;0.005;0.16;0.14;0.4];
for i=1:31
M05=i-1;
M06=31-i;
[x(:,i),fval(:,i),exitflag(i),output(:,i),Jacobian(i,:,

]=fsolve(@fafun,x0,options);
end
x
fval
exitflag
output
end

function f=fafun(x)
global M05 M06
M01=30;
M02=30;
M03=5;
M04=5;
M11=M01/60;
M12=M02/56;
M13=M03/102;
M14=M04/40.3;
M15=M05/62;
M16=M06/78;
M21=M11/(M11+M12+M13+M14+M15+M16);
M22=M12/(M11+M12+M13+M14+M15+M16);
M23=M13/(M11+M12+M13+M14+M15+M16);
M24=M14/(M11+M12+M13+M14+M15+M16);
M25=M15/(M11+M12+M13+M14+M15+M16);
M26=M16/(M11+M12+M13+M14+M15+M16);
N31=0.5*M21/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N32=0.6875*M22/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N33=0.3542*M23/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N34=0.4583*M24/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N35=0.6736*M25/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
N36=0.7444*M26/(0.5*M21+0.6875*M22+0.3542*M23+0.4583*M24+0.6736*M25+0.7444*M26);
T=1623.15;
N11=(243.2+0.031*T)/1000;
N12=(791-0.0935*T)/1000;
N13=(1024-0.177*T)/1000;
N14=(1770-0.636*T)/1000;
N15=(438-0.116*T)/1000;
N16=(1604.6-0.72*T)/1000;
N0=6.02*10^23;
R=8.314;
N21=R*T/( N0^(1/3)*(27.516*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N22=R*T/( N0^(1/3)*(20.7*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N23=R*T/(N0^(1/3)*(28.3*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N24=R*T/(N0^(1/3)*(16.1*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N25=R*T/(N0^(1/3)*(33.0*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
N26=R*T/(N0^(1/3)*(31.3*[1+1*10^(-4)*(T-1773)]*10^-6)^(2/3));
f=[N31*exp((x(7)-N11)/N21)-x(1);
N32*exp((x(7)-N12)/N22)-x(2);
N33*exp((x(7)-N13)/N23)-x(3);
N34*exp((x(7)-N14)/N24)-x(4);
N35*exp((x(7)-N15)/N25)-x(5);
N36*exp((x(7)-N16)/N26)-x(6);
N31*exp((x(7)-N11)/N21)+N32*exp((x(7)-N12)/N22)+N33*exp((x(7)-N13)/N23)+N34*exp((x(7)-N14)/N24)+N35*exp((x(7)-N15)/N25)+N36*exp((x(7)-N16)/N26)-1];
end

贊一下

回復(fù)此樓

showmethemoney

7樓2013-05-02 15:05:46

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 3209的煙屁的主題更新

返回列表

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

[考研] 311求調(diào)劑一志愿合肥工業(yè)大學(xué) +5	秋二十二 2026-03-30	5/250	2026-03-31 06:55 by 無際的草原
[考研] 一志愿鄭大材料工程290求調(diào)劑 +12	Youth_ 2026-03-30	12/600	2026-03-31 03:34 by 蒙奇奇521
[考研] 一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑 +25	稻妻小編 2026-03-26	25/1250	2026-03-30 20:11 by 滴滴上岸呀
[考研] 327求調(diào)劑 +5	小卡不卡. 2026-03-29	5/250	2026-03-30 19:30 by Wang200018
[考研] 332求調(diào)劑 +6	@MZB382400 2026-03-28	6/300	2026-03-30 16:57 by 無際的草原
[考研] 311求調(diào)劑 +10	lin0039 2026-03-26	10/500	2026-03-30 10:26 by herarysara
[考研] 279求調(diào)劑 +4	蝶舞輕繞 2026-03-29	4/200	2026-03-29 09:45 by laoshidan
[考研] 調(diào)劑求院校招收 +6	鶴鯨鴿 2026-03-28	6/300	2026-03-29 08:15 by fmesaito
[考研] 299求調(diào)劑 +7	嗯嗯嗯嗯2 2026-03-27	7/350	2026-03-28 13:09 by 唐沐兒
[考研] 283求調(diào)劑 +7	A child 2026-03-28	7/350	2026-03-28 12:05 by zllcz
[考研] 286求調(diào)劑 +4	丟掉懶惰 2026-03-27	7/350	2026-03-28 08:07 by baoball
[考研] 279 分求調(diào)劑 +4	睡個好覺_16 2026-03-24	4/200	2026-03-27 15:05 by 醉在風里
[考研] 考研調(diào)劑 +10	呼呼？~+123456 2026-03-24	10/500	2026-03-27 11:46 by wangjy2002
[考研] 求調(diào)劑 +3	劉柯@ 2026-03-24	4/200	2026-03-27 11:28 by shangxh
[考研] 一志愿鄭大085600，310分求調(diào)劑 +5	李瀟可 2026-03-26	5/250	2026-03-27 11:14 by 不吃魚的貓
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6	kongweilin 2026-03-26	8/400	2026-03-27 10:18 by FF_16
[考研] 調(diào)劑求收留 +7	果然有我 2026-03-26	7/350	2026-03-27 00:26 by wxiongid
[考研] 290分調(diào)劑求助 +3	吉祥止止陳 2026-03-25	3/150	2026-03-25 19:58 by barlinike
[考研] 【2026考研調(diào)劑】制藥工程 284分求相關(guān)專業(yè)調(diào)劑名額 +4	袁奐奐 2026-03-25	8/400	2026-03-25 14:32 by lbsjt
[考研] 一志愿武理085500機械專業(yè)總分300求調(diào)劑 +3	an10101 2026-03-24	7/350	2026-03-25 00:00 by 山鬼0-

24小時熱門版塊排行榜

[求助] 解帶參數(shù)方程，參數(shù)在變

» 猜你喜歡

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

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

[求助] 解帶參數(shù)方程，參數(shù)在變

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