完整的维纳滤波器Matlab源程序说课材料.pdfVIP

完 整 的 维 纳 滤 波 器 Ma t l a b 源 程 序 精品文档 clear;clc; %输入信号 A=1; %信号的幅值 f=1000; %信号的频率 fs=10^5; %采样频率 t=(0:999); %采样点 Mlag=100; %相关函数长度变量 x=A*cos(2*pi*f*t/fs); %输入正弦波信号 xmean=mean(x); %正弦波信号均值 xvar=var(x,1); %正弦波信号方差 xn=awgn(x,5); %给正弦波信号加入信噪比为 20dB 的高斯白噪声 figure(1) plot(t,xn) %绘制输入信号图像 title( 输入信号图像 ) xlabel(x 轴单位： t/s,color,b) ylabel(y 轴单位： f/HZ,color,b) 收集于网络，如有侵权请联系管理员删除 精品文档 xnmean=mean(xn) %计算输入信号均值 xnms=mean(xn.^2) %计算输入信号均方值 xnvar=var(xn,1) %计算输入信号方差 Rxn=xcorr(xn,Mlag,biased); %计 算输入信号自相关函数 figure(2) subplot(221) plot((- Mlag:Mlag),Rxn) %绘制自相关函数图像 title( 输入信号自相关函数图像 ) [f,xi]=ksdensity(xn); %计算输入信号的概率密度， f 为样本点 xi 处的概率密 度 subplot(222) plot(xi,f) %绘制概率密度图像 title( 输入信号概率密度图像 ) X=fft(xn); %计算输入信号序列的快速离散 傅里叶变换 Px=X.*conj(X)/600; %计算信号频谱 subplot(223) semilogy(t,Px) %绘制在半对数坐标系下频谱图像 title( 输入信号在半对数坐标系下频谱图像 ) 收集于网络，如有侵权请联系管理员删除 精品文档 xlabel(x 轴单位： w/rad,color,b) ylabel(y 轴单位： w/HZ,color,b) pxx=periodogram(xn); %计算输入信号的功率谱密度 subplot(224) semilogy(pxx) %绘制在半对数坐标系下功率谱密度图像 title( 输入信号在半对