匹配追蹤原子的MATLAB Wigner圖

Question

使用MATLAB我將Matching Pursuit應用於近似信號。我的問題是，我很難想象所選原子的時頻表示。我正在嘗試製作類似於以下圖像的Wigner圖（source）。

我都看着小波工具箱，信號處理工具箱以及開源的時頻工具箱，但我可能只是使用了錯誤的參數，因爲我與信號處理經驗是相當有限。

例

使用this data我的目標是從上面複製的情節。

% fit the signal using MP 
itermax = 50; 
signal = load('signal.txt'); 
dict = wmpdictionary(length(signal)); 
[signal_fit, r, coeff, iopt, qual, X] = wmpalg('OMP', signal, dict, ... 
               'itermax', itermax); 

% wigner plot of the simulated signal 
tfrwv(signal_fit) % wigner-ville function from time-frequency toolbox 

% wigner plot of each atom 
atoms = full(dict(:, iopt)) % selected atoms 
for i = 1:itermax 
    tfrwv(atoms(:, i)) 
end 

不幸的是，沒有一個結果圖接近目標可視化。請注意，在這個例子中，我使用tfrwv和標準參數，這些標準參數隨着它打開的GUI進行調整。

我非常感謝您的幫助。

更新

我想我現在已經明白，一個需要使用的Gabor原子得到斑點與形狀類似拉伸高斯。不幸的是，在Signal Processing Toolbox的預定義字典中沒有Gabor函數。然而，this question幫我實現所需的字典，這樣，我得到原子看起來頗爲相似，例如：

由於我的情節接近，但還不完善，還存在兩個問題打開：

• 我們在第一個例子中看到的所有斑點都可以單獨由Gabor原子建模，還是我需要另一個函數字典？
• 我該如何將獨立的imagesc圖組合成一個可視化圖？

Answer 1

要回答你的第二個問題「我怎樣才能結合indidividual於imagesc地塊成一個單一的可視化？

如果您有多個要疊加並使用imagesc進行顯示的2d矩陣，我建議採用元素最大值。

例如，我生成兩個具有不同均值和方差的高斯分佈的31x31網格。

function F = generate2dGauss(mu, Sigma) 
    x1 = -3:.2:3; x2 = -3:.2:3; 
    [X1,X2] = meshgrid(x1,x2); 
    F = mvnpdf([X1(:) X2(:)],mu,Sigma); 
    F = reshape(F,length(x2),length(x1)); 
end 

F1 = generate2dGauss([1 1], [.25 .3; .3 1]); 
F2 = generate2dGauss([-1 -1], [.1 .1; .1 1]); 

我可以積小樣方爲它們在你的榜樣，

figure; 
subplot(1,2,1); 
title('Atom 1'); 
imagesc(F1); 

subplot(1,2,2); 
title('Atom 2'); 
imagesc(F2); 

或者，我可以將每個元素最多兩個網格的繪製。

figure; 
title('Both Atoms'); 
imagesc(max(F1, F2)); 

您也可以嘗試用元素方面的手段，資金等，但根據你給的例子，我覺得最大會給你最乾淨的期待結果。不同的功能

可能的利弊：

1. 最大將工作最好的，如果你的原子總是零值的背景和沒有負值。如果背景爲零值，但原子也包含負值，則負值可能被其他原子的背景掩蓋。如果你的原子的重疊，更高的值當然會占主導地位。
2. 意味着你的峯值會變得更低，但是在原子之間有重疊的地方可能更直觀。
3. 總和會使重疊的區域更大。
4. 如果您有非零背景，您也可以嘗試使用邏輯索引。你必須做出一些關於在重疊區域做什麼的決定，但是它可以很容易地濾除背景。

Answer 2

問：如何將獨立圖像c圖組合成單個可視化圖？ A.使用子圖繪製多個圖，在圖中查找下面的樣本和2乘以2的圖。更改公式中的代碼

x = linspace(-5,5); 
y1 = sin(x); 
subplot(2,2,1) 
plot(x,y1) 
title('First subplot') 

y2 = sin(2*x); 
subplot(2,2,2) 
plot(x,y2) 
title('Second subplot') 

y3 = sin(4*x); 
subplot(2,2,3) 
plot(x,y3) 
title('Third subplot') 

y4 = sin(6*x); 
subplot(2,2,4) 
plot(x,y4) 
title('Fourth subplot')