查看: 615|回复: 6|关注: 0

[已解决] 二重积分解出来还是int(int)

[复制链接]

新手

7 麦片

财富积分


050


3

主题

11

帖子

0

最佳答案
二重积分解出来还是没有解析解,得到的结果还含有int(int),请问该如何解决,
  1. syms psi yl=0.2/2;%矩形长的一半
  2. d=0.1;%矩形的宽度
  3. t=0;%两者中心线之间的水平距离
  4. h=0.2;%两者之间的间距
  5. r=0.1;%圆的半径
  6. u0=4*pi*10^-7;%真空磁导率
  7. M=u0./(4*pi)*int(int(r*cos(psi)./sqrt((r*cos(psi)-d)^2+(r*sin(psi)+t-y)^2+h^2),psi,0,2*pi),y,-l,l)
复制代码


doubleintigrate.m

314 Bytes, 下载次数: 2

MATLAB 基础讨论
版块优秀回答者

入门

287 麦片

财富积分


50500


1

主题

257

帖子

59

最佳答案
  • 关注者: 9
发表于 2018-5-17 10:08:35 | 显示全部楼层 |此回复为最佳答案
  1. syms psi y
  2. l=0.2/2;%矩形长的一半
  3. d=0.1;%矩形的宽度
  4. t=0;%两者中心线之间的水平距离
  5. h=0.2;%两者之间的间距
  6. r=0.1;%圆的半径
  7. u0=4*pi*10^-7;%真空磁导率
  8. fun=@(psi,y) r*cos(psi)./sqrt((r*cos(psi)-d).^2+(r*sin(psi)+t-y).^2+h^2);
  9. M=u0/(4*pi)*integral2(fun,0,2*pi,-l,l);
复制代码

  1. M =

  2.    4.2225e-09
复制代码

新手

7 麦片

财富积分


050


3

主题

11

帖子

0

最佳答案
 楼主| 发表于 2018-5-29 16:01:32 | 显示全部楼层

今天才看到,太感谢啦,我看看你编的

新手

7 麦片

财富积分


050


3

主题

11

帖子

0

最佳答案
 楼主| 发表于 2018-5-30 16:44:26 | 显示全部楼层
  1. syms psi y
  2. l=0.2/2;%矩形长的一半
  3. d=0.1;%矩形的宽度
  4. t=0:0.001:1;%两者中心线之间的水平距离
  5. h=0.2;%两者之间的间距
  6. r=0.1;%圆的半径
  7. u0=4*pi*10^-7;%真空磁导率
  8. fun=@(psi,y) r.*cos(psi)./sqrt((r.*cos(psi)-d).^2+(r.*sin(psi)+t-y).^2+h.^2);
  9. M1=u0./(4*pi).*integral2(fun,0,2*pi,-l,l);
  10. M=4.*M1;
  11. plot(t,M)
复制代码

新手

7 麦片

财富积分


050


3

主题

11

帖子

0

最佳答案
 楼主| 发表于 2018-5-30 16:47:46 | 显示全部楼层

想请教一下,如果t变化,做二维图为什么会报错

MATLAB 基础讨论
版块优秀回答者

入门

287 麦片

财富积分


50500


1

主题

257

帖子

59

最佳答案
  • 关注者: 9
发表于 2018-5-31 18:12:26 | 显示全部楼层
Dawnshy 发表于 2018-5-30 16:47
想请教一下,如果t变化,做二维图为什么会报错
  1. syms psi y t

  2. l=0.2/2;%矩形长的一半
  3. d=0.1;%矩形的宽度
  4. h=0.2;%两者之间的间距
  5. r=0.1;%圆的半径
  6. u0=4*pi*10^-7;%真空磁导率
  7. i=1;
  8. for t=0:0.001:1
  9.     fun=@(psi,y) r*cos(psi)./sqrt((r*cos(psi)-d).^2+(r*sin(psi)+t-y).^2+h^2);
  10.     M(i)=u0/(4*pi)*integral2(fun,0,2*pi,-l,l);
  11.     i=i+1;
  12. end
  13. t=0:0.001:1;
  14. plot(t,M)
复制代码

新手

7 麦片

财富积分


050


3

主题

11

帖子

0

最佳答案
 楼主| 发表于 2018-6-1 09:14:10 | 显示全部楼层

谢谢大神,我在琢磨琢磨
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

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

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