带2参数的数值积分

yxlnbu

回复 15 # qibbxxt 的帖子

谢谢，我是初学matlab。刚开始是用c的，对于单元数组之类的还是不太熟悉，谢谢您的指教。

qibbxxt

回复 16 # yxlnbu 的帖子

对于输入或者输出不多的情况，也可以不用单元数组

1. rr = linspace(40.00 , 45.00 , 100) ;          %产生r=40到45的数组
   
2. zz = linspace(-4.5 , 4.5 , 6) ;               %产生每匝的z位移从-4.5到4.5
   
3. f = @(r,z)quad(@(x) 2*(1600 -40*r*cos(x))./(r^2-80*r*cos(x)+1600+z^2).^1.5,0,pi);
   
4. [z,r]=meshgrid(zz,rr);
   
5. b = arrayfun(f,r,z);
   
6. b=[b,sum(b,2),rr'];

复制代码