Сгенерируйте импульс с 1000 выборками и тональный сигнал с 1000 выборками с нормированной частотой $\pi /2$. Вычислите распределение Вигнера-Вилля суммы двух сигналов.

x = zeros(1001,1);
x(500) = 10;

y = sin(pi*(0:1000)/2)';

[d,f,t] = wvd(x+y);

Постройте график распределения Вигнер-Вилле.

imagesc(t,f,d)
axis xy
colorbar

Воспроизведите результат путем вызова wvd без выходных аргументов.

wvd(x+y)

Сгенерируйте сигнал, состоящий из синусоиды 200 Гц, дискретизированной на 1 кГц в течение 1,5 секунд.

fs = 1000;
t = (0:1/fs:1.5)';
x = cos(2*pi*t*200);

Вычислите распределение сигнала по Вигнеру-Вилле.

wvd(x,fs)

Добавьте к сигналу щебет, частота которого изменяется синусоидально между 250 Гц и 450 Гц. Преобразуйте сигнал в расписание MATLAB ®. Вычислите распределение Вигнера-Вилле.

x = x + vco(cos(2*pi*t),[250 450],fs);
xt = timetable(seconds(t),x);

wvd(xt)

Установите для нуля распределительных элементов с амплитудой менее 0.

wvd(xt,'MinThreshold',0)

Сгенерируйте сигнал, дискретизированный с частотой 1 кГц в течение 1 секунды. Одним из компонентов сигнала является щебет, который увеличивается в частоте квадратично с 100 Гц до 400 Гц во время измерения. Другой компонент сигнала является щебетом, которое уменьшается в частоте линейно с 350 Гц до 50 Гц за тот же промежуток времени.

Сохраните сигнал в расписании.

fs = 1000;
t = 0:1/fs:1;

x = chirp(t,100,1,400,'quadratic') + chirp(t,350,1,50);

Вычислите распределение сигнала по Вигнеру-Вилле.

wvd(x,fs)

Вычислите сглаженное псевдо распределение Wigner-Ville сигнала. Задайте 501 частотные и 502 моменты времени.

wvd(x,fs,'smoothedPseudo','NumFrequencyPoints',501,'NumTimePoints',502)

Увеличьте количество временных точек, чтобы квадратичный щебет стал видимым.

wvd(x,fs,'smoothedPseudo','NumFrequencyPoints',501,'NumTimePoints',522)

Увеличьте частотные и временные точки, чтобы получить более резкое изображение.

wvd(x,fs,'smoothedPseudo','NumFrequencyPoints',1000,'NumTimePoints',1502)

Сгенерируйте двухкомпонентный сигнал, дискретизированный на частоте 3 кГц в течение 1 секунды. Первый компонент является квадратичным щебетом, частота которого увеличивается с 300 Гц до 1300 Гц во время измерения. Второй компонент является щебетом с синусоидально изменяющимся содержимым частоты. Сигнал встроен в белый Гауссов шум. Выразите время между последовательными выборками как duration скаляр.

fs = 3000;
t = 0:1/fs:1-1/fs;
dt = seconds(t(2)-t(1));

x1 = chirp(t,300,t(end),1300,'quadratic');
x2 = exp(2j*pi*100*cos(2*pi*2*t));

x = x1 + x2 + randn(size(t))/10;

Вычислите и постройте график сглаженного псевдосигнала Wigner Ville сигнала. Отобразите распределение во времени с помощью 601-выборочного окна Хэмминга и по частоте с помощью 305-образного прямоугольного окна. Используйте 600 точек частоты для отображения. Установите, чтобы нуль эти компоненты распределения с амплитудой меньше $-50$.

wvd(x,dt,'smoothedPseudo',hamming(601),rectwin(305), ...
    'NumFrequencyPoints',600,'MinThreshold',-50)

Сгенерируйте сигнал, состоящий из четырех атомов Гауссова. Каждый атом состоит из синусоиды, модулированной Гауссовым. Синусоиды имеют частоты 100 Гц и 400 Гц. Гауссы центрированы со скоростью 150 миллисекунд и 350 миллисекунд и имеют отклонение $0.{01}^2$. Все атомы имеют единичную амплитуду. Дискретизация сигнала производится на частоте 1 кГц в течение полсекунды.

fs = 1000;
t = (0:1/fs:0.5)';

f1 = 100;
f2 = 400;

mu1 = 0.15;
mu2 = 0.35;

gaussFun = @(A,x,mu,f) exp(-(x-mu).^2/(2*0.01^2)).*sin(2*pi*f.*x)*A';

s = gaussFun([1 1 1 1],t,[mu1 mu1 mu2 mu2],[f1 f2 f1 f2]);

Вычислите и отобразите распределение сигнала по Вигнеру-Вилле. Интерференционные условия, которые могут иметь отрицательные значения, появляются на полпути между каждой парой автоматических членов.

wvd(s,fs)

Вычислите и отобразите сглаженное псевдо распределение Wigner-Ville сигнала. Сглаживание во времени и частоте ослабляет условия интерференции.

wvd(s,fs,'SmoothedPseudo')