数字信号处理实验讲义
信号时频分析-讲义-WVD
Wigner-Ville 分布Wigner-Ville 分布可以看作是一大类分布的原型，它们和短时傅立叶变换谱有着本质的不同。
它首先由Wigner 提出，用于量子力学领域问题的研究，后由Ville 引入到信号分析。
因为在计算中，信号需要用到两次，因此Wigner-Ville 分布被称为一种二次型分布。
基本定义及计算Wigner-Ville 分布可由信号x (t )本身或它的频谱)(ωX 定义为如下两种等价方式ττ+τ-=ωτω-+∞∞-⎰d )e 21()21(π21)(i t x t x ,t WVD *x , (2.1.1)τθ+ωθ-ω=ωθ+∞∞-⎰d )e 21()21(π21)(i t *x X X ,t WVD . (2.1.2) 其中*表示复数共轭。
要证明上面两式是等价的，只需将信号写成它的频谱形式，然后将其代入到(2.1.1)式，即可得到(2.1.2)式。
式(2.1.1)中，)2/()2/(*ττ+-t x t x 称为信号的瞬时相关函数，因此Wigner-Ville 分布实质上是对信号的瞬时相关函数的傅立叶变换，它的结果能够反映信号的时频特征。
例2.1.1 对于信号)π400sin()(t t x = )10(≤≤t (2.1.3)其采样频率为1000 Hz 。
图2.1.1是其Wigner-Ville 分布，频率轴划分区间数为512。
图中清楚显示，该信号在整个时间段上，只含有一个频率为200Hz 的分量。
需要说明的是，图中显示的是Wigner-Ville 分布的绝对值，后面所有图中，如果没有特别注明，都默认显示的是绝对值。
图2.1.1 信号(2.1.3)的Wi gn er-Vi ll e 分布W i g n e r -V il l e 分布5000.2 0.4 0.6 10.20.40.60.8例2.1.2)π50π100sin()(2t t t x += )20(≤≤t , (2.1.4)这是一个线性调频信号。
数字信号处理讲义-2抽取与内插滤波器
L：内插的倍数 M：抽样的倍数。
例：离散信号x[k]是由抽样频率为10Hz,试求出抽样频
率为15Hz的序列y[k]。
f=0.35;N=40;
fs=10;fs1=15;
k=0:N-1;t=k/fs;
k1=0:N*1.5-1;t1=k1/fs1;
x=cos(2*pi*f*t);
xr=cos(2*pi*f*t1);
l=3时，要求的阻带为[(6p0.4p)/4,(6p+0.4p)/4]=[1.4p,1.6p]
综上所述，抽取滤波器阻带为 [0.4p,0.6p]，[0.9p, p] 7
2021/3/18
x = firls(511,[0 0.1 0.1 0.5 0.5 1],[1 1 1 0 0 0]); f=[0.1 0.4 0.6 0.9];a=[1 0 0];dev=[0.01 0.001 0.001]; %设计滤波器 [N,fo,ao,w] = remezord(f,a,dev); h = remez(N,fo,ao,w); xd=filter(h,[1],x); y=xd(1:4:end); w=linspace(0,pi,512); mag=freqz(h,[1],w); subplot(3,1,1);plot(w/pi,20*log10(abs(mag))); xlabel('Normalized frequency'); ylabel('Gain, dB'); X=freqz(x,[1],w);Y=freqz(y,[1],w); subplot(3,1,2);plot(w/pi,(abs(X))); xlabel('Normalized frequency'); ylabel('magnitude');title('Spectrum of x'); Subplot(3,1,3);plot(w/pi,abs(Y)); xlabel('Normalized frequency'); ylabel('magnitude');title('Spectrum of y');
数字信号处理实验
数字信号处理实验实验一信号、系统及系统响应1、实验目的认真复习采样理论、离散信号与系统、线性卷积、序列的z 变换及性质等有关内容；掌握离散时间序列的产生与基本运算，理解离散时间系统的时域特性与差分方程的求解方法，掌握离散信号的绘图方法；熟悉序列的z 变换及性质，理解理想采样前后信号频谱的变化。
2、实验内容a. 产生长度为500 的在[0,1]之间均匀分布的随机序列，产生长度为500 的均值为0 单位方差的高斯分布序列。
b. 线性时不变系统单位脉冲响应为h(n)=(0.9)nu(n)，当系统输入为x(n)=R10(n)时，求系统的零状态响应，并绘制波形图。
c. 描述系统的差分方程为：y(n)-y(n-1)+0.9y(n-2)=x(n)，其中x(n)为激励，y(n)为响应。
计算并绘制n=20,30,40,50,60,70,80,90,100 时的系统单位脉冲响应h(n)；计算并绘制n=20,30,40,50,60,70,80,90,100 时的系统单位阶跃响应s(n)；由h(n)表征的这个系统是稳定系统吗？d. 序列x(n)=(0.8)nu(n)，求DTFT[x(n)]，并画出它幅度、相位，实部、虚部的波形图。
观察它是否具有周期性？e. 线性时不变系统的差分方程为y(n)=0.7y(n-1)+x(n)，求系统的频率响应H(ejω)，如果系统输入为x(n)=cos(0.05πn)u(n)，求系统的稳态响应并绘图。
f. 设连续时间信号x(t)=e-1000|t|，计算并绘制它的傅立叶变换；如果用采样频率为每秒5000 样本对x(t)进行采样得到x1(n)，计算并绘制X1(ejω)，用x1(n)重建连续信号x(t)，并对结果进行讨论；如果用采样频率为每秒1000 样本对x(t)进行采样得到x2(n)，计算并绘制X2(ejω)，用x2(n)重建连续信号x(t)，并对结果进行讨论。
加深对采样定理的理解。
g. 设X1(z)=z+2+3z-1，X2(z)=2z2+4z+3+5z-1，用卷积方法计算X1(z)X2(z)。
数字信号处理EXPIV型教学实验系统实验四常规实验exp4_常规
第四章常规实验指导实验一常用指令实验一、实验目的1、了解DSP开发系统的组成和结构；2、熟悉DSP开发系统的连接；3、熟悉CCS的开发界面；4、熟悉C54X系列的寻址系统；5、熟悉常用C54X系列指令的用法。
二、实验设备计算机，CCS 2.0版软件，DSP仿真器，实验箱。
三、实验步骤与内容1、系统连接进行DSP实验之前，先必须连接好仿真器、实验箱及计算机，连接方法如下所示：2、上电复位在硬件安装完成后，确认安装正确、各实验部件及电源连接正常后，接通仿真器电源，启动计算机，此时，仿真器上的“红色小灯”应点亮，否则DSP开发系统有问题。
3、运行CCS程序待计算机启动成功后，实验箱后面220V输入电源开关置“ON”，实验箱上电，启动CCS，此时仿真器上的“绿色小灯”应点亮，并且CCS正常启动，表明系统连接正常；否则仿真器的连接、JTAG接口或CCS相关设置存在问题，掉电，检查仿真器的连接、JTAG 接口连接，或检查CCS相关设置是否正确。
注：如在此出现问题，可能是系统没有正常复位或连接错误，应重新检查系统硬件并复位；也可能是软件安装或设置有问题，应尝试调整软件系统设置，具体仿真器和仿真软件CCS的应用方法参见第三章。
●成功运行程序后，首先应熟悉CCS的用户界面●学会CCS环境下程序编写、调试、编译、装载，学习如何使用观察窗口等。
4、修改样例程序，尝试DSP其他的指令。
注：实验系统连接及CCS相关设置是以后所有实验的基础，在以下实验中这部分内容将不再复述。
5、填写实验报告。
6、样例程序实验操作说明仿真口选择开关K9拨到右侧，即仿真器选择连接右边的CPU：CPU2；启动CCS 2.0，在Project Open菜单打开exp01_cpu2目录下面的工程文件“exp01.pjt”注意：实验程序所在的目录不能包含中文，目录不能过深，如果想重新编译程序，去掉所有文件的只读属性。
用下拉菜单中Project/Open，打开“exp01.pjt”，双击“Source”，可查看源程序在File Load Program菜单下加载exp01_cpu2\debug目录下的exp01.out文件：加载完毕，单击“Run”运行程序；实验结果：可见指示灯D1定频率闪烁；单击“Halt”暂停程序运行，则指示灯停止闪烁，如再单击“Run”，则指示灯D1又开始闪烁；注：指示灯D1在CPLD单元的右上方关闭所有窗口，本实验完毕。
数字信号处理实验(1-7)原始实验内容文档(含代码)
实验要求1.每个实验进行之前须充分预习准备，实验完成后一周内提交实验报告；2.填写实验报告时，分为实验题目、实验目的、实验内容、实验结果、实验小结五项；3.实验报告要求：实验题目、实验目的、实验内容、实验结果四项都可打印；但每次实验的实验内容中的重要代码（或关键函数）后面要用手工解释其作用。
实验小结必须手写！（针对以前同学书写实验报告时候抄写代码太费时间的现象，本期实验报告进行以上改革）。
实验一信号、系统及系统响应实验目的：1. 掌握使用MATLAB进行函数、子程序、文件编辑等基本操作；2. 编写一些数字信号处理中常用序列的3. 掌握函数调用的方法。
实验内容：1.在数字信号处理的基本理论和MATLAB信号处理工具箱函数的基础上，可以自己编写一些子程序以便调用。
（1）单位抽样序列δ(n-n0)的生成函数impseq.m（2）单位阶跃序列u(n-n0)的生成函数stepseq.m（3）两个信号相加的生成函数sigadd.m（4）两个信号相乘的生成函数sigmult.m（5）序列移位y(n)=x(n-n0)的生成函数sigshift.m（6）序列翻褶y(n)=x(-n)生成函数sigfold.m（7）奇偶综合函数evenodd.m（8）求卷积和2.产生系列序列，并绘出离散图。
（1） x1(n)=3δ(n-2)-δ(n+4) -5≤n≤5（2） x3(n)=cos(0.04πn)+0.2w(n) 0≤n≤50其中：w(n)是均值为0，方差为1 的白噪声序列。
3.设线性移不变系统的抽样响应h(n)=(0.9)^n*u(n),输入序列x(n)=u(n)-u(n-10),求系统的输出y(n).实验二 系统响应及系统稳定性1.实验目的（1）掌握 求系统响应的方法。
（2）掌握时域离散系统的时域特性。
（3）分析、观察及检验系统的稳定性。
2.实验原理与方法在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。
精品文档-数字信号处理实验(MATLAB版)刘舒帆-第13章
y(n)＝x((n＋m)N)RN(n) x(n)左移m位的过程可由以下步骤获得：
(1)将x(n)以N为周期进行周期延拓，得到
＝
x((n)N)； (2)将
左移m位，得到
； ~x(n)
(3)取 y(n)。
~x(n)
~x(n
m)
的主值序列，~x(得n到mx)(n)循环移位序列
10
有限长序列的移位也称为循环移位，原因是将x(n)左移m 位时，移出的m位又依次从右端进入主值区。下面举例说明。
11
例13-2 已知有限长序列x(n)＝［1，2，3，4，5，6］，
求x(n)左移2位成为新的向量y(n)，并画出循环移位的中间过
程。
解 MATLAB程序如下：
xn＝［1，2，3，4，5，6］；
%建立xn序列
Nx＝length(xn)；nx＝0：Nx－1；
nx1＝－Nx：2*Nx－1；%设立周期延拓的范围
subplot(4，1，2)，stem(nx1，x1)；%画出x1 subplot(4，1，3)，stem(ny1，y1)；%画出y1 subplot(4，1，4)，stem(ny1，RN1.*y1)； %画出y1的 主值部分 运行结果如图13-2所示。
13 图13-2 例13-2有限长序列的循环移位
x1＝xn(mod(nx1，Nx)＋1)；%建立周期延拓序列
ny1＝nx1－2；y1＝x1；%将x1左移2位，得到y1
12
RN＝(nx1>＝0)&(nx1<Nx)；%在x1的位置向量nx1上设置 主值窗
RN1＝(ny1>＝0)&(ny1<Nx)；%在y1的位置向量ny1上设置 主值窗
subplot(4，1，1)，stem(nx1，RN.*x1)；%画出x1的主 值部分
数字信号处理实验指导吴镇扬
实验一快速Fourier变换(FFT)及其应用一、实验目的1．在理论学习的基础上，通过本实验，加深对FFT的理解，熟悉FFT子程序。
2．熟悉应用FFT对典型信号进行频谱分析的方法。
3. 了解应用FFT进行信号频谱分析过程中可能出现的问题以便在实际中正确应用FFT。
4．熟悉应用FFT实现两个序列的线性卷积的方法。
5．初步了解用周期图法作随机信号谱分析的方法。
返回页首二、实验原理与方法在各种信号序列中，有限长序列信号处理占有很重要地位，对有限长序列，我们可以使用离散Fouier变换(DFT)。
这一变换不但可以很好的反映序列的频谱特性，而且易于用快速算法在计算机上实现，当序列x(n)的长度为N时，它的DFT定义为：反变换为：有限长序列的DFT是其Z变换在单位圆上的等距采样，或者说是序列Fourier 变换的等距采样，因此可以用于序列的谱分析。
FFT并不是与DFT不同的另一种变换，而是为了减少DFT运算次数的一种快速算法。
它是对变换式进行一次次分解，使其成为若干小点数的组合，从而减少运算量。
常用的FFT是以2为基数的，其长度。
它的效率高，程序简单，使用非常方便，当要变换的序列长度不等于2的整数次方时，为了使用以2为基数的FFT，可以用末位补零的方法，使其长度延长至2的整数次方。
(一)、在运用DFT进行频谱分析的过程中可能产生三种误差：(1)混叠序列的频谱时被采样信号的周期延拓，当采样速率不满足Nyquist定理时，就会发生频谱混叠，使得采样后的信号序列频谱不能真实的反映原信号的频谱。
避免混叠现象的唯一方法是保证采样速率足够高，使频谱混叠现象不致出现，即在确定采样频率之前，必须对频谱的性质有所了解，在一般情况下，为了保证高于折叠频率的分量不会出现，在采样前，先用低通模拟滤波器对信号进行滤波。
(2)泄漏实际中我们往往用截短的序列来近似很长的甚至是无限长的序列，这样可以使用较短的DFT来对信号进行频谱分析，这种截短等价于给原信号序列乘以一个矩形窗函数，也相当于在频域将信号的频谱和矩形窗函数的频谱卷积，所得的频谱是原序列频谱的扩展。
数字信号处理讲义第7章滤波器的设计方法
第７章滤波器的设计方法教学目的1．掌握由连续时间滤波器设计离散时间IIR滤波器的方法，包括冲激响应不变法，双线性变换法等；2．了解常用的窗函数，掌握低通IIR滤波器的频率变换法、用窗函数法设计FIR滤波器的方法；3．掌握FIR滤波器的逼近原理与设计方法。
教学重点与难点重点：本章是本课程的重中之重，滤波器的设计是核心内容之一。
1．连续时间滤波器设计离散时间IIR滤波器的方法，包括冲激响应不变法，双线性变换法等；2．常用的窗函数，掌握低通IIR滤波器的频率变换法、用窗函数法设计FIR滤波器的方法；3．掌握FIR滤波器的逼近原理与设计方法。
难点：1.冲激响应不变法，双线性变换法2.用窗函数法设计FIR滤波器FIR滤波器的逼近原理与设计方法基本概念7.0.1 选频滤波器的分类数字滤波器是数字信号处理的重要基础。
在对信号的过滤、检测与参数的估计等处理中, 数字滤波器是使用最广泛的线性系统。
数字滤波器是对数字信号实现滤波的线性时不变系统。
它将输入的数字序列通过特定运算转变为输出的数字序列。
因此，数字滤波器本质上是一台完成特定运算的数字计算机。
我们已经知道，一个输入序列x(n)，通过一个单位脉冲响应为h(n)的线性时不变系统后，其输出响应y(n)为∑∞-)(y))()()(n(nn=m*=xmhnhx将上式两边经过傅里叶变换，可得式中，Y (e j ω)、X (e j ω)分别为输出序列和输入序列的频谱函数， H (ejω)是系统的频率响应函数。
可以看出，输入序列的频谱X (e j ω)经过滤波后，变为X (e j ω)H (e j ω)。
如果|H (e j ω)|的值在某些频率上是比较小的，则输入信号中的这些频率分量在输出信号中将被抑制掉。
因此，只要按照输入信号频谱的特点和处理信号的目的，适当选择H (ej ω)，使得滤波后的X (e j ω)H (e j ω)符合人们的要求，这就是数字滤波器的滤波原理。
和模拟滤波器一样，线性数字滤波器按照频率响应的通带特性可划分为低通、高通、带通和带阻几种形式。
数字信号处理实验-音频信号处理
图1 任务一程序流程图1、音频信号采集道，只取第一个声道进行处理,接着使用sound函数以fs频率进行音频回放。
2、音频信号频域分析以采样间隔T划分时域并绘制出signal信号的时域波形；调用fft函数,对signal 进行快速傅里叶变换，用abs函数取傅里叶变换后结果的幅值进行幅频分析，绘制出频谱图。
在绘制频谱图时由于考虑到快速傅里叶变换的对称性，只取序列的前半部分进行观察分析。
3、音频信号分解为了实现音频信号的分解及合成，先对原信号的频谱图进行观察分析，发现原信号的主要能量集中在三个主要频率上，于是考虑用这三频率的正弦信号合成原信号。
为了求得这三个频率，先调用findpeaks函数找到频谱图上的各个局部极大值peak及其对应的位置locs,然后用sort对峰值点进行排序，找到最大的三个值，接着用find 函数找到这三个最大值在locs中的位置，也就知道了对应的频率。
这里有一个问题就是最小的峰值频率并不是在sort排序后的第三位而是在第四位，需要有一个调整；确定了主要谱线后，使用text函数进行峰值标注；4、音频信号合成接着将这三个谱线还原回时域正弦信号，幅度的比例等于对应频率上的幅度比例然后然后叠加，得到合成后的信号，绘制出时域波形，与原信号波形进行比较,接着对两个正弦信号进行fft,绘制出他们的频谱，然后对合成的信号进行fft,做出频谱图和原信号的频谱图进行比较.5、音频信号回放用sound函数进行原信号和合成信号的回放，比较差异。
实验内容二：任意音频信号的时域和频域分析及数字滤波器设计通过对任务具体内容的分析，可以建立出任务二程序框图如下，之后将对编程思想及思路进行介绍：图2任务二程序流程图1、音频信号采样自己录音频并另存为”ding.wav”后,先用audioread函数读取音频信号得到采样序列signal及对应采样频率fs,由于获取的音频信号是双声道，只取第一个声道进行处理。
2、时域采样使用audioread函数得到的采样序列signal及采样频率fs为过采样状态，此时我们对signal再进行等间隔采样，达到减少采样点数和降低采样频率的效果，进而实现合理采样状态signal2、fs2和欠采样状态signal1、fs1;使用sound函数分别对这两种采样状态进行回放。
数字信号处理实验一
数字信号处理实验一实验目的：掌握利用Matlab产生各种离散时间信号，实现信号的相加、相乘及卷积运算实验函数：参考课本77-19页，注意式(2.11.1)的表达与各matlab子函数间的关系。
1、stem(x,y) % 绘制以x为横轴，y为纵轴的离散序列图形2、[h ,t] = impz(b, a) % 求解数字系统的冲激响应h，取样点数为缺省值[h, t] = impz(b, a, n) % 求解数字系统的冲激响应h，取样点数为nimpz(b, a) % 在当前窗口用stem(t, h)函数出图3、[h ,t] = dstep(b, a) % 求解数字系统的阶跃响应h，取样点数为缺省值[h, t] = dstep (b, a, n) % 求解数字系统的阶跃响应h，取样点数为ndstep (b, a) % 在当前窗口用stairs(t, h)函数出图4、y = filter(b,a,x) % 在已知系统差分方程或转移函数的情况下求系统输出实验原理：一、常用的时域离散信号及其程序1、产生单位抽样函数δ(n)n1 = -5;n2 = 5;n0 = 0;n = n1:n2;x = [n==n0]; % x在n=n0时为1，其余为0stem(n,x,'filled'); %filled:序列圆心处用实心圆表示axis([n1,n2,0,1.1*max(x)])title('单位抽样序列')xlabel('time(n)')ylabel('Amplitude:x(n)')2、产生单位阶跃序列u(n)n1 = -2;n2 = 8;n0 = 0;n = n1:n2;x = [n>=n0]; % x在n>=n0时为1，其余为0stem(n,x,'filled');axis([n1,n2,0,1.1*max(x)])title('单位阶跃序列')xlabel('time(n)')ylabel('Amplitude:x(n)')3、复指数序列复指数序列的表示式为()(),00,0j n e n x n n σω+⎧≥⎪=⎨<⎪⎩，当0ω=时，()x n 为实指数序列；当0σ=时，()x n 为虚指数序列，即()()cos sin j n e n j n ωωω=+，即其实部为余弦序列，虚部为正弦序列。