fgp824.m

matlab数字信号处理代码

%《数字信号处理教程——MATLAB释义与实现》第八章图8.2.4生成程序fgp824
% 切贝雪夫I型滤波器的幅特性函数
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
N=input('N= ');                 % 输入滤波器阶数
epsilon=0.2;                    % 设定epsilon
x1=[0:0.01:1];                  % 分段设定横坐标x1,x2
x2=[1:0.01:1.1];
x=[x1,x2];                      % 合成x
y1=cos(N*acos(x1));             % 分段计算切贝雪夫多项式的值y1,y2
y2=cosh(N*acosh(x2));
y=[y1,y2];                      % 合成y
M=1./(1+epsilon^2*y.^2);        % 求切贝雪夫滤波器的幅特性
plot(x,M)                       % 画切贝雪夫滤波器的幅特性图
set(gcf,'color','w')            % 置图形背景色为白