RC一阶滤波器matlab滤波器节选

%---------------------------滤波器部分------------------------------------
%-----------------------滤波器本身特性（波特图）----------------------------
fc=15*10^3;%截止频率
w=2*pi*f;%角频率
Para=(1/(2*pi*fc))*1i;
for n=1:length(f)
   S(n)=abs(1/(1+Para*w(n)));%滤波器对于不同频率幅值衰减系数
   P(n)=angle(1/(1+Para*w(n)))*180/pi;%滤波器对于不同频率相移系数
end
figure(2);
subplot(2,1,1);
plot(f,S,'r');%幅值曲线
title('幅值衰减特性');
xlabel('频率/hz');
ylabel('增益倍数');
subplot(2,1,2);
plot(f,P,'blue');%相位曲线
title('相位特性');
xlabel('频率/hz');
ylabel('相位');
Func=tf(1,[1/(2*pi*fc),1]);%系统的传递函数
figure(3);
bode(Func);%系统的波特图
title('幅频特性');