维纳滤波法matlab代码,完整的维纳滤波器Matlab源程序

维纳滤波器 fir滤波器程序

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) %输入正弦波信号 正弦波信号均值 正弦波信号方给正弦波信号加绘制输入计算输入信号均值 计算输入信号均方值 计算输入信号方计算输入信号自相关函绘制自相关函数图计算输入信号绘制概计算输入计算信号频谱 绘制在半对数