matlab各种功率谱画法及比较

matlab各种功率谱画法及比较 matlab各种功率谱画法及比较

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%用不同的方法绘制功率谱并比较%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%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