数文字信号处理MATLAB实验.docxVIP

实验一 熟悉 MATLAB 环境 一、实验目的 (1)熟悉 MATLAB 的主要操作命令。 (2)学会简单的矩阵输入和数据读写。 (3)掌握简单的绘图命令。 (4)用 MATLAB 编程并学会创建函数。 (5)观察离散系统的频率响应。 二、实验内容 认真阅读本章附录，在 MATLAB 环境下重新做一遍附录中的例 子，体会各条命令的含义。在熟悉了 MATLAB 基本命令的基础上， 完成以下实验。 上机实验内容： (1)数组的加、减、乘、除和乘方运算。输入 A=[1 2 3 4]，B=[3 4 5 6]， C=A+B ，D=A-B ，E=A.*B ，F=A./B ，G=A.^B 并用 stem 语句画出 A、B、C、D、E、F、G。 (2)用MATLAB 实现以下序列。 a）x(n)=0.8n 0≤n≤15 b）x(n)=e(0.2+3j)n 0≤n≤15 c）x(n)=3cos(0.125πn+0.2π)+2sin(0.25πn+0.1π) 0≤n≤15 d）将 c)中的 x(n)扩展为以 16为周期的函数 x16(n)=x(n+16),绘出四个周期。 e）将 c)中的 x(n)扩展为以 10为周期的函数 x10(n)=x(n+10),绘出四个周期。 (3)x(n)=[1,-1,3,5]，产生并绘出下列序列的样本。 a）x1(n)=2x(n+2)-x(n-1)-2x(n) 5 b） x2 (n) nx(n k) k 1 (4)绘出下列时间函数的图形， 对x轴、y轴以及图形上方均须加上适当 的标注。 a) x(t)=sin(2πt) 0≤t≤10s b) x(t)=cos(100πt)sin(πt) 0≤t≤4s (5)编写函数 stepshift(n0,n1,n2)实现 u(n-n0),n1n0n2,绘出该函数的图 形，起点为 n1，终点为 n2。 (6)给定一因果系统 H(z) (1 2z-1 1)/(1 - 0.67z -1 0.9z -2 ) 求出并绘制 H(z) 的幅频响应与相频响应。 (7)计算序列｛ 8 -2 -1 2 3｝和序列｛ 2 3 -1 -3｝的离散卷积，并作图表 示卷积结果。 (8)求以下差分方程所描述系统的单位脉冲响应 h(n), 0≤n50 y(n)+0.1y(n-1)-0.06y(n-2)=x(n)-2x(n-1) 三、思考题 (1)对于有限长序列，如何用 MATLAB 计算其 DTFT ？ (2)对于由两个子系统级联或并联的系统， 如何用 MATLAB 计算它们的 幅频响应与相频响应？ 四、实验报告要求 (1)简述实验目的及原理。 (2)按实验步骤附上实验程序。 (3)按实验步骤附上有关离散系统的频率特性曲线。 (4)简要回答思考题。 五、与本实验有关的 MATLAB 函数 x=sin（2*pi*f/fs*n ）;生成频率为 f，采样频率为 fs的正弦信号，式 中， n=[0 1 2 ⋯ N] 。 sum（X）； 于向量 X， 算 X各元素的和。 于矩 X， 算 X 各列元素之和 成的行向量。 plot（t，y）；画出以向量 t 坐 的向量 y（行或列）的曲 。向 t和向量 y具有相同的 数。命令 plot（s1,t1,s2,t2,s3,t3）;将在同一 上画出分 以 t1,t2,t3 坐 的向量 s1,s2,s3的曲 。 xlabel(‘samples’);在x 上加上 注。ylabel(‘amplitude’);在y 上加上 注。 title(‘sinusoidal signal’);在 的上部加上 。 y=conv（h，x）； 算向量 h和x的卷 ， 果放在 y中。 y=filter （b，a，x）；以向量 b和a 参数的 波器 入信号向 x 行 波 理。 h=impz （b，a，N）； 算 b和a 参数的 N 点 波器脉冲响 。 [H,f]=freqz （b，a，N,Fs）； 定以 Hz 位的采 率 Fs， 算以（ b，a） 参数的 波器 N点 率向量 f和N点复 率向量 H。 命令用于 制 波器的幅 和相 响 。如果省略左 的 [H,f] ， 命 令将直接 出 波器的 数幅 和相 响 。 实验二 信号的采样与重建 一、 实验目的 （1）在学习本章内容的基础上，通过实验加强本章内容的有关信号 采样与重建的基本概念，熟悉相关 MATLAB 函数。 2）通过观察采样信号的混叠现象，进一步理解奈奎斯特采样频率的意义。 3）通过实验，了解数字信号采样率转换过程中的频谱特性。 4）对实际的音频文件做内插和提取操作，体会低通滤波器在内插和提取中的作用。 二、 实验内容 认真阅读相关的 MATLAB 函数帮助文件。再熟悉 MATLAB 函数的基础上，完成以下试验。 上机实验内