实验二-信号的运算

1、信号与系统实验二实验二 信号的运算、信号 f (t) (t 1) (t 1) (t 1) 。请编写程序实现些列运算并画出波形。1、f(2t) f(t) (t)解：>>t=li nspace(-5,5,1001);% 时间 t 的取值范围 自定>>y=(t/2+1).*stepfu n(t/2,-0.5)-(t/2+1).*stepfu n(t/2,0.5)+(t+1).*stepfu n( t,-1).*stepfu n( t,0)-(t+1).*stepfun(t,1).*stepfun(t,O);%式中的 stepfun 为阶跃函数2、f (t 2) (t 2)解：

2、>> t=linspace(-6,6,1001);%时间 t 的取值范围>> y=(t-1).*stepfun(t,1).*stepfun(t,2)-(t-1).*stepfun(t,3).*stepfun(t,2);%式中 stepfun 为阶跃函数>> plot(t,y);xlabel('t');ylabel('y');3、f (t 2) (t)解：>> t=linspace(-6,6,1001);%时间 t 的取值范围>> y=(t-1).*stepfun(t,1).*stepfun(t,0)-(

3、t-1).*stepfun(t,3).*stepfun(t,0);%式中的 stepfun 为阶跃函数>> plot(t,y);xlabel('t');ylabel('y');4、f(2 t)解：>> t=linspace(-10,10,1001);%时间 t 的取值范围>> y=(3-t)*stepfun(-t,-3)-(3-t)*stepfun(-t,-1);%式中的 stepfun 为阶跃函数>> plot(t,y);xlabel('t');ylabel('y');叮U.S0

4、4112Q0i-O.B-10 电 -6-4-? o5、f( 2 t) ( t)%stepfun为阶跃函数解：>> t=linspace(-10,10,1001);%时间 t 的取值范围>> y=(-t-1).*stepfu n(-t,1).*stepfu n(-t,0)-(-t-1).*stepfu n(-t,3).*stepfu n(-t,0);>> plot(t,y);xlabel('t');ylabel('y');二、编写一个利用 MATLAB提供的数值积分计算函数计算两个连续函数的卷积积分的函数 程序，用其计算以下函数

5、在 010内的卷积积分的结果，并与理论计算出的结果相比拟，观 测其计算结果e (t 1) (t 1)1. fi(t)( t 1) 2 (t 1),f2(t)解：>>dt=0.01;>> k1=0:dt:10;>> f仁stepfu n(-k1,-1)+2*stepfu n(k1,1);>> k2=k1;>> f2=exp(-(k2+1).*stepfu n(k2,1);>> f=dt*co nv(f1,f2);>> k0=k1(1)+k2(1);>> k3=le ngth(f1)+le ngth(f

6、2)-2;>>k=k0:dt:k0+k3*dt;>> subplot(2,2,1);>> plot(k1,f1);title('f1(t)');xlabel('t');>> subplot(2,2,2);>> plot(k2,f2);title('f2(t)');xlabel('t');>> subplot(2,2,3);>> plot(k,f);>> h=get(gca,'positi on');>> h(

7、3)=2.5*h(3);set(gca,'positio n',h);>> title('f(t)=f1(t)*f2(t)');xlabel('t');%生成信号f1%生成信号f2%计算卷积结果f%计算序列f非零样值的起点位置%计算卷积和f的非零样值的宽度%确定卷积和f的非零样值的时间向量%在图1绘f2(t)时域波形图%在图2绘f2(t)时域波形图%画卷积f(t)的时域波形%将第三个示图横坐标范围扩至原来的倍fl (I)2. fi(t)sint (t), f2(t)(t 1)解：>> dt=0.01;>>k1=

8、0:dt:10;>>f1=si n( k1).*stepfu n(k1,0);>>k2=k1;>> f2=stepfu n(k2,-1);>>f=dt*co nv(f1,f2);>>k0=k1(1)+k2(1);>>k3=le ngth(f1)+le ngth(f2)-2;>>k=k0:dt:k0+k3*dt;>> subplot(2,2,1);>>plot(k1,f1);title('f1(t)');xlabel('t');>>subplot

9、(2,2,2);>>plot(k2,f2);title('f2(t)');xlabel('t');>>subplot(2,2,3);>> plot(k,f);%画卷积f(t)的时域波形>>h=get(gca,'positi on');>>h(3)=2.5*h(3);set(gca,'positio n',h);>>title('f(t)=f1(t)*f2(t)');xlabel('t');%生成信号f1%生成信号f2%计算卷积结

10、果f%计算序列f非零样值的起点位置%计算卷积和f的非零样值的宽度%确定卷积和f的非零样值的时间向量%在图1绘f2(t)时域波形图%在图2绘f2(t)时域波形图%将第三个示图横坐标范围扩至原来的倍1.51C.5W°0%生成信号f1%生成信号f2%计算卷积结果f%计算序列f非零样值的起点位置%计算卷积和f的非零样值的宽度%确定卷积和f的非零样值的时间向量%在图1绘f2(t)时域波形图%在图2绘f2(t)时域波形图%画卷积f(t)的时域波形3. fi(t) 2 (t) (t 1),f2(t) sin( t) (t) (t 1)解：>> dt=0.01;>> k1=0

11、:dt:10;>> f1=2*stepfu n( k1,0)-stepfu n(k1,1);>> k2=k1;>> f2=si n( pi*k2).*(stepfu n( k2,0)-stepfu n( k2,1);>> f=dt*co nv(f1,f2);>>k0=k1(1)+k2(1);>>k3=le ngth(f1)+le ngth(f2)-2;>> k=k0:dt:k0+k3*dt;>>subplot(2,2,1);>> plot(k1,f1);title('f1(t)&

12、#39;);xlabel('t');>>subplot(2,2,2);>> plot(k2,f2);title('f2(t)');xlabel('t');>>subplot(2,2,3);>> plot(k,f);>>h=get(gca,'positi on');>> h(3)=2.5*h(3);set(gca,'positio n',h);>>title('f(t)=f1(t)*f2(t)');xlabel('t');fi(t)co)