24小时热门版块排行榜    

查看: 1298  |  回复: 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的回帖

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的回帖

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的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

可能版本不一样吧,我把数据传给你吧,x取小点
我的是2011

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : Book1.xlsx
  • 2012-06-13 17:58:08, 25.08 K
showmethemoney
6楼2012-06-13 17:58:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 watertxf 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +6 司空. 2026-03-18 6/300 2026-03-18 23:03 by 星空星月
[考研] 274求调剂 +5 S.H1 2026-03-18 5/250 2026-03-18 21:27 by guosr9609
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 286求调剂 +6 lemonzzn 2026-03-16 9/450 2026-03-18 15:31 by vgtyfty
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 9/450 2026-03-18 15:05 by stone_128
[考研] 331求调剂(0703有机化学 +7 ZY-05 2026-03-13 8/400 2026-03-18 14:13 by 007_lilei
[考研] 材料,纺织,生物(0856、0710),化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 321求调剂(食品/专硕) +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
信息提示
请填处理意见