完整的维纳滤波器Matlab源程序学术2009-11-20 20:31:58 阅读308 评论0 字号:大中小
完整的维纳滤波器Matlab源程序
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) %绘制在半对数坐标系下频