24小时热门版块排行榜    

查看: 1517  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 科研孤儿太难了 +14 我4大白菜 2026-08-20 15/750 2026-08-20 22:33 by 袁向阳007
[基金申请] 让我中一个面上吧! +8 大萍1987 2026-08-20 10/500 2026-08-20 21:48 by 大萍1987
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 22/1100 2026-08-20 17:03 by yaoyewhu2008
[基金申请] 放榜前的不淡定 40+3 snowwithsea 2026-08-19 12/600 2026-08-20 15:56 by lfy8008
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见