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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 怎样用matlab或者是mathcad求解含有变参数的四次方程？
81/1返回列表
查看: 2958  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +31 医学老男孩 2026-08-13 69/3450 2026-08-16 00:16 by mo_mo_mo_mo
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
关闭
请填处理意见
关闭