getailijingxing.m

来自「利用平稳随机信号各态历经性」· M 代码 · 共 36 行

%随机信号各态历经性应用
%利用各态历经性获取信号的功率谱

clear all;
Fs=1000;                                 %采样频率
%产生含有噪声的序列
var=sqrt(1/exp(1.0));
n=0:1/Fs:1;
N=length(n);                             %为总采样点数1001
e=var*randn(1,N);
w0=100*pi;
w1=50*pi;
xn=exp(j*w0*n-j*pi)+exp(j*w1*n-j*0.7*pi)+e;
%绘制信号波形
subplot(2,1,1)
plot(n,abs(xn))
xlabel('n')
title('x(n)=exp(j*w0*n-j*pi)+exp(j*w1*n-j*0.7*pi)+e(n)')
%计算序列的DFT
nfft=1024;
Xk=fft(xn,nfft);
%计算序列的PSD
Pxx=abs(Xk).^2/N;
%绘制功率谱图形
%index=0:round(nfft/2-1);
%k=index*N/nfft;
k=(0:1023)/1024*Fs;
%plot_Pxx=10*log10(Pxx(index+1));
plot_Pxx=10*log10(Pxx);
subplot(2,1,2)
plot(k,plot_Pxx)
xlabel('频率(Hz)')
ylabel('公式直接计算的功率谱')