维纳滤波器

用matlab编写的维纳滤波器，实现了维纳滤波器。

clc;
clear;
pi=3.1415926;
N=1024;
Nfft=128;
f1s=0.2;
f2s=0.25;
for n=1:N
x1（n）=sin（2*pi*f1s*n）;
x2（n）=sin（2*pi*f2s*n）;
end
x=x1+x2;
y=awgn（x10）; 

%Levenson递推法数据长度N=1024，
%模型阶数分别为15、20、30时信号的功率谱
[Pxx1f]=pyulear（y15Nfft1）;
Pxx1=10*log10（Pxx1）;
[Pxx2f]=pyulear（y20Nfft1）; 
Pxx2=10*log10（Pxx2）;
[Pxx3f]=pyulear（y30Nfft1）;
Pxx3=10*log10（Pxx3）;
hold on
subplot（331）;plot（fPxx1‘m‘）;axis（[0 0.5 -15 15]）;
title（‘order=15N=1024‘）;grid on;
subplot（332）;plot（fPxx2‘m‘）;axis（[0 0.5 -15 15]）;
title（‘order=20N=1024‘）;grid on;
subplot（333）;plot（fPxx3‘m‘）;axis（[0 0.5 -15 15]）;
title（‘order=30N=1024‘）;grid on;

%Levenson递推法数据长度N=512，
%模型阶数分别为15、20、30时信号的功率谱
N=512;
for n=1:N
x1（n）=sin（2*pi*f1s*n）;
x2（n）=sin（2*pi*f2s*n）;
end
x=x1+x2;
y=awgn（x10）; 
[Pxx4f]=pyulear（y15Nfft1）;
Pxx4=10*log10（Pxx4）;
[Pxx5f]=pyulear（y20Nfft1）; 
Pxx5=10*log10（Pxx5）;
[Pxx6f]=pyulear（y30Nfft1）;
Pxx6=10*log10（Pxx6）;
hold on
subplot（334）;plot（fPxx4‘m‘）;axis（[0 0.5 -15 15]）;
title（‘order=15N=512‘）;grid on;
subplot（335）;plot（fPxx5‘m‘）;axis（[0 0.5 -15 15]）;
title（‘order=20N=512‘）;grid on;
subplot（336）;plot（fPxx6‘m‘）;axis（[0 0.5 -15 15]）;
title（‘order=30N=512‘）;grid on;

%Levenson递推法数据长度N=256，
%模型阶数分别为15、20、30时信号的功率谱
N=256;
for n=1:N
x1（n）=sin（2*pi*f1s*n）;
x2（n）=sin（2*pi*f2s*n）;
end
x=x1+x2;
y=awgn（x10）; 
[Pxx7f]=pyulear（y15Nfft1）;
Pxx7=10*log10（Pxx7）;
[Pxx8f]=pyulear（y20Nfft1）; 
Pxx8=10*log10（Pxx8）;
[Pxx9f]=pyulear（y30Nfft1）;
Pxx9=10*log10（Pxx9）;
hold on
subplot（337）;plot（fPxx7‘m‘）;axis（[0 0.5 -15 15]）;
title（‘order=15N=256‘）;grid on;
subplot（338）;plot（fPxx8‘m‘）;axis（[0 0.5 -15 15]）;
title（‘order=20N=256‘）;grid on;
subplot（339）;plot（fPxx9‘m‘）;axis（[0 0.5 -15 15]）;
title（‘order=30N=256‘）;grid on;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-08-29 16:12  维纳滤波器\
     文件          55  2018-08-27 11:33  维纳滤波器\【源码使用必读】.url
     目录           0  2018-08-29 16:07  维纳滤波器\维纳滤波器设计\
     文件        1905  2009-12-01 16:07  维纳滤波器\维纳滤波器设计\a.m
     文件        1868  2009-12-01 16:09  维纳滤波器\维纳滤波器设计\b.m
     文件      188416  2012-11-20 13:56  维纳滤波器\维纳滤波器设计\维纳滤波器设计.doc