信号与系统上机实验matlab第一第二次实验课.doc

实验一 典型连续时间信号描述及运算 1、单边指数衰减 t=0:1:500; f1=200*exp(-1.*t./62.5); f2=200*exp(-1.*t./125.0); f3=200*exp(-1.*t./250.0); plot(t,f1,-,t,f2,--,t,f3,-.) 2、正弦信号 t=-250:1:250; pi=3.1415926 f1=150*sin(2*pi*t/100); f2=150*sin(2*pi*t/200); f3=150*sin(2*pi*t/200+pi/5); plot(t,f1,-,t,f2,--,t,f3,-) 3、正弦衰减信号 t=0:1:500; f1=200*sin(((2*pi*t)/100)).*exp(-t/250); plot(t,f1) 4、钟形信号 t=-250:1:250; f1=400*exp((-t.*t)./(100*100)); f2=400*exp((-t.*t)./(150*150)); f3=400*exp((-t.*t)./(250*250)); plot(t,f1,-,t,f2,--,t,f3,-.) 5、符号信号 t=-5:0.001:5; f1=sign(t); plot(t,f1) 6、阶跃信号 t=-5:0.001:5; f1=0.5+0.5*sign(t); plot(t,f1) 二、连续时间信号的运算 已知及信号，用MATLAB绘出满足下列要求的信号波形。 （1） （2） （3） （4） （5） (1): syms t f1=sym((-t+4)*(heaviside(t)-heaviside(t-4))); subplot(1,2,1); ezplot(f1); y1=subs(f1,t,-t); f3=f1+y1; subplot(1,2,2); ezplot(f3); (2) syms t f1=sym((-t+4)*(heaviside(t)-heaviside(t-4))); subplot(1,2,1); ezplot(f1); y2=subs(f1,-2*t); f4=-(y2+f1); subplot(1,2,2); ezplot(f4); (3) syms t f2=sym(sin(2*pi*t)); subplot(1,3,1); ezplot(f2); subplot(1,3,2); ezplot(f3); f5=f2*f3; subplot(1,3,3); (4) syms t subplot(1,3,1); ezplot(f1); subplot(1,3,2); ezplot(f2); f6=f1*f2; subplot(1,3,3); ezplot(f6) (5) syms t y3=subs(f6,t-2); f7=y3*f2; subplot(1,1,1); ezplot(f7) 程序随堂测试：以f(t)为三角信号为例，求f(2t) , f(2-2t) t=-3:0.001:3; ft=tripuls(t,4,0.5); subplot(3,1,1); plot(t,ft); grid on; title (f(t)); ft1= tripuls(2*t,4,0.5); subplot(3,1,2); plot(t,ft1); grid on; title (f(2t)); ft2= tripuls(2-2*t,4,0.5); subplot(3,1,3); plot(t,ft2); grid on; 实验二 线性系统时域分析 、连续时间信号卷积 利用MATLAB实现下述两个信号的卷积积分： p=0.01; k1=-0.5:p:1; f1=1; k2=0:p:2; f2=0.5*k2; [f,k]=sconv(f1,f2,k1,k2,p) 二、线性系统时域分析 连续系统的冲激响应、阶跃响应及MATLAB实现 （3）实验内容 求解如下LTI系统的冲激响应、阶跃响应 a=[2 1 8]; b=[