1、1实验一 连续时间信号分析一、实验目的(一)掌握使用 Matlab 表示连续时间信号1、学会运用 Matlab 表示常用连续时间信号的方法2、观察并熟悉常用信号的波形和特性(二)掌握使用 Matlab 进行连续时间信号的相关运算1、学会运用 Matlab 进行连续时间信号的时移、反褶和尺度变换2、学会运用 Matlab 进行连续时间信号微分、积分运算3、学会运用 Matlab 进行连续时间信号相加、相乘运算4、学会运用 Matlab 进行连续时间信号卷积运算二、实验条件一台电脑、winXP 系统、matlab7.0 软件三、实验内容1、利用 Matlab 命令画出下列连续信号的波形图。(1)
2、)4/3t(2cos代码:clear all;close all;clc;K=2;a=3;t=0:0.01:3;ft=K*cos(a*t+pi/4);plot(t,ft),grid onaxis(-5,5,-2.2,2.2)title(2cos(3t+4)2-5 -4 -3 -2 -1 0 1 2 3 4 5-2-1.5-1-0.500.511.522cos(3t+4)(2) )t(ue(-3 -2 -1 0 1 2 3-3-2-10123 与 与 与 与 与 与 与 与 与 与 与 与3代码:函数文件:function f=uCT(t)f=(t=0);命令文件:clear all;close
3、 all;clc;a=-1;t=-5:0.01:5;ft=(2-exp(a*t).*uCT(t);%y=2-exp(a*t);%plot(t,y),grid onplot(t,ft),grid onaxis(-3,3,-3,3);title(指数信号与阶跃信号的乘积)(3) )2()ut(cos1t4-3 -2 -1 0 1 2 3-3-2-10123 与 与 与 与 与 与 与 与 与 与 与 与代码:函数文件:function f=uCT(t)f=(t=0);代码文件:clear all;close all;clc;t=-5:0.01:5;ft=(uCT(t)-uCT(t-2).*(1+c
4、os(pi*t);%(1+cos(pi*t)*plot(t,ft),grid onaxis(-3,3,-3,3);title(余弦信号与阶跃信号的乘积)52、利用 Matlab 命令画出复信号 的实部、虚部、模和辐角。)4/t(j2e)t(f解:图像:-5 -4 -3 -2 -1 0 1 2 3 4 5-5-4-3-2-1012345 f(t)=2exp(j*(t+/4)代码:clear all; close all;clc;a=1;t=-10:0.01:10;ft=a*exp(j*(t+pi/4);r=real(ft);im=imag(ft);ab=abs(ft);6an=angle(ft)
5、;plot(t,r,r,t,im,b,t,ab,g,t,an,m),grid onaxis(-5,5,-5,5);title(f(t)=2exp(j*(t+/4)3、已知信号的波形(课本 P13 例题 1-1) ,画出 的2332tftftftf ,波形图。解:-2 0 2 4-1012 f1(t)-2 0 2 4-1012 f2(t)-2 0 2 4-1012 f3(t)-2 0 2 4-1012 f4(t)代码:clear all;close all;clc;t=-5:0.01:5;y1=funct1(t-2);y2=funct1(3*t);y3=funct1(-t);y4=funct1(
6、-3*t-2);subplot(2,2,1);plot(t,y1);7axis(-3,5,-1,2);title(f1(t),grid on;subplot(2,2,2);plot(t,y2);axis(-3,5,-1,2);title(f2(t),grid on;subplot(2,2,3);plot(t,y3);axis(-3,5,-1,2);title(f3(t),grid on;subplot(2,2,4);plot(t,y4);axis(-3,5,-1,2);title(f4(t),grid on;4、使用微分命令求 关于变量 x 的一阶导数;使用积分命令计算不定积xsinly分 ,
7、定积分 。dxax25 dex102解:代码:clear all;close all;clc;syms x;y=x*sin(x)*log(x);dy=diff(y,x,1)z=x5-2*x*x+sqrt(x)/2;lz=int(z,x)w=x*exp(x)/(1+x)2;lw=int(w,x,0,1)dy=sin(x)*log(x)+x*cos(x)*log(x)+sin(x);lz=1/6*x6-2/3*x3+1/3*x(3/2)lw=1/2*exp(1)-18运行结果:dy =sin(x)*log(x)+x*cos(x)*log(x)+sin(x)lz =1/6*x6-2/3*x3+1/3
8、*x(3/2)lw =1/2*exp(1)-1 lz =1/6*x6-2/3*x3+1/3*x(3/2)lw =0.35915、已知 ,使用命令画出两信号和及两信号乘积的波形图。ttftf 8sin,sin21其中, Hz2解:9-2 -1.5 -1 -0.5 0 0.5 1 1.5 2-2-1.5-1-0.500.511.52ty代码:clear all;close all;clc;t=-2:0.01:2;y1=sin(2*pi*t);y2=sin(16*pi*t);y3=y1.*y2;plot(t,y1,r,t,y2,b,t,y3,g),grid onaxis(-2,2,-2,2);xlabel(t),ylabel(y)6、用 Matlab 命令绘出下列信号的卷积积分 的时域波形图。)(21tf10解:-4 -2 0 2 4 6 8 10 12-2-10123456 f=f1(t)*f2(t)代码:clear all;close all;clc;t=-2:0.01:4;dt=0.01;f1=uCT(t)-uCT(t-2);f2=uCT(t)+uCT(t-1)-uCT(t-2);f=conv(f1,f2)*dt;n=length(f);tt=(0:n-1)*dt-2;plot(tt,f),grid onaxis(-4,12,-2,6);title(f=f1(t)*f2(t))