load?sumsin;
x=sumsin;

figure（1）;
%?subplot（611）;
plot（x）;
ylabel（‘x‘）;
title（‘原始信號‘）
%?title（‘原始信號以及各層近似信號‘）;

figure（2）;
%?對信號進行傅里葉變換
f=fft（x2048）;
f=abs（f）;????
subplot（211）;
plot（x）;
subplot（212）;
plot（f）;
title（‘原始信號以及傅里葉變換信號‘）;

figure（3）;
%?title（‘各層近似信號‘）
%?使用db3小波進行5層分解
[cl]=wavedec（x5‘db3‘）;
for?i=1:5
%?對分解的第5層到第1層的低頻系數分別進行重構
????a=wrcoef（‘a‘cl‘db3‘6-i）;
????subplot（51i）;
????plot（a）;?
????title（[‘第‘num2str（i）‘層近似信號‘]）
????ylabel（[‘a‘num2str（5-i）]）;
end

figure（4）;
subplot（611）
plot（x）;
title（‘原始信號‘）
ylabel（‘x‘）
for?i=1:5
%?對分解的第5層到第1層的高頻系數進行重構
????d=wrcoef（‘d‘cl‘db3‘6-i）;
????subplot（61i+1）;
????plot（d）;
????title（[‘第‘num2str（i）‘層高頻系數重構信號‘]）
????ylabel（[‘d‘num2str（6-i）]）;
end

%?使用小波分析一個由三個不同頻率的正弦信號疊加的信號，
%?看是否能將這三個正弦信號區分開來，
%?結果證明小波分析可以很好地識別某一頻率區間的信號。

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????435??2020-05-13?16:55??temp\readme.txt

?????文件???????1025??2020-05-13?16:50??temp\wavelet_final.m

?????目錄??????????0??2020-05-13?20:27??temp

-----------?---------??----------?-----??----

?????????????????1460????????????????????3