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