用MATLAB求幅频,相频特性

用MATLAB求幅频，相频特性

%对于有理分式H(jw),MATLAB提供freqs函数处理方法。调用格式为H=freqs(b,a,w)

%其中，b为分子多项式的系数，a为分母多项式系数，w为需计算的频率特性

%函数的取样点数。程序如下(低通滤波器为例)：

clear all

clc

w=linspace(0,5,200);

b=[1];

a=[1 2 2 1];

H=freqs(b,a,w);

subplot(2,1,1);

plot(w,abs(H));

set(gca,'xtick',[0 1 2 3 4 5]);

set(gca,'ytick',[0 0.4 0.707 1]);grid

xlabel('\omega');

ylabel('|H(j\omega|');

subplot(2,1,2);

plot(w,angle(H));

set(gca,'xtick',[0 1 2 3 4 5]);grid

xlabel('\omega');

ylabel('\phi(\omega)');

用MATLAB求幅频，相频特性

%对于有理分式H(jw),MATLAB提供freqs函数处理方法。调用格式为H=freqs(b,a,w)

%其中，b为分子多项式的系数，a为分母多项式系数，w为需计算的频率特性

%函数的取样点数。程序如下(低通滤波器为例)：

clear all

clc

w=linspace(0,5,200);

b=[1];

a=[1 2 2 1];

H=freqs(b,a,w);

subplot(2,1,1);

plot(w,abs(H));

set(gca,'xtick',[0 1 2 3 4 5]);

set(gca,'ytick',[0 0.4 0.707 1]);grid

xlabel('\omega');

ylabel('|H(j\omega|');

subplot(2,1,2);

plot(w,angle(H));

set(gca,'xtick',[0 1 2 3 4 5]);grid

xlabel('\omega');

ylabel('\phi(\omega)');