" Arrays have incompatible sizes for this operation" when computing nested integral

chenjie ruan

2021 9 4

0 个回答

2 次查看（30 天）

0 个投票

function output = tsp_up(Z)     
ap = active_p(Z); % A constant
disp(ap);
theta = const.alpha / 2;
fun = @(x,y) 1 ./ (1 + const.gamma_B ^(-1) .* (x.^(theta)).* (y.^(- theta))) .* exp(- ((Z - const.gamma_B).* const.Noise .* (x.^(theta))) ./ (const.gamma_B .* const.P));
tsp_r = @(y) y .* exp( -pi * ap * const.lambda_BS .* integral(@(x)fun(x,y),0,Inf) - pi .* y.^2 .* const.lambda_BS);
result = integral(tsp_r,0.01,100);
output = 2 * pi * const.lambda_BS * result;
end
y = tsp_up(5.1);

Here I want to compute the nested Integral with variable x(from 0 to Inf) and y(from 0.01 to 100). Everything else is constant. The " Arrays have incompatible sizes for this operation" error has occured in line 5(fun) . I can't see what it has anything to do with the array size. Does anyone know where I am wrong? Any kind of advice is sincerely appreciated!

I try to modify the domain of x and y to be (0.01 to 100) both, this time the error disappear but I totally don't understand what is happening.