| 5 | 1/1 | 返回列表 |
| 查看: 3828 | 回復(fù): 7 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
[求助]
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 |
|
對于你昨天的那個(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)主

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 22408 359分調(diào)劑 +4 | Qshers 2026-03-27 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 285求調(diào)劑 +6 | AZMK 2026-03-29 | 9/450 |
|
|
[考研]
|
鐘llll 2026-03-26 | 5/250 |
|
|
[考研] 289求調(diào)劑 +16 | 新時(shí)代材料 2026-03-27 | 16/800 |
|
|
[考研] 085602 化學(xué)工程專碩 340分求調(diào)劑 +4 | qianbai11 2026-03-29 | 4/200 |
|
|
[考研] 295求調(diào)劑 +5 | wei-5 2026-03-26 | 5/250 |
|
|
[考研] 調(diào)劑310 +12 | 溫柔的晚安 2026-03-25 | 13/650 |
|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
[考研] 279求調(diào)劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
[考研] 330分求調(diào)劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑,求導(dǎo)師收 +9 | 天天好運(yùn)來上岸?/a> 2026-03-24 | 10/500 |
|
|
[考研] 一志愿太原理工安全工程300分,求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|
[考研] 一志愿上海理工能源動(dòng)力(085800)310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 考研調(diào)劑 +10 | 呼呼?~+123456 2026-03-24 | 10/500 |
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
[考研] 281求調(diào)劑 +6 | Koxui 2026-03-24 | 7/350 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|