dwt2_chirp.m

来自「数字信号处理的视频分析中短时傅立叶变化STFT和连续小波变化CWT和离散小波变换」· M 代码 · 共 61 行

chirp;

x = chirp';
sampling_rate = 1;

clf;
subplot(4,2,1)
plot(x);
xlabel('Time (seconds)');
title('Oringinal Signal');
axis([0 256 -1 1]);

subplot(4,2,2)
scale_level =  1;
y = DWT2(x, sampling_rate, scale_level);
xlabel('Time (seconds)');
ylabel('Frequency (Hz)');
title('Scale Level = 1');

subplot(4,2,3)
scale_level =  2;
y = DWT2(x, sampling_rate, scale_level);
xlabel('Time (seconds)');
ylabel('Frequency (Hz)');
title('Scale Level = 2');

subplot(4,2,4)
scale_level =  3;
y = DWT2(x, sampling_rate, scale_level);
xlabel('Time (seconds)');
ylabel('Frequency (Hz)');
title('Scale Level = 3');

subplot(4,2,5)
scale_level =  4;
y = DWT2(x, sampling_rate, scale_level);
xlabel('Time (seconds)');
ylabel('Frequency (Hz)');
title('Scale Level = 4');

subplot(4,2,6)
scale_level =  5;
y = DWT2(x, sampling_rate, scale_level);
xlabel('Time (seconds)');
ylabel('Frequency (Hz)');
title('Scale Level = 5');

subplot(4,2,7)
scale_level =  6;
y = DWT2(x, sampling_rate, scale_level);
xlabel('Time (seconds)');
ylabel('Frequency (Hz)');
title('Scale Level = 6');

subplot(4,2,8)
scale_level =  7;
y = DWT2(x, sampling_rate, scale_level);
xlabel('Time (seconds)');
ylabel('Frequency (Hz)');
title('Scale Level = 7');