数字信号处理B实验指导书.doc

- PAGE 11 - 实验1 矩阵、序列及系统的时域表示与运算 一、实验目的 掌握用MATLAB表示离散时间信号和系统，以及它们的运算和显示。 二、实验内容与要求 用MATLAB产生并画出下列序列的样本。 a． b． c． d．，式中是在【-1，1】之间均匀分布的随机序列，你如何表征这个序列？ e．为周期的，画出5个周期。 2若线性时不变系统的单位样值响应为，输入序列，求系统的输出，并画出其波形图。(思考：你可以用几种方法来实现？) 三、实验所用部分函数如下 1.单位冲激序列（信号）生成函数impseq [x,n] = impseq(n0,n1,n2) 2.阶跃序列（信号）生成函数stepseq [x,n] = stepseq(n0,n1,n2) 3.序列（信号）相加函数sigadd [y,n] = sigadd(x1,n1,x2,n2) 以上为MATLAB没有，需外加入的函数（将相应函数拷贝到自己当前目录下） 4. 正（余）弦生成函数sin、cos y = sin(x) ,y = cos(x) (注意:x以弧度为单位) 5. 随机序列生成函数rand，用法如： Y = rand (n) 生成n×n阶的均匀分布随机阵； Y = rand (m, n) 生成m×n阶的随机阵； rand 返回在[0，1]区间上的一个随机数； 将上面的rand写成randn则可以生成均值为0、方差为1的正态分布的随机变量。 6．全1矩阵生成函数ones (m, n) ：生成m×n阶全1矩阵 7．全0矩阵生成函数zeros (m, n) ：生成m×n阶全0矩阵 8．离散序列绘图函数stem stem (y) 以1、2、3…为横坐标，?y为纵坐标画杆形图； stem(x, y) 以x为横坐标，?y为纵坐标画杆形图(x与y数据个数必须一致)； stem (…,’fill’) 选项’fill’指定杆顶为实心，若无此选项则默认空心。 stem(...,LineSpec) 参数LineSpec指定杆形图的线形、数据标志符号及颜色，具体用法可查看MATLAB帮助 9．线性坐标平面绘图函数plot 用法与stem?类似，具体用法可查看MATLAB帮助 以上为MATLAB内置函数（在此仅为同学复习MATLAB提供） 四、在MATLAB程序中变量赋值注意问题 在MATLAB?中，对变量赋值时其维数可以按需要动态地改变，这样虽然方便程序设计但同时容易出错。另外，频繁分配变量空间会大大降低程序的执行速度，因而应该尽量避免不必要的矩阵、向量维数的改变。通常先用zeros()函数给变量分配足够大小的空间，再对变量进行赋值。 例：依次执行下面的语句 tic %开始计时 for i=1:10000 c (i) = i; %每次都重新分配空间 end toc %读取计时时间 tic %开始计时 d=zeros(1,10000); %预先分配空间 for i=1:10000 d (i) = i; %直接赋值，不必重新分配空间 end toc %读取计时时间 运行结果如下： elapsed_time = 1．1560 elapsed_time = 0．0470 从结果可以看出，第2种赋值方法所用的时间比第1种方法所用时间少得多（以上是在主频为2．66GHZ的机器上运行的结果）。 实验2 离散信号、系统的频域表示 一、实验目的 考察抽样间隔对信号频谱的影响； 掌握用FFT做谱分析的方法。 二、实验内容与要求 1、用DFT/FFT对模拟信号做傅里叶分析 以频率fs 对以下信号抽样N点 xa(t) = cos (a t) + cos (b t) + cos (c t) 相应的参数是 a = 2*pi*6500, b = 2*pi *7000, c = 2*pi*9000 fs = 32000，N = 16 对这N点序列作N点DFT，观察其幅频特性，如果 X = fft (x) w是频率坐标向量，你可以考虑用stem (w, abs(X)), plot (w, abs(X)), plot (w, abs(X), *)来显示，然后确定用哪种显示方式。 接下来，对x作M = 256点FFT，这意味着在x后补了M-N个0，