matlab中频响函数

自己编写的matlab里面的freqz函数

function [ H ] = my_freqz（ baNfs ）
%UNtitleD2 Summary of this function goes here
%   Detailed explanation goes he
%%计算滤波器系数的长度
num = length（b）;
%%定义横坐标的数组
w = 0:2*pi/N:2*pi;
%%数字变量进行离散化
z = exp （ j * w）;
X = zeros（size（z））;
Y = zeros（size（z））;
%%计算分子、分母的数组
for i = 1: num
    X  = X + a（i） * z.^（-（i-1））;
    Y = Y + b（i） * z.^（-（i-1））;
end
%%注意这里都是用的点运算
H = Y ./ X;
w = w / pi * fs;
%%画出图像来
plot（w（1:N/2）20*log10（abs（H（1:N/2））））;

end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         518  2018-06-28 11:34  my_freqz.m