信号与系统实验_2资料.ppt

信号与系统实验（二）;目 录;例：已知某线性时不变离散系统的差分方程为： 试用MATLAB命令绘出当激励信号为 时，该系统的零状态响应。;例：已知某线性时不变离散系统的差分方程为： 试用MATLAB命令绘出，该系统的单位取样响应。;例：已知某线性时不变离散系统的单位取样响应为： 试用MATLAB求当激励信号为 时系统的零状态响应。;;例：描述某线性时不变系统的方程为： 试求系统的冲激响应和阶跃响应。;连续系统的时域分析;利用卷积求连续系统响应;解：;周期信号的分解;周期信号的合成;二、非周期信号的傅里叶变换;非周期信号的频谱图;syms t w ft=sym((t+2)*Heaviside(t+2)-2*t*Heaviside(t)+(t-2)*Heaviside(t-2)); Fw=simple(fourier(ft)); ezplot(abs(Fw),[-4*pi,4*pi]),axis([-4*pi,4*pi,-0.5,4.5]);grid on;替换符号表达式函数： subs(f) f为符号表达式;fourier( )和ifourier( )函数的局限性;傅里叶变换的数值法求解;dt=0.01; t=-4:dt:4; ft=(t+2).*Heaviside(t+2)-2*t.*Heaviside(t)+(t-2).*Heaviside(t-2); N=2000; k=-N:N; w=pi*k/(N*dt); F=dt*ft*exp(-j*t*w); plot(w,abs(F)) axis([-4*pi,4*pi,-0.5,4.5]) xlabel(w),ylabel(Fw) title(amplitude spectrum),grid on;傅里叶变换的数值法求解;求系统频率响应函数： freqs(b,a,w) b和a分别为系统函数分子和分母多 项式的系数向量。 w为频率范围;连续时间系统的频率特性;解： syms w t f=sym(heaviside(t)-heaviside(t-1)); F=fourier(f);F=simple(F); H1=sym(1/(i*w+1));H2=sym(10/(i*w+10)); Yw1=simple(H1*F);Yw2=simple(H2*F); yt1=simple(ifourier(Yw1));yt2=simple(ifourier(Yw2)); figure(); Subplot(3,3,1),ezplot(f,[-1,5]);axis([-1,5,-0.5,1.5]) subplot(3,3,3),ezplot(abs(F),[-8*pi,8*pi]);axis([-8*pi,8*pi,-0.5,1.5]) subplot(3,3,4),ezplot(abs(H1),[-8*pi,8*pi]);axis([-8*pi,8*pi,-0.5,1.5]) subplot(3,3,7),ezplot(abs(H2),[-8*pi,8*pi]);axis([-8*pi,8*pi,-0.5,1.5]) subplot(3,3,5),ezplot(abs(Yw1),[-8*pi,8*pi]);axis([-8*pi,8*pi,-0.5,1.5]) subplot(3,3,8),ezplot(abs(Yw2),[-8*pi,8*pi]);axis([-8*pi,8*pi,-0.5,1.5]) subplot(3,3,6),ezplot(yt1,[-1,5]);axis([-1,5,-0.5,1.5]) subplot(3,3,9),ezplot(yt2,[-1,5]);axis([-1,5,-0.5,1.5]);低通滤波器的响应;连续时间系统的稳态响应;系统的稳态响应;第 6 章 复频域分析;二、连续时间系统的拉普拉斯变换分析法;求代数方程的根函数： roots(b) b为方程多项式系数向量;画零极点图函数：pzmap(b,a) b为分子多项式向量 a为分母多项式向量 pzmap(sys);第 7 章 抽样定理;Ts=1; dt=0.01;t1=-4:dt:4; ft=((1+cos(t1))/2).*[Heaviside(t1+pi)-Heaviside(t1-pi)]; subplot(221),plot(t1,ft),grid on,axis([-4 4 -0.1 1.1]) xlabel(‘Time(sec)’),ylabel(‘f(t)’),title(‘升余弦脉冲信号) N=500;k=-N:N