| 5 | 1/1 | 返回列表 |
| 查看: 3829 | 回復(fù): 7 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
[求助]
matlab解一元三次方程,得到的都是復(fù)數(shù)根。
|
||||
|
我想問一下,matlab解一元三次方程,得到的都是復(fù)數(shù)根。但是虛部前面的系數(shù)很小很小了,這個可以把虛部忽略嗎,直接取實部作為根。解如下: v = 28.786359708932167 - 0.000000000000449i -0.056673985902386 - 0.000000000229294i -0.000315349578379 + 0.000000000229743i |

新蟲 (初入文壇)
主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗: +1059 |
|
對于你昨天的那個方程,是典型的非線性方程,我修改了程序,用fsolve函數(shù)來解,程序在下面。先新建一個m文件,再把以下全部程序直接復(fù)制進m文件中,運行即可,不要在主窗口中運行,用m文件來寫程序,比直接在主窗口中操作方便,也算是一種好習(xí)慣吧。 由計算結(jié)果可知,解更符合實際,而且沒有出現(xiàn)虛數(shù)。 function solution_nonl clear all;clc global aa bb tt pp %-------------根據(jù)計算結(jié)果,修改迭代初值------------------------------------ xx0=[1]; %迭代初值,如果計算結(jié)果不符合實際,可修改,一般由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; 計算結(jié)果: 28.7864 28.6871 28.5380 28.3969 28.1413 28.0084 27.8591 |

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

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

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 26年申博 +3 | staryer 2026-03-30 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 296求調(diào)劑 +4 | 汪?! 2026-03-31 | 7/350 |
|
|
[考研] 320分,材料與化工專業(yè),求調(diào)劑 +14 | 一定上岸aaa 2026-03-27 | 18/900 |
|
|
[考研] 290求調(diào)劑 +5 | dfffsar 2026-03-29 | 5/250 |
|
|
[考研] 303分 0807學(xué)碩求調(diào)劑 +3 | TYC3632 2026-04-01 | 3/150 |
|
|
[考研] 324求調(diào)劑 +10 | hanamiko 2026-03-26 | 12/600 |
|
|
[考研] 267求調(diào)劑 +13 | uiybh 2026-03-31 | 13/650 |
|
|
[考研] 考研材料工程351分調(diào)劑 +5 | 整個好的 2026-03-31 | 5/250 |
|
|
[考研] 333求調(diào)劑 +4 | 阿科逸 2026-03-31 | 4/200 |
|
|
[基金申請] 面上5B能上會嗎? +8 | redcom 2026-03-29 | 8/400 |
|
|
[考研] 354求調(diào)劑 +3 | lxb598 2026-03-31 | 4/200 |
|
|
[考研] 287求調(diào)劑 +17 | land xuxu 2026-03-26 | 17/850 |
|
|
[有機交流]
10+3
|
kaobao456 2026-03-29 | 4/200 |
|
|
[考研] 0703 化學(xué) 求調(diào)劑,一志愿山東大學(xué) 342 分 +7 | Shern—- 2026-03-28 | 7/350 |
|
|
[考研] 085600,材料與化工321分求調(diào)劑 +10 | 大饞小子 2026-03-28 | 10/500 |
|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +4 | IceGuo 2026-03-26 | 5/250 |
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 求調(diào)劑 +7 | 爭取九點睡 2026-03-28 | 8/400 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|