24小时热门版块排行榜    

查看: 1304  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +5 吃吃吃才有意义 2026-03-19 5/250 2026-03-19 16:18 by 30660438
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 10/500 2026-03-19 16:05 by 余麟余
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +11 yangfz 2026-03-17 11/550 2026-03-19 15:06 by houyaoxu
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 7/350 2026-03-19 14:52 by 功夫疯狂
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +3 @taotao 2026-03-19 3/150 2026-03-19 14:07 by peike
[考研] 287求调剂 +3 晨昏线与星海 2026-03-19 4/200 2026-03-19 12:32 by peike
[考研] 0703化学调剂,求各位老师收留 +10 秋有木北 2026-03-14 10/500 2026-03-19 05:52 by anny19840123
[考研] 一志愿华中科技大学,080502,354分求调剂 +4 守候夕阳CF 2026-03-18 4/200 2026-03-18 22:16 by li123456789.
[考研] 354求调剂 +4 Tyoumou 2026-03-18 7/350 2026-03-18 21:45 by Tyoumou
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 268求调剂 +6 简单点0 2026-03-17 6/300 2026-03-18 09:04 by 无际的草原
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 材料与化工专硕调剂 +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by 勇敢太监王公公
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
信息提示
请填处理意见