| 8 | 1/1 | 返回列表 |
| 查看: 3797 | 回復(fù): 7 | |||||
[求助]
matlab解一元三次方程,得到的都是復(fù)數(shù)根。
|
|
我想問一下,matlab解一元三次方程,得到的都是復(fù)數(shù)根。但是虛部前面的系數(shù)很小很小了,這個(gè)可以把虛部忽略嗎,直接取實(shí)部作為根。解如下: v = 28.786359708932167 - 0.000000000000449i -0.056673985902386 - 0.000000000229294i -0.000315349578379 + 0.000000000229743i |

主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |
|
對(duì)于你昨天的那個(gè)方程,是典型的非線性方程,我修改了程序,用fsolve函數(shù)來解,程序在下面。先新建一個(gè)m文件,再把以下全部程序直接復(fù)制進(jìn)m文件中,運(yùn)行即可,不要在主窗口中運(yùn)行,用m文件來寫程序,比直接在主窗口中操作方便,也算是一種好習(xí)慣吧。 由計(jì)算結(jié)果可知,解更符合實(shí)際,而且沒有出現(xiàn)虛數(shù)。 function solution_nonl clear all;clc global aa bb tt pp %-------------根據(jù)計(jì)算結(jié)果,修改迭代初值------------------------------------ xx0=[1]; %迭代初值,如果計(jì)算結(jié)果不符合實(shí)際,可修改,一般由pV=nRT給出初始的V,這里我直接取1了 t =[350.1500 348.9500 347.1500 345.4500 342.3500 340.7500 338.9500]; p =[101.3300 101.3300 101.3300 101.3300 101.3300 101.3300 101.3300 101.3300]; a =[0.9098 0.9365 0.9763 1.0223 1.0549 1.0968 1.1340 1.1773]; b =[0.0573 0.0565 0.0551 0.0535 0.0523 0.0507 0.0491 0.0470]; v=[]; for n=1:7 aa=a(n); bb=b(n); tt=t(n); pp=p(n); x(n)=fsolve(@nonlin,xx0); end disp(x) function f= nonlin(x) global aa bb tt pp f=8.314*tt/(x-bb)-aa/(x*(x+bb))-pp; 計(jì)算結(jié)果: 28.7864 28.6871 28.5380 28.3969 28.1413 28.0084 27.8591 |

主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |

木蟲 (知名作家)
迷糊領(lǐng)域榮譽(yù)領(lǐng)主


木蟲 (知名作家)
|
銅蟲 (小有名氣)
新蟲 (初入文壇)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 266求調(diào)劑 +3 | 陽陽哇塞 2026-03-14 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +7 | yangfz 2026-03-17 | 7/350 |
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +3 | 脫穎而出 2026-03-16 | 6/300 |
|
|
[考研] 302求調(diào)劑 +9 | 負(fù)心者當(dāng)誅 2026-03-11 | 9/450 |
|
|
[考研] 085600材料與化工 +4 | 安全上岸! 2026-03-16 | 4/200 |
|
|
[考研] 211本,11408一志愿中科院277分,曾在中科院自動(dòng)化所實(shí)習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
[考研] 一志愿985,本科211,0817化學(xué)工程與技術(shù)319求調(diào)劑 +5 | Liwangman 2026-03-15 | 5/250 |
|
|
[考研] 311求調(diào)劑 +5 | 26研0 2026-03-15 | 5/250 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧, 2026-03-12 | 4/200 |
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
[考研] 290求調(diào)劑 +9 | ADT 2026-03-11 | 9/450 |
|
|
[考研] 一志愿西南交大,材料專碩317求調(diào)劑 +5 | lx8568 2026-03-11 | 5/250 |
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
[考研] 277求調(diào)劑 +4 | anchor17 2026-03-12 | 4/200 |
|