四川大学信号与系统实验2(题目二)

信号与系统第二次实验-1-信号与系统第二次实验（题目二）2、给定一个连续LTI系统，其微分方程为：1101000100100ytytytxtxt输入信号：123sin2sin2sin2xtftftft(其中1f=0.06Hz，2f=5Hz，3f=32Hz)理论计算：1)计算系统的幅度响应，判断该系统是哪一类频率选择性滤波器？（低通、高通、带通、带阻）该系统是带通滤波器。）（）（）（）（)100)(10(1100)100j)(10j()1j(100j)100j)(10j()1j(1001000j100100j100jjj2222HXYH2)求系统的输出信号。3-10)]28.5664(sin919.2)59.110(sin019.9)71.1712.0(sin064.1[)(21)(ytttdejYttj仿真分析：程序如下：ts=0;%起始时刻te=50;%终止时刻fs=100;%采样频率t=ts:1/fs:te;N=1000;x=sin(2*pi*0.06*t)+sin(2*pi*5*t)+sin(2*pi*32*t);%信号描述X=fft(x,N)/N;fx=(-(N-1)/2:(N-1)/2)/N*fs;%建立频率轴num=[11];%分子系数den=[11101000];%分母系数sys=tf(num,den);%建立系统syssubplot(334);plot(t,x);xlabel('时间');ylabel('时域');title('x(t)时域波形');mag=abs(fftshift(X));theta=angle(fftshift(X))/pi*180;%%%%%%修正幅度谱%%%%%%delta_1=[1,2*ones(1,length(fftshift(X))-1)];mag=mag.*delta_1;%修正后的幅度谱%%%%%%修正相位谱%%%%%%delta_2=(mag0.1);%将幅值为0的频率分量的相位置为0”theta=theta.*delta_2;%修正后的相位谱subplot(335);plot(fx,mag);xlabel('频率/Hz');ylabel('幅值');title('x(t)幅频谱');subplot(336);bar(fx,theta,3);xlabel('频率/Hz');ylabel('相位');title('x(t)相频谱');y=lsim(sys,x,t);%将输入信号x加在信号与系统第二次实验-2-系统sys上，得到输出信号ysubplot(337);plot(t,y);xlabel('时间');ylabel('时域');title('y(t)时域波形');Y=fft(y,N)/N;mag=abs(fftshift(Y))';theta=(angle(fftshift(Y))/pi*180)';%%%%%%修正幅度谱%%%%%%delta_1=[1,2*ones(1,length(fftshift(Y))-1)];mag=mag.*delta_1;%修正后的幅度谱%%%%%%修正相位谱%%%%%%delta_2=(mag0.0005);%将幅值为0的频率分量的相位置为0”theta=theta.*delta_2;%修正后的相位谱subplot(338);plot(fx,mag);xlabel('频率/Hz');ylabel('幅值');title('y(t)幅频谱');subplot(339);bar(fx,theta,3);xlabel('频率/Hz');ylabel('相位');title('y(t)相频谱');%%%%%画bode图%%%%%w=0:0.1:10000;W=logspace(-1,2,100);[mag,pha]=bode(num,den,w);magdB=20*log10(mag);set(gcf,'color','w')%设置背景颜色为白色subplot(321);semilogx(w,magdB);gridon;xlabel('\omega');ylabel('20log_1_0|H(j\omega)|');title('伯德图（a）模');subplot(322);semilogx(w,pha);gridon;set(gca,'yticklabel',{'-\pi/2','-\pi/4','0','\pi/4','\pi/2'})xlabel('\omega');ylabel('argH(j\omega)');title('伯德图（b）相位');1)绘制系统的频率响应Bode图；频率响应Bode图如下：信号与系统第二次实验-3-2)绘制输入、输出信号的时域波形和频谱，分析和解释滤波过程；输入信号x(t)时域波形及频谱图如下：输出信号y(t)时域波形及频谱图如下：信号与系统第二次实验-4-分析和解释滤波过程：根据计算和bode图可以看出可看做一个带通滤波器无明显衰减的幅值有明显衰减，、系统后号的叠加，当信号输入输入信号是三个频率信递减）（时递增）（时基本不变）（时，与总之，十倍频程下降这一因式，产生的转折频率对应于述上升抵消其产生的频程下降与上这一因式的转折频率对应于十倍频程的上升开始产生这一因式，由的转折频率对应于）（z5z32z06.0j,100j,101j1001010/dB20)100/j1/(1100,)10/j1/(110/dB201)j1(1)j1)(100/j11)(10/j11)(101(j231HfHfHfHHHH3)对输入信号的采样时长改变时，频谱怎样变化？为什么？采样频率为50、200时频谱如下：由图可看出，采样时长变大，采样点减少，分散度变大；采样时间变小，采样点增多，分散度变小。采样时长改变，频谱幅值不变。但由于采样频率为50时未大于信号中最高频率32的2倍，出现了一个相位倒置。原因：采样时间的长短影响采样点的多少，从而影响频谱分散程度，但是不会影响频谱幅值。