武大电气信号与自动控制综合实验报告.pdfVIP

信号与自动控制 综合实验 （信号与系统部分） 学院：电气工程学院 姓名：吕知彼 班级：2013 级3 班 学号：2013301610345 WeChat:lvzhibi 实验 1 连续时间信号的表示及可视化 (1)、(2)： f (t) = δ (t)以及f (t) =ε(t)； clc;clear; t=-2:0.001:2; u=(t=0); subplot(211); plot(t,u); axis([-2 2 -0.1 1.1]) title(单位阶跃函数ε（t）); y=diff(u)/0.001; % 相邻元素逐个求差，除以dt = 0.001得到导数值 y=[0 y]; % 前面再加个零保持行向量y的大小不变 subplot(212); plot(t,y); axis([-2 2 -100 1100]) title(单位冲击函数δ（t）); 实验结果分析：在matlab中，单位阶跃函数为y=heaviside(t)，而单位冲 激函数为y=dirac(t)，但是用这两个函数来画图并没有太大的意义，因为前者 在0处的值为0.5，后者在0处的值为无穷大，即dirac （0）=Inf，若画出图来则 横坐标为0处没有图像，故采用u=(t=0)的方法来定义阶跃函数并对其求导得到 冲击函数，得到上面所示图像。 (3) f (t) = （分别取a 0及a 0） clc;clear; t=-2:0.001:2; y1=exp(t); subplot(211); plot(t,y1);title(y1=e^t);axis([-2,2,0,8]); y2=exp(-t); subplot(212); plot(t,y2);title(y2=e^{-t});axis([-2,2,0,8]); 实验结果分析：除非是符号运算，MATLAB 本身不能进行连续系统的运算。 用 MATLAB 分析连续时间信号，可以通过时间间隔取足够小的离散时间信号的数 值计算方法来实现。为了画出连续的指数函数图像，上图中的时间信号间隔取 0.001，相对于区间长度来说足够小。 (4)f (t) = R(t) （分别取a 0 及 a 0） clc;clear; t=-2:0.001:2; y=heaviside(t+1)-heaviside(t-1); plot(t,y); title(y=R(t));axis([-2,2,-0.2,1.2]);xlable(t);ylable(y); 实验结果分析：同（3）。 (5)f (t) = Sa(ωt)； clc;clear; w1=1;t=-20:0.001:20; y1=sin(w1*t)./(w1*t); subplot(211); plot(t,y1);title(y1=Sa(t));axis([-20,20,-0.2,1.2]); w2=2*pi;t=-20:0.001:20; y2=sin(w2*t)./(w2*t); subplot(212); plot(t,y2);title(y2=Sa(2πt));axis([-20,20,-0.2,1.2]); (6) f (t) = Sin(2 πft) （分别画出不同周期个数的波形）； clc;clear; t=-10:0.0001:10; f=50;y=sin(2*pi*f*t);subplot(211); plot(t,y); axis([-0.04,0.04,-1.2,1.2]);title(y=sin(2 π*50*t)); f=100;y=sin(2*pi*f*t);subplot(212); plot(t,y); axis([-0.04,0.04,-1.2,1.2]);title(y=sin(2 π*100*t)); 实验结果分析：同（3）。 实验 2 离散时间信号的表示及可视化 (1)、(2):f (n) =δ (n)以及f (n) = ε (n) clc;clear; n=-10:1:10; x=(n=0); subplot(211); stem(n,x); axis([-10 10 -0.1 1.1]) title(单位阶跃函数ε（n）); y=diff(x)/0.001; y=[0 y]; subplot