Why do I get an error when I ran the integral?

Question

Hadeel Obaid 2024-2-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2079061-why-do-i-get-an-error-when-i-ran-the-integral

评论： Walter Roberson 2024-3-23

Hi everyone, I got an error when I run following the code, I think because of the symbolic piecewise functions or the nested integrals, or maybe something else. Could anyone help here, please?

clear all; clear; clc;
hU = 100;
NU = 10;
x0 = 150;
a1 = 11.95;
b1 = 0.135;
drm = sqrt(hU^2 + x0^2);
DML = 0.57 * drm^1.6;
% AL = zeros(size(50:10:300));
j = 1;
Rs = 50:10:300;
for i = 50:10:300
   
    xl = sqrt(i^2 - x0^2 + hU^2);
    xh = sqrt(i^2 + x0^2 + hU^2);
     syms x r
    
    f(x) = piecewise(x >= hU & x <= xl, (2 * x) / i^2, x > xl & x <= xh, (2 * x) / (pi * i^2) * acos(((x^2) - hU^2 + x0^2 - i^2) / (2 * x0 * sqrt(x^2 - hU^2))));
    f(r) = piecewise(r >= hU & r <= xl, (2 * r) / i^2, r > xl & r <= xh, (2 * r) / (pi * i^2) * acos(((r^2) - hU^2 + x0^2 - i^2) / (2 * x0 * sqrt(r^2 - hU^2))));
    
  
    DLM(r) = r^1.6;
    
    PL(r) = 1 / (1 + a1 * exp(-b1 * ((180/pi)*asin(hU / r) - a1)));
    PN(r) = 1 - PL(r);
    
    PL(x) = 1 / (1 + a1 * exp(-b1 * ((180/pi)*asin(hU / x) - a1)));
    PN(x) = 1 - PL(x);
    
    
    AL(j) = NU * vpa(int(PL(r) * f(r) * (int(PL(x) * f(x), r, xh) + int(PN(x) * f(x), DLM(r), xh))^((NU)-1), hU, DML), 6);
    
    j = j + 1;
  
end
plot(Rs, AL, 'b--o')
grid on
xlabel('Avg. Cell radius in meters')
ylabel('Association Probability')
legend('LOS Association Prob.')

5 个评论
显示 3更早的评论隐藏 3更早的评论

Walter Roberson 2024-2-7

编辑：Walter Roberson 2024-2-7

在 MATLAB Online 中打开

clear all; clear; clc;

hU = 100;

NU = 10;

x0 = 150;

a1 = 11.95;

b1 = 0.135;

drm = sqrt(hU^2 + x0^2);

DML = 0.57 * drm^1.6;

% AL = zeros(size(50:10:300));

j = 1;

%Rs = 50:10:300;

Rs = 50:10:100;

for i = Rs

xl = sqrt(i^2 - x0^2 + hU^2);

xh = sqrt(i^2 + x0^2 + hU^2);

syms x r

f(x) = piecewise(x >= hU & x <= xl, (2 * x) / i^2, x > xl & x <= xh, (2 * x) / (pi * i^2) * acos(((x^2) - hU^2 + x0^2 - i^2) / (2 * x0 * sqrt(x^2 - hU^2))));

f(r) = piecewise(r >= hU & r <= xl, (2 * r) / i^2, r > xl & r <= xh, (2 * r) / (pi * i^2) * acos(((r^2) - hU^2 + x0^2 - i^2) / (2 * x0 * sqrt(r^2 - hU^2))));

DLM(r) = r^1.6;

PL(r) = 1 / (1 + a1 * exp(-b1 * ((180/pi)*asin(hU / r) - a1)));

PN(r) = 1 - PL(r);

PL(x) = 1 / (1 + a1 * exp(-b1 * ((180/pi)*asin(hU / x) - a1)));

PN(x) = 1 - PL(x);

AL(j) = NU * vpa(int(PL(r) * f(r) * (int(PL(x) * f(x), r, xl, xh) + int(PN(x) * f(x), DLM(r), xh))^((NU)-1), hU, DML), 6);

AL(j)

j = j + 1;

end

Warning: Unable to check whether the integrand exists everywhere on the integration interval.

ans =

NaN

ans =

NaN

ans =

NaN

ans =

NaN

ans =

NaN

Warning: Unable to check whether the integrand exists everywhere on the integration interval.

ans =

NaN

dAL = double(AL);

plot(Rs, dAL, 'b--o')

grid on

xlabel('Avg. Cell radius in meters')

ylabel('Association Probability')

legend('LOS Association Prob.')

Hadeel Obaid 2024-2-7

Below is the error as @Walter Roberson posted:

Warning: Unable to check whether the integrand exists everywhere on the integration interval.

Walter Roberson 2024-2-8

Several of your int() do not indicate which variable to integrate over, and are guessing incorrectly.

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2024-2-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2079061-why-do-i-get-an-error-when-i-ran-the-integral#answer_1404096

在 MATLAB Online 中打开

    f(x) = piecewise(x >= hU & x <= xl, (2 * x) / i^2, x > xl & x <= xh, (2 * x) / (pi * i^2) * acos(((x^2) - hU^2 + x0^2 - i^2) / (2 * x0 * sqrt(x^2 - hU^2))));

There is no clause for x > xh or x < hU

5 个评论
显示 3更早的评论隐藏 3更早的评论

Wafaa 2024-3-23

编辑：Walter Roberson 2024-3-23

在 MATLAB Online 中打开

Warning: Unable to check whether the integrand

exists everywhere on the integration interval.

> In symengine

In sym/int (line 162)

In pp (line 31)

Warning: Unable to check whether the integrand

exists everywhere on the integration interval.

R(x,y)=piecewise(y <= x & x <= 1 & 0 <= y, y + 1, x < y & 0 <= x & y <= 1, x + 1);

%define xi
xi=zeros(1,m);
for i=1:m
    xi(i)=(i-1)/(m-1);   
end
%define psi(i)=L(R(x,y)) where y=xi
Lpsi=sym(zeros(1,m));
psi=sym(zeros(1,m));
for i=1:m
    % psi(i)=ff(xi(i))-(R(x,xi(i))^3)*(x/2)+diff(R(x,xi(i))^3)*((x/3)-(x^2/2))+diff(R(x,xi(i))^3,2)*((x/4)-(2*x^2/4)+(x^3/3));%int(K(x,y)*(R(x,y))^3,y,0,xi(i))+int(K(x,y)*(R(x,y))^3,y,xi(i),1);  
    psi(i)=ff(xi(i))-int(K(x,y)*(R(x,y))^3,y,0,1);
    Lpsi(i)=ff(x)-int(K(x,y)*subs(psi(i),x,y)^3,y,0,1);
    % Lpsi(i)=psi(i)^3*(x/2)+diff(psi(i)^3)*((x/3)-(x^2/2))+diff(psi(i)^3,2)*((x/4)-(2*x^2/4)+(x^3/3));
end

please help me

Walter Roberson 2024-3-23

在 MATLAB Online 中打开

your xl is mostly complex-valued, which is outside the valid piecewise() range, so you mostly get NaN

hU = 100;

NU = 10;

x0 = 150;

a1 = 11.95;

b1 = 0.135;

drm = sqrt(hU^2 + x0^2);

DML = 0.57 * drm^1.6;

% AL = zeros(size(50:10:300));

j = 1;

%Rs = 50:10:300;

Rs = 50:10:100;

for i = Rs

xl = sqrt(i^2 - x0^2 + hU^2)

xh = sqrt(i^2 + x0^2 + hU^2)

syms x r

f(x) = piecewise(x >= hU & x <= xl, (2 * x) / i^2, x > xl & x <= xh, (2 * x) / (pi * i^2) * acos(((x^2) - hU^2 + x0^2 - i^2) / (2 * x0 * sqrt(x^2 - hU^2))))

%f(r) = piecewise(r >= hU & r <= xl, (2 * r) / i^2, r > xl & r <= xh, (2 * r) / (pi * i^2) * acos(((r^2) - hU^2 + x0^2 - i^2) / (2 * x0 * sqrt(r^2 - hU^2))));

DLM(r) = r^1.6;

%PL(r) = 1 / (1 + a1 * exp(-b1 * ((180/pi)*asin(hU / r) - a1)));

%PN(r) = 1 - PL(r);

PL(x) = 1 / (1 + a1 * exp(-b1 * ((180/pi)*asin(hU / x) - a1)));

PN(x) = 1 - PL(x);

inner3 = PL(x) * f(x)

inner1 = int(inner3, x, xl, xh)

inner4 = PN(x) * f(x)

inner2 = int(inner4, x, DLM(r), xh)

outer = int(PL(r) * f(r) * (inner1 + inner2)^((NU)-1), r, hU, DML)

AL(j) = NU * vpa(outer, 6);

AL(j)

j = j + 1;

end

xl = 0.0000e+00 + 1.0000e+02i

xh = 187.0829

f(x) =

inner3 =

inner1 =

NaN

inner4 =

Warning: Unable to check whether the integrand exists everywhere on the integration interval.

inner2 =

outer =

NaN

ans =

NaN

xl = 0.0000 +94.3398i

xh = 190

f(x) =

NaN

inner3 =

NaN

inner1 =

NaN

inner4 =

NaN

inner2 =

NaN

outer =

NaN

ans =

NaN

xl = 0.0000 +87.1780i

xh = 193.3908

f(x) =

NaN

inner3 =

NaN

inner1 =

NaN

inner4 =

NaN

inner2 =

NaN

outer =

NaN

ans =

NaN

xl = 0.0000 +78.1025i

xh = 197.2308

f(x) =

NaN

inner3 =

NaN

inner1 =

NaN

inner4 =

NaN

inner2 =

NaN

outer =

NaN

ans =

NaN

xl = 0.0000 +66.3325i

xh = 201.4944

f(x) =

NaN

inner3 =

NaN

inner1 =

NaN

inner4 =

NaN

inner2 =

NaN

outer =

NaN

ans =

NaN

xl = 0.0000 +50.0000i

xh = 206.1553

f(x) =

inner3 =

inner1 =

NaN

inner4 =

Warning: Unable to check whether the integrand exists everywhere on the integration interval.

inner2 =

outer =

NaN

ans =

NaN

dAL = double(AL);

plot(Rs, dAL, 'b--o')

grid on

xlabel('Avg. Cell radius in meters')

ylabel('Association Probability')

legend('LOS Association Prob.')

请先登录，再进行评论。

Why do I get an error when I ran the integral?

5 个评论
显示 3更早的评论隐藏 3更早的评论

回答（1 个）

5 个评论
显示 3更早的评论隐藏 3更早的评论

另请参阅

类别

标签

Community Treasure Hunt

Why do I get an error when I ran the integral?

5 个评论 显示 3更早的评论隐藏 3更早的评论

回答（1 个）

5 个评论 显示 3更早的评论隐藏 3更早的评论

另请参阅

类别

标签

Community Treasure Hunt

WeChat

5 个评论
显示 3更早的评论隐藏 3更早的评论

5 个评论
显示 3更早的评论隐藏 3更早的评论