2

Googleでこれに近いものを見つけることができなかったので、私の質問自体に欠陥があるのではないかと心配しています.

さまざまな固定ダイナミック レンジで値 (Z) のマトリックスを表示したいと考えています。この場合、0dB、10dB、...、40dB に固定されます。

私の現在のアプローチは、Zmag = abs(Z).^2, Zn = normalized(Zmag), Zdb = 10*log10(1+Zn) を見つけることです。

異なるダイナミック レンジ (たとえば 10dB) を表示するには、Zdb を見つける前に 'Zn(Zn<0.1)=0.1' を含めます。20dB についても同じことを行いますが、対象の値が 0.01 に変更されます。

次に、Zn のカラー メッシュ プロットを作成し、XY (上、3D 視点から) プロットを表示して、imagesc(Zn) が与えるものと同様のものを確認します。意図は、ダイナミック レンジを大きくすると、より詳細なプロットが表示されるようにすることです (この場合、最大値と最小値の間により多くの色が表示されます)。

私の現在の方法は、10dB の場合 (そうあるべきだと思います) 実行されています: 10dB のダイナミック レンジ メッシュ 40dB と比較: 40dB のダイナミック レンジ メッシュ プロット

しかし、0、20、30、および 40dB のプロットの違いはわかりません。値が 0dB から 40dB まで徐々に増加すると予想されます。

-ディラン

編集: ここにいくつかのサンプル コードがあります。これは実際のコードの抜粋ですが、実行する必要があります。

%% Constants
fnum = 1;
Fc = 1/16;
taup = 128;
taumin = 1;
taumax = 512;
taux = taumin:taumax;

%% Signal
l = 1:16; %Signal length
s = sin(2*pi*Fc*l); %Original Signal
sig = zeros([1 taup+512]);
sig(taup:taup+size(l,2)-1) = s;

[mfr,fdy] = MatchedFilterResponse(sig,taup,l);
Z = mfr;

slices = true;
%full dynamic range
name = 'Short Tone Ping results with 0dB range';
Zmag = abs(Z).^2;
Zn = normalizeMat(Zmag);
Zdb = 10*log10(1+Zn);
fnum = plotSurfaces(taux,fdy,Zdb,fnum,name,slices);

slices = false;
%40dB dynamic range
name = 'Short Tone Ping results with 40dB range';
Z40mag = Zmag;
Z40n = normalizeMat(Z40mag);
Z40n(Z40n<0.0001) = 0.0001;
Z40db = 10*log10(1+Z40n);
fnum = plotSurfaces(taux,fdy,Z40db,fnum,name,slices);

function [mfr,fdy] = MatchedFilterResponse(sig,taup,l)
    Fdmin = -1/16;
    Fdmax = 1/16;
    Fdinc = (0.125)/(255);
    fdy = linspace(Fdmin,Fdmax,256);
    i = 0;
    for tau = 1:512
        i = i+1;
        j = 0;
        for Fd = Fdmin:Fdinc:Fdmax
            j = j+1;
            a = sig(l+taup-1);
            b = sig(l+tau).*exp(1i*2*pi*Fd*l);
            mfr(j,i) = sum(a.*b);
        end
    end
return
end

function [fnum] = plotSurfaces(taux,fdy,z,fnum,name,slices)
    fid = figure(fnum);
    axes1 = axes('Parent',fid);
    grid(axes1,'on');
    hold(axes1,'all');
    msh = mesh(taux,fdy,z,'Parent',axes1);
    xlabel ('Delay - seconds');
    ylabel ('Frequency offset from center frequency - Cycles/sample');
    zlabel ('Ambiguity function (Normalized Magnitude-Squared)','Visible','off');
    fname = strcat(name,' (Ambiguity Function z(\tau;F_d))');
    title(fname);
    ax = axis;
    axis([50 200 ax(3) ax(4)])
    cb = colorbar('peer',axes1);
    set(get(cb,'ylabel'),'String','Magnitude-Squared (dB)');
    hold off;
    fnum = fnum + 1;
    return
end
4

1 に答える 1

1

ダイナミックレンジを圧縮/拡張する方法はいくつかあり、下からのしきい値処理は1つにすぎません。よりスムーズな関数を使用して、上からしきい値を設定することも、ダイナミックレンジの「よりソフトな」圧縮を実行することもできます。

予想される段階的な増加は、Znの値の分布に依存し、使用するデータとコードの例がないと、コードの原因や問題を見つけるのは困難です（バグである場合）。

とにかく、データを表示することが目的である場合（これ以上の分析を行うのではなく）、データ自体ではなく、カラーマップのダイナミックレンジを圧縮することをお勧めします。これを行う最も速い方法は、カラーバーを右クリックして[インタラクティブカラーマップシフト]を選択することです。これにより、ダイナミックレンジを変更するためにマウスをコルバー上で直接使用できます。

プログラムで実行する場合は、値の分布に応じていくつかのカスタムカラーマップを作成します。たとえば（簡単にするために、赤の値のみを使用しています）

data = 10 .^ randn(100,100);                      % large dynamic range data
ncols = 128;                                      % number of colors in the colormap
R1 = [linspace(0,1, 10)'; ones(ncols - 10, 1); ]; % colormap for the small values
R2 = [zeros(ncols - 10, 1); linspace(0,1,10)'];   % colormap for the large values
G = zeros(ncols, 1);
B = zeros(ncols, 1);
cmap1 = [R1, G, B];
cmap2 = [R2, G, B];

figure
subplot(1,2,1)
rgbplot(cmap1) % plot the colormap values
subplot(1,2,2)
imagesc (data) % plot the data
colormap(cmap1)
colorbar

figure
subplot(1,2,1)
rgbplot(cmap2)
subplot(1,2,2)
imagesc (data)
colormap(cmap2)
colorbar
于 2011-11-01T08:08:37.230 に答える