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

watertxf

铁虫 (初入文坛)

[求助] 用MATLAB编程求复数四次方程

方程如下
y^4-2i*(10^4)*y^3-y^2*((2.4*10^10)^2+((4.8096*10^10)^2-(1*10^9)^2)/x)+2i*10^4*(4.8096*10^10)^2*y/x+(4.8096*10^10*2.4*10^10)^2/x=0
x的取值为0~0.2
求y的四个根的实部和虚部分别随x的变化关系图
得出类似附件中的图形

求高手帮忙,非常感谢!
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 2000PhysRevE.62.bmp
  • 2012-06-12 16:08:10, 1.83 M

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

watertxf

铁虫 (初入文坛)

引用回帖:
3楼: Originally posted by csgt0 at 2012-06-13 17:25:13
如果y实部虚部对x画图,就用
plot(x,real(y0))
figure
plot(x,imag(y0))

但是结果好像跟你给的图不一样

我运行上面的程序怎么会出现下面的问题:
Warning: The argument for the %s format specifier must be of type char (a string).
> In solve>getEqns at 182
  In solve at 67
  In paser20001 at 6
??? Error using ==> solve>getEqns at 182
'  ' is not a valid expression or equation.

Error in ==> solve at 67
[eqns,vars] = getEqns(varargin{:});

Error in ==> paser20001 at 6
t=solve(f,y,'MaxDegree', 4)
5楼2012-06-13 17:49:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与,应助指数 +1
watertxf: 金币+70, ★★★★★最佳答案, 非常感谢! 2012-06-14 09:12:29
syms x y
f=y^4-2i*(10^4)*y^3-y^2*((2.4*10^10)^2+((4.8096*10^10)^2-(1*10^9)^2)/x)+...
    2i*10^4*(4.8096*10^10)^2*y/x+(4.8096*10^10*2.4*10^10)^2/x;
t=solve(f,y,'MaxDegree', 4)
x=0:0.01:0.2;
y0=subs(t);
还不知道你要拿什么画图
计算就用上面的就可以了
showmethemoney
2楼2012-06-13 17:15:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

如果y实部虚部对x画图,就用
plot(x,real(y0))
figure
plot(x,imag(y0))

但是结果好像跟你给的图不一样
showmethemoney
3楼2012-06-13 17:25:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

watertxf

铁虫 (初入文坛)

引用回帖:
2楼: Originally posted by csgt0 at 2012-06-13 17:15:11
syms x y
f=y^4-2i*(10^4)*y^3-y^2*((2.4*10^10)^2+((4.8096*10^10)^2-(1*10^9)^2)/x)+...
    2i*10^4*(4.8096*10^10)^2*y/x+(4.8096*10^10*2.4*10^10)^2/x;
t=solve(f,y,'MaxDegree', 4)
x=0:0.01:0.2;
y0=s ...

非常感谢您的回复!就是画解出的四个y的实部和虚部分别随x的变化图。按照上边您给的程序运行,有点错误,我看不怎么懂,您在帮忙看看吧?非常感谢!
4楼2012-06-13 17:35:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-29 14:24 by etmYJ6d2rquH
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-29 13:56 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +3 tianxiaoxian 2026-09-28 6/300 2026-09-28 22:44 by tianxiaoxian
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见