[未答复] 绘制积分图形

[复制链接]
twinsmiki 发表于 2017-10-27 08:53:25
本帖最后由 twinsmiki 于 2017-10-27 09:22 编辑

       如图1所示,绘制该函数的积分图形,式中a=1;R=0.2,r在[a-R,a+R]之间。本人用matlab绘制出来的图形一直是不正确的,根据错误提示,Imaginary parts of complex X and/or Y arguments ignored,进行更改程序,但是得到的图像还是不正确的,正确的图像应为图2所示,横坐标在[0.8,1.2]之间。请大家帮忙看看这个程序是哪里出了错,谢谢!!!
       程序如下:F=@(r,theta) 2.*sqrt(1-(r.*cos(theta)-0.2).^2-(r.*sin(theta)).^2);
                      b=@(r) acos((2.5).*(r.^2-0.96)./r);%积分上限
                      Q=@(r) quadgk(@(theta) F(r,theta),0,b(r));
                      B=0.8:0.05:1.2;%r的变化范围
                      Qr=arrayfun(@(r) Q(r),B);
                      plot(B,real(Qr));hold on;
                      plot(B,imag(Qr))


图1

图1

图2

图2

1 条回复


twinsmiki 发表于 2017-10-27 15:27:19
是本人之前计算错误了,上面的内容是可以相对应的,没有问题的!!!
您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
热门教程
站长推荐
快速回复 返回顶部 返回列表