Тень и вычислить определенную область

Question

Я попытался изменить код так, чтобы только первая область была затушевана серым цветом. Как установить горизонтальную линию таким образом, чтобы она отображалась только в области, которую я хочу затенять?

Кроме того, я хочу рассчитать площадь ОДНОГО региона. Как мне это достичь? Я знаю, что это trapz, но я не уверен, как установить границы. Благодаря!

x = 0:.01:4*pi; %// x data 
y = sin(x);  %// y data 
level = 0.5;  %// level 
plot(x, y) 
hold on 
area(x, max(y, level), level, 'EdgeColor', 'none', 'FaceColor', [.7 .7 .7]) 

Curve: -

Answer 1

вы можете попробовать и этот простой вариант:

x = 0:.01:4*pi; %// x data 
y = sin(x);  %// y data 
level = 0.5;  %// level 
lineStart = find(y>=level,1); 
lineEnd = find(y(lineStart:end)<=level,1)+lineStart; 
plot(x,y) 
hold all 
area(x(lineStart:lineEnd),y(lineStart:lineEnd),... 
    level,'EdgeColor', 'none', 'FaceColor', [.7 .7 .7],'ShowBaseLine','off') 
line([x(lineStart),x(lineEnd)],[level level ]) 
hold off 

без определения областей, представляющих интерес а-априорная:

И не забудьте hold off. ..

Калориметрическое пространство: A = trapz(x(lineStart:lineEnd),y(lineStart:lineEnd))

Answer 2

Вы можете ограничить диапазон вашей оси х в area участке в интересующем диапазоне, например, от 0 до 4, а затем вычислить результирующие значения функции в этом диапазоне. Для базовой линии вы можете скрыть ее в команде area и добавить ее вручную, используя команду line.

x = 0:.01:4*pi; %// x data 
y = sin(x);  %// y data 
level = 0.5;  %// level 
plot(x, y) 
hold on 

x_interest = 0:.01:4; 
y_interest = sin(x_interest); 
area(x_interest, max(y_interest, level), level, ... 
    'EdgeColor', 'none', 'FaceColor', [.7 .7 .7], ... 
    'ShowBaseLine', 'off'); 
line([ min(x_interest) max(x_interest) ], [ level level ])