资源描述:
《信号分析与处理实验报告基于MATLAB.doc》由会员上传分享,免费在线阅读,更多相关内容在应用文档-天天文库。
1、《信号分析与处理》实验报告武汉工程大学电气信息学院专业班级14通信3班实验时间12-16周周二上午学生学号1404200529实验地点4B315学生姓名赵晶指导教师许楠实验项目信号分析与处理实验(基于MATLAB)实验类别基础实验实验学时10学时实验目的及要求一、掌握连续信号的MATLAB表示方法及用MATLAB描绘二维图像二、掌握用MATLAB对连续信号进行基本运算和时域变换的方法;三、掌握两个连续时间信号卷积的计算方法和MATLAB编程技术。四、掌握LTI系统的微分方程描述方法及其MATLAB编程的求解方法。五、掌握周期函数的傅里叶级数计算方法和编程技术,掌握用MATLAB
2、进行傅里叶正反变换的方法。六、掌握系统频率响应特性的计算方法和特性曲线的绘制方法;掌握用MATLAB语言进行系统频响特性分析的方法成绩评定表类别评分标准分值得分合计上机表现按时出勤、遵守纪律认真完成各项实验内容30分报告质量程序代码规范、功能正确填写内容完整、体现收获70分说明:评阅教师:日期:年月日《信号分析与处理》实验报告实验内容实验一:连续时间信号的时域表示一、实验内容1、参考示例程序,绘制信号的图形,t取-1到4,步长值设为0.01。2、产生一个指数为的复指数函数,绘出函数的实部、虚部、幅度和相位的波形,t取0到20,步长值设为0.1。二、实验方法与步骤1、绘制信号的图
3、形程序代码如下:clearallt=-1:0.01:4;%信号时间样本点向量f=exp(-2*t).*cos(3*pi*t).*(heaviside(t)-heaviside(t-3));%函数描述plot(t,f)%命令绘制信号的时域波形title('信号1')%添加标题xlabel('t')%添加横坐标说明ylabel('f')%添加纵坐标说明axis([-1,4,-0.7,1.1])%对横纵坐标进行限定gridon2、绘制指数为的复指数函数的实部、虚部、幅度和相位的波形程序代码如下:clearallt=0:0.1:20;%信号时间样本点向量《信号分析与处理》实验报告z=(
4、-0.1+(pi/4)*1i)*t;%函数描述f=exp(z);%定义指数信号fr=real(f);%描述函数实部fi=imag(f);%描述函数虚部fa=abs(f);%描述函数幅度fg=angle(f);%描述函数相位subplot(2,2,1)%将当前窗口分成2行2列个子窗口,并在第1个子窗口绘图plot(t,fr)title('实部')xlabel('t')axis([-0.5,20.5,-0.8,1.2])gridonsubplot(2,2,2)%将当前窗口分成2行2列个子窗口,并在第2个子窗口绘图plot(t,fi)title('虚部')xlabel('t')axi
5、s([-0.5,20.5,-0.8,1])gridonsubplot(2,2,3)%将当前窗口分成2行2列个子窗口,并在第3个子窗口绘图plot(t,fa)《信号分析与处理》实验报告title('幅度')xlabel('t')axis([-0.5,20.5,0,1.1])gridonsubplot(2,2,4)%将当前窗口分成2行2列个子窗口,并在第4个子窗口绘图plot(t,fg)title('相位')xlabel('t')axis([-0.5,20.5,-3.5,3.5])gridon三、实验数据与结果分析1、《信号分析与处理》实验报告2、四、思考:1、为什么图二中t=0处
6、曲线是间断的,如何使其成为连续的曲线?因为axis函数对纵坐标的的上边界限定过小,使图形在边界处不能完整的显示。《信号分析与处理》实验报告实验内容实验二:连续时间信号的时域分析一、实验内容1、已知,,绘制x(t)和y(t)的图形,t取-3到5,步长值设为0.01。2、根据符号函数和单位阶跃函数的关系,利用符号函数sign实现单位阶跃函数。要求图形窗口的横坐标范围为-5~5,纵坐标范围为-1.5~1.5。3、任意定义一个有限长时间信号y1(t),根据式2.1产生一个周期信号,绘制y1(t)和y(t)的图形。二、实验方法与步骤1、绘制x(t)和y(t)的图形程序代码如下:clear
7、t=-3:0.01:5;x=exp(-0.5*t).*heaviside(t);%描述因果信号x(t)y=exp(-0.5*1.5*(t+2)).*heaviside(1.5*(t+2));%描述因果信号x(1.5t+3)subplot(2,1,1)plot(t,x)title('x(t)')xlabel('t')ylabel('x')axis([-3,5,-0.1,1.1])gridon《信号分析与处理》实验报告subplot(2,1,2)plot(t,y)title('y(t)')xla