测试信号处理实验五
信号分析与处理实验报告
华北电力大学实验报告||实验名称FFT的软件实现实验（Matlab）IIR数字滤波器的设计课程名称信号分析与处理||专业班级：电气化1308 学生姓名：袁拉麻加学号： 2 成绩：指导教师：杨光实验日期： 2015-12-17快速傅里叶变换实验一、实验目的及要求通过编写程序，深入理解快速傅里叶变换算法（FFT）的含义，完成FFT和IFFT算法的软件实现。
二、实验内容利用时间抽取算法，编写基2点的快速傅立叶变换（FFT）程序；并在FFT程序基础上编写快速傅里叶反变换（IFFT）的程序。
三：实验要求1、FFT和IFFT子程序相对独立、具有一般性，并加详细注释；2、验证例6-4，并能得到正确结果。
3、理解应用离散傅里叶变换（DFT）分析连续时间信号频谱的数学物理基础。
四、实验原理：a.算法原理1、程序输入序列的元素数目必须为2的整数次幂，即N=2M,整个运算需要M 级蝶形运算；2、输入序列应该按二进制的码位倒置排列，输出序列按自然序列排列；3、每个蝶形运算的输出数据军官占用其他输入数据的存储单元，实现“即位运算”；4、每一级包括N/2个基本蝶形运算，共有M*N/2个基本蝶形运算；5、第L级中有N/2L个群，群与群的间隔为2L。
6、处于同一级的各个群的系数W分布相同，第L级的群中有2L-1个系数；7、处于第L级的群的系数是(p=1,2,3，…….,2L-1)而对于第L级的蝶形运算，两个输入数据的间隔为2L-1。
b.码位倒置程序流程图开始检测A序列长度nk=0j=1x1(j)=bitget(k,j);j=j+1Yj<m?Nx1=num2str(x1);y(k+1)=bin2dec(x1);clear x1k=k+1c.蝶形运算程序流程图五、程序代码与实验结果a.FFT程序：%%clear all;close all;clc;%输入数据%A=input('输入x(n)序列','s');A=str2num(A);% A=[1,2,-1,4]; %测试数据%%%%校验序列,%n=length(A);m=log2(n);if (fix(m)~=m)disp('输入序列长度错误，请重新输入!');A=input('输入x(n)序列','s');A=str2num(A);elsedisp('输入正确，请运行下一步')end%%%码位倒置%for k=0:n-1for j=1:m %取M位的二进制数%x1(j)=bitget(k,j); %倒取出二进制数%endx1=num2str(x1); %将数字序列转化为字符串%y(k+1)=bin2dec(x1); %二进制序列转化为十进制数%clear x1endfor k=1:nB(k)=A(y(k)+1); %时间抽取序列%endclear A%%%计算%for L=1:m %分解为M级进行运算%LE=2^L; %第L级群间隔为2^L%LE1=2^(L-1); %第L级中共有2^(L-1)个Wn乘数，进行运算蝶运算的两数序号相隔LE1%W=1;W1=exp(-1i*pi/LE1);for R=1:LE1 %针对第R个Wn系数进行一轮蝶运算，共进行LE1次%for P=R:LE:n %每个蝶的大小为LE% Q=P+LE1;T=B(Q)*W;B(Q)=B(P)-T;B(P)=B(P)+T;endW=W*W1;endendB %输出X(k)%%%验证结果：例6-4b.IFFT程序：%%clear all;close all;clc;%输入数据%A=input('输入X(k)序列','s');A=str2num(A);% A=[6,2+2i,-6,2-2i]; %测试数据%%%%校验序列,%n=length(A);m=log2(n);if (fix(m)~=m)disp('输入序列长度错误，请重新输入!');A=input('输入x(n)序列','s');A=str2num(A);elsedisp('输入正确，请运行下一步')end%%%码位倒置%for k=0:n-1for j=1:m %取M位的二进制数%x1(j)=bitget(k,j); %倒取出二进制数%endx1=num2str(x1); %将数字序列转化为字符串%y(k+1)=bin2dec(x1); %二进制序列转化为十进制数%clear x1endfor k=1:nB(k)=A(y(k)+1); %时间抽取序列%endclear A%%%计算%for L=1:m %分解为M级进行运算%LE=2^L; %第L级群间隔为2^L%LE1=2^(L-1); %第L级中共有2^(L-1)个Wn乘数，进行运算蝶运算的两数序号相隔LE1%W=1;W1=exp(-1i*pi/LE1);for R=1:LE1 %针对第R个Wn系数进行一轮蝶运算，共进行LE1次%for P=R:LE:n %每个蝶的大小为LE%Q=P+LE1;T=B(Q)*W;B(Q)=B(P)-T;B(P)=B(P)+T;endW=W*W1;endendB=conj(B); %取共轭%B=B/n %输出x(n)%验证结果：六、实验心得与结论本次实验借助于Matlab软件，我避开了用C平台进行复杂的复数运算，在一定程度上简化了程序，并添加了简单的检错代码，码位倒置我通过查阅资料，使用了一些函数，涉及到十-二进制转换，数字-文本转换，二-文本转换，相对较复杂，蝶运算我参考了书上了流程图，做些许改动就能直接实现。
实验五 信号的采样与恢复
F
0
t
(a) 连续信号的频谱
m
0
m
f s t
Fs
t
0 Ts
1 TS
s
0
m
m
s
（b） 高抽样频率时的抽样信号及频谱（不混叠）
f s t
Fs
1 TS
0 Ts
t
s
0
m
பைடு நூலகம்
m
s
（c） 低抽样频率时的抽样信号及频谱（混叠） 图 5-2 抽样过程中出现的两种情况 4、点频抽样还原实验采用分立方式，对2kHz正弦波进行抽样和还原，首先2kHz的 方波经过截止频率为2.56kHz低通滤波器得到2kHz的正弦波，然后用可调窄脉冲对正弦 波进行抽样得到抽样信号，抽样信号经低通滤波器后还原出正弦波。 考虑下面的正弦信号： x(t ) cos(
s
2
2
t)
该信号在采样周期 2
s 整数倍点上
的值都是零；因此在这个采样频率下所产生的信号全是零。当这个零输入加到理想低通 滤波器上时，所得输出当然也都是零。 5、为了实现对连续信号的抽样和抽样信号的复原，除选用足够高的抽样频率外， 常采用前置低通滤波器来防止原信号频谱宽而造成抽样后信号频谱的混叠。但这也会造 成失真。原始的语音信号带宽为 40Hz 到 10000Hz,但实际中传输的语音信号的带宽为 300Hz 到 3400Hz,并不影响我们的听觉效果,因此本实验加了前置滤波器。 6、语音抽样还原实验采用集成方式，本实验采用PCM编译码器TP3067专用大规模集 成电路，它是CMOS工艺制造的单片PCM A律编译码器．片内带有输入输出话路滤波器．它 把编译码器（Codec）和滤波器(Filter)集成在一个芯片上。 脉冲编码调制（PCM)就是把一个时间连续、取值连续的模拟信号变换成时间离散、 取值离散的数字信号后在信道中进行传输。而脉冲编码调制就是对模拟信号先进行抽样 后，再对样值的幅度进行量化、编码的过程。话音信号先经过防混叠低通滤波器，得到 限带信号（300Hz～3400Hz），进行脉冲抽样，变成 8kHz 重复频率的抽样信号（即离散 的脉冲调幅 PAM 信号），然后将幅度连续的 PAM 信号用“四舍五入”办法量化为有限个 幅度取值的信号，再经编码，转换成二进制码。对于电话，CCITT（国际电话与电报顾 问委员会 International Telephone and Telegraph Consultative Committee）规定 8 抽样率为 8kHz，每抽样值编 8 位码，即共有 2 =256 个量化值，因而每话路 PCM 编码后 的标准数码率是 64kb/s。
数字信号处理实验五
实验五：抽样定理一、实验目的1、了解用MA TLAB 语言进行时域、频域抽样及信号重建的方法。
2、进一步加深对时域、频域抽样定理的基本原理的理解。
3、观察信号抽样与恢复的图形，掌握采样频率的确定方法和内插公式的编程方法。
二、实验内容及步骤1、阅读并输入实验原理中介绍的例题程序，观察输出的数据和图形，结合基本原理理解每一条语句的含义。
2、已知一个连续时间信号f(t)=sinc(t)，取最高有限带宽频率f m =1Hz 。
（1）分别显示原连续信号波形和F s =f m 、F s =2f m 、F s =3f m 三种情况下抽样信号的波形；dt=0.1;f0=1;T0=1/f0; fm=1;Tm=1/fm; t=-2:dt:2; f=sinc(t);subplot(4,1,1);plot(t,f);axis([min(t),max(t),1.1*min(f),1.1*max(f)]); title('原连续信号和抽样信号'); for i=1:3;fs=i*fm;Ts=1/fs; n=-2:Ts:2; f=sinc(n);subplot(4,1,i+1);stem(n,f,'filled');axis([min(n),max(n),1.1*min(f),1.1*max(f)]);课程名称 数字信号处理 实验成绩 指导教师实 验 报 告院系 信息工程学院 班级 学号 姓名 日期end-2-1.5-1-0.50.511.5200.51原连续信号和抽样信号（2）求解原连续信号和抽样信号的幅度谱； dt=0.1;f0=1;T0=1/f0; fm=1;Tm=1/fm; t=-2:dt:2; N=length(t); f=sinc(t); wm=2*pi*fm; k=0:N-1; w1=k*wm/N; F1=f*exp(-j*t'*w1)*dt;subplot(4,1,1);plot(w1/(2*pi),abs(F1));axis([0,max(4*fm),1.1*min(abs(F1)),1.1*max(abs(F1))]); for i=1:3;if i<=2 c=0;else c=1;end fs=(i+c)*fm;Ts=1/fs; n=-2:Ts:2; N=length(n); f=sinc(n); wm=2*pi*fs; k=0:N-1; w=k*wm/N; F=f*exp(-j*n'*w)*Ts;subplot(4,1,i+1);plot(w/(2*pi),abs(F));axis([0,max(4*fm),1.1*min(abs(F)),1.1*max(abs(F))]); end00.511.522.533.540.20.40.60.811.200.511.522.533.54012（3）用时域卷积的方法（内插公式）重建信号。
数字信号管理方案计划实验报告实验五
物理与电子信息工程学院实验报告实验课程名称：数字信号处理实验名称：FIR数字滤波器设计与软件实现班级：1012341姓名：严娅学号：101234153成绩：_______实验时间：2012年12月20 日一、实验目的（1）掌握用窗函数法设计FIR 数字滤波器的原理和方法。
（2）掌握用等波纹最佳逼近法设计FIR 数字滤波器的原理和方法。
（3）掌握FIR 滤波器的快速卷积实现原理。
（4）学会调用MATLAB 函数设计与实现FIR 滤波器。
二、实验原理1、用窗函数法设计FIR 数字滤波器的原理和方法。
如果所希望的滤波器的理想频率响应函数为 )(ωj d e H ，则其对应的单位脉冲响应为)(n h d =π21ωωωππd e e H j j d )(⎰- （2-1）窗函数设计法的基本原理是用有限长单位脉冲响应序列)(n h 逼近)(n h d 。
由于)(n h d 往往是无限长序列，且是非因果的，所以用窗函数)(n ω将)(n h d 截断，并进行加权处理，得到：)(n h ＝)(n h d )(n ω (2-2))(n h 就作为实际设计的FIR 数字滤波器的单位脉冲响应序列，其频率响应函数)(ωj d e H 为：)(ωj d e H ＝∑-=-1)(N n j e n h ω (2-3) 式中，N 为所选窗函数)(n ω的长度。
由第七章可知，用窗函数法设计的滤波器性能取决于窗函数)(n ω的类型及窗口长度N 的取值。
设计过程中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N 。
各种类型的窗函数可达到的阻带最小衰减和过渡带宽度见第七章。
这样选定窗函数类型和长度N 后，求出单位脉冲响应)(n h ＝)(n h d ·)(n ω，并按式（2-3）求出)(ωj e H 。
)(ωj e H 是否满足要求，要进行验算。
一般在)(n h 尾部加零使长度满足于2的整数次幂，以便用FFT 计算)(ωj e H 。
《测试信号分析报告与处理》实验报告材料
《测试信号分析与处理》实验一差分方程、卷积、z变换一、实验目的通过该实验熟悉 matlab软件的基本操作指令，掌握matlab软件的使用方法，掌握数字信号处理中的基本原理、方法以及matlab函数的调用。
二、实验设备1、微型计算机1台；2、matlab软件1套三、实验原理Matlab 软件是由mathworks公司于1984年推出的一套科学计算软件，分为总包和若干个工具箱，其中包含用于信号分析与处理的sptool工具箱和用于滤波器设计的fdatool工具箱。
它具有强大的矩阵计算和数据可视化能力，是广泛应用于信号分析与处理中的功能强大且使用简单方便的成熟软件。
Matlab软件中已有大量的关于数字信号处理的运算函数可供调用，本实验主要是针对数字信号处理中的差分方程、卷积、z变换等基本运算的matlab函数的熟悉和应用。
差分方程（difference equation)可用来描述线性时不变、因果数字滤波器。
用x表示滤波器的输入，用y表示滤波器的输出。
a0y[n]+a1y[n-1]+…+a N y[n-N]=b0x[n]+b1x[n-1]+…+b M x[n-M] (1)ak,bk 为权系数，称为滤波器系数。
N为所需过去输出的个数，M 为所需输入的个数卷积是滤波器另一种实现方法。
y[n]= ∑x[k] h[n-k] = x[n]*h[n] (2) 等式定义了数字卷积，*是卷积运算符。
输出y[n] 取决于输入x[n] 和系统的脉冲响应h[n]。
传输函数H(z)是滤波器的第三种实现方法。
H(z)=输出/输入= Y(z)/X(z) （3）即分别对滤波器的输入和输出信号求z变换，二者的比值就是数字滤波器的传输函数。
序列x[n]的z变换定义为X (z)=∑x[n]z-n (4) 把序列x[n] 的z 变换记为Z{x[n]} = X(z)。
由X(z) 计算x[n] 进行z 的逆变换x[n] = Z-1{X(z)}。
数字信号处理实验五谱分析
用FFT 对信号作频谱分析1．实验目的学习用FFT 对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析 误差及其原因，以便正确应用FFT 。
2. 实验原理用FFT 对信号作频谱分析是学习数字信号处理的重要内容。
经常需要进行谱分析的信号是模拟信号和时域离散信号。
对信号进行谱分析的重要问题是频谱分辨率D 和分析误差。
频谱分辨率直接和FFT 的变换区间N 有关，因为FFT 能够实现的频率分辨率是N /2π，因此要求D N ≤/2π。
可以根据此式选择FFT 的变换区间N 。
误差主要来自于用FFT 作频谱分析时，得到的是离散谱，而信号（周期信号除外）是连续谱，只有当N 较大时离散谱的包络才能逼近于连续谱，因此N 要适当选择大一些。
周期信号的频谱是离散谱，只有用整数倍周期的长度作FFT ，得到的离散谱才能代表周期信号的频谱。
如果不知道信号周期，可以尽量选择信号的观察时间长一些。
对模拟信号进行谱分析时，首先要按照采样定理将其变成时域离散信号。
如果是模拟周期信号，也应该选取整数倍周期的长度，经过采样后形成周期序列，按照周期序列的谱分析进行。
3．实验步骤及内容（1）对以下序列进行谱分析。
⎪⎩⎪⎨⎧≤≤-≤≤-=⎪⎩⎪⎨⎧≤≤-≤≤+==其它nn n n n n x 其它n n n n n n x n R n x ,074,330,4)(,074,830,1)()()(3241 选择FFT 的变换区间N 为8和16 两种情况进行频谱分析。
分别打印其幅频特性曲线。
并进行对比、分析和讨论。
（2）对以下周期序列进行谱分析。
4()cos 4x n n π=5()cos(/4)cos(/8)x n n n ππ=+选择FFT 的变换区间N 为8和16 两种情况分别对以上序列进行频谱分析。
分别打印其幅频特性曲线。
并进行对比、分析和讨论。
（3）对模拟周期信号进行谱分析6()cos8cos16cos20x t t t t πππ=++选择 采样频率Hz F s 64=，变换区间N=16,32,64 三种情况进行谱分析。
测试技术与信号处理试验报告
测试技术与信号处理实验报告学院：班级：学号：姓名：指导老师：实验二金属箔式应变片——半桥性能实验实验目的：比较半桥与单臂电桥的不同性能，了解其特点。
二、基本原理：不同受力方向的两片应变片接入电桥作为邻边，电桥输出灵敏度提高，非线性得到改善。
当两片应变片阻值和应变量相同时，其桥路输出电压U02=EK/ε2。
三、需用器件与单元：应变式传感器实验模板、应变式传感器、砝码、直流电压表、±15V电源、±4V电源、万用表(自备)。
四、实验步骤：1、根据图(2-1)应变式传感器的插头插入应变传感器模块（Ti）上。
传感器中各应变片就接入了模板。
的左上方的R1、R2、R3、R4。
加热丝贴在应变传感器上，用时插入+5V直流电源，可用万用表进行测量判别，R1=R2=R3=R4=350Ω，加热丝阻值为50Ω左右。
图2-12、接入模板电源±15V(从主控台引入)，检查无误后，合上主控台电源开关，将实验模板调节增益电位器R w3顺时针调节大致到中间位置，再进行差动放大器调零，方法为将差放的正、负输入端与地短接，V o1与Vi2连接，输出V o2与主控台面板上直流电压表输入端+相连，调节实验模板上调零电位器RW4使直流电压表显示为零(直流电压表的切换开关打到2V档)。
关闭主控台电源。
3、根据图2-2接线。
R1、R2为实验模板左上方的应变片，注意R2应和R1受力状态相反，即将传感器中两片受力相反(一片受拉、一片受压)的电阻应变片作为电桥的相邻边。
接入桥路电源±4V，调节电桥调零电位器R w1进行桥路调零。
图 2-24\在电子称上放置一只砝码，读取直流电压表数值，依次增加砝码和读取相应的直流电压表值，直到200g砝码加完。
记下实验结果填入表2-1，关闭电源。
根据表2-1计算系统灵敏度S：S=Δu/ΔW(Δu输出电压变化量；ΔW重量变化量)；计算非线性误差：δf1=Δm/y F·S×100%式中Δm为输出电压值(多次测量时为平均值)与拟合直线的最大偏差；y F·S为满量程输出平均值。
测试信号实验——模拟滤波及数字滤波报告
信号的调理与滤波器设计实验报告一、实验目的掌握模拟滤波器的设计方法和实现过程；掌握数字滤波器的设计方法和实现过程。
二、实验原理在信号传感和传输过程中，由于热噪声、漏电流和电源干扰等因素的影响，不可避免地会有干扰信号叠加到有用信号上，当这种干扰信号非常强时，将严重影响有用信号的识别和利用，因而，通常都有必要对这些干扰信号进行滤波处理。
干扰信号按照频谱分布可分为低频、中频和高频信号，因而，滤波器也相应设计成高通、带通、低通和带阻等形式，具体的滤波器原理和设计方法可参考模拟电子技术和其它相关资料。
在本实验中，要求在对干扰信号频谱分析的基础上，确定滤波器的形式，设计滤波器的截止频率和具体的RC参数，实现对干扰信号的抑制，通过对滤波后信号的时频域分析，评估滤波效果。
三、实验仪器1、电子称1台2、万用表1个3、采集卡1块4、面包板1块5、计算机1台6、信号发生器1台7、Labview软件1套8、运算放大器若干片9、电阻、电容等若干四、实验内容和步骤1、数字滤波器设计：①将电子称、电源、万用表、噪声发生器、采集卡和计算机连接，构成一个完整的测试系统；②利用Labview软件对采集到的信号进行频谱分析，判断干扰信号的频谱分布特征；③根据干扰信号的频谱分布特征进行滤波器的设计，并在面包板上实现；④利用Labview软件对加入滤波器的采集信号进行频谱分析，判断滤波后的干扰信号被抑制的情况，并评价滤波器的功效，如果滤波效果不好，分析具体原因，进一步改进滤波器，直至滤波效果达到预期要求；⑤改变干扰噪声的频率，比较滤波效果，并重新设计滤波器，重复2～4步骤。
2、模拟滤波器设计：①将信号发生器的噪声信号叠加到表示电子称输出的信号上；②将叠加了噪声的信号连接到数据采集卡的接口板上；③利用labview将信号采集到计算机中；④分析信号的频谱，得到信号的幅度谱；⑤根据信号特点提出滤波器设计参数、截止频率；⑥设计出滤波器的传递函数；⑦根据滤波器传递函数设计电路，完成电路的搭接；⑧将滤波器的输出送到采集卡，用计算机程序求出重物重量。
数字信号处理实验1--5含代码
数字信号处理实验1--5含代码实验一离散时间信号的时域分析 1. 在MATLAB中利用逻辑关系式n,,0来实现序列，显示范围。
(产生如下，，,n,nn,n,n012图所示的单位脉冲信号的函数为impseq(n0,n1,n2)，程序如示例所示),3,n,10并利用impseq函数实现序列:; ，，，，，，yn,2,n,3，,n,6，，xn1nnnn120源代码:impseq.mfunction y=impseq(n0,n1,n2)n=[n1:n2]y=[(n-n0)==0]exp01-1.mfunction impseq(n0,n1,n2)n=-3:1:10y=2*impseq(3,-3,10)+impseq(6,-3,10);stem(n,y)n,,0，，2. 在MATLAB中利用逻辑关系式来实现序列，显示范围。
(自己编写un,nn,n,n012产生单位阶跃信号的函数，函数命名为stepseq(n0,n1,n2)) 并利用编写的stepseq函数实现序列: ，，，，，，yn,un，2，un,2,5,n,10源代码:stepseq.mfunction y=stepseq(n0,n1,n2)n=n1:1:n2y=[(n-n0)>=0]exp01-2.mfunction stepseq(n0,n1,n2)n=-5:1:20y=stepseq(-2,-5,20)+stepseq(2,-5,20)stem(n,y)3. 在MATLAB中利用数组运算符“.^”来实现一个实指数序列。
如: n ，，，，xn,0.30,n,15源代码:n=0:1:15;x=0.3.^nstem(n,x)4. 在MATLAB中调用函数sin或cos产生正余弦序列，如:π,, ，，，，xn,3sin0.4πn，，5cos0.3πn0,n,20,,5,,源代码:n=0:1:20x=11*sin(0.3*pi*n+pi/5)+5*cos(0.3*pi*n)stem(n,x)思考题:1.在MATLAB环境下产生单位脉冲序列和单位阶跃序列各有几种方法,如何使用,2.在MATLAB环境下进行序列的相乘运算时应注意什么问题,实验二离散时间系统的时域分析1. 在MATLAB中利用内部函数conv来计算两个有限长序列的卷积。
信号分析与处理实验五
实验五：FIR数字滤波器设计与软件实现一、实验目的（略）二、实验内容及步骤（略）三、相关程序1、信号产生函数xtg清单Function xt=xtg(N) %信号x(t)产生函数，并显示信号的幅频特性曲线%xt=xtg产生一个长度为N，有加性高频噪声的单频调幅信号xt，采样频率Fs=1khz%载波频率fc=Fs/10=100Hz，调制正弦波频率f0=Fc/10=10HzN=2000;Fs=1000;T=1/Fs;Tp=N*T;T=0:T:(N-1)*T;fc=Fs/10;f0=fc/10; %载波频率fc=Fs/10，单频调制信号频率为f0=fc/10mt=cos(2*pi*f0*t); %产生单频正弦波调制信号mt，频率为f0ct=cos(2*pi*fc*t); %产生载波正弦波信号ct，频率为fcxt=mt.*ct; %相乘产生单频调制信号xtnt=2*rand(1,N)-1; %产生随机噪声nt%=========设计高通滤波器hn，用于滤除噪声nt中的低频成分，生成高频噪声===== fp=150;fs=200;rp=0.1;rs=70; %滤波器指标fb=[fp,fs];m=[0,1]; %计算remenzord函数所需参数f，m，devdev=[10^(-rs/20),(10^(rp/20)-1)/ (10^(rp/20)+1)];[n,fo,mo,w]=remenzord(fb,m,dev,Fs); %确定remenz函数所需参数hn=remenz(n,fo,mo,w); %调用remenz函数进行设计，用于滤除噪声nt中的低频成分yt=filter(hn,1,10*nt); %滤除随机噪声中低频成分，生成高频噪声yt%=================================================================== xt=xt+yt; %噪声加信号fst=fft(xt,N);k=0:N-1;f=k/Tp;subplot(3,1,1);plot(t,xt);grid;xlabel([0,Tp/5,min(xt),max(xt)]);title(‘(a)信号加噪声波形’)subplot(3,1,2);plot(f,abs(fst)/max(abs(fst)));grid;title(‘(b)信号加噪声的频谱’)axis([0,Fs/2,0,1.2]);xlabel(‘f/Hz’);ylabel(‘幅度’)2、实验程序清单%实验五程序exp5.m%FIR数字滤波器设计与软件实现clear all;close all;%=====调用xtg产生信号xt，xt长度N=1000，并显示xt及其频谱===========N=1000;xt=xtg(N);fp=120;fs=150;rp=0.2;rs=60; %输入给定指标Fs=1000; T=1/Fs;Tp=N*T; k=0:N-1;f=k/Tp; t=0:T:(N-1)*T;%(1)用窗函数法设计滤波器wc=(fp+fs)/Fs; %理想低通滤波器截止频率（关于pi归一化）B=2*pi*(fs-fp)/Fs; %过渡带宽度指标Nb=ceil(11*pi/B); %blackman窗的长度Nhn=fir(Nb-1,wc,blackman(Nb));Hw=abs(fft(hn,1024)); %求设计的滤波器频率特性ywt=fftfilt(hn,xt,N); %调用函数fftfilt对xt滤波subplot(3,1,1);plot(f,Hw/max(Hw)); grid;title(‘(a)低通滤波器幅频特性’) axis([0,Fs/2,-120,20]);xlabel(‘f/Hz’);ylabel(‘幅度’)subplot(3,1,2);plot(t,ywt); grid;title(‘(b)滤除噪声后的信号波形’)axis([0,0.5,-1,1]);xlabel(‘t/s’);ylabel(‘yw(t)’)%=======(2)用等波纹逼近法设计滤波器======================== fb=[fp,fs];m=[1,0]; %确定remenzord函数所需参数f,m,devdev=[(10^(rp/20)-1)/(10^(rp/20)+1),10^(-rs/20)];[Ne,fo,mo,W]=remenzord(fb,m,dev,fs); %确定remenz函数所需参数hn=remenz(Ne,fo,mo,W); %调用remenz函数进行设计Hw=abs(fft(hn,1024)); %求设计的滤波器频率特性yet=fftfilt(hn,xt,N); %调用函数fftfilt对xt滤波figure(2)subplot(3,1,1);plot(f,Hw/max(Hw)); grid;title(‘(a)低通滤波器幅频特性’) axis([0,Fs/2,-80,10]);xlabel(‘f/Hz’);ylabel(‘幅度’)subplot(3,1,2);plot(t,yet); grid;title(‘(b)滤除噪声后的信号波形’)axis([0,0.5,-1,1]);xlabel(‘t/s’);ylabel(‘ye(t)’)。
数字信号处理实验报告实验五
物理与电子信息工程学院实验报告实验课程名称：数字信号处理实验名称：FIR数字滤波器设计与软件实现班级：1012341姓名：严娅学号：101234153成绩：_______实验时间：2012年12月20 日一、实验目的（1）掌握用窗函数法设计FIR 数字滤波器的原理和方法。
（2）掌握用等波纹最佳逼近法设计FIR 数字滤波器的原理和方法。
（3）掌握FIR 滤波器的快速卷积实现原理。
（4）学会调用MATLAB 函数设计与实现FIR 滤波器。
二、实验原理1、用窗函数法设计FIR 数字滤波器的原理和方法。
如果所希望的滤波器的理想频率响应函数为 )(ωj d e H ，则其对应的单位脉冲响应为)(n h d =π21ωωωππd e e H j j d )(⎰- （2-1）窗函数设计法的基本原理是用有限长单位脉冲响应序列)(n h 逼近)(n h d 。
由于)(n h d 往往是无限长序列，且是非因果的，所以用窗函数)(n ω将)(n h d 截断，并进行加权处理，得到：)(n h ＝)(n h d )(n ω (2-2))(n h 就作为实际设计的FIR 数字滤波器的单位脉冲响应序列，其频率响应函数)(ωj d e H 为：)(ωj d e H ＝∑-=-10)(N n j e n h ω (2-3) 式中，N 为所选窗函数)(n ω的长度。
由第七章可知，用窗函数法设计的滤波器性能取决于窗函数)(n ω的类型及窗口长度N 的取值。
设计过程中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N 。
各种类型的窗函数可达到的阻带最小衰减和过渡带宽度见第七章。
这样选定窗函数类型和长度N 后，求出单位脉冲响应)(n h ＝)(n h d ·)(n ω，并按式（2-3）求出)(ωj e H 。
)(ωj e H 是否满足要求，要进行验算。
一般在)(n h 尾部加零使长度满足于2的整数次幂，以便用FFT 计算)(ωj e H 。
数字信号处理实验五
实验报告实验名称____利用DFT分析离散信号频谱课程名称____数字信号处理________院系部:电气与电子工程专业班级：信息1002学生姓名：王萌学号: 11012000219同组人:实验台号:指导教师：范杰清成绩：实验日期:华北电力大学一、实验目的应用离散傅里叶变换(DFT)，分析离散信号x [k ]的频谱。
深刻理解DFT 分析离散信号频谱的原理，掌握改善分析过程中产生的误差的方法。
二、实验原理根据信号傅里叶变换建立的时域与频域之间的对应关系，可以得到有限长序列的离散傅里叶变换(DFT)与四种确定信号傅里叶变换的之间的关系，实现由DFT 分析其频谱。
Matlab 中提供了fft 函数，FFT 是DFT 的快速算法X=fft(x)：用于计算序列x 的离散傅里叶变换（DFT ）X=fft(x,n)：对序列x 补零或截短至n 点的离散傅里叶变换。
当x 的长度小于n 时，在x 的尾部补零使x 的长度达到n 点； 当x 的长度大于n 时，将x 截短使x 的长度成n 点；x=ifft(X)和x=ifft(X ，n)是相应的离散傅里叶反变换。
fftshift(x)将fft 计算输出的零频移到输出的中心位置。
利用DFT 计算离散周期信号 的频谱 分析步骤为：(1) 确定离散周期序列 ][~k x 的基本周期N ；(2) 利用fft 函数求其一个周期的DFT ，得到X [m ];(3) ][][~m X m X。
][~k x利用DFT 计算离散非周期信号x [k ] 的频谱分析步骤为：(1) 确定序列的长度M 及窗函数的类型。
当序列为无限长时，需要根据能量分布，进行截短。
(2) 确定作FFT 的点数N ；根据频域取样定理.为使时域波形 不产生混叠，必须取。
(3) 使用fft 函数作N 点FFT 计算X [m ]。
三、实验内容1、利用FFT 分析信号31,1,0 ),8π3cos(][ ==k k k x的频谱；(1) 确定DFT 计算的参数；(2) 进行理论值与计算值比较，讨论信号频谱分析过程中误差 原因及改善方法。
硕士信号处理实验报告(3篇)
第1篇一、实验背景随着信息技术的飞速发展，数字信号处理（DSP）技术已成为通信、图像处理、语音识别等领域的重要工具。
本实验旨在通过一系列实验，加深对数字信号处理基本原理和方法的理解，提高实际应用能力。
二、实验目的1. 理解数字信号处理的基本概念和原理。
2. 掌握常用信号处理算法的MATLAB实现。
3. 培养分析和解决实际问题的能力。
三、实验内容本实验共分为五个部分，具体如下：1. 离散时间信号的基本操作（1）实验目的：熟悉离散时间信号的基本操作，如加法、减法、乘法、除法、延时、翻转等。
（2）实验步骤：- 使用MATLAB生成两个离散时间信号。
- 对信号进行基本操作，如加法、减法、乘法、除法、延时、翻转等。
- 观察并分析操作结果。
2. 离散时间系统的时域分析（1）实验目的：掌握离散时间系统的时域分析方法，如单位脉冲响应、零状态响应、零输入响应等。
（2）实验步骤：- 使用MATLAB设计一个离散时间系统。
- 计算系统的单位脉冲响应、零状态响应和零输入响应。
- 分析系统特性。
（1）实验目的：掌握离散时间信号的频域分析方法，如快速傅里叶变换（FFT）、离散傅里叶变换（DFT）等。
（2）实验步骤：- 使用MATLAB生成一个离散时间信号。
- 对信号进行FFT和DFT变换。
- 分析信号频谱。
4. 数字滤波器的设计与实现（1）实验目的：掌握数字滤波器的设计与实现方法，如巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器等。
（2）实验步骤：- 使用MATLAB设计一个低通滤波器。
- 使用窗函数法实现滤波器。
- 对滤波器进行性能分析。
5. 信号处理在实际应用中的案例分析（1）实验目的：了解信号处理在实际应用中的案例分析，如语音信号处理、图像处理等。
（2）实验步骤：- 选择一个信号处理应用案例。
- 分析案例中使用的信号处理方法。
- 总结案例中的经验和教训。
四、实验结果与分析1. 离散时间信号的基本操作实验结果表明，离散时间信号的基本操作简单易懂，通过MATLAB可以实现各种操作，方便快捷。
数字信号处理实验报告_五个实验
实验一 信号、系统及系统响应一、 实验目的1、熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解；2、熟悉时域离散系统的时域特性；3、利用卷积方法观察分析系统的时域特性；4、掌握序列傅立叶变换的计算机实现方法，利用序列的傅立叶变换对连续信号、离散信号及系统响应进行频域分析。
二、 实验原理及方法采样是连续信号数字处理的第一个关键环节。
对采样过程的研究不仅可以了解采样前后信号时域和频域特性发生变化以及信号信息不丢失的条件，而且可以加深对傅立叶变换、Z 变换和序列傅立叶变换之间关系式的理解。
对一个连续信号)(t x a 进行理想采样的过程可用下式表示：)()()(^t p t t x x 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 jwae 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)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]);endk=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=yesinput1=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);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)]'); end end end子函数:离散傅立叶变换及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、认真复习采样理论，离散信号与系统，线性卷积，序列的傅立叶变换及性质等有关内容，阅读本实验原理与方法。
数字信号处理实验五用DFT(FFT)对信号进行频谱分析
开课学院及实验室：电子楼3172018年 4月 29 日3()x n ：用14()()x n R n =以8为周期进行周期性延拓形成地周期序列.(1> 分别以变换区间N ＝8,16,32,对14()()x n R n =进行DFT(FFT>,画出相应地幅频特性曲线；(2> 分别以变换区间N ＝4,8,16,对x 2(n >分别进行DFT(FFT>,画出相应地幅频特性曲线； (3> 对x 3(n >进行频谱分析,并选择变换区间,画出幅频特性曲线.<二）连续信号 1. 实验信号：1()()x t R t τ=选择 1.5ms τ=,式中()R t τ地波形以及幅度特性如图7.1所示.2()sin(2/8)x t ft ππ=+式中频率f 自己选择.3()cos8cos16cos 20x t t t t πππ=++2. 分别对三种模拟信号选择采样频率和采样点数.对1()x t ()R t τ=,选择采样频率4s f kHz =,8kHz ,16kHz ,采样点数用τ.s f 计算.对2()sin(2/8)x t ft ππ=+,周期1/T f =,频率f 自己选择,采样频率4s f f =,观测时间0.5p T T =,T ,2T ,采样点数用p s T f 计算.图5.1 R(t>地波形及其幅度特性对3()cos8cos16cos 20x t t t t πππ=++,选择采用频率64s f Hz =,采样点数为16,32,64. 3. 分别对它们转换成序列,按顺序用123(),(),()x n x n x n 表示.4. 分别对它们进行FFT.如果采样点数不满足2地整数幂,可以通过序列尾部加0满足.5. 计算幅度特性并进行打印.五、实验过程原始记录<数据、图表、计算等）(一> 离散信号%14()()x n R n = n=0:1:10。
实验五典型电信号的观察与测量
实验五：典型电信号的观察与测量一、实验目的：1、熟悉函数信号发生器的各种旋钮、开关的作用及其使用方法。
2、初步掌握用示波器观察电信号波形，定量测出正弦信号和脉冲信号的波形参数。
二、实验原理：1、正弦交流信号和方波脉冲信号是常用的电激励信号，可以由函数信号发生器提供2、电子示波器是一种信号波形观察和测量仪器，可定量测出电信号的波形参数。
三、实验仪器：1、双踪示波器1台2、函数信号发生器1台四、实验内容与照片1、双踪示波器的自检2、信号的观测（1）幅值1v频率50Hz(2)幅值1v频率1.5kHz(3) 幅值1v频率20kHz信号发生器读数项目测定正弦信号频率的测定50Hz 1500Hz 20000Hz示波器“t/div”值10ms 0.2ms 0.02ms一个周期占有的格数2 3.2 2.2信号周期(s) 0.02 0.00064 0.000044计算所得频率（Hz）50 1562.5 22727.27在这里，“t/div”表示横向每个所代表的时间。
可以看出，经计算后的频率与实际频率基本一致，但是仍有一定误差，误差应该主要来源于观察一个周期占了几个格子时产生。
而且信号从函数信号发生器到示波器的导线接口处信号可能会略有变化。
（4）频率1kHz 有效值0.3v（5）频率1kHz 有效值1v（6）频率1kHz 有效值3v信号发生器读数正弦信号幅值的测定项目测定0.3v 1v 3v示波器“v/div”值0.48v 1v 5v峰峰值波形格数 2 2.96 1.72峰值0.48v 1.48v 4.3v计算所得有效值0.339v 1.047v 3.041v 在这里“t/div”代表纵向每一格代表多少电压。
可以看到，这里在误差允许范围内计算值等于理论值。
但是仍有一定误差，误差应该主要来源于观察波峰与波谷之间占了几个格子时产生。
而且信号从函数信号发生器到示波器的导线接口处信号可能会略有变化。
3.方波脉冲信号的测定（1）幅值3.0v 频率300Hz(2) 幅值3.0v 频率3kHz(3) 幅值3.0v 频率30kHz信号发生器读数方波脉冲信号频率的测定在这里，“t/div”表示横向每个所代表的时间。
数字信号处理--实验五-用DFT(FFT)对信号进行频谱分析
学生实验报告开课学院及实验室：电子楼3172013年4月29日、实验目的学习DFT 的基本性质及对时域离散信号进行频谱分析的方法，进一步加深对频域概念和数字频率的理解，掌握 MATLAB 函数中FFT 函数的应用。
二、实验原理离散傅里叶变换(DFT)对有限长时域离散信号的频谱进行等间隔采样，频域函数被离散化了, 便于信号的计算机处理。
设x(n)是一个长度为 M 的有限长序列，x(n)的N 点傅立叶变换:X(k)N 1j 三 knDFT[x(n)]N x(n)e N0 k N 1n 0其中WNe.2 jN，它的反变换定义为:1X(n)NkN 1nkX(k)W N0 令z W N k，X(zz WN k则有:N 1x( n)Wj kn 0可以得到，X(k)X(Z)Z WN kZ W N*是Z 平面单位圆上幅角为2kN 的点，就是将单位圆进行N 等分以后第 K 个点。
所以, X(K)是Z 变换在单位圆上的等距采样，或者说是序列傅立叶变换的等距采样。
时域采样在满足Nyquist 定理时，就不会发生频谱混叠。
DFT 是对序列傅立叶变换的等距采样，因此可以用于序列的频谱分析。
如果用FFT 对模拟信号进行谱分析，首先要把模拟信号转换成数字信号，转换时要求知道模拟 信号的最高截至频率，以便选择满足采样定理的采样频率。
般选择采样频率是模拟信号中最高频率的3~4倍。
另外要选择对模拟信号的观测时间，如果采样频率和观测时间确定，则采样点数也确定 了。
这里观测时间和对模拟信号进行谱分析的分辨率有关,最小的观测时间和分辨率成倒数关系。
最小的采样点数用教材相关公式确定。
要求选择的采样点数和观测时间大于它的最小值。
如果要进行谱分析的模拟信号是周期信号,最好选择观测时间是信号周期的整数倍。
如果不知道■ 厂1*1IE向i1A I1f Ii i 0r 1 疋0Jfb-4W0 70000图5.1 R(t)的波形及其幅度特性xn=[on es(1,4),zeros(1,7)];%输入时域序列向量 xn=R4( n)%计算xn 的8点DFTXk16=fft(x n,16);%计算xn 的16点DFTXk32=fft(x n,32); %计算xn 的32点DFTk=0:7;wk=2*k/8;对 x 3(t) cos8 t cos16 t cos20 t ，选择采用频率 f s 64Hz ，采样点数为 16 , 32 , 64。
实验五无失真传输系统
无失真传输系统一、实验目的1、理解无失真传输的概念2、理解无失真传输的条件二、实验内容1、观察信号在无失真系统中的波形2、观察信号在无失真系统中的波形三、实验仪器1、信号与系统实验箱一台2、系统频域与复域的分析模块一块3、20M双踪示波器一台四、实验原理1、什么是无失真传输无失真传输是指响应信号与激励信号相比，只是大小与出现的时间不同，而无波形上的变化。
设激励信号为f(t)，响应为y(t)，可知y(t)=Kf(t-t0)。
2、如何实现无失真传输为满足y(t)=Kf(t-t0) （1）这一条件，其频域上的关系应满足（2）从信号与系统的观点看，信号在通过系统时，系统相当于一个频谱变换器，由系统的无失真条件可知，只有当系统频谱（3）时系统才能对输入信号做到无失真传输，由（3）式可知，为使信号在通过线性系统时产生任何失真，系统频率复频特性必须在信号的全部频带内为一常数，相频特性应为一通过原点的直线。
3、实际无失真系统实际电路中，可使R1、R2、C1、C2中的之一固定，第四个元件函数可调，以满足R1C1=R2C2，本实验采用电阻可调。
五、实验测试点说明1、测试点分别为：“输入”：模拟信号的输入。
“输出”：模拟信号经过系统后的输出。
“GND”：与实验箱的地相连。
2、调节点分别为：“失真调节”：调节此电仪器，可以观察信号失真的过程。
六、实验步骤（1）将“系统频域与复域分析模块”插到实验箱上。
（2）将上述模块上的电源接入插孔，用导线与实验箱上的电源输出插孔对应相接。
（3）将“常用信号分类与观察模块”上的输出插口与“系统频域与复域分析模块”上的无失真传输系统的输入插口相接。
（4）示波器面板上的两个“VOLTS/DIV”旋钮均打到“5”，“TIME/DIV”旋钮打到“10ms”位置，按下“×10MAG”按钮。
（5）示波器面板上的“MODE”转换开关打到“CH1 ”,其余旋钮及开关均保持平常的测试位置。
（6）打开示波器电源，打开实验箱电源，按下“常用信号分类与观察模块”区中的三个电源按钮SP1、SP2及S5 ,把实验箱最左上角的模块上的“+12V，-12V，GND”和“复域分析模块”上的相对应的“+12V，-12V，GND”用导线连接上，可以看到实验箱左上角的4 个二极管以及“复域分析模块”上的2 个二极管发光。
数字信号处理--实验五
实验五FIR数字滤波器的设计04011344 王晨一、实验目的(1) 掌握用窗函数法、频率采样法及优化设计法设计FIR滤波器的原理及方法，熟悉相应的MATLAB编程。
(2) 熟悉线性相位FIR滤波器的幅频特性和相频特性。
(3) 了解各种不同窗函数对滤波器性能的影响。
二、实验原理①线性相位实系数FIR滤波器按其N值奇偶和h(n)的奇偶对称性分为四种：1、h(n)为偶对称，N为奇数;H(ejω)的幅值关于ω=0，π，2π成偶对称。
2、h(n)为偶对称，N为偶数;H(ejω)的幅值关于ω=π成奇对称，不适合作高通。
3、h(n)为奇对称，N为奇数;H(ejω)的幅值关于ω=0，π，2π成奇对称，不适合作高通和低通。
4、h(n)为奇对称，N为偶数;H(ejω) ω=0、2π＝0，不适合作低通。
② 窗口法窗函数法设计线性相位FIR 滤波器步骤：➢ 确定数字滤波器的性能要求：临界频率k {}ω，滤波器单位脉冲响应长度N ； ➢ 根据性能要求，合理选择单位脉冲响应(n)h 的奇偶对称性，从而确定理想频率响应j (e )d H ω的幅频特性和相频特性；➢ 求理想单位脉冲响应(n)d h ，在实际计算中，可对j (e )d H ω按M(M 远大于N)点等距离采样，并对其求IDFT 得(n)M h ，用(n)M h 代替(n)d h ；➢ 选择适当的窗函数(n)ω，根据d (n)h (n)(n)h ω=求所需设计的FIR 滤波器单位脉冲响应；➢ 求j (e )H ω，分析其幅频特性，若不满足要求，可适当改变窗函数形式或长度N ，重复上述设计过程，以得到满意的结果。
窗函数的傅式变换j (e )W ω的主瓣决定了j (e )H ω过渡带宽。
j (e )W ω的旁瓣大小和多少决定了j (e )H ω在通带和阻带范围内波动幅度，常用的几种窗函数有：（1） 矩形窗(Rectangle Window) N (n)R (n)ω=（2） 汉宁(Hanning)窗，又称升余弦窗（3） 汉明(Hamming)窗，又称改进的升余弦窗（4） 布莱克曼(Blankman)窗，又称二阶升余弦窗（5） 凯塞(Kaiser)窗其中：β是一个可选参数，用来选择主瓣宽度和旁瓣衰减之间的交换关系，一般说来，β越大，过渡带越宽，阻带越小衰减也越大。
信号处理综合实验报告(3篇)
第1篇一、实验目的1. 深入理解信号处理的基本原理和方法。
2. 掌握信号处理在各个领域的应用，如语音信号处理、图像处理等。
3. 熟悉实验设备的使用，提高实际操作能力。
4. 培养团队协作和问题解决能力。
二、实验内容本次实验主要分为以下几个部分：1. 语音信号处理（1）采集语音信号：使用麦克风采集一段语音信号，并将其转换为数字信号。
（2）频谱分析：对采集到的语音信号进行频谱分析，观察其频谱特性。
（3）噪声消除：设计并实现噪声消除算法，对含噪语音信号进行处理，提高信号质量。
（4）语音增强：设计并实现语音增强算法，提高语音信号的清晰度。
2. 图像处理（1）图像采集：使用摄像头采集一幅图像，并将其转换为数字图像。
（2）图像增强：对采集到的图像进行增强处理，如对比度增强、亮度增强等。
（3）图像滤波：设计并实现图像滤波算法，去除图像中的噪声。
（4）图像分割：设计并实现图像分割算法，将图像中的不同区域分离出来。
3. 信号处理算法实现（1）傅里叶变换：实现离散傅里叶变换（DFT）和快速傅里叶变换（FFT）算法，对信号进行频谱分析。
（2）小波变换：实现离散小波变换（DWT）算法，对信号进行时频分析。
（3）滤波器设计：设计并实现低通滤波器、高通滤波器、带通滤波器等，对信号进行滤波处理。
三、实验原理1. 语音信号处理（1）语音信号采集：通过麦克风将声音信号转换为电信号，再通过模数转换器（ADC）转换为数字信号。
（2）频谱分析：利用傅里叶变换将时域信号转换为频域信号，分析信号的频谱特性。
（3）噪声消除：采用噪声消除算法，如维纳滤波、谱减法等，去除信号中的噪声。
（4）语音增强：利用语音增强算法，如谱峰增强、长时能量增强等，提高语音信号的清晰度。
2. 图像处理（1）图像采集：通过摄像头将光信号转换为电信号，再通过模数转换器（ADC）转换为数字图像。
（2）图像增强：通过调整图像的亮度、对比度等参数，提高图像的可视效果。
（3）图像滤波：利用滤波器去除图像中的噪声，如均值滤波、中值滤波、高斯滤波等。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
机电工程学院实验报告
课程名称：测试信号与分析Array实验项目名称：离散傅立叶变换DFT 实验时间：2020.04.07
班级：测控172 姓名：梁宇学号：201711501218
实验目的:
１.运用MATLAB计算有限长序列的DFT和IDFT。
实验环境:
计算机、Matlab软件
实验内容及过程:
１. 构造离散傅立叶的MATLAB程序正变换，画出xn=R4（x）的N=8、N=16、和N=64的离散傅里叶变换图n=[1,1,1,1];
程序如下：
>>
clear;clc;close all;
n=[1 1 1 1];
H=8;
J=16;
W=64;
y=fft(n,H);N=0:H-1;
subplot(3,1,1);stem(N,y);
title('8点DFT LY-201711501218');
y2=fft(n,J);N=0:J-1;
subplot(3,1,2);stem(N,y2);title('16点DFT');
y3=fft(n,W);N=0:W-1;
subplot(3,1,3);stem(N,y3);title('64点DFT');
结果如下:
2. 画出的()sin(2/)sin(2*2/)x n n N n N ππ=+的波形图，求其离散傅立叶变换，画出幅频曲线和相频曲线，并用离散傅立叶变换出来的结果做离散傅立叶反变换，画出反变换的结果曲线，并与
()sin(2/)sin(2*2/)x n n N n N ππ=+的波形图对比。
程序如下：
>>
clear;clc;close all;
n=[1:8];N=8;
a=sin(2*n*pi/N)+sin(4*n*pi/N);
x=fft(a);
am=abs(x);
sng=angle(x);
y=ifft(x);
subplot(2,2,1),plot(a);title('波形图 LY-201711501218');
subplot(2,2,2),plot(am);title('幅频响应');
subplot(2,2,3),plot(sng);title('相频响应');
subplot(2,2,4),plot(y);title('反变换结果曲线');
结果如下:
实验结果及分析：
１.讨论实验程序中的.^代表什么含义？
答：.^是矩阵中的每个元素都求平方
２.讨论实验程序N对分析结果的影响
答：N值越大会使结果曲线越倾向于平滑。
实验心得：
通过这次设计实践。
我学会了卷积运算计算系统的输出序列，加深对离散系统的差分方程、冲激响应和卷积分析方法的理解，对matlab的工作原理和使用方法也有了更深刻的理解。
在对理论的运用中，提高了我对matlab软件的运用熟练程度，在没有做实验以前，我们所知道的都是理论上的知识，而通过了这次的实验让我们把理论上的知识通过软件做出图来，让我们更加直观的认识离散傅立叶变换DFT。
附 录：
离散傅立叶变换DFT 的定义
一个有限长度的序列x(n)（0≤n<N-1）, 它的 DFT X(k) 可以通过在ω轴（πω20<≤）上对)
(ωj e X 均匀采样得到 2/2/()()
()j j kn N k N n X k X e x n e ωπωπ∞-==-∞==∑ 10-≤≤N k
可以看到)(k X 也是频域上的有限长序列，长度为N 。
序列)(k X 称为序列x(n)的N 点DFT 。
N 称为DFT 变换区间长度。
通常表示
N
j N e W /2π-=
可将定义式表示为
∑∞-∞==
n kn W n x k X )()( 10-≤≤N k
X(k)的离散傅里叶逆变换(IDFT)为
∑∞-∞=-=n kn W k X N n x )(1)( 10-≤≤N n。