《2022年数字信号处理MATLAB实验.pdf》由会员分享，可在线阅读，更多相关《2022年数字信号处理MATLAB实验.pdf（12页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、数字信号处理MATLAB 实验实验一熟悉 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语句画出
2、 A、B、C、D、E、F、G。(2)用MATLAB 实现以下序列。a)x(n)=0、8n0n15 b)x(n)=e(0、2+3j)n0n15 c)x(n)=3cos(0 、125 n+0、2 )+2sin(0、25 n+0、1 ) 0n15 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, 产生并绘出下列序列的样本。精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - -
3、- - - - - - -第 1 页，共 12 页 - - - - - - - - - - 数字信号处理MATLAB 实验a)x1(n)=2x(n+2)-x(n-1)-2x(n) b)51k2)kn(nx(n) x(4)绘出下列时间函数的图形,对x轴、 y轴以及图形上方均须加上适当的标注。a) x(t)=sin(2 t) 0t10s b) x(t)=cos(100 t)sin( t) 0t4s (5)编写函数 stepshift(n0,n1,n2) 实现u(n-n0),n1n0n2,绘出该函数的图形,起点为 n1,终点为n2。(6)给定一因果系统)0.9z0.67z-1)/(1z2(1H(z)
4、-2-1-1求出并绘制 H(z)的幅频响应与相频响应。(7)计算序列 8 -2 -1 2 3 与序列 2 3 -1 -3 的离散卷积 ,并作图表示卷积结果。(8)求以下差分方程所描述系统的单位脉冲响应h(n), 0n50 y(n)+0、1y(n-1)-0 、06y(n-2)=x(n)-2x(n-1) 三、思考题(1)对于有限长序列 ,如何用 MATLAB 计算其DTFT ？(2)对于由两个子系统级联或并联的系统,如何用MATLAB 计算它们的幅频响应与相频响应？四、实验报告要求(1)简述实验目的及原理。(2)按实验步骤附上实验程序。(3)按实验步骤附上有关离散系统的频率特性曲线。(4)简要回答
5、思考题。五、与本实验有关的MATLAB 函数x=sin(2*pi*f/fs*n); 生成频率为 f,采样频率为 fs的正弦信号 ,式中,n=0 1 2 N。精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 2 页，共 12 页 - - - - - - - - - - 数字信号处理MATLAB 实验sum(X);对于向量 X,计算X各元素的与。对于矩阵 X,计算X各列元素之与组成的行向量。plot(t,y);画出以向量 t为坐标的向量 y(行或列)的曲线。向量t与向量y具有相同的维数。命令 plot(s1,t
6、1,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,
7、a)为参数的滤波器 N点频率向量 f与N点复频率向量 H。该命令用于绘制滤波器的幅频与相频响应。 如果省略左边的 H,f,该命令将直接绘出滤波器的对数幅频与相频响应。实验二信号的采样与重建一、实验目的（1）在学习本章内容的基础上,通过实验加强本章内容的有关信号采样与重建的基本概念 ,熟悉相关 MATLAB 函数。（2）通过观察采样信号的混叠现象,进一步理解奈奎斯特采样频率的意义。（3）通过实验 ,了解数字信号采样率转换过程中的频谱特性。（4）对实际的音频文件做内插与提取操作,体会低通滤波器在内插与提取中的作用。精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载
8、 名师归纳 - - - - - - - - - -第 3 页，共 12 页 - - - - - - - - - - 数字信号处理MATLAB 实验二、实验内容认真阅读相关的MATLAB 函数帮助文件。再熟悉 MATLAB 函数的基础上 ,完成以下试验。上机实验内容 : （1）一信号就是三个正弦信号的与,正弦信号的频率分别为50、500、1000Hz, 该信号以 8kHz 采样。用适当数量的样本画出该信号。（2）一信号就是三个正弦信号的与,正弦信号的频率分别为50、500、1000Hz, 该信号以 800Hz采样。用适当数量的样本画出该信号,并讨论信号的混叠情况。（3）令,其中 f/fs=1/1
9、6, 即每个周期内有16 个点。试利用 MATLAB 编程实现 : 1 作 M=4 倍的抽取 ,使每个周期变成 4 点。 2 作 L=3 倍的插值 ,使每个周期变成48点。（4）输入信号 x(n)为归一化频率分别为f1=0、04,f2=0、3 的正弦信号相加而成 ,N=50,内插因子为 5,抽取因子为 3,给出按有理因子 5/3 做采样率变换的输入输出波形。（5）常见的音频文件采样率为44、1khz。请找一个wav 格式、采样率为 44、1khz 的音频文件 ,用 MATLAB 编写程序 ,把它转换为采样率为48、32、22、05、16 与 8khz 的音频文件 ,用播放器分别进行播放 ,比较
10、音质的变化 ,并解释原因。（6）请找一个wav 格式、采样率为11、025khz 的音频文件 ,用精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 4 页，共 12 页 - - - - - - - - - - 数字信号处理MATLAB 实验MATLAB 编写程序 ,把它转换为采样率为44、1khz 的音频文件,采用两种方法: 1 直接插零 ; 2 插零后滤波 (可直接利用MATLAB 相关函数 ), 用播放器分别进行播放 ,比较音质的变化,并解释原因。三、思考题（1）试说明对于周期信号 ,应当如何采样 ,才
11、能保证周期扩展后与原信号保持一致。（2）模拟抗混叠滤波器的指标就是如何确定的,欠采样的情况下就是否需要模拟抗混叠滤波器？（3）抽取就是否会造成信号频谱成分的缺失？为什么还要这样做？四、实验报告及要求（1）简述实验目的及原理。（2）按实验步骤附上试验程序。（3）按实验内容附上有关离散信号的波形或关键样本,对音频信号给出测听的结果。（4）简要回答思考题。五、与本实验相关的MATLAB 函数y=decimate(x,M); 对信号 x 按整数 M 作抽取 ,抽取前后作抗混叠低通滤波 ,结果放在 y中。y=interp(x,L); 对信号 x 按整数 L 插零,然后作抗镜像低通滤波 ,结果放在 y 中
12、。精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 5 页，共 12 页 - - - - - - - - - - 数字信号处理MATLAB 实验y=resample(x,L,M); 对信号 x 按有理因子 L/M 作采样率转换 ,结果放在 y 中。y,fs,bits=wavread( );读取音频波形文件 ,这里 y 就是声音的数组,fs就是声音的采样频率 ,bits就是采样样本的二进制位数。wavwrite(y,fs,bits, );将声音的数组 y写入音频波形文件 ,fs就是声音的采样频率 bits 就
13、是采样样本的二进制位数。实验三快速傅立叶变换及其应用一、实验目的(1)在理论学习的基础上,通过本实验 ,加深对FFT 的理解 ,熟悉MATLAB 中函数。(2)熟悉应用 FFT 对典型信号进行频谱分析的方法。(3)了解应用 FFT 进行信号频谱分析过程中可能出现的问题,以便在实际中正确应用FFT。(4)熟悉应用 FFT 实现两个序列的线性卷积与相关。二、实验内容实验中用到的信号序列: 高斯序列衰减正弦序列精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 6 页，共 12 页 - - - - - - - -
14、- - 数字信号处理MATLAB 实验三角波序列反三角波序列上机实验内容 : (1)观察高斯序列的时域与幅频特性,固定信号 xa(n)中参数 p=8,改变 q 的值,使 q 分别等于 2、4、8,观察她们的时域与频域特性,了解当 q取不同值时 ,对信号时域频域特性的影响;固定 q=8,改变 p,使 p 分别等于 8,13,14, 观察参数 p 变化对信号序列的时域与频域特性的影响,注意当 p=多少时会发生明显的泄漏现象,混叠就是否也随时出现？记录实验中观察到的现象 ,绘出相应的时域序列与幅频特性曲线。(2)观察衰减正弦序列xb(n)的时域与幅频特性。=0、1,f=0、0625,检查谱峰出现位置
15、就是否正确,注意频谱的形状 ,绘出幅频特性曲线 ,改变 f,使 f 分别等于 0、4375 与 0、5625, 观察这两种情况下 ,频谱的形状与谱峰出现位置 ,有无混叠与泄露现象？说明产生现象的原因。(3)观察三角波与反三角波序列的时域与幅频特性,用 N=8 点 FFT分析信号序列 xc(n)与 xd(n)的幅频特性 ,观察两者的序列情况与频谱曲线有什么异同？绘出两序列及其幅频特性曲线。在 xc(n)与 xd(n)末尾补零 ,用 N=32 点 FFT 分析这两个信号的幅频特性,观察幅频特性发生了什么变化？这些变化说明了什么？精品资料 - - - 欢迎下载 - - - - - - - - - -
16、 - 欢迎下载 名师归纳 - - - - - - - - - -第 7 页，共 12 页 - - - - - - - - - - 数字信号处理MATLAB 实验(4)一个连续信号含有两个频率分量,经采样得已知 N=16,f 分别为 1/16 与 1/64,观察其频谱 ;当 N=128 时,f 不变,其结果有何不同 ,为什么？(5)用 FFT 分别计算 xa(n) (p=8,q=2)与 xb(n) (a=0、1,f=0、0625) 的16点循环卷积与线性卷积。(6)产生一 512 点的随机序列xe(n), 并用 xc(n)与 xe(n)做线性卷积 ,观察卷积前后 xe(n)频谱的变化。要求将xe
17、(n)分成 8 段,分别采用重叠相加法与重叠保留法。(7) 用 FFT 分别计算)2, 8)(qpnxa与 xb(n) (a=0、1,f=0、0625)的 16点循环相关与线性相关 ,问一共有多少种结果 ,她们之间有何异同点。(8) 用 FFT 分别计算)2, 8)(qpnxa与 xb(n) (a=0、1,f=0、0625)的自相关函数。三、思考题(1)实验中的信号序列xc(n)与 xd(n),在单位圆上的z 变换频谱与会相同不？如果不同 ,说明哪一个低频分量更多一些,为什么？(2)对一个有限长序列进行DFT 等价于将该序列周期延拓后进行DFS展开,因为 DFS也只就是取其中一个周期来运算,所
18、以 FFT 在一定条 件 下 也 可 以 用 于 分 析 周 期 信 号 序 列 。 如 果 时 正 弦 序 列,f=0、1 用 16点 FFT 来作 DFS运算,得到的频谱就是信号本精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 8 页，共 12 页 - - - - - - - - - - 数字信号处理MATLAB 实验身的真实谱不？为什么？四、实验报告要求(1)简述实验目的及原理。(2)按实验步骤附上实验信号序列与幅频特性曲线,分析所得到的图形,说明参数改变对时域与频域的影响。(3)总结实验中的主要结
19、论。(4)简要回答简答题。五、与本实验相关的MATLAB 函数y=exp(X);对向量 X 的各元素做指数运算 ,结果为一向量。conj(X); 对向量 X 的各元素做复共轭运算 ,即将虚部改变符号。real(X);对向量 X 的各元素取其实部。v=randn(size(X); 生成同 X 具有相同维数的正态分布的随机矩阵,通常用于生成高斯白噪声。fft(X,N); 计算 X 的 N 点 FFT,如果 X 的长度小于 N,则将在 X 后补零。反之 ,如果 X 的长度大于 N,则对 X 进行截取。ifft(X); 计算 X 的 N 点 IFFT。fftshift(Y);如果 Y为向量 ,该命令将
20、 Y分成左右两部分并交换位置。实验四 IIR数字滤波器设计一 、实验目的(1) 掌握双线性变换法及脉冲响应不变法设计IIR数字滤波器的具体设计方法及其原理 , 熟悉用双线性变换法及脉冲响应不变法设计低通、高通与带通IIR 数字滤波器的计算机编程。(2) 观察双线性变换及脉冲响应不变法设计的滤波器的频域特性, 了解双线精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 9 页，共 12 页 - - - - - - - - - - 数字信号处理MATLAB 实验性变换法及脉冲响应不变法的特点。(3) 熟悉巴特沃思
21、滤波器、切比雪夫滤波器与椭圆滤波器的频率特性。二 、实验内容实验中有关变量的定义 : fc 通带边界频率fr 阻带边界频率通带波动At 最小阻带衰减fs 采样频率T 采样周期上机实验内容 : (1)kHzfc3.0,dB8.0,kHzfr2 .0,dBAt20,msT1, 设计一切比雪夫高通滤波器 , 观察其通带损耗与阻带衰减就是否满足要求。(2)kHzfc2.0,dB1,kHzfr3. 0,dBAt25,msT1, 分别用脉冲响应不变法及双线性变换法设计一巴特沃思数字低通滤波器, 观察所设计数字滤波器的幅频特性曲线 , 记录带宽与衰减量 , 检查就是否满足要求。 比较这两种方法的优缺点。(3
22、) 利用双线性变换法分别设计满足下列指标的巴特沃思型、切比雪夫型与椭 圆 型 数 字 低 通 滤 波 器 , 并 作 图 验 证 设 计 结 果 : f =1、 2kHz, 0、5dB,fr=2kHz,At40dB,fs=8kHz。 比较这三种滤波器的阶数。(4) 分别用脉冲响应不变法与双线性变换法设计一巴特沃思型数字带通滤波器, 已知 fs=30kHz, 其等效的模拟滤波器指标为3dB,2kHzf 3kHz;At 5dB,f 6kHz;At 20dB,f 1、5kHz。(5) 利用双线性变换法设计满足下列指标的切比雪夫型数字带阻滤波器, 并作图验证设计结果 ; 当 1kHzf 2kHz;At
23、 18dB;当 f 500Hz以及 f 3kHz,3dB;采样频率 fs=10kHz。三、思考题(1) 双线性变换法中与之间的关系就是非线性的, 在实验中您注意到这精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 10 页，共 12 页 - - - - - - - - - - 数字信号处理MATLAB 实验种非线性变换了不？从哪几种数字滤波器的幅频特性曲线中可以观察到这种非线性关系？(2) 能否利用公式H(z)=H(s),s=lnz/T完成脉冲响应不变法的数字滤波设计？为什么？四、实验报告及要求(1) 简述
24、实验目的及原理。(2) 按实验步骤附上所设计滤波器系统函数H(z) 及相应的幅频特性曲线 , 定性分析她们的性能 , 判断设计就是否满足要求。(3) 谈谈双线性变换法的特点 , 简述用双线性变换法设计滤波器的全过程。(4) 简要回答思考题。实验五 FIR数字滤波器的设计一、实验目的(1) 掌握用窗函数法 , 频率采样法及优化设计法设计FIR 滤波器的原理及方法, 熟悉响应的计算机编程 ; (2) 熟悉线性相位 FIR 滤波器的幅频特性与相频特性; (3) 了解各种不同窗函数对滤波器性能的影响。二、实验内容(1) N=45,计算并画出矩形窗、汉明窗、布莱克曼窗的归化幅度谱, 并比较各自的特点。(
25、2) N=15,带通滤波器的两个通带边界分别为1=0、3 , 2=0、5 。用汉宁窗设计线性相位带通滤波器, 观察它的实际 3dB与 20dB带宽。 N=45,重复着一设计, 观察幅度与相位特性的变化, 注意长度 N变化的影响。(3) 分别改用矩形窗与布莱克曼窗, 设计 (2) 中的带通滤波器 , 观测并记录窗函数对滤波器幅独特性的影响, 比较三种窗的特点。(4) 用 Kaiser窗设计一专用线性相位滤波器,N=40, Hd(ejw) 如图 , 当0=4,6,10 时, 分别设计 , 比较它们的幅频与相频特性, 注意 0取不同值时的影响。精品资料 - - - 欢迎下载 - - - - - -
26、- - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 11 页，共 12 页 - - - - - - - - - - 数字信号处理MATLAB 实验(5) 用 频 率 采样法 设计 (4) 中的 滤 波 器 , 过渡 带 分别 设 一 个过 渡点, 令H(k)=0、5。比较两种不同方法的结果。(6) 用雷米兹 (Remez)交替算法设计 (4) 中的滤波器 , 并比较 (4) 、(5) 、(6) 三种不同方法的结果。(7) 利用雷米兹交替算法 , 设计一个线性相位高通FIR 数字滤波器 , 其指标为 :通带边界频率 fc=800Hz, 阻带边界频率 fr=500Hz
27、, 通带波动 =1dB阻带最小衰减At=40dB,采样频率 fs=5kHz。三、思考题(1) 定性地说明用本实验程序设计的FIR 滤波器的3dB 截至频率在什么位置？它等于理想频率响应的截至频率不？(2) 如果没有给定h(n) 长度 N,而就是给定了通带边缘截至频率与阻带临界频率, 以及相应的衰减 , 能根据这些条件用窗函数设计线性相位FIR 低通滤波器不？四、实验报告及要求(1) 简述实验目的及要求。(2) 按实验步骤附上所设计滤波器的h(n) 及相应幅频与相频特性曲线, 比较她们的性能 , 说明不同的设计方法对滤波器性能的影响。(3) 总结三种不同的设计方法的特点, 归纳设计中的主要公式。(4) 简要回答思考题。精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 12 页，共 12 页 - - - - - - - - - -