离散时间信号与系统实验报告.docVIP

实验报告一 学院：物理与机电工程学院 系：电子科学系 专业：无线电 年级：09级 姓名：刘奇 学号：19820092203790 组别：4 实验时间：2012年02月29日晚上 指导老师签字： 实验一 离散时间信号和系统 一． 实验目的： 掌握时域离散信号的产生及基本时域处理的Matlab实现方法；掌握离散信号卷积和相关的Matlab实现；线性时不变系统及性质验证；进一步理解和巩固理论知识，提高分析和解决实际问题的能力。 二． 实验原理： 1．时域离散信号的产生及时域处理 2．卷积与相关 3．线性时不变系统 三．实验内容： 1．用Matlab产生并画出（用stem函数）下列序列的样本： 1) 2) 其中u(n)是阶跃信号。 2．已知序列： 其中R(n)均匀分布的随机序列。 1) 在四个子图中画出其幅度、相位、实部和虚部。 2) 求该序列的能量和功率。 3．卷积与相关： 已知序列 经过系统h(n)后得到输出y(n)。系统函数为： 1) 求输出y(n)及y(n)的定义范围； 2) 求x(n)和y(n)的互相关。 4．系统线性性质验证： 一个特定的线性时不变系统，描述它的差分方程为： ， 1) 确定系统的稳定性； 2) 在之间求解并画出系统的脉冲响应，从脉冲响应确定系统的稳定性； 3) 如果此系统的输入为。在间求出y(n)的响应。 四. 实验结果＆结果分析： 1.解答： 命令窗口输出 代码 a.m文件代码： x1=[]; y=0; for n=0:25 for m=0:10 y=y+(m+1).*[impseq(2*m,0,25)-impseq(2*m+1,0,25)]; end x1=[x1,y(n+1)]; end n=[0:25]; subplot(2,1,1);stem(n,x1);title(1.(1)的序列); xlabel(n);ylabel(x1(n)); n=[-25:25]; s1=n.*n.*[stepseq(-5,-25,25)-stepseq(6,-25,25)]; s2=10*impseq(0,-25,25); s3=20*0.5.^n.*[stepseq(4,-25,25)-stepseq(10,-25,25)]; x2=s1+s2+s3; subplot(2,1,2);stem(n,x2);title(1.(2)的序列); xlabel(n);ylabel(x2(n)); impseq.m文件代码： function [x,n]=impseq(n0,n1,n2) n=[n1:n2]; x=[(n-n0)==0]; end stepseq.m文件代码： function[x,n]=stepseq(n0,n1,n2) n=[n1:n2]; x=[(n-n0)=0]; end 1.(1)的序列和(2)的序列图像如下： 2.解答： 命令窗口输出 代码 Ex = 47.8188 Px = 2.2771 b.m文件代码： n=-10:10; R=rand(1,21); s=0.1-0.3*pi*1i; x3=exp(s*n)+cos(2*n+0.2*pi)+R/10; subplot(2,2,1);stem(n,real(x3));title(实部);xlabel(n); subplot(2,2,2);stem(n,imag(x3));title(虚部);xlabel(n); subplot(2,2,3);stem(n,abs(x3));title(振幅);xlabel(n); subplot(2,2,4);stem(n,(180/pi)*angle(x3));title(相位);xlabel(n); Ex=sum(abs(x3).^2) Px=Ex/21 其幅度、相位、实部和虚部图如下： 3.解答： 命令窗口输出 代码 y = Columns 1 through 7 -6.4000 -11.2000 -4.0000 6.2000 12.5000 4.3000 -13.5000 Columns 8 through 14 -20.3000 -18.8000 -9.3000 10.4000 19.1000 19.8000 11.2000 ny = Columns 1 through 11 -1 0 1 2 3 4 5 6 7 8 9 Columns 12 through 14 10 11 12 c.m文件代码： x=[3.2,4,0,-1.5,-0.3,2,5,6,7.1,5.6]; n