 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1934 | 回復: 6 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
[求助]
解帶參數(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)齻�量，所以想請教怎么辦？ |
回復此樓» 猜你喜歡
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復
-
北工大宋曉艷教授團隊招收2026年9月入學博士【學博和專博名額均有】！
已經(jīng)有1人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有51人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復
-
計算與實驗科學工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學換能器分會S8
已經(jīng)有0人回復
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有4人回復
-
歡迎調(diào)劑到贛南師范大學智能制造與未來能源學院
已經(jīng)有23人回復
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有4人回復
-
歡迎調(diào)劑到贛南師范大學智能制造與未來能源學院
已經(jīng)有6人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 微分方程組參數(shù)擬合的問題（多參數(shù)動力系統(tǒng)參數(shù)估計）
已經(jīng)有11人回復
- 如何用得到的吸附數(shù)據(jù)進行Langmuir-Freundlich方程擬合，而得出相關(guān)參數(shù)的具體值
已經(jīng)有5人回復
- 請問這個4階的非線性偏微分方程組在PDE中怎么解？等待高手解答。
已經(jīng)有17人回復
- 求高版本1stopt，解決一個問題，微分方程組參數(shù)擬合
已經(jīng)有3人回復
- 動力學方程參數(shù)估計方法
已經(jīng)有14人回復
- 用matlab解帶參數(shù)的一元三次方程--求解
已經(jīng)有3人回復
- 拜求，這種帶參數(shù)的方程怎么解比較好？
已經(jīng)有4人回復
- origin8.0擬合方程中參數(shù)怎么求
已經(jīng)有5人回復
- 一階微分方程求解，并優(yōu)化參數(shù)值
已經(jīng)有14人回復
- matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進行數(shù)值比較
已經(jīng)有4人回復
- matlab擬合方程參數(shù)時初值的選擇
已經(jīng)有15人回復
- matlab解微分方程
已經(jīng)有10人回復
- matlab-常微分方程參數(shù)估計
已經(jīng)有12人回復
- 求回歸方程參數(shù)估計的最優(yōu)化算法matlab代碼
已經(jīng)有9人回復
- 求各位前輩幫忙！matlab解復數(shù)方程，出現(xiàn)圖中的警告，怎來修改呢？謝謝了啊
已經(jīng)有5人回復
- 請教Matlab高手解多元非線性方程組中的參數(shù)
已經(jīng)有6人回復
- 【求助】matlab 求解微分方程中的未知參數(shù)
已經(jīng)有20人回復
- 求非線性方程中的參數(shù)問題
已經(jīng)有13人回復
- 【求助】如何從Murnaghan方程的結(jié)果中得到優(yōu)化的晶格參數(shù)？
已經(jīng)有11人回復
- 【求助】不僅控制方程里含有未知參數(shù)，邊界條件里也有
已經(jīng)有5人回復
3樓2013-04-28 11:30:48
csgt0
榮譽版主 (著名寫手)
彩色掛圖
-
專家經(jīng)驗: +2 - 計算強帖: 2
- 應助: 367 (碩士)
- 貴賓: 0.457
- 金幣: 7807
- 散金: 8897
- 紅花: 20
- 沙發(fā): 7
- 帖子: 2809
- 在線: 515.8小時
- 蟲號: 449046
- 注冊: 2007-11-02
- 專業(yè): 金屬非晶態(tài)、準晶和納米晶
- 管轄: 計算模擬
2樓2013-04-27 17:15:11
|
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ù)點，這樣一個一個來太慢了 |
4樓2013-04-28 16:13:46
tracuer2011
木蟲 (小有名氣)
- 應助: 2 (幼兒園)
- 金幣: 1745.2
- 紅花: 1
- 帖子: 223
- 在線: 93.2小時
- 蟲號: 1260856
- 注冊: 2011-04-10
- 性別: GG
- 專業(yè): 光學
★
csgt0: 金幣+1, 謝謝 2013-05-02 09:10:58
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 %%類似于這樣行不行 |
5樓2013-04-30 00:09:22
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑：085600材料與化工，考材科基，總分319 +8 | 678lucky 2026-03-31 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-29 | 3/150 |
|
|
[考研] 291求調(diào)劑 +12 | Y-cap 2026-03-29 | 15/750 |
|
|
|
[考研] 材料與化工調(diào)劑一志愿大連海事085600，349 +6 | 吃的不少 2026-03-30 | 6/300 |
|
|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|
|
[考研] 材料工程專碩求調(diào)劑 +8 | hyl3153942 2026-03-29 | 8/400 |
|
|
|
[考研] 295材料工程專碩求調(diào)劑 +10 | 1428151015 2026-03-27 | 10/500 |
|
|
|
[考研] 348求調(diào)劑 +6 | 小懶蟲不懶了 2026-03-28 | 6/300 |
|
|
|
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷 +8 | cccchenso 2026-03-29 | 8/400 |
|
|
|
[考研] 一志愿北京理工大學本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
|
[考研] 2026年華南師范大學歡迎化學，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學子加入！ +3 | llss0711 2026-03-28 | 6/300 |
|
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 283求調(diào)劑 +7 | A child 2026-03-28 | 7/350 |
|
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
|
[考研] 070300化學求調(diào)劑 +4 | 起個名咋這么難 2026-03-27 | 4/200 |
|
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
|
[考研]
 085602化學工程求調(diào)劑。
+4
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
|
[考研] 機械學碩總分317求調(diào)劑�。。�！ +4 | Acaciad 2026-03-25 | 4/200 |
|