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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 怎样用matlab或者是mathcad求解含有变参数的四次方程？
81/1返回列表
查看: 2466  |  回复: 7
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

watertxf

铁虫 (初入文坛)

[求助] 怎样用matlab或者是mathcad求解含有变参数的四次方程？

怎样用matlab或者是mathcad求解含有变参数的四次方程？方程中含有虚数，要分别求出方程的四个根的实部和虚部与a的关系，画出关系曲线。a的取值范围为[0 0.2]方程如下

我写的用循环计算，算不出结果
clear
clc
syms x
a=0.001:0.001:0.2;
for i=1:200
    f=a(i)*x^4-5.76*10^20*a(i)*x^2-2.273746798*10^21*x^2+1.332417724416*10^42-20000i*a*x^3+4.626450432*10^25i*x;
    result=solve(f)
end
出现提示如下：
result =

[ empty sym ]

Warning: 200 equations in 1 variables.
> In F:\Program Files\MATLAB\R2013a\toolbox\symbolic\symbolic\symengine.p>symengine at 56
  In mupadengine.mupadengine>mupadengine.evalin at 97
  In mupadengine.mupadengine>mupadengine.feval at 150
  In solve at 170
  In wenxian2000 at 8
Warning: Explicit solution could not be found.
> In solve at 179
  In wenxian2000 at 8
怎样用matlab或者是mathcad求解含有变参数的四次方程？
OWJ))%RZ~M1SBS{OH0F}@%6.jpg

[ Last edited by watertxf on 2013-10-17 at 10:20 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2013-10-17 10:18:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jerkwin

专家顾问 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
用roots函数就可以了
赞一下 回复此楼
2楼2013-10-17 10:25:00
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

watertxf

铁虫 (初入文坛)
引用回帖:
2楼: Originally posted by jerkwin at 2013-10-17 10:25:00
用roots函数就可以了

谢谢！可以具体说一下吗？
赞一下 回复此楼
3楼2013-10-17 11:06:39
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baobiao007

木虫 (职业作家)

中国特色

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
watertxf: 金币+10 2013-10-28 08:25:21
f=a(i)*x^4-5.76*10^20*a(i)*x^2-2.273746798*10^21*x^2+1.332417724416*10^42-20000i*a*x^3+4.626450432*10^25i*x;

仔细看这段代码 a*x^3,  知道问题了吗？
赞一下 回复此楼
我同意叔本华的观点，人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活，逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦
4楼2013-10-17 11:07:05
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
CODE:
clear
clear
clc
syms x
a=0.001:0.001:0.2;
for i=1:200
    f=a(i)*x^4-5.76*10^20*a(i)*x^2-2.273746798*10^21*x^2+1.332417724416*10^42-20000i*a(i)*x^3+4.626450432*10^25i*x;
    result=solve(f)
end

把方程写对就可以了。
赞一下(1人) 回复此楼

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

watertxf
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
5楼2013-10-17 13:48:19
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

watertxf

铁虫 (初入文坛)
送红花一朵
引用回帖:
5楼: Originally posted by 月只蓝 at 2013-10-17 13:48:19
clear
clear
clc
syms x
a=0.001:0.001:0.2;
for i=1:200
    f=a(i)*x^4-5.76*10^20*a(i)*x^2-2.273746798*10^21*x^2+1.332417724416*10^42-20000i*a(i)*x^3+4.626450432*10^25i*x;
    result=solve(f)
...

太感谢了！算出来的结果怎么样提取出每个根的实部和虚部呢？
赞一下 回复此楼
6楼2013-10-17 15:27:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
6楼: Originally posted by watertxf at 2013-10-17 15:27:26
太感谢了！算出来的结果怎么样提取出每个根的实部和虚部呢？...

real(x)
imag(x)
回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
7楼2013-10-17 15:31:33
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

watertxf

铁虫 (初入文坛)
引用回帖:
7楼: Originally posted by 月只蓝 at 2013-10-17 15:31:33
real(x)
imag(x)...

它有四个根，怎么样分别提取出来，并画出，每一个根的实部和虚部随a的变化曲线？就是需要画出八条曲线。谢谢！
赞一下 回复此楼
8楼2013-10-17 15:38:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 watertxf 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭