| 查看: 3807 | 回复: 7 | |||||
[求助]
matlab解一元三次方程，得到的都是复数根。
|
|
我想问一下，matlab解一元三次方程，得到的都是复数根。但是虚部前面的系数很小很小了，这个可以把虚部忽略吗，直接取实部作为根。解如下： v = 28.786359708932167 - 0.000000000000449i -0.056673985902386 - 0.000000000229294i -0.000315349578379 + 0.000000000229743i |
回复此楼» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
用matlab解带参数的一元三次方程--求解
已经有3人回复
- Matlab,解复数方程
已经有8人回复
- 求各位前辈帮忙！matlab解复数方程，出现图中的警告，怎来修改呢？谢谢了啊
已经有5人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
dutudou: 金币+5, ★★★★★最佳答案, 真心感谢啊！！！ 2013-07-05 11:37:40
感谢参与，应助指数 +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 |
2楼2013-07-05 09:23:16
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
3楼2013-07-05 09:25:44
tgvictor99
木虫 (知名作家)
迷糊领域荣誉领主
- 应助: 40 (小学生)
- 金币: 2667.4
- 散金: 1907
- 红花: 49
- 沙发: 1
- 帖子: 5238
- 在线: 923.8小时
- 虫号: 2524003
- 注册: 2013-06-27
- 专业: 流体力学
dutudou4楼2013-07-05 11:12:48
5楼2013-07-05 11:37:57
tgvictor99
木虫 (知名作家)
 |
6楼2013-07-05 15:02:49
7楼2013-07-06 15:50:48
8楼2015-11-03 09:52:58
40