在matlab中把超过阈值的区域用线标出来

问题提出

本问题来源于“问答”上看到的一个问题，感觉挺有趣的，就来试一下。

想法

对于已知的数据，我们可以遍历一遍，然后获得大于阈值的数据所在的位置，然后通过判断位置是否相邻可以获得以闭区间表示的区域，然后将各区域画在图上就可以了。

代码实现

 
x = 0.5:0.5:50;
y = 2*rand(1,100);
v = 1;
plotAreaOverV(x, y, v);
 
function plotAreaOverV(x, y, v)
lineY=0;%黄色线所在位置
 
plot(x,y);
hold on
plot(x, v*ones(size(x,2)));%画出阈值线
 
areaOverX = find(y>v);%获得大于阈值的数据所在位置
startp = areaOverX(1);
startx = x(startp);
endp = areaOverX(1);
endx = x(endp);
 
for i=2:size(areaOverX,2)
    if areaOverX(i)-endp <= 1 %前后数据相邻
        endp = areaOverX(i);
        endx = x(endp);
    else
        hold on;
        dstartx = (x(max(startp-1,1)) - startx)/2; %起始点向前移动半个单位
        dendx = (x(min(endp+1,size(x,2))) - endx)/2;%结束点向后移动半个单位
        plot([startx+dstartx endx+dendx],[lineY lineY],Color='y',LineWidth=5);%画出前一段连续数据
        startp = areaOverX(i);
        startx = x(startp);
        endp = areaOverX(i);
        endx = x(endp);
    end
end
 
%画出最后一段线
dstartx = x(max(startp-1,1)) - startx;
dendx = x(min(startp+1,size(x,2))) - endx;
plot([startx+dstartx/2 endx+dendx/2],[lineY lineY],Color='y',LineWidth=5);
 
end

总结

结果看起来还是可以的，不过边界点我这里是向左右延伸了半个单位，处理地比较粗糙，可以改进为根据阈值在前后两点连线上的x值。