%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%用不同的方法繪制功率譜并比較%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%2017/11/13????????????????%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc;clear;close?all;
%%?定義信號及基礎變量
Fs?=?5000;????????????????????%采樣頻率5000Hz
T?=?1/Fs;?????????????????????%采樣周期
L?=?1000;?????????????????????%信號長度
t?=?（0:L-1）*T;????????????????%時間向量
xn0?=?10*sin（2*pi*100*t）;?????%原始信號?
xn??=?xn0?+?randn（size（t））;???%加噪聲信號
%%?信號時域圖
figure（1）
f1?=?subplot（211）;
plot（f1txn0）;
axis（[00.01-1515]）;
title（‘原始信號波形‘）;
f2?=?subplot（212）;
plot（f2txn）;
axis（[0?0.01?-15?15]）;
title（‘加高斯噪聲信號‘）;
%%?fft畫功率譜
figure（‘Numbertitle‘?‘off‘?‘Name‘?‘fft畫功率譜‘）;
subplot（211）;
Pxx=10*log10（（abs（fft（xn01024））/512）.^2）;
f?=?Fs/2*linspace（01513）;
figure（2）
plot（fPxx（1:513））;
ylim（[-100100]）;
title（‘原始信號功率譜‘）;
grid
subplot（212）;
Pxx=10*log10（（abs（fft（xn1024））/512）.^2）;
f?=?Fs/2*linspace（01513）;
plot（fPxx（1:513））;
ylim（[-100100]）;
grid
title（‘噪聲信號功率譜‘）;
%%?周期圖法與平均周期圖法
figure（‘Numbertitle‘?‘off‘?‘Name‘?‘周期圖法‘）;
Hrect?=?spectrum.periodogram（）;
subplot（311）;psd（Hrectxn‘Fs‘Fs‘NFFT‘1024）;
title（‘周期圖法\_—矩形窗‘）;
Hrect?=?spectrum.periodogram（‘hamming‘）;
subplot（312）;psd（Hrectxn‘Fs‘Fs‘NFFT‘1024）;
title（‘周期圖法-漢明窗‘）;
Hrect?=?spectrum.welch（‘rectangular‘10075）;
subplot（313）;psd（Hrectxn‘Fs‘Fs‘NFFT‘512）;
title（‘平均周期圖法‘）;

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

?????文件???????1559??2018-09-13?10:04??matlab各種功率譜畫法及比較\a.m

?????文件?????????55??2018-08-27?11:33??matlab各種功率譜畫法及比較\【源碼使用必讀】.url

?????目錄??????????0??2018-09-13?12:52??matlab各種功率譜畫法及比較

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

?????????????????1614????????????????????3