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