matlab绘制函数图像
(2011-06-02 14:46:40)
标签：
杂谈
前两天在一个偶然机会下画了这个函数图像，命令如下：
plot(linspace(-2,2,100),(linspace(-2,2,100).^(1/3))
发现一个很有意思的现象：负半轴的图像处于正半轴，命令行提示：Warning: Imaginary parts of complex X
and/or Y arguments
ignored.且图形画出来并不对称，经查证提示的错误是指算出的数值为复数，但在绘图时自动省略虚部。这样，就把一个典型的奇函数画成了一个“四不象”！思考后个人估计出现复数的根本原因如下：在MATLAB点的分数次方计算中，底数不是以我们正常认为的实数做法出现的，例如：初中生都知道(-8)^(1/3)=-2,但用MATLAB算出的结果却是1+1.7321i。
说明在MATLAB中-8不是-8，而是8(cos(π)+isin(π))，这样，在作开三次根号时就有2(cos(π/3)+isin(π/3))＝1+1.7321i。
当时的解决办法如下：
x=linspace(-2,2,200);
x1=x([1:length(x)/2]);
x2=x([length(x)/2+1:end]);
y=[-(-x1).^(1/3)’;x2.^(1/3)’];
line(x,y,’color’,’r’)
呵呵…办法比较老土！但在MATLAB群中交流这个问题时，朋友们给出了更好的建议：第一位朋友建议对上述方法作如下更改：
x=linspace(-2,2,100); y=sign(x).*(abs(x).^(1/3));
plot(x,y,’r’)
这样看起来就好多啦！然后另外一位朋友发现MATLAB7.0中已经自带一个“nthroot”命令来解决这个问题(据另一网友说在6.5中未发现这个命令，我还没有确实这个问题)，这样这个问题就圆满解决了
分享：
喜欢
0
赠金笔
加载中，请稍候……
评论加载中，请稍候…
发评论
登录名： 密码： 找回密码 注册记住登录状态
昵 称：
评论并转载此博文
发评论
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。
发布者：全栈程序员-站长，转载请注明出处：https://javaforall.net/217013.html原文链接：https://javaforall.net
