查看: 380|回复: 2|关注: 0

[已解决] 使用integral2求二重积分并绘制图像

[复制链接]

新手

9 麦片

财富积分


050


5

主题

9

帖子

0

最佳答案
本人初学matlab,请大神帮帮忙看一下~
感谢~

问题如下:求二重积分,并求F关于h的图像
  1. clear;
  2. clc;
  3. u=4*pi*10.^(-7);
  4. n1=6;
  5. n2=6;
  6. r1=10;
  7. r2=10;
  8. t=0;
  9. c=0;
  10. fun=@(a,h,b)(cos(a)*cos(b)+sin(a)*sin(b)*cos(c))./((r1*cos(a)-r2*cos(b)).^2+(r1*sin(a)*cos(c)+t-r2*sin(b)).^2+(-r1*sin(a)*sin(c)+h).^2).^0.5;
  11. h=-20:2:20;
  12. F=(u*r1*r2*n1*n2/4*pi)*arrayfun(@(h)integral2(@(a,b)fun(a,h,b),0,2*pi,0,2*pi),h);
  13. plot(h,F);
复制代码


拜托大神看看错误出在哪里,万分感谢~

错误提示

错误提示

求解二重积分

求解二重积分

论坛优秀回答者

权威

3372 麦片

财富积分



2

主题

3640

帖子

758

最佳答案
  • 关注者: 154
发表于 2019-2-10 16:21:13 | 显示全部楼层 |此回复为最佳答案
fun=@(a,h,b)(cos(a).*cos(b)+sin(a).*sin(b)*cos(c))./((r1*cos(a)-r2*cos(b)).^2+(r1*sin(a)*cos(c)+t-r2*sin(b)).^2+(-r1*sin(a)*sin(c)+h).^2).^0.5;

新手

9 麦片

财富积分


050


5

主题

9

帖子

0

最佳答案
 楼主| 发表于 2019-2-10 20:56:24 | 显示全部楼层
maple1314168 发表于 2019-2-10 16:21
fun=@(a,h,b)(cos(a).*cos(b)+sin(a).*sin(b)*cos(c))./((r1*cos(a)-r2*cos(b)).^2+(r1*sin(a)*cos(c)+t-r2 ...

:lol感谢大神,改过已经可以做出来,真的还需要认真学习呢
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /4 下一条

快速回复 返回顶部 返回列表