幅频响应xlabel('频率/Hz');ylabel('振幅/dB');gridon;title('幅频图');subplot(2,1,2),plot(w/2/pi,unwrap(angle(H)));%绘制相频响应xlabel('频率/Hz');ylabel('相位/^o');gridon;title('相频图')%--------------------------------------------------------------------------figure(2)dt=1/10000;%模拟信号采样间隔f1=100;f2=1500;f3=2900;%输入信号的三个频率成分t=0:dt:0.04;%给定模拟时间段x=sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t)+0.5*sin(2*pi*f3*t);%输入信号h=[tf(bt,at)];%滤波器在MATLAB系统中的表示[y,t1]=lsim(h,x,t);%模拟输出subplot(2,1,1),plot(t,x),title('输入信号')%绘出输入信号subplot(2,1,2),plot(t1,y)%绘制输出信号title('输出信号'),xlabel('时间/s')%--------------------------------------------------------------------------figure(3)Fy=abs(fft(y));w=(-(length(Fy)-1)/2:(length(Fy)-1)/2)/(length(Fy)*dt);plot(w((length(Fy)-1)/2+1:(length(Fy)-1)),Fy(1:(length(Fy)-1)/2)*2/length(Fy))title('输出信号频谱')xlabel('频率/Hz')
全文预览
