ch3example1a.m

来自「清华大学出版社 邵玉斌编写的《通信系统建模与仿真实例分析》一书的所有MATLAB」· M 代码 · 共 16 行

% ch3example1A.m
clear;
f_p=2400; f_s=5000; R_p=3; R_s=25;     % 设计要求指标
[n, fn]=buttord(f_p,f_s,R_p,R_s, 's'); % 计算阶数和截止频率
Wn=2*pi*fn;                            % 转换为角频率
[b,a]=butter(n, Wn, 's');              % 计算H(s)
f=0:100:10000;                         % 计算频率点和频率范围
s=j*2*pi*f;                            % s=jw=j*2*pi*f
H_s=polyval(b,s)./polyval(a,s);        % 计算相应频率点处H(s)的值
figure(1);
subplot(2,1,1); plot(f, 20*log10(abs(H_s))); % 幅频特性
axis([0 10000 -40 1]);
xlabel('频率 Hz');ylabel('幅度 dB');
subplot(2,1,2); plot(f, angle(H_s));         % 相频特性
xlabel('频率 Hz');ylabel('相角 rad');
figure(2); freqs(b,a);  % 也可用指令freqs直接画出H(s)的频率响应曲线。