数字信号实验6报告

合集下载

数字信号实验报告材料 (全)

数字信号处理实验报告实验一：用 FFT 做谱分析一、实验目的1、进一步加深 DFT 算法原理和基本性质的理解。

2、熟悉 FFT 算法原理和 FFT 子程序的应用。

3、学习用FFT 对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用 FFT 。

二、实验原理用FFT 对信号作频谱分析是学习数字信号处理的重要内容。

经常需要进行谱分析的信号是模拟信号和时域离散信号。

对信号进行谱分析的重要问题是频谱分辨率D 和分析误差。

频谱分辨率直接和FFT 的变换区间N 有关，因为FFT 能够实现的频率分辨率是2π/N ≤D 。

可以根据此时选择FFT 的变换区间N 。

误差主要来自于用FFT 作频谱分析时，得到的是离散谱，而信号（周期信号除外）是连续谱，只有当N 较大时离散谱的包络才能逼近于连续谱，因此N 要适当选择大一些。

周期信号的频谱是离散谱，只有用整数倍周期的长度作FFT ，得到的离散谱才能代表周期信号的频谱。

如果不知道信号周期，可以尽量选择信号的观察时间长一些。

对模拟信号的频谱时，首先要按照采样定理将其变成时域离散信号。

如果是模拟周期信号，也应该选取整数倍周期的长度，经过采样后形成周期序列，按照周期序列的谱分析进行。

三、实验内容和步骤对以下典型信号进行谱分析：⎪⎩⎪⎨⎧≤≤-≤≤-=⎪⎩⎪⎨⎧≤≤-≤≤+==其它nn n n n n x 其它nn n n n n x n R n x ,074,330,4)(,074,830,1)()()(32414()cos4x n n π=5()cos(/4)cos(/8)x n n n ππ=+6()cos8cos16cos20x t t t t πππ=++对于以上信号，x1(n)~x5(n) 选择FFT 的变换区间N 为8和16 两种情况进行频谱分析。

分别打印其幅频特性曲线。

并进行对比、分析和讨论;；x6(t)为模拟周期信号，选择采样频率Hz F s 64=，变换区间N=16,32,64 三种情况进行谱分析。

数字信号处理实验报告_五个实验

实验一信号、系统及系统响应一、实验目的1、熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解；2、熟悉时域离散系统的时域特性；3、利用卷积方法观察分析系统的时域特性；4、掌握序列傅立叶变换的计算机实现方法，利用序列的傅立叶变换对连续信号、离散信号及系统响应进行频域分析。

二、实验原理及方法采样是连续信号数字处理的第一个关键环节。

对采样过程的研究不仅可以了解采样前后信号时域和频域特性发生变化以及信号信息不丢失的条件，而且可以加深对傅立叶变换、Z 变换和序列傅立叶变换之间关系式的理解。

对一个连续信号)(t x a 进行理想采样的过程可用下式表示：)()()(^t p t t xx aa=其中)(^t x a 为)(t x a 的理想采样，p(t)为周期脉冲，即∑∞-∞=-=m nT t t p )()(δ)(^t x a的傅立叶变换为)]([1)(^s m a m j X T j a XΩ-Ω=Ω∑∞-∞=上式表明^)(Ωj Xa为)(Ωj Xa的周期延拓。

其延拓周期为采样角频率（T /2π=Ω）。

只有满足采样定理时，才不会发生频率混叠失真。

在实验时可以用序列的傅立叶变换来计算^)(Ωj X a 。

公式如下：Tw jw ae X j X Ω==Ω|)()(^离散信号和系统在时域均可用序列来表示。

为了在实验中观察分析各种序列的频域特性，通常对)(jw e X 在[0，2π]上进行M 点采样来观察分析。

对长度为N 的有限长序列x(n)，有：n jw N n jw k ke m x eX--=∑=)()(1其中，k Mk πω2=，k=0,1,……M-1 时域离散线性非移变系统的输入/输出关系为 ∑∞-∞=-==m m n h m x n h n x n y )()()(*)()(上述卷积运算也可在频域实现)()()(ωωωj j j e H e X eY =三、实验程序s=yesinput(Please Select The Step Of Experiment:\n 一.(1时域采样序列分析 s=str2num(s); close all;Xb=impseq(0,0,1); Ha=stepseq(1,1,10);Hb=impseq(0,0,3)+2.5*impseq(1,0,3)+2.2*impseq(2,0,3)+impseq(3,0,3); i=0;while(s);%时域采样序列分析 if(s==1) l=1; k=0;while(1)if(k==0)A=yesinput('please input the Amplitude:\n',...444.128,[100,1000]); a=yesinput('please input the Attenuation Coefficient:\n',...222.144,[100,600]); w=yesinput('please input the Angle Frequence(rad/s):\n',...222.144,[100,600]); end k=k+1;fs=yesinput('please input the sample frequence:\n',...1000,[100,1200]); Xa=FF(A,a,w,fs); i=i+1;string+['fs=',num2str(fs)]; figure(i)DFT(Xa,50,string); 1=yesinput 1=str2num(1); end%系统和响应分析else if(s==2)kk=str2num(kk);while(kk)if(kk==1)m=conv(Xb,Hb);N=5;i=i+1;figure(i)string=('hb(n)');Hs=DFT(Hb,4,string);i=i+1;figure(i)string('xb(n)');DFT(Xb,2,string);string=('y(n)=xb(n)*hb(n)');else if (kk==2)m=conv(Ha,Ha);N=19;string=('y(n)=ha(n)*(ha(n)');else if (kk==3)Xc=stepseq(1,1,5);m=conv(Xc,Ha);N=14;string=('y(n)=xc(n)*ha(n)');endendendi=i+1;figure(i)DFT(m,N,string);kk=yesinputkk=str2num(kk);end卷积定理的验证else if(s==3)A=1;a=0.5;w=2,0734;fs=1;Xal=FF(A,a,w,fs);i=i+1;figure(i)string=('The xal(n)(A=1,a=0.4,T=1)'); [Xa,w]DFT(Xal,50,string);i=i+1;figure(i)string =('hb(n)');Hs=DFT(Hb,4,string);Ys=Xs.*Hs;y=conv(Xal,Hb);N=53;i=i+1;figure(i)string=('y(n)=xa(n)*hb(n)');[yy,w]=DFT(y,N,string);i=i+1;figure(i)subplot(2,2,1)plot(w/pi,abs(yy));axis([-2 2 0 2]);xlabel('w/pi');ylabel('|Ys(jw)|');title(FT[x(n)*h(n)]');subplot(2,2,3)plot(w/pi,abs(Ys));axis([-2 2 0 2]);xlabel('w/pi');ylabel('|Ys(jw)|');title('FT[xs(n)].FT[h(n)]');endendend子函数:离散傅立叶变换及X(n),FT[x(n)]的绘图函数function[c,l]=DFT(x,N,str)n=0:N-1;k=-200:200;w=(pi/100)*k;l=w;c=x*Xc=stepseq(1,1,5);子函数：产生信号function c=FF(A,a,w,fs)n=o:50-1;c=A*exp((-a)*n/fs).*sin(w*n/fs).*stepseq(0,0,49); 子函数：产生脉冲信号function [x,n]=impseq(n0,n1,n2)n=[n1:n2];x=[(n-n0)==0];子函数：产生矩形框信号function [x,n]=stepseq(n0,n1,n2) n=[n1:n2];x=[(n-n0>=0)];四、实验内容及步骤1、认真复习采样理论，离散信号与系统，线性卷积，序列的傅立叶变换及性质等有关内容，阅读本实验原理与方法。

数字信号分析实训报告

一、实训目的本次实训旨在通过实际操作，让学生掌握数字信号分析的基本原理和方法，提高学生对数字信号处理技术的理解和应用能力。

通过本次实训，使学生能够：1. 理解数字信号的基本概念、特性及分类；2. 掌握数字信号的采样、量化、编码等基本过程；3. 掌握离散傅里叶变换（DFT）及其快速算法（FFT）；4. 熟悉数字滤波器的设计与实现；5. 能够运用所学知识分析和处理实际信号。

二、实训内容1. 数字信号基本概念及特性（1）数字信号的定义及分类；（2）模拟信号与数字信号的转换过程；（3）数字信号的采样定理及奈奎斯特准则；（4）数字信号的量化及编码；（5）数字信号的时域分析。

2. 离散傅里叶变换（DFT）及其快速算法（FFT）（1）DFT的定义及性质；（2）DFT的计算过程；（3）FFT算法的基本原理及实现方法；（4）DFT与FFT在信号处理中的应用。

3. 数字滤波器的设计与实现（1）滤波器的基本概念及分类；（2）低通、高通、带通、带阻滤波器的设计方法；（3）无限脉冲响应（IIR）滤波器与有限脉冲响应（FIR）滤波器的特点及设计方法；（4）数字滤波器在信号处理中的应用。

4. 实际信号分析（1）选取实际信号进行采集；（2）对采集到的信号进行预处理；（3）运用所学知识对信号进行时域分析、频域分析及滤波处理；（4）对分析结果进行总结。

三、实训过程1. 准备工作（1）熟悉实训设备，了解数字信号分析仪器的功能及操作方法；（2）复习数字信号处理相关理论知识；（3）明确实训目标，制定实训计划。

2. 实训操作（1）数字信号基本概念及特性通过实验，观察不同采样频率下的信号波形，验证采样定理；分析不同量化位数对信号质量的影响；对采集到的信号进行时域分析，了解信号的特性。

（2）离散傅里叶变换（DFT）及其快速算法（FFT）运用DFT算法对信号进行频域分析，观察信号频谱特性；通过FFT算法提高计算效率，实现快速频域分析。

（3）数字滤波器的设计与实现根据实际需求，设计合适的数字滤波器，对信号进行滤波处理；比较IIR滤波器与FIR滤波器的性能差异。

西安交通大学数字信号处理实验报告

数字信号处理实验报告班级：硕姓名：学号：实验1 常见离散信号的MATLAB 产生和图形显示实验目的：加深对常用离散信号的理解；实验内容：（1）单位抽样序列clc;x=zeros(1,11); x(1)=1; n=0:1:10;stem(n,x, 'fill'); title('单位抽样序列'); xlabel('n'); ylabel('x[n]')延迟5个单位：clc;x=zeros(1,11); x(6)=1; n=0:1:10;stem(n,x, 'fill'); title('单位抽样序列'); xlabel('n'); ylabel('x[n]')nx [n ]（2）单位阶跃序列clc;x=[zeros(1,5),ones(1,6)]; n=-5:1:5;stem(n,x,'fill'); title('单位阶跃序列'); xlabel('n'); ylabel('x[n]');nx [n ]（3）正弦序列clc; N=50; n=0:1:N-1; A=1; f=1; Fs=50; fai=pi;x=A*sin(2*pi*f*n/Fs+fai); stem(n,x,'fill'); title('正弦序列'); xlabel('n'); ylabel('x[n]'); axis([0 50 -1 1]);nx [n ]（4）复正弦序列clc; N=50; n=0:1:N-1; w=2*pi/50; x=exp(j*w*n); subplot(2,1,1); stem(n,real(x)); title('复正弦序列实部'); xlabel('n');ylabel('real(x[n])'); axis([0 50 -1 1]); subplot(2,1,2); stem(n,imag(x)); title('复正弦序列虚部'); xlabel('n');ylabel('imag(x[n])'); axis([0 50 -1 1]);nx [n ]（5）指数序列clc; N=10; n=0:1:N-1; a=0.5; x=a.^n;stem(n,x,'fill'); title('指数序列'); xlabel('n'); ylabel('x[n]'); axis([0 10 0 1]);nr e a l (x [n ])ni m a g (x [n ])（6）复指数序列性质讨论：0(j )()enx n σω+=将复指数表示成实部与虚部为00()e cos j sin n n x n n e n σσωω=+1.当σ=0时，它的实部和虚部都是正弦序列。

数字信号处理实验报告

一、实验目的1. 理解数字信号处理的基本概念和原理。

2. 掌握离散时间信号的基本运算和变换方法。

3. 熟悉数字滤波器的设计和实现。

4. 培养实验操作能力和数据分析能力。

二、实验原理数字信号处理（Digital Signal Processing，DSP）是利用计算机对信号进行采样、量化、处理和分析的一种技术。

本实验主要涉及以下内容：1. 离散时间信号：离散时间信号是指时间上离散的信号，通常用序列表示。

2. 离散时间系统的时域分析：分析离散时间系统的时域特性，如稳定性、因果性、线性等。

3. 离散时间信号的变换：包括离散时间傅里叶变换（DTFT）、离散傅里叶变换（DFT）和快速傅里叶变换（FFT）等。

4. 数字滤波器：设计、实现和分析数字滤波器，如低通、高通、带通、带阻滤波器等。

三、实验内容1. 离散时间信号的时域运算（1）实验目的：掌握离散时间信号的时域运算方法。

（2）实验步骤：a. 使用MATLAB生成两个离散时间信号；b. 进行时域运算，如加、减、乘、除等；c. 绘制运算结果的时域波形图。

2. 离散时间信号的变换（1）实验目的：掌握离散时间信号的变换方法。

（2）实验步骤：a. 使用MATLAB生成一个离散时间信号；b. 进行DTFT、DFT和FFT变换；c. 绘制变换结果的频域波形图。

3. 数字滤波器的设计和实现（1）实验目的：掌握数字滤波器的设计和实现方法。

（2）实验步骤：a. 设计一个低通滤波器，如巴特沃斯滤波器、切比雪夫滤波器等；b. 使用MATLAB实现滤波器；c. 使用MATLAB对滤波器进行时域和频域分析。

4. 数字滤波器的应用（1）实验目的：掌握数字滤波器的应用。

（2）实验步骤：a. 采集一段语音信号；b. 使用数字滤波器对语音信号进行降噪处理；c. 比较降噪前后的语音信号，分析滤波器的效果。

四、实验结果与分析1. 离散时间信号的时域运算实验结果显示，通过MATLAB可以方便地进行离散时间信号的时域运算，并绘制出运算结果的时域波形图。

数字信号处理实验报告

数字信号处理实验报告一、实验目的本次数字信号处理实验的主要目的是通过实际操作和观察，深入理解数字信号处理的基本概念和方法，掌握数字信号的采集、处理和分析技术，并能够运用所学知识解决实际问题。

二、实验设备与环境1、计算机一台，安装有 MATLAB 软件。

2、数据采集卡。

三、实验原理1、数字信号的表示与采样数字信号是在时间和幅度上都离散的信号，可以用数字序列来表示。

在采样过程中，根据奈奎斯特采样定理，为了能够准确地恢复原始信号，采样频率必须大于信号最高频率的两倍。

2、离散傅里叶变换（DFT）DFT 是将时域离散信号变换到频域的一种方法。

通过 DFT，可以得到信号的频谱特性，从而分析信号的频率成分。

3、数字滤波器数字滤波器是对数字信号进行滤波处理的系统，分为有限冲激响应（FIR）滤波器和无限冲激响应（IIR）滤波器。

FIR 滤波器具有线性相位特性，而 IIR 滤波器则在性能和实现复杂度上有一定的优势。

四、实验内容与步骤1、信号的采集与生成使用数据采集卡采集一段音频信号，或者在 MATLAB 中生成一个模拟信号，如正弦波、方波等。

2、信号的采样与重构对采集或生成的信号进行采样，然后通过插值算法重构原始信号，观察采样频率对重构信号质量的影响。

3、离散傅里叶变换对采样后的信号进行DFT 变换，得到其频谱，并分析频谱的特点。

4、数字滤波器的设计与实现（1）设计一个低通 FIR 滤波器，截止频率为给定值，观察滤波前后信号的频谱变化。

（2）设计一个高通 IIR 滤波器，截止频率为给定值，比较滤波前后信号的时域和频域特性。

五、实验结果与分析1、信号的采集与生成成功采集到一段音频信号，并在MATLAB 中生成了各种模拟信号，如正弦波、方波等。

通过观察这些信号的时域波形，对不同类型信号的特点有了直观的认识。

2、信号的采样与重构当采样频率足够高时，重构的信号能够较好地恢复原始信号的形状；当采样频率低于奈奎斯特频率时，重构信号出现了失真和混叠现象。

数字信号处理实验报告完整版[5篇模版]

数字信号处理实验报告完整版[5篇模版]第一篇：数字信号处理实验报告完整版实验 1利用 T DFT 分析信号频谱一、实验目的1.加深对 DFT 原理的理解。

2.应用 DFT 分析信号的频谱。

3.深刻理解利用DFT 分析信号频谱的原理，分析实现过程中出现的现象及解决方法。

二、实验设备与环境计算机、MATLAB 软件环境三、实验基础理论T 1.DFT 与与 T DTFT 的关系有限长序列的离散时间傅里叶变换在频率区间的N 个等间隔分布的点上的 N 个取样值可以由下式表示：212 /0()|()()0 1Nj knjNk NkX e x n e X k k Nπωωπ--====≤≤-∑由上式可知，序列的 N 点 DFT ,实际上就是序列的 DTFT 在 N 个等间隔频率点上样本。

2.利用 T DFT 求求 DTFT方法 1 1：由恢复出的方法如下：由图 2.1 所示流程可知：101()()()Nj j n kn j nNn n kX e x n e X k W eNωωω∞∞----=-∞=-∞=⎡⎤==⎢⎥⎣⎦∑∑∑由上式可以得到：IDFT DTFT第二篇：数字信号处理实验报告JIANGSUUNIVERSITY OF TECHNOLOGY数字信号处理实验报告学院名称：电气信息工程学院专业：班级：姓名：学号：指导老师：张维玺（教授）2013年12月20日实验一离散时间信号的产生一、实验目的数字信号处理系统中的信号都是以离散时间形态存在的，所以对离散时间信号的研究是数字信号的基本所在。

而要研究离散时间信号，首先需要产生出各种离散时间信号。

使用MATLAB软件可以很方便地产生各种常见的离散时间信号，而且它还具有强大绘图功能，便于用户直观地处理输出结果。

通过本实验，学生将学习如何用MATLAB产生一些常见的离散时间信号，实现信号的卷积运算，并通过MATLAB中的绘图工具对产生的信号进行观察，加深对常用离散信号和信号卷积和运算的理解。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

广州大学学生实验报告
开课学院及实验室：电子楼317EDA 2012年5月13日
学院机械与电气
工程学院
年级、
专业、
班
信工101 姓名齐同远学号1007100031
实验课
程名称
数字信号处理实验成绩
实验项目名称IIR滤波器的设计
指导教
师
张承云
一、实验目的
设计IIR滤波器，实现对存在加性干扰的时域离散信号进行滤波。

二、实验内容和步骤
已知带加性干扰的信号
()
x n表示，()()()
s
x n x n n
η
=+
，式中
()
s
x n
是有用的信号，是一个0~0.2πrad
的带限信号。

()n
η
是一个干扰信号，其频谱分布在0.3πrad以上。

要求设计一个巴特沃斯IIR数字滤
波器对信号
()
x n进行滤波，将干扰()n
η
滤除。

要求在
()
s
x n
所在的通带内滤波器幅度平坦，在0.2πrad
处幅度衰减不大于1dB，在噪声所在的0.3πrad以上的频带内滤波器幅度衰减大于等于40dB。

程序：
%-------------------------------------IIR滤波器设计------------------------------------- T = 2;%采样周期T取2
wp = 0.2 * pi;%设计的通带数字边界频率
ws = 0.3 * pi;%设计的阻带数字边界频率
Wp = 2 / T * tan(wp /2);%通带边界频率数字滤波器指标转换为模拟滤波器指标
Ws = 2 / T * tan(ws /2);%阻带边界频率数字滤波器指标转换为模拟滤波器指标
Ap = 1;%确定通带最大衰减
As = 10;%确定阻带最小衰减
E = sqrt(10^(Ap / 10) - 1);%计算波纹幅度参数E
A = 10^(As / 20);%计算波纹幅度参数A
N = ceil(log10(E / sqrt(A^2 - 1)) / log10(Wp / Ws));%计算选用滤波器阶数
Wc = Ws / (A^2 - 1)^(1 / (2 * N));%计算3dB截止频率
syms p s z ;%定义中间变量
Gp = 1 / ((p^2 + 0.7654 * p + 1) * (p^2 + 1.8478 * p + 1));%进行测试后得知阶数为4，查表得
多项式系数
Has = subs(Gp,p,s / Wc);%去归一化
Hz = subs(Has,s,(1 - z^-1) / (1 + z^-1));%双线性变换法转换
%---------------------------采用Matlab内置函数设计IIR滤波器--------------------------- Matlab_wp = 0.2;%对π归一化的设计的通带数字边界频率
Matlab_ws = 0.3;%对π归一化的设计的阻带数字边界频率
rp = 1;%确定通带最大衰减r
rs = 10;%确定阻带最小衰减
[Matlab_N,Matlab_wc] = buttord(Matlab_wp,Matlab_ws,rp,rs);%计算数字滤波器N和3dB截止频率[B,A] = butter(Matlab_N,Matlab_wc);%设计数字滤波器
Matlab_Hz = (B(1) + B(2)*z^-1 + B(3)*z^-2 + B(4)*z^-3) / (A(1) + A(2) * z^-1 + A(3) * z^-2 + A(4) * z^-3); %设计数字滤波器
figure('Name','自己编写的双线性变换法得出的滤波器');ezplot(abs(Hz),[0,3]); %画图
axis tight;%调整画图范围
xlabel('归一化角频率');%设置X轴
ylabel('幅度');%设置Y轴
title('自己编写的双线性变换法得出的滤波器');%设置标题
figure('Name','Matlab的双线性变换法得出的滤波器');ezplot(abs(Matlab_Hz),[0,3]); %画图
axis tight;% 调整画图范围
xlabel('归一化角频率');%设置X轴
ylabel('幅度');%设置Y轴
title('Matlab的双线性变换法得出的滤波器');%设置标题
明显，得到的两个滤波器的幅频特性不同，原因是：采样周期T的选取不同。

如果把参数T = 2;这一句改成T = 100;观察结果。

这样，两者的形状就接近了，但是幅度还是有差距，成比例关系。

数字信号实验6报告

数字信号实验报告材料 (全)

数字信号处理实验报告_五个实验

数字信号分析实训报告

最新数字信号处理实验报告

西安交通大学数字信号处理实验报告

数字信号处理实验报告

数字信号处理实验报告

数字信号处理实验报告完整版[5篇模版]