资源描述:
《上机实验一 连续时间信号的时域分析》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、上机实验一连续时间信号的时域分析一、实验目的1、掌握连续时间信号时域运算的基本方法2、掌握相关函数的调用格式及应用3、掌握连续信号的基本运算4、掌握利用计算机进行卷积运算的原理和方法5、熟悉连续信号卷积运算函数conv的应用二、实验内容1、验证性实验(1)相加实现两个连续信号的相加Matlab程序:clearall;>>t=0:0.0001:3;>>b=3;>>t0=1;>>u=stepfun(t,t0);>>n=length(t);>>fori=1:n;u(i)=b*u(i)*(t(i)-t0);end>>y=sin(2*pi*t
2、);>>f=y+u;>>plot(t,f);>>xlabel('时间(t)');ylabel('幅值f(t)');title('连续信号的相加');(2)相乘实现两个连续信号的相乘Matlab程序:>>clearall;>>t=0:0.0001:5;>>b=3;>>t0=1;u=stepfun(t,t0);>>n=length(t);>>fori=1:n;u(i)=b*u(i)*(t(i)-t0);end>>y=sin(2*pi*t);>>f=y.*u;>>plot(t,f);>>xlabel('时间(t)');ylabel('幅值
3、f(t)');title('连续信号的相乘');(3)移位实现连续信号的移位Matlab程序:>>clearall;>>t=0:0.0001:2;>>y=sin(2*pi*t);>>y1=sin(2*pi*(t-0.2));>>plot(t,y,'-',t,y1,'--');>>ylabel('f(t)');xlabel('t');title('信号的移位');(1)翻转信号的翻转就是将信号的波形以纵轴为对称轴翻转180。将信号f(t)中的自变量t替换为-t即可得到其翻转信号。Matlab程序:>>clearall;>>t=0:0.
4、02:1;>>t1=-1:0.02:0;>>g1=3*t;g2=3*(-t1);>>gridon;>>plot(t,g1,'--',t1,g2);>>xlabel('t');ylabel('g(t)');title('信号的反折');2、程序设计实验(1)已知信号f1(t)=(-t+4)[U(t)-U(t-4)],f2(t)=sin(2πt),用Matlab绘出信号的时域波形。a.f3(t)=f1(-t)+f1(t)b.f4(t)=-[f1(-t)+f1(t)]c.f5(t)=f2(t)*f3(t)d.f6(t)=f1(t)*f2(
5、t)Matlab程序:symstf1=sym('(-t+4)*(heaviside(t)-heaviside(t-4))');f2=sin(2*pi*t);f3=f1+subs(f1,t,-t);f4=-f3;f5=f2*f3;f6=f1*f2;subplot(3,2,1);ezplot(f1,[-1,5]);title('f1(t)');subplot(3,2,2);ezplot(f2,[-5,5]);title('f2(t)');subplot(3,2,3);ezplot(f3,[-5,5]);title('f3(t)');su
6、bplot(3,2,4);ezplot(f4,[-5,5]);title('f4(t)');subplot(3,2,5);ezplot(f5,[-5,5]);title('f5(t)');subplot(3,2,6);ezplot(f6,[-5,5]);title('f6(t)');(3)证明卷积满足如下结论:a.f1(t)*f2(t)=f2(t)*f1(t);b.f1(t)*[f2(t)+f3(t)]=f1(t)*f2(t)+f1(t)*f3(t);Matlab程序:>>clearall;a=1000;t1=-5:1/a:5;f1
7、=stepfun(t1,0);f2=stepfun(t1,-1/a)-stepfun(t1,1/a);f3=stepfun(t1,0)-stepfun(t1,4);y1=conv(f1,f2);y2=conv(f2,f1);g1=conv(f1,(f2+f3));g2=conv(f1,f2)+conv(f1,f3);r=2*length(t1)-1;t=-10:1/a:10;subplot(2,2,1);plot(t,y1);axis([-5,5,0,1.2]);title('y1(t)');subplot(2,2,2);plot(
8、t,y2);axis([-5,5,0,1.2]);title('y2(t)');subplot(2,2,3);plot(t,g1);axis([-5,5,0,1.2]);title('g1(t)');subplot(2,2,4);pl