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

23 次查看（过去 30 天）
先生
先生 2022-10-18
编辑： 先生 2022-10-18
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)
】

请先登录，再回答此问题。

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Statistics and Machine Learning Toolbox 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!