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

主管區(qū)長(zhǎng) (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |
|
對(duì)于你昨天的那個(gè)方程,是典型的非線性方程,我修改了程序,用fsolve函數(shù)來(lái)解,程序在下面。先新建一個(gè)m文件,再把以下全部程序直接復(fù)制進(jìn)m文件中,運(yùn)行即可,不要在主窗口中運(yùn)行,用m文件來(lái)寫程序,比直接在主窗口中操作方便,也算是一種好習(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ū)長(zhǎng) (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |

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


木蟲 (知名作家)
|
銅蟲 (小有名氣)
新蟲 (初入文壇)
| 8 | 1/1 | 返回列表 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
|
百度網(wǎng)盤 |
360云盤 |
千易網(wǎng)盤 |
華為網(wǎng)盤
在新窗口頁(yè)面中打開自己喜歡的網(wǎng)盤網(wǎng)站,將文件上傳后,然后將下載鏈接復(fù)制到帖子內(nèi)容中就可以了。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +7 | yangfz 2026-03-17 | 7/350 |
|
|---|---|---|---|---|
|
[基金申請(qǐng)] 被我言中:新模板不強(qiáng)調(diào)格式了,假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 274求調(diào)劑0856材料化工 +13 | z2839474511 2026-03-11 | 14/700 |
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +3 | 困于星晨 2026-03-17 | 3/150 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 一志愿985,本科211,0817化學(xué)工程與技術(shù)319求調(diào)劑 +5 | Liwangman 2026-03-15 | 5/250 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研]
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 四川大學(xué)085601材料工程專碩 初試294求調(diào)劑 +4 | 祝我們好在冬天 2026-03-11 | 4/200 |
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
[考研] 工科278分求調(diào)劑 +5 | 周慢熱啊 2026-03-12 | 7/350 |
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過(guò) 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-12 | 3/150 |
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|