| 7 | 1/1 | 返回列表 |
| 查看: 1058 | 回復(fù): 6 | ||
xwl196銀蟲 (小有名氣)
|
[求助]
matlab求解含中間變量的非線性方程組 已有1人參與
|
|
方程組如附件所示,其中: T=328.2; ρmc=1.0254; Eii=0.432*195.2; C=1.23; Eis=-1250; ρads為中間變量 求大神指點:1.如何求解想x和y;2.或者如何從方程(1)中把ρads表示成關(guān)于x的表達式。 1.jpg |
鐵桿木蟲 (職業(yè)作家)
銀蟲 (小有名氣)
|
ρb=x=[[0.0116 0.0240 0.0502 0.0804 0.1170 0.1608 0.2152 0.2972 0.4123 0.5257 0.6100]; 所以問題轉(zhuǎn)化為在方程(1)中用matlab的solve求解ρads,然后將ρads帶入方程(2)中求解nex就好了! matlab程序: >> T=328.2; rhomc=1.0254; Eii=0.432*195.2; C=1.23; Eis=-1250; x=[0.0116 0.0240 0.0502 0.0804 0.1170 0.1608 0.2152 0.2972 0.4123 0.5257 0.6100]; >> syms rhoads >> [rhoads]=solve(log(rhoads*(rhomc-x)/(x*(rhomc-rhoads)))+(7*rhoads-8*x)/rhomc*Eii/T+Eis/T==0) 運行后出現(xiàn)的提示是: ??? Error using ==> mupadmex Error in MuPAD command: Array sizes must match. Error in ==> sym.plus at 14 X = mupadmex('symobj::zip',A.s,B.s,'_plus'); 沒找到問題出現(xiàn)在哪? |
鐵桿木蟲 (職業(yè)作家)
|
Solve命令是求符號解,不一定有的。可以是數(shù)值解。 下面是用1stOpt求解,很簡單: Ñ-»·3£êy x rhoads nex 0.0116 0.257851015773962 0.590771892728639 0.024 0.378121213925542 0.849559378054255 0.0502 0.520370658449439 1.1279693970993 0.0804 0.618855300376744 1.29178860827656 0.117 0.699976313181814 1.3985973575456 0.1608 0.769244989810605 1.45969833716997 0.2152 0.831497991791605 1.47853818978677 0.2972 0.896153280504329 1.43692712116311 0.4123 0.952359132344292 1.2956363034591 0.5257 0.984696573619486 1.10116205490924 0.61 0.999729281534294 0.934985403329788 |
銀蟲 (小有名氣)
|
你用的1stopt是哪個版本的,我用1stopt1.5運行的結(jié)果如下: "model nex" ====== 結(jié)果 ====== 迭代數(shù): 16 計算用時(時:分:秒:毫秒): 00:00:00:50 計算中止原因: 達到收斂判定標準 優(yōu)化算法: 麥夸特法(Levenberg-Marquardt) + 通用全局優(yōu)化法 函數(shù)表達式 ln(rhoads*(1.0254-x)/(x*(1.0254-rhoads)))+(7*rhoads-8*x)/1.0254*84.3264/328.2+(-1250)/328.2 目標函數(shù)值: 0 rhoads: 0.485072231445505 x: 0.0421283782771957 ====== 計算結(jié)束 ====== 沒有得到你運行的結(jié)果,是不是我用的1stopt版本太低的原因? |
鐵桿木蟲 (職業(yè)作家)
銀蟲 (小有名氣)
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +3 | sbdksD 2026-03-19 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +6 | sllhht 2026-03-19 | 7/350 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分,一志愿四川大學(xué),誠心求調(diào)劑 +5 | 吃吃吃才有意義 2026-03-19 | 5/250 |
|
|
[考研] 一志愿福大288有機化學(xué),求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +3 | Ncdx123456 2026-03-19 | 3/150 |
|
|
[考研] 281求調(diào)劑(0805) +9 | 煙汐憶海 2026-03-16 | 19/950 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 0817調(diào)劑 +3 | 沒有答案_ 2026-03-14 | 3/150 |
|
|
[考研] 材料專碩英一數(shù)二306 +5 | z1z2z3879 2026-03-18 | 5/250 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識可樂 2026-03-14 | 5/250 |
|