24小时热门版块排行榜    

查看: 3806  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

dutudou

铜虫 (小有名气)

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

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


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

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

我就是我!!!
已阅   回复此楼   关注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的回帖
查看全部 8 个回答

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +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的量级,工程上完全可以忽略掉。

» 本帖已获得的红花(最新10朵)

每个人的人生都是时代的产物。每个人的血统、命运和梦想组成了现在这个时代。
4楼2013-07-05 11:12:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 一志愿深大,0703化学,总分302,求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 307求调剂 +3 余意卿 2026-03-18 3/150 2026-03-21 17:31 by ColorlessPI
[考研] 265求调剂 +12 梁梁校校 2026-03-19 14/700 2026-03-21 13:38 by lature00
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 265求调剂 +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[考研] 265求调剂 +9 梁梁校校 2026-03-17 9/450 2026-03-21 02:17 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 材料专业求调剂 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 材料学求调剂 +4 Stella_Yao 2026-03-20 4/200 2026-03-20 20:28 by ms629
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 广西大学家禽遗传育种课题组2026年硕士招生(接收计算机专业调剂) +3 123阿标 2026-03-17 3/150 2026-03-20 15:58 by 飞行琦
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
信息提示
请填处理意见