如何在matlab中绘制这个？

Question

我有以下代码

xr=randi([1 150],1,20)
    z=numel(xr);
    N=10;  %Window length
    gAll=zeros(1,z-N+1);
    for n=0:z-N;
        x=xr(1+n:N+n)
      d=max(x);
        m=numel(x);
        y=zeros(d,1);
        p=zeros(d,d);
        for k=1:m-1
        y(x(k))=y(x(k))+1; 
        p(x(k),x(k+1))=p(x(k),x(k+1))+1;
        end
        p=bsxfun(@rdivide,p,y);
        p(isnan(p)) = 0; 
        j=prod(p(p~=0));
        [~,~,idx] = unique(x);
        q=prod(hist(idx,1:max(idx))/numel(x));
        s=log(j);
        l=log(q);
        g=s+l
    gAll(n+1)=g;
    end
    plot(gAll)

我想要一个图，对于 gAll =-22 的阈值线，阈值线上方的图形应为红色，阈值线下方的图形应为蓝色，但图形应与这两种不同颜色连续连接，如何去做吧。

Answer 1

您可以通过将值设置为NaN 来屏蔽图形区域。我将创建数据gAll 的高分辨率插值，然后创建两个副本，一个带有gAll>-22 掩码，一个带有gAll<-22 掩码，并将它们绘制在相同的轴上。

您可以使用interp1(x, gAll, x_HR) 将gAll_HR 设为x_HR 上的矢量gAll 的高分辨率版本，然后使用逻辑索引将值替换为NaN：

gAll_low = gAll_HR;
gAll_low(gAll_HR>=-22) = NaN;
gAll_high = gAll_HR;
gAll_high(gAll_HR<-22) = NaN;

plot(x_HR, gAll_low, 'b-', x_HR, gAll_high, 'r-')