Як я можу побудувати графік частоти проти величини вейвлет-перетворення?

Я веду безперервне вейвлетське перетворення Morlet. У мене є wscalogramсигнал, і тепер я хочу накреслити частоту частоти, як на наступному малюнку. Але я не знаю, як це зробити:

Я використовував scal2freqфункцію MATLAB для перетворення масштабів у псевдочастоти. Крім того, у мене є деякі частоти в моєму сигналі, які мають велике коефіцієнт демпфування (4%), тому вони не добре помітні в сюжеті. Як я можу перебільшити ці режими з сильним затуханням?

Я використовую MATLAB, ось мої коди:

% Import the text4.txt to matlab workspace. and save it under name "data"
t=linspace(0,30,301);
Fs=ceil(inv(t(2)-t(1)));
x=data(:,4); % use x=data(:,3),x=data(:,5) too. first column is time,second is refrence
wname = 'morl';
scales = 1:1:256;
coefs = cwt(x,scales,wname,'lvlabs');
freq = scal2frq(scales,wname,1/Fs);
surf(t,freq,abs(coefs));shading('interp');
axis tight; xlabel('Seconds'); ylabel('Pseudo-Frequency (Hz)');
axis([0 30 0 1 0 60])
xlabel('Time'); ylabel('Frequency');
figure;
sc=wscalogram('image',coefs,'scales',freq,'ydata',x);
hold on 

Мій поточний сюжет виглядає приблизно так:

Трохи змінити свій код, але серйозних змін не буде, і в будь-якому випадку я отримую правильні результати. Використовуйте цей код шаблону тут, і ви не повинні бачити жодних проблем. Я отримую правильні результати.

clear all; 
t=linspace(0,30,301); 
Fs = (inv(t(2)-t(1))); 
x=randn(100,1);  
wname = 'morl'; 
scales = 1:1:256; 
chefs = cwt(x,scales,wname,'lvlabs');

freq = scal2frq(scales,wname,1/Fs);

figure; 
coefsSquared = abs(coefs).^2; 
imagesc(coefsSquared); 
grid off;

%Pick one of the columsn to plot: 
figure; 
plot(coefsSquared(:,47))

Щоб отримати фігуру з правого боку першого зображення, яке ви поставили, це просто спектр часових відрізків. Іншими словами, один стовпець вашої матриці частота-час.

Наприклад, ми можемо отримати таку матрицю квадратичного коефіцієнта за частотою та часом:

Припустимо, тоді ми хочемо подивитися спектр часового відрізка в стовпці 47. Тоді це дано тут:

Тут ви бачите потужність коефіцієнтів проти частоти для часового відрізка, що відповідає стовпцю 47.

Безперервне перетворення вейвлетів (CWT) - метод для аналізу масштабу в часі . Так, ви читаєте це правильно, масштаб , а не частота . Однак можна віднести масштаби до частот, і навіть досить легко. Оскільки ви користувач MATLAB, ви, ймовірно, захочете використовувати цю функцію , яка виконує такі дії:

F = scal2frq (A, 'wname', DELTA) повертає псевдочастоти, що відповідають масштабам, заданим A, вейвлет-функція 'wname' (див. Хвильову функцію для отримання додаткової інформації) та період вибірки DELTA.

Як ви можете прочитати, частоти - це не реальні частоти в сигналі, а наближення. Якщо ви хочете реальної частоти, вам слід скористатися короткочасним перетворенням Фур'є (STFT).

Отже, коли слід використовувати CWT? Як загальний принцип, вейвлети добре працюють, коли сигнали є тимчасовими (тобто мають швидкі зміни).

EDIT: ха-ха, я не знаю, що я тут роблю. вибачте :-D повністю пропустив пункт вашого питання. Просто прочитайте заголовок….

Не впевнений, чи потрібна допомога ...

Якщо ви хочете частоту проти величини, скористайтеся командою "mesh / surf", щоб побудувати коефіцієнти.

Потім змініть орієнтацію сюжету, щоб у вас був результат, якого ви прагнете. Для цього див. Команду "view".