exa1_45.m

来自「新版Matlab 7辅助信号处理技术与应用」· M 代码 · 共 22 行

% exa1-45_transfer_function_reponse.m for example 1-45
% transfer function reponse

clear all;
fs=1000; %sample frequency unit in Hz
b=[0.3 0.2]; %sytem parameter
a=[1 -0.4 -0.7]; %sytem parameter
[h, f]=freqz(b, a, 256, fs);
mag=abs(h); %compute magnitude
ph=angle(h); %compute phase
ph=ph*180/pi; %unit in degrees
subplot(211)
plot(f,mag)
grid
xlabel('frequency(Hz)');
ylabel('magnitude');
subplot(212)
plot(f,ph)
grid
xlabel('frequency(Hz)');
ylabel('phase');