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

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

主管區(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)寫(xiě)程序,比直接在主窗口中操作方便,也算是一種好習(xí)慣吧。 由計(jì)算結(jié)果可知,解更符合實(shí)際,而且沒(méi)有出現(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 |


| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 316求調(diào)劑 +7 | 梁茜雯 2026-03-19 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 350求調(diào)劑 +6 | weudhdk 2026-03-19 | 6/300 |
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分,求調(diào)劑 +4 | jiajunX 2026-03-22 | 4/200 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-17 | 10/500 |
|
|
[考研] 336化工調(diào)劑 +3 | 王大坦1 2026-03-23 | 4/200 |
|
|
[考研] 291求調(diào)劑 +5 | 孅華 2026-03-22 | 5/250 |
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
[考研] 324求調(diào)劑 +6 | lucky呀呀呀鴨 2026-03-20 | 6/300 |
|
|
[考研] 303求調(diào)劑 +5 | 安憶靈 2026-03-22 | 6/300 |
|
|
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考! +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
[考研] 一志愿中國(guó)石油大學(xué)(華東) 本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-17 | 3/150 |
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
[考研] 一志愿華南師大 070300(化學(xué))304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|