求二重积分【错误使用 .^ 矩阵维度必须一致。】

clc,clear

m = 1;

p =0.7;%可调

A1 = 0.242665422341411;

A2 = 1.743150000000000;

yita1=exp(-2);

yita2=0.2;

alpha=0;%B分布时，alpha>0

beta=2;

if alpha>0

yita_min=0;

else

yita_min=exp(-beta);

end

yita_max=1;

if alpha>0

B=quadl(@(u) (u.^(alpha-1)) .* ( (1-u).^(beta-1) ), 0, 1);

f_p_d_f = @(yita) 1/B .* ( yita.^(alpha-1) ) .* ( (1-yita).^(beta-1) );

else

f_p_d_f = @(yita) 2/-beta^2 .* log(yita) ./ yita;

end

f = @(x) 1./m .* ( ...

integral(@(yita) (m./x).^(A2./yita +1).* A2./yita.*f_p_d_f(yita) , yita_min,yita1) + ...

integral(@(yita) (m./x).^(A1./yita +1).* A1./yita.*f_p_d_f(yita) , yita1,yita2) + ...

integral(@(yita) (m./x).^(A2./yita +1).* A2./yita.*f_p_d_f(yita) , yita2,yita_max) ...

);

G_func =@(x,y) 1./(2.*.2.*m) .* abs(x-y) .* f(x) .* f(y);

x_min = 1;

x_max = Inf;

y_min = 1;

y_max = Inf;

G = integral2(G_func, x_min, x_max, y_min, y_max)

【

错误使用 .^

矩阵维度必须一致。

出错 Gini_quad3>@(yita)(m./x).^(A1./yita+1).*A1./yita.*(2/-beta^2.*log(yita)./yita)

出错 integralCalc/iterateScalarValued (line 314)

fx = FUN(t);

出错 integralCalc/vadapt (line 132)

[q,errbnd] = iterateScalarValued(u,tinterval,pathlen);

出错 integralCalc (line 75)

[q,errbnd] = vadapt(@AtoBInvTransform,interval);

出错 integral (line 88)

Q = integralCalc(fun,a,b,opstruct);

出错

Gini_quad3>@(x)1./m.*(integral(@(yita)(m./x).^(A2./yita+1).*A2./yita.*(2/-beta^2.*log(yita)./yita),yita_min,yita1)+integral(@(yita)(m./x).^(A1./yita+1).*A1./yita.*(2/-beta^2.*log(yita)./yita),yita1,yita2)+integral(@(yita)(m./x).^(A2./yita+1).*A2./yita.*(2/-beta^2.*log(yita)./yita),yita2,yita_max))

出错 integralCalc/iterateScalarValued (line 323)

fx = FUN(t).*w;

出错 integralCalc/vadapt (line 132)

[q,errbnd] = iterateScalarValued(u,tinterval,pathlen);

出错 integralCalc (line 83)

[q,errbnd] = vadapt(@AToInfInvTransform,interval);

出错 integral (line 88)

Q = integralCalc(fun,a,b,opstruct);

出错 Gini_quad3 (line 39)

test_int = integral(f, 1 ,Inf)

】