2020-06-05MATLAB 使用匿名函数进行三重积分求解

2020-06-06  本文已影响0人  锅炉工的自我修养

需求,求解含有两个自由参数的三重积分

https://blog.csdn.net/qq_34374664/article/details/79186465

int 的用法,以及二重积分

f=@(x,y,X,Y,Z) Y.*cos(Z)./sqrt(x.^2+Y.^2+(x-X).^2-2.*y.*Y.*cos(Z));
Q1=@(x,y) integral3(@(X,Y,Z) f(x,y,X,Y,Z),3.47,3.55,0.3,0.6,0,2*pi,'Method','tiled');
delta_x=0.05;
delta_y=0.05;
x=[0.1:delta_x:3.1];
len_x=length(x);
y=[0.1:delta_y:0.3];
len_y=length(y);
[x0,y0]=meshgrid([0.1:delta_x:3.1],[0.1:delta_y:0.3]);

for i= 1:len_x
    for j= 1:len_y
        val(j,i)=Q1(x(i),y(j));
    end
end
% val_every(:,:,k)=val;
pcolor(x0,y0,val);
shading interp
figure
plot(val(1,:),'r-')

refs: https://ww2.mathworks.cn/help/matlab/ref/integral3.html

上一篇 下一篇

猜你喜欢

热点阅读