Làm thế nào tôi có thể vẽ tần số so với cường độ của biến đổi wavelet?

Tôi đang chạy Morlet biến đổi wavelet liên tục. Tôi đã nhận được wscalogramtín hiệu và bây giờ tôi muốn vẽ biểu đồ tần số như hình ảnh sau đây, nhưng tôi không biết làm thế nào để làm điều đó:

Tôi đã sử dụng scal2freq hàm MATLAB để chuyển đổi tỷ lệ thành tần số giả. Ngoài ra tôi có một số tần số trong tín hiệu của mình có tỷ lệ giảm chấn lớn (4%) nên chúng không thể nhìn thấy rõ trong âm mưu. Làm thế nào tôi có thể phóng đại các chế độ ẩm ướt này?

Tôi đang sử dụng MATLAB, đây là mã của tôi:

% 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 

Cốt truyện hiện tại của tôi trông như thế này:

Sửa đổi mã của bạn một chút, nhưng không có thay đổi lớn và tôi cũng nhận được kết quả chính xác. Sử dụng mã mẫu này ở đây và bạn sẽ không thấy bất kỳ vấn đề nào. Tôi nhận được kết quả chính xác.

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

Để có được hình ở phía bên phải của hình ảnh đầu tiên bạn đưa lên, đó chỉ đơn giản là phổ cắt thời gian. Nói cách khác, một cột trong ma trận tần số thời gian của bạn.

Ví dụ: chúng ta có thể nhận được ma trận hệ số bình phương sau, theo tần suất so với thời gian:

Sau đó, chúng ta hãy giả sử rằng chúng ta muốn xem phổ cho lát cắt thời gian ở cột 47. Sau đó, điều này được đưa ra ở đây:

Tại đây, bạn sẽ thấy sức mạnh của các hệ số, so với tần số, cho lát cắt thời gian tương ứng với cột 47.

Biến đổi wavelet liên tục (CWT) là một phương pháp để phân tích quy mô thời gian . Có, bạn đọc nó một cách chính xác, tỷ lệ , không phải tần số . Tuy nhiên, có thể ánh xạ tỷ lệ theo tần số, và thậm chí khá dễ dàng. Vì bạn là người dùng MATLAB, nên có thể bạn sẽ muốn sử dụng chức năng này , thực hiện như sau:

F = scal2frq (A, 'wname', DELTA) trả về tần số giả tương ứng với các thang đo được đưa ra bởi A, hàm sóng con 'wname' (xem wavefun để biết thêm thông tin) và thời gian lấy mẫu DELTA.

Như bạn có thể đọc, tần số không phải là tần số thực trong tín hiệu, mà là xấp xỉ. Nếu bạn muốn tần số thực, bạn nên sử dụng biến đổi Fourier thời gian ngắn (STFT).

Vậy, khi nào nên sử dụng CWT? Theo nguyên tắc chung, wavelet hoạt động tốt khi các tín hiệu là thời (tức là có sự thay đổi nhanh chóng).

EDIT: haha, tôi không biết tôi đang làm gì ở đây. xin lỗi :-D hoàn toàn bỏ lỡ điểm câu hỏi của bạn Chỉ cần đọc tiêu đề ....

Không chắc chắn nếu bạn vẫn cần trợ giúp ...

Nếu bạn muốn tần số so với cường độ, hãy sử dụng lệnh "lưới / lướt" để vẽ các hệ số.

Sau đó thay đổi hướng của cốt truyện để bạn có kết quả mà bạn tìm kiếm. Đối với điều này xem lệnh "xem".