Question

MATLAB定积分求解，如下图所示，求各位大神指教

Answer 1

因e和a分别为在一定范围里取值的变量，积分的结果应该是关于这两个变量的二元函数，画出图像来是一个曲面。

 

参考代码：

K1=0.0229;
E = linspace(-1,1,15);
A = linspace(0,1,10);
m = length(E);
n = length(A);
J = zeros(m,n);
for i = 1 : m
    for j = 1 : n
        e = E(i);
        a = A(j);
        f = @(B)(1-e^2*sind(B).^2)./((K2-K1)*pi*a^2*cosd(B).^2+(1-e^2*sind(B).^2)*K2);
        J(i,j) = quadgk(f,0,90);
    end
end
mesh(A,E,J)
xlabel a; ylabel e
view(-35,10)