zyulewalk.m

来自「用MATLAB实现数字信号处理的程序」· M 代码 · 共 9 行

%Yulewalk函数：设计递归数字滤波器
m=[0 0 1 1 0 0 1 1 0 0];
f=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 1];
[b,a]=yulewalk(10,f,m);
[H,W]=freqz(b,a,128);
plot(f,m,'r:',W/pi,abs(H),'b-');grid;
xlabel('频率/Hz'),ylabel('幅度');
legend('期望曲线','实际曲线');