これをMATLABでどのようにプロットするのですか？私は胆嚢= -22の閾値ラインのようにAプロットしたい

Question

Iコード

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) 

を以下た、閾値ライン上のグラフは、閾値線の下赤色グラフであるべきである青色べきであるが、グラフがすべきこれらの2つの異なる色との連続継ぎ目で、どのように行うか。

Answer 1

NaNに値を設定することによって、グラフの領域をマスクすることができます。データgAllの高解像度補間を作成し、gAll>-22をマスクし、gAll<-22をマスクしたものを2つ作成し、それらを同じ軸にプロットします。

あなたは、論理インデックスを使用してNaNで値を置き換え、interp1(x, gAll, x_HR)でx_HRにgAll_HRあなたのベクトルgAllの高解像度版を作ることができる：

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-')