myfreqre.m

matlab例子

function [H]=myfreqre(b,a,w)
%由差分方程求频率响应函数
%[H]=freqresp(b,a,w)
%H=在w频率点上求得的频率响应数组
%b=分子系数数组
%a=分母系数数组
%w=频率点位置向量
m=0:length(b)-1;l=0:length(a)-1;
num=b*exp(-j*m'*w);
den=a*exp(-j*l'*w);
H=num./den;