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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 用MATLAB编程求复数四次方程
61/1返回列表
查看: 1238  |  回复: 5
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-06-12 16:08:51
已阅   回复此楼   关注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 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭