I am not able to integrate heaviside(y-f) for x= 0 to pi/2 and y = 0 to pi/2 (only syms integration "int"). After running the code, output is shown in below and i am not able to get the numerical answer. Can you give some ideas on how to proceed?

2 次查看（过去 30 天）

CHINTALAPUDI GNANA SURYA NARAYANA SAI 2020-8-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/579165-i-am-not-able-to-integrate-heaviside-y-f-for-x-0-to-pi-2-and-y-0-to-pi-2-only-syms-integration

评论： hosein Javan 2020-8-14

syms x y
f = sin(x);
integ = int( int(heaviside(y-f), x, 0, pi/2),y, 0, pi/2)

Result:

integ =

int(int(heaviside(y - sin(x)), x, 0, pi/2), y, 0, pi/2)

4 个评论
显示 2更早的评论隐藏 2更早的评论

Walter Roberson 2020-8-14

For your purposes, is heaviside(0) equal to 0, or is it equal to 1, or is it equal to 1/2, or is it undefined ?

CHINTALAPUDI GNANA SURYA NARAYANA SAI 2020-8-14

编辑：CHINTALAPUDI GNANA SURYA NARAYANA SAI 2020-8-14

yes, it is heaviside(=0 or >0) equal to 1 and heaviside(<0) equal to 0.

请先登录，再进行评论。

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

采纳的回答

Walter Roberson 2020-8-13

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/579165-i-am-not-able-to-integrate-heaviside-y-f-for-x-0-to-pi-2-and-y-0-to-pi-2-only-syms-integration#answer_479580

在 MATLAB Online 中打开

Reverse the order of integration.

>> integ = int( int(heaviside(y-f), y, 0, pi/2),x, 0, pi/2)
 
integ =
 
pi^2/4 - 1

更多回答（1 个）

hosein Javan 2020-8-13

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/579165-i-am-not-able-to-integrate-heaviside-y-f-for-x-0-to-pi-2-and-y-0-to-pi-2-only-syms-integration#answer_479556

在 MATLAB Online 中打开

if you need numerical integration, some functions cannot be easily evaluated by symbolic. instead use "integral". in here we have a double integral. the integral calculation method is 'iterated' which means on each iteration it is integrating on one dimension and then the next. it takes a few seconds to accomplish though.

fun = @(x,y) double( y > sin(x) )
q = integral2(fun, 0, pi/2, 0, pi/2, 'Method', 'iterated', 'Reltol', 1e-8)
% answer in command window:
q =
    1.4674

here's your function plot. blue part is zero. yellow part equals one.