应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab解一元三次方程，得到的都是复数根。
81/1返回列表
查看: 3849  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

dutudou

铜虫 (小有名气)

[求助] matlab解一元三次方程，得到的都是复数根。

我想问一下，matlab解一元三次方程，得到的都是复数根。但是虚部前面的系数很小很小了，这个可以把虚部忽略吗，直接取实部作为根。解如下：
v =


28.786359708932167 - 0.000000000000449i
-0.056673985902386 - 0.000000000229294i
-0.000315349578379 + 0.000000000229743i
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
我就是我！！！
1楼 2013-07-05 08:28:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
dutudou: 金币+5, ★★★★★最佳答案, 真心感谢啊！！！ 2013-07-05 11:37:40
对于你昨天的那个方程，是典型的非线性方程，我修改了程序，用fsolve函数来解，程序在下面。先新建一个m文件，再把以下全部程序直接复制进m文件中，运行即可，不要在主窗口中运行，用m文件来写程序，比直接在主窗口中操作方便，也算是一种好习惯吧。
由计算结果可知，解更符合实际，而且没有出现虚数。

function solution_nonl
clear all;clc
global aa bb tt pp

%-------------根据计算结果，修改迭代初值------------------------------------
xx0=[1];  %迭代初值，如果计算结果不符合实际，可修改，一般由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;

计算结果：

28.7864   28.6871   28.5380   28.3969   28.1413   28.0084   27.8591
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
2楼2013-07-05 09:23:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
解精确到小数点后16位的结果：

28.786359708135375  28.687090969223252  28.537987804689500  28.396887298198173  28.141321545916146  28.008427048788107  27.859124010004489
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
3楼2013-07-05 09:25:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tgvictor99

木虫 (知名作家)

迷糊领域荣誉领主
看你需要几位有效数字。虚部已经到了e-9的量级，工程上完全可以忽略掉。
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

dutudou
每个人的人生都是时代的产物。每个人的血统、命运和梦想组成了现在这个时代。
4楼2013-07-05 11:12:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dutudou

铜虫 (小有名气)
送红花一朵
引用回帖:
4楼: Originally posted by tgvictor99 at 2013-07-05 11:12:48
看你需要几位有效数字。虚部已经到了e-9的量级，工程上完全可以忽略掉。

谢谢啦，
回复此楼
我就是我！！！
5楼2013-07-05 11:37:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tgvictor99

木虫 (知名作家)
引用回帖:
5楼: Originally posted by dutudou at 2013-07-05 11:37:57
谢谢啦，...

6楼2013-07-05 15:02:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

戏佳人

铜虫 (小有名气)
读不懂啊，
回复此楼
7楼2013-07-06 15:50:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

DICKYldg

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by 月只蓝 at 2013-07-05 09:23:16
对于你昨天的那个方程，是典型的非线性方程，我修改了程序，用fsolve函数来解，程序在下面。先新建一个m文件，再把以下全部程序直接复制进m文件中，运行即可，不要在主窗口中运行，用m文件来写程序，比直接在主窗口 ...

你好。我也是解方程时出现复数解，是由于有非整数次方出现复数运算，这个如何来避免呢？
赞一下 回复此楼
8楼2015-11-03 09:52:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dutudou 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 070300化学学硕311分求调剂 +17 梁富贵险中求 2026-04-04 19/950 2026-04-11 00:29 by wangjihu
[考研] 广东省 085601 329分求调剂 +13 Eddieddd 2026-04-10 13/650 2026-04-10 23:27 by 314126402
[考研] 085404 298分求调剂 +10 呼啦呼啦呼呼呼 2026-04-10 11/550 2026-04-10 16:44 by wangy0907
[考研] 085601初试330分找调剂 +10 流心奶黄包l 2026-04-09 10/500 2026-04-10 08:14 by Sammy2
[考研] 生物学求调剂 一志愿沪9，326分 +7 刘墨墨 2026-04-06 7/350 2026-04-10 08:11 by kangsm
[考研] 297求调剂 +8 Kwgyz 2026-04-09 8/400 2026-04-09 23:22 by may_新宇
[考研] 生物学调剂，一志愿西南大学348，Top期刊一区二作、二区三作，三等奖学金三次 +4 candyyyi 2026-04-09 4/200 2026-04-09 18:39 by l_paradox
[考研] 材料调剂 +10 18815505510 2026-04-09 11/550 2026-04-09 17:07 by 544594351
[考研] 一志愿电子科技大学085600材料与化工 329分求调剂 +14 Naiko 2026-04-04 14/700 2026-04-09 16:56 by luoyongfeng
[考研] 求调剂 +3 猪肉墩粉条cc 2026-04-08 4/200 2026-04-09 10:05 by 猪肉墩粉条cc
[考研] 求调剂 +7 chenxrlkx 2026-04-05 9/450 2026-04-09 09:04 by wj165256
[考博] 材料方向考博，求推荐 +3 言语aaa 2026-04-05 4/200 2026-04-08 22:22 by nxgogo
[考研] 一志愿华东理工085601材料工程303分求调剂 +15 a1708 2026-04-06 15/750 2026-04-08 16:23 by luoyongfeng
[考研] 298求调剂 +4 残荷新柳 2026-04-07 4/200 2026-04-07 23:02 by lbsjt
[考研] 328求调剂 +4 ghhh88888 2026-04-06 5/250 2026-04-07 14:45 by ghhh88888
[考研] 285求调剂 +5 mapmath 2026-04-06 6/300 2026-04-06 17:18 by 蓝云思雨
[考研] 复试调剂 +5 asdasdassda 2026-04-05 5/250 2026-04-06 09:32 by dongzh2009
[考研] 男生，一志愿沪9生物学071000，初试308求调剂 +3 刘墨墨 2026-04-04 3/150 2026-04-05 08:26 by barlinike
[考研] 290求调剂 +7 luoziheng 2026-04-04 7/350 2026-04-04 23:17 by lqwchd
[考研] 可跨专业调剂 +3 周的得地 2026-04-04 6/300 2026-04-04 22:21 by barlinike
信息提示
关闭
请填处理意见
关闭