数字信号 实验二
数字信号处理 实验报告 实验二 应用快速傅立叶变换对信号进行频谱分析
数字信号处理实验报告实验二应用快速傅立叶变换对信号进行频谱分析2011年12月7日一、实验目的1、通过本实验，进一步加深对DFT 算法原理合基本性质的理解，熟悉FFT 算法 原理和FFT 子程序的应用。
2、掌握应用FFT 对信号进行频谱分析的方法。
3、通过本实验进一步掌握频域采样定理。
4、了解应用FFT 进行信号频谱分析过程中可能出现的问题，以便在实际中正确应用FFT 。
二、实验原理与方法1、一个连续时间信号)(t x a 的频谱可以用它的傅立叶变换表示()()j t a a X j x t e dt +∞-Ω-∞Ω=⎰2、对信号进行理想采样，得到采样序列()()a x n x nT =3、以T 为采样周期，对)(n x 进行Z 变换()()n X z x n z +∞--∞=∑4、当ωj ez =时，得到序列傅立叶变换SFT()()j j n X e x n e ωω+∞--∞=∑5、ω为数字角频率sT F ωΩ=Ω=6、已经知道：12()[()]j a m X e X j T T Tωωπ+∞-∞=-∑ ( 2-6 )7、序列的频谱是原模拟信号的周期延拓，即可以通过分析序列的频谱，得到相应连续信号的频谱。
（信号为有限带宽，采样满足Nyquist 定理）8、无线长序列可以用有限长序列来逼近，对于有限长序列可以使用离散傅立叶变换（DFT ）。
可以很好的反映序列的频域特性，且易于快速算法在计算机上实现。
当序列()x n 的长度为N 时，它的离散傅里叶变换为：1()[()]()N knN n X k DFT x n x n W-===∑ 其中2jNN W eπ-=，它的反变换定义为：101()[()]()N knN k x n IDFT X k X k W N --===∑比较Z 变换式 ( 2-3 ) 和DFT 式 ( 2-7 )，令kN z W -=则1()()[()]|kNN nkN N Z W X z x n W DFT x n ---====∑ 因此有()()|kNz W X k X z -==k N W -是Z 平面单位圆上幅角为2kNπω=的点，也即是将单位圆N 等分后的第k 点。
数字信号处理第二版(实验二) 信号的采样与重建
实验二信号的采样与重建一，实验目的（1）通过观察采样信号的混叠现象，进一步理解奈奎斯特采样频率的意义。
（2）通过实验，了解数字信号采样转换过程中的频率特征。
（3）对实际的音频文件作内插和抽取操作，体会低通滤波器在内插和抽取中的作用。
二，实验内容（1）采样混叠，对一个模拟信号Va(t)进行等间采样，采样频率为200HZ，得到离散时间信号V(n).Va（t）由频率为30Hz,150Hz,170Hz,250Hz,330Hz的5个正弦信号的加权和构成。
Va(t)=6cos(60pi*t)+3sin(300pi*t)+2cos(340pi*t)+4cos(500pi*t )+10sin(660pi*t)观察采样后信号的混叠效应。
程序：clear,close all,t=0:0.1:20;Ts=1/2;n=0:Ts:20;V=8*cos(0.3*pi*t)+5*cos(0.5*pi*t+0.6435)-10*sin(0.7*pi*t);Vn=8*cos(0.3*pi*n)+5*cos(0.5*pi*n+0.6435)-10*sin(0.7*pi*n);subplot(221)plot(t,V),grid on,subplot(222)stem(n,Vn,'.'),gridon,05101520-40-200204005101520-40-2002040（2）输入信号X(n)为归一化频率f1=0.043，f2=0.31的两个正弦信号相加而成，N=100，按因子M=2作抽取：（1）不适用低通滤波器；（2）使用低通滤波器。
分别显示输入输出序列在时域和频域中的特性。
程序：clear;N=100; M=2;f1=0.043; f2=0.31; n=0:N-1;x=sin(2*pi*f1*n)+sin(2*pi*f2*n); y1=x(1:2:100);y2=decimate(x,M,'fir'); figure(1);stem(n,x(1:N));title('input sequence'); xlabel('n');ylabel('fudu'); figure(2); n=0:N/2-1; stem(n,y1);title('output sequence without LP'); xlabel('n');ylabel('fudu'); figure(3); m=0:N/M-1;stem(m,y2(1:N/M));title('output sequence with LP'); xlabel('n');ylabel('fudu'); figure(4);[h,w]=freqz(x);plot(w(1:512),abs(h(1:512)));title('frequency spectrum of the input sequence'); xlabel('w');ylabel('fudu'); figure(5);[h,w]=freqz(y1);plot(w(1:512),abs(h(1:512)));title('frequency spectrum of the output sequence without LP'); xlabel('w');ylabel('fudu'); figure(6);[h,w]=freqz(y2);plot(w(1:512),abs(h(1:512)));title('frequency spectrum of the output sequence without LP'); xlabel('w');ylabel('fudu');0102030405060708090100-2-1.5-1-0.500.511.52input sequencenf u d u05101520253035404550-2-1.5-1-0.500.511.52output sequence without LPnf u d u05101520253035404550-1.5-1-0.50.511.5output sequence with LPnf u d u0.511.522.533.505101520253035404550frequency spectrum of the input sequencewf u d u00.51 1.52 2.53 3.551015202530frequency spectrum of the output sequence without LPwf u d u00.51 1.52 2.53 3.5510152025frequency spectrum of the output sequence without LPwf u d u（3）输入信号X(n)为归一化频率f1=0.043，f2=0.31的两个正弦信号相加而成，长度N=50，内插因子为2.（1）不适用低通滤波器；（2）使用低通滤波器。
实验2 数字信号处理的基本作图工具
实验二数字信号处理的基本作图工具一、实验目的1. 进一步加深对离散时间系统的理解。
2. 学习在MATLAB中怎样表示离散时间信号。
3. 熟悉离散时间信号的作图。
二、实验步骤1. 复习离散时间系统的有关内容。
2. 复习MATLAB的基本语法。
3. 按实验内容熟悉stem。
4. 编写程序。
5. 输出结果，总结结论，按要求写出实验报告。
三、实验内容1．掌握stem函数STEM(Y) plots the data sequence Y as stems from the x axis terminated with circles for the data value.STEM(X,Y) plots the data sequence Y at the values specified in X.例：t=[0:0.1:2]; x=cos(pi*t+0.6); stem(t,x);xn=[4,2,2,3,6,7]; stem(xn);思考：STEM(Y)与STEM(X,Y)有什么不同？STEM与PLOT函数有什么不同？2．掌握subplot函数H = SUBPLOT(m,n,p), or SUBPLOT(mnp), breaks the Figure window into an m-by-n matrix of small axes, selects the p-th axes for the current plot, and returns the axis handle. The axes are counted along the top row of the Figure window, then the second row, etc.例：n1=0:3;x1=[1,1,1,1];subplot(221);stem(n1,x1);title('x1序列');n2=0:7;x2=[1,2,3,4,4,3,2,1];subplot(222);stem(n2,x2);title('x2序列');n3=0:7;x3=[4,3,2,1,1,2,3,4];subplot(223);stem(n3,x3);title('x3序列');n4=0:7;x41=cos((pi/4)*n4);subplot(224);stem(n4,x41);title('x4序列');思考：subplot是怎样分配各个作图分区的顺序号的？3．信号的运算]0,1.0,4.0,7.0,1[)(1=n x ，]9.0,7.0,5.0,3.0,1.0[)(2=n x ，请作出)()(21n x n x +，)()(21n x n x 的图形。
“数字信号处理”实验报告二
实验报告课程名称：数字信号处理实验任课教师：杨鉴实验名称：离散时间系统的时域分析年级、专业：2015级通信工程学号：**********姓名：***日期：2017 年10 月9 日云南大学信息学院一、实验1.通过MATLAB仿真一些简单的离散时间系统，并研究他们的时域特性。
2.掌握卷积在MATLAB的算法并理解滤波的概念。
二、实验内容1. 假定另一个系统为y[n]=x[n]x[n-1],修改程序P2.3，计算这个系统的输出序列y1[n],y2[n]和y[n]。
比较y[n]和yt[n]。
这两个序列是否相等？该系统是线性系统吗？2. 考虑另一个系统：y[n]=nx[n]+x[n-1],修改程序P2.4，以仿真上面的系统并确定该系统是否为时不变系统。
3.修改程序P2.7，计算长度为15的序列h[n]和长度为10的序列x[n]的卷积，重做问题Q2.28。
h[n]和x[n]的样本值你自己给定。
4.修改程序P2.9，将输入序列改变成扫频正弦序列（长度为301、最低频率为0、最高频率为0.5）。
那个滤波器能更好的抑制输入信号x[n]的高频分量？三、主要算法与程序Q2.11：clf;n = 0:40;a = 2;b = -3;f1=0.1;f2=0.4;x11=[0 cos(2*pi*f1*n) 0];x12=[0 0 cos(2*pi*f1*n)];x21=[0 cos(2*pi*f2*n) 0];x22=[0 0 cos(2*pi*f2*n)];x = a*x11 + b*x21;y1 = x11.*x12;y2 = x21.*x22;xd = a*x12+b*x22;y = x.*xd;yt = a*y1 + b*y2;d = y - yt; % Compute the difference output d[n]% Plot the outputs and the difference signalsubplot(3,1,1)stem([0 n 0],y);ylabel('Amplitude');title('Output Due to Weighted Input: a \cdot x_{1}[n] + b \cdot x_{2}[n]');subplot(3,1,2)stem([0 n 0],yt);ylabel('Amplitude');title('Weighted Output: a \cdot y_{1}[n] + b \cdot y_{2}[n]'); subplot(3,1,3)stem([0 n 0],d);xlabel('Time index n');ylabel('Amplitude');title('Difference Signal');Q2.17：clf;n = 0:40; D = 10;a = 3.0;b = -2;x = a*cos(2*pi*0.1*n) + b*cos(2*pi*0.4*n);xd = [zeros(1,D) x];nd=0:length(xd)-1;y=(n.*x)+[0 x(1:40)];yd=(nd.*xd)+[0 xd(1:length(xd)-1)];d = y - yd(1+D:41+D);subplot(3,1,1)stem(n,y);ylabel('振幅');title('输出 y[n]'); grid;subplot(3,1,2)stem(n,yd(1:41));ylabel('振幅');title('由于延时输入 x[n-10]的输出'); grid;subplot(3,1,3)stem(n,d);xlabel('时间序号 n'); ylabel('振幅');title('差值信号');grid;Q2.29：clf;h = [3 2 1 -2 1 0 -4 0 3 1 5 4 0 3 5]; % impulse responsex = [1 -2 3 -4 3 2 1 5 6 1]; % input sequencey = conv(h,x);n = 0:23;subplot(2,1,1);stem(n,y);xlabel('时间序号n');ylabel('振幅');title('用卷积得到的输出'); grid;x1 = [x zeros(1,14)];y1 = filter(h,1,x1);subplot(2,1,2);stem(n,y1);xlabel('时间序号 n'); ylabel('振幅');title('由滤波生成的输出'); grid;Q2.35：f=w/2pi=(2a*n+b)/2pi=[0,0.5],所以b=0，a*n<=0.5*pi,当n=300时，a取pi/600。
数字信号处理实验报告一二
数字信号处理课程实验报告实验一 离散时间信号和系统响应一. 实验目的1. 熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解2. 掌握时域离散系统的时域特性3. 利用卷积方法观察分析系统的时域特性4. 掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对离散信号及系统响应进行频域分析二、实验原理1. 采样是连续信号数字化处理的第一个关键环节。
对采样过程的研究不仅可以了解采样前后信号时域和频域特性的变化以及信号信息不丢失的条件，而且可以加深对离散傅里叶变换、Z 变换和序列傅里叶变换之间关系式的理解。
对连续信号()a x t 以T 为采样间隔进行时域等间隔理想采样，形成采样信号： 式中()p t 为周期冲激脉冲，()a x t 为()a x t 的理想采样。
()a x t 的傅里叶变换为()a X j Ω：上式表明将连续信号()a x t 采样后其频谱将变为周期的，周期为Ωs=2π/T 。
也即采样信号的频谱()a X j Ω是原连续信号xa(t)的频谱Xa(jΩ)在频率轴上以Ωs 为周期，周期延拓而成的。
因此，若对连续信号()a x t 进行采样，要保证采样频率fs ≥2fm ，fm 为信号的最高频率，才可能由采样信号无失真地恢复出原模拟信号ˆ()()()a a xt x t p t =1()()*()21()n a a a s X j X j P j X j jn T π∞=-∞Ω=ΩΩ=Ω-Ω∑()()n P t t nT δ∞=-∞=-∑计算机实现时，利用计算机计算上式并不方便，因此我们利用采样序列的傅里叶变换来实现，即而()()j j n n X e x n e ωω∞-=-∞=∑为采样序列的傅里叶变换2. 时域中，描述系统特性的方法是差分方程和单位脉冲响应，频域中可用系统函数描述系统特性。
已知输入信号，可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应。
(完整版)数字信号处理实验二
y = filter(num,den,x,ic);
yt = a*y1 + b*y2;
d = y - yt;
subplot(3,1,1)
stem(n,y);
ylabel('振幅');
title('加权输入: a \cdot x_{1}[n] + b \cdot x_{2}[n]的输出');
subplot(3,1,2)
%扫频信号通过2.1系统：
clf;
n = 0:100;
s1 = cos(2*pi*0.05*n);
s2 = cos(2*pi*0.47*n);
a = pi/2/100;
b = 0;
arg = a*n.*n + b*n;
x = cos(arg);
M = input('滤波器所需的长度=');
num = ones(1,M);
三、实验器材及软件
1.微型计算机1台
2. MATLAB 7.0软件
四、实验原理
1.三点平滑滤波器是一个线性时不变的有限冲激响应系统，将输出延时一个抽样周期，可得到三点平滑滤波器的因果表达式，生成的滤波器表示为
归纳上式可得
此式表示了一个因果M点平滑FIR滤波器。
2.对线性离散时间系统，若y1[n]和y2[n]分别是输入序列x1[n]和x2[n]的响应，则输入
plot(n, y);
axis([0, 100, -2, 2]);
xlabel('时间序号 n'); ylabel('振幅');
数字信号处理实验二
实验二信号的分析与处理综合实验一、实验目的目的：综合运用数字信号处理的理论知识进行信号的采样，重构，频谱分析和滤波器的设计，通过理论推导得出相应结论，再利用Matlab作为编程工具进行计算机实现，从而加深对所学知识的理解，建立概念。
二、基本要求1．掌握数字信号处理的基本概念、基本理论和基本方法；2．学会MA TLAB的使用，掌握MA TLAB的程序设计方法；3．掌握用MA TLAB设计简单实验验证采样定理的方法；4．掌握在Windows环境下语音信号采集的方法；5．学会用MA TLAB对信号进行频谱分析；6．掌握MA TLAB设计FIR和IIR数字滤波器的方法；三、实验内容1．利用简单正弦信号设计实验验证采样定理：（1）Matlab产生离散信号的方法，作图的方法，以及基本运算操作（2）对连续正弦信号以不同的采样频率作采样（3）对采样前后信号进行傅立叶变换，并画频谱图（4）分析采样前后频谱的有变化，验证采样定理。
掌握画频谱图的方法，深刻理解采样频率，信号频率，采样点数，频率分辨率等概念2．真实语音信号的采样重构：录制一段自己的语音信号，并对录制的信号进行采样；画出采样前后语音信号的时域波形和频谱图；对降采样后的信号进行插值重构，滤波，恢复原信号。
（1）语音信号的采集（2）降采样的实现（改变了信号的采样率）（3）以不同采样率采样后，语音信号的频谱分析（4）采样前后声音的变化（5）对降采样后的信号进行插值重构，滤波，恢复原信号3．带噪声语音信号的频谱分析（1）设计一频率已知的噪声信号，与实验2中原始语音信号相加，构造带噪声信号（2）画出原始语音信号和加噪声后信号，以及它们的频谱图（3）利用频谱图分析噪声信号和原语音信号的不同特性4．对带噪声语音信号滤波去噪：给定滤波器性能指标，采样窗函数法或双线性变换设计滤波器，并画出滤波器的频率响应；然后用自己设计的滤波器对采样的语音信号进行滤波，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化；回放语音信号；（1）分析带噪声信号频谱，找出噪声所在的频率段（2）利用matlab中已有的滤波器滤波（3）根据语音信号特点，自己设计滤波器滤波（4）比较各种滤波器性能（至少四种），选择一种合适的滤波器将噪声信号滤除（5）回放语音信号，比较滤波前后声音的变化四、实验原理参考《数字信号处理》教材《数字信号处理的MATLAB实现》万永革编著五、主要实验仪器及材料微型计算机、Matlab。
数字信号实验报告 (全)
三、实验内容和步骤
对以下典型信号进行谱分析：
x1 (n) R4 (n) n 1, x 2 (n) 8 n, 0 , 4 n, x3 (n) n 3, 0, 0n3 4n7
其它n
0n3 4n7
其它n
nห้องสมุดไป่ตู้
x4 ( n) cos
4
用 FFT 对信号作频谱分析是学习数字信号处理的重要内容。 经常需要进行谱 分析的信号是模拟信号和时域离散信号。 对信号进行谱分析的重要问题是频谱分 辨率 D 和分析误差。 频谱分辨率直接和 FFT 的变换区间 N 有关， 因为 FFT 能够实 现的频率分辨率是 2π /N≤D。可以根据此时选择 FFT 的变换区间 N。误差主要 来自于用 FFT 作频谱分析时，得到的是离散谱，而信号（周期信号除外）是连续 谱，只有当 N 较大时离散谱的包络才能逼近于连续谱，因此 N 要适当选择大一 些。 周期信号的频谱是离散谱，只有用整数倍周期的长度作 FFT，得到的离散谱 才能代表周期信号的频谱。 如果不知道信号周期，可以尽量选择信号的观察时间 长一些。 对模拟信号的频谱时， 首先要按照采样定理将其变成时域离散信号。如果是 模拟周期信号，也应该选取整数倍周期的长度，经过采样后形成周期序列，按照 周期序列的谱分析进行。
3 k 1
1
0.7051 z 2
0.0007378 1 z 1 1 1.0106 z 1 0.3583 z 2 1 0.9044 z 1 0.2155 z 2
6
H k z
（2.1）
式中：
H k z
A 1 2 z 1 z 2 ，k 1， 2， 3 1 Bk z 1 Ck z 2
数字信号处理实验二时域采样和频域采样
数字信号处理实验二时域采样和频域采样数字信号处理是一门研究信号的数字化表示、处理和传输的学科。
在数字信号处理中，时域采样和频域采样是两种常用的信号分析方法。
下面我们将对这两种采样方法进行详细介绍和比较。
一、时域采样时域采样是数字信号处理中最基本的采样方法之一。
它通过对连续时间信号进行离散时间采样，将连续时间信号转换为离散时间信号。
时域采样的基本原理是，如果一个连续时间信号f(t)在采样时刻t=kT(k=0,1,2,)上的值f(kT)能够被准确地测量，则可以通过这些采样值重建出原始信号。
时域采样的优点是简单易行，适用于大多数信号的采样。
但是，时域采样也存在一些缺点。
首先，如果信号中含有高于采样率的频率成分，这些高频成分将会被混叠到低频部分，导致信号失真。
这种现象被称为混叠效应。
其次，时域采样需要大量的采样数据才能准确地重建出原始信号，这会占用大量的存储空间和计算资源。
二、频域采样频域采样是一种在频域上对信号进行采样的方法。
它通过对信号进行傅里叶变换，将信号转换到频域，然后对频域中的信号进行采样。
频域采样的基本原理是，如果一个离散时间信号f(n)的傅里叶变换在频域上有有限的带宽，那么频域上的信号可以被认为是无穷多个离散的冲激函数的线性组合。
通过对这些冲激函数的幅度和相位进行采样，可以得到频域采样值。
相比时域采样，频域采样具有一些优点。
首先，频域采样可以避免混叠效应，因为高频成分在频域中可以被准确地表示和处理。
其次，频域采样只需要采样信号的幅度和相位信息，而不必存储大量的采样数据，可以节省存储空间和计算资源。
此外，频域采样还可以用于对信号进行压缩和编码，以便于信号的传输和存储。
然而，频域采样也存在一些缺点。
首先，傅里叶变换需要将信号从时域转换到频域，这需要使用复杂的数学运算和计算。
其次，频域采样的结果通常需要经过逆傅里叶变换才能得到原始信号的离散时间表示，这同样需要复杂的数学运算和计算。
此外，频域采样的结果可能存在频率混叠和泄漏现象，这会影响到重建出的原始信号的质量。
数字信号处理实验二
实验二 快速傅里叶变换(FFT)及其应用一、思考题(1) 实验中的信号序列()c x n 和()d x n 在单位圆上的z 变换频谱()()c j j d X e X e ωω和会相同吗?如果不同，说出哪一个低频分量更多一些，为什么?答：设j Z r e ω=⨯ ()()n n G z g n z ∞-=-∞=⨯∑因为为单位圆，故r=1.因为()()j j n n G e g n eωω∞-=-∞=⨯∑，故3723456704()(8)23432j j n j n j j j j j j j c n n X e nen e e e e e e e e ωωωωωωωωωω---------===+-=++++++∑∑7235670()(4)43223j j n j j j j j j d n X e n ee e e e e e ωωωωωωωω-------==-=+++---∑比拟可知频谱不相同，()c X n 的低频分量多。
(2) 对一个有限长序列进行DFT 等价于将该序列周期延拓后进行DFS 展开，因为DFS 也只是取其中一个周期来运算，所以FFT 在一定条件下也可以用以分析周期信号序列。
如果实正弦信号()sin(2),0.1x n fn f π== 用16点FFT 来做DFS 运算，得到的频谱是信号本身的真实谱吗?为什么?答：针对原来未经采样的连续时间信号来说，FFT 做出来的永远不会是信号本身的真实频谱，只能够是无限接近。
FFT 频谱泄露问题是一定会存在的，因为毕竟采样率再高，也不能完全到达原来的连续时间信号准确。
原题的采样率是1/10，就是将2*pi 分成10份，即每个正弦波周期进行10次采样，这样的采样率很低，而最后你只截取16个点来做分析，泄露一般会挺严重，看到的频谱，应该是一个上头尖，下面慢慢变宽的尖锥形，而纯粹的正弦波的理想频谱应该是在某频点只有一个尖峰。
二． 实验原理：〔1〕混叠：采样序列的频谱是被采样信号频谱的周期延拓，当采样频率不满足奈奎斯特采样定理的时候，就会发生混叠，使得刺痒后的序列信号的频谱不能真实的反映原采样信号的频谱。
数字信号处理实验二
实验二 快速傅里叶变换(FFT)及其应用一、思考题(1) 实验中的信号序列()c x n 和()d x n 在单位圆上的z 变换频谱()()c j j d X e X e ωω和会相同吗如果不同，说出哪一个低频分量更多一些，为什么答：设j Z r e ω=⨯ ()()n n G z g n z ∞-=-∞=⨯∑因为为单位圆，故r=1.因为()()j j n n G e g n eωω∞-=-∞=⨯∑，故3723456704()(8)23432j j n j n j j j j j j j c n n X e nen e e e e e e e e ωωωωωωωωωω---------===+-=++++++∑∑ 7235670()(4)43223j j n j j j j j j d n X e n ee e e e e e ωωωωωωωω-------==-=+++---∑比较可知频谱不相同，()c X n 的低频分量多。
(2) 对一个有限长序列进行DFT 等价于将该序列周期延拓后进行DFS 展开，因为DFS 也只是取其中一个周期来运算，所以FFT 在一定条件下也可以用以分析周期信号序列。
如果实正弦信号()sin(2),0.1x n fn f π== 用16点FFT 来做DFS 运算，得到的频谱是信号本身的真实谱吗为什么答：针对原来未经采样的连续时间信号来说，FFT 做出来的永远不会是信号本身的真实频谱，只能够是无限接近。
FFT 频谱泄露问题是一定会存在的，因为毕竟采样率再高，也不能完全达到原来的连续时间信号准确。
原题的采样率是1/10，就是将2*pi 分成10份，即每个正弦波周期进行10次采样，这样的采样率很低，而最后你只截取16个点来做分析，泄露一般会挺严重，看到的频谱，应该是一个上头尖，下面慢慢变宽的尖锥形，而纯正的正弦波的理想频谱应该是在某频点只有一个尖峰。
二．?实验原理：?（1）混叠：采样序列的频谱是被采样信号频谱的周期延拓，当采样频率不满足奈奎斯特采样定理的时候，就会发生混叠，使得刺痒后的序列信号的频谱不能真实的反映原采样信号的频谱。
数字信号实验二离散系统的时域分析
实验报告学院：计信学院专业：网络工程班级：网络091 姓名学号实验组实验时间2012-5-10 指导教师成绩实验项目名称实验2 离散系统的时域分析实验目的1、熟悉并掌握离散系统的差分方程表示法；2、加深对冲激响应和卷积分析方法的理解。
实验要求1、在MATLAB中，熟悉利用函数实现差分方程的仿真；2、在MATLAB中，熟悉用函数计算卷积，用求系统冲激响应的过程。
实验原理在时域中，离散时间系统对输入信号或者延迟信号进行运算处理，生成具有所需特性的输出信号，具体框图如下:][nx][nyDiscrete-timesystme其输入、输出关系可用以下差分方程描述：∑=∑=-=-MkkNkkknxpknyd][][输入信号分解为冲激信号,∑-=∞-∞=m m n m x n x ][][][δ记系统单位冲激响应][][n h n →δ，则系统响应为如下的卷积计算式：∑∞-∞=-=*=m m n h m x n h n x n y ][][][][][当N k dk,...2,1,0==时，h[n]是有限长度的（n ：[0，M]），称系统为FIR 系统；反之，称系统为IIR 系统。
实验步骤1. 在MATLAB 环境中编写程序。
2. 运行程序。
3. 得出实验结果。
实验内容1、以下程序中分别使用conv 和filter 函数计算h 和x 的卷积y 和y1，运行程序，并分析y 和y1是否有差别，为什么要使用x[n]补零后的x1来产生y1；具体分析当h[n]有i 个值，x[n]有j 个值，使用filter 完成卷积功能，需要如何补零？% Program P2_7 clf;h = [3 2 1 -2 1 0 -4 0 3]; %impulse response x = [1 -2 3 -4 3 2 1]; %input sequence y = conv(h,x);n = 0:14; subplot(2,1,1); stem(n,y);xlabel('Time index n'); ylabel('Amplitude'); title('Output Obtained by Convolution'); grid; x1 = [x zeros(1,8)]; y1 = filter(h,1,x1);subplot(2,1,2);stem(n,y1);xlabel('Time index n'); ylabel('Amplitude');title('Output Generated by Filtering'); grid;2、编制程序求解下列两个系统的单位冲激响应和阶跃响应，并绘出其图形。
数字信号处理实验02
实验二2019年11月16日一、实验目的（1）通过实验进一步理解卷积定理，了解卷积的过程。
（2）掌握应用线性卷积求解离散时间系统响应的基本方法。
（3）加深对离散LSI系统时域特性的认识。
（4）掌握MATLAB求解离散时间系统相应的基本方法。
（5）了解MATLAB中求解系统响应的子函数及其应用方法。
（6）加深对离散系统变换域分析——z变换的理解。
（7）掌握进行z变换和z反变换的基本方法，了解部分分式法在z反变换中的应用。
（8）掌握使用MATLAB语言进行z变换和z反变换的常用子函数。
（9）了解离散系统的基本描述模型。
（10）掌握各种模型相互间的关系及转换方法。
（11）熟悉MATLAB中进行离散系统模型间转换的常用子函数。
二、实验用到的MATLAB函数（1）conv 进行两个序列间的卷积运算。
（2）sum 求各元素之和。
（3）hold 控制当前图形是否刷新的双向切换开关。
（4）pause 暂停执行文件。
（5）dlsim 求解离散系统的响应。
（6）ztrans 返回无限长序列函数x[n]的z变换。
（7）iztrans 求函数X(z)的z反变换x[n]。
（8）syms 定义多个符号变量。
（9）residuez 有理多项式的部分分式展开。
（10）tf2zp 将系统传递函数（tf）模型转换为系统函数的零-极点增益（zpk）模型。
（11）zp2tf 将系统函数的零-极点增益（zpk）模型转换为系统传递函数（tf）模型。
（12）tf2sos 将系统传递函数（tf）模型转换为系统函数的二次分式（sos）模型。
（13）sos2tf 将系统函数的二次分式（sos）模型转换为系统传递函数（tf）模型。
（14）sos2zp 将系统函数的二次分式（sos）模型转换为系统函数的零-极点增益（zpk）模型。
（15）zp2sos 将系统函数的零-极点增益（zpk）模型转换为系统函数的二次分式（sos）模型。
（16）ss2tf 将系统状态空间（ss）模型转换为系统传递函数（tf）模型。
数字信号实验二离散时间傅里叶变换
一、实验项目名称离散时间傅里叶变数二、实验目的理解数值计算在离散时间傅里叶变换（DTFT)中的作用.没下载券联系企鹅2417677728给你传原文件三、实验内容与步骤1、脉冲信号的DTFT设矩形脉冲r[n]由下式定义r[n]=1 , 00 ,n L≤<⎧⎨⎩其它a.证明r[n]的DTFT可由下面的数学表示式得出R(e jω)=1sin()21sin()2Lωω·e-jω(L-1)/2该变换的第一项时常具有与DTFT相关的特殊形式，称为混叠sinc函数：a s i n c(ω,L)=1sin()21 sin()2Lωωb. 使用dtft函数计算12点脉冲信号的DTFT。
绘出在区间-π≤ω<π上对ω的DTFT。
把实部和虚部分开绘出，但是注意这些图不是很有用。
另绘出DTFT的幅度（参见M A TLAB 中的abs函数）。
选择频率样本的数量是脉冲长度的5到10倍，以使绘出的图看上去平滑。
用不同数量的频率样本做试验。
绘图时，注意要正确地标注频率坐标轴的变量ω。
c. 注意asinc函数零点的位置是规则分布的。
对奇数长脉冲，比如L=15的脉冲重复进行DTFT计算并绘出幅度；同样再次检验零点位置，注意峰值高度。
d. 对asinc函数零点的间距与asinc函数的直流值，确定出通用规则。
2、验证频移特性x(n)=cos(n*π/2),y(n)=exp(j*n*π/4)*x(n) 0≤n≤10绘出x(n)和y(n)的幅频和相频特性并比较3、系统分析一个LTI系统的差分方程为：y(n)=0.8y(n-1)+x(n)a、求H(e jω)b、求出并画出输入信号为x(n)=cos(0.5πn)u(n)的响应y(n)4、指数信号对于信号x[n] = (0.9)nu[n]，使用freqz函数计算其DTFT X(e jω)。
a. 对ω在区间-π≤ω<π上绘出幅度与相位特性。
这需要从freqz返回的[X, W]向量的移位。
数字信号处理--实验二_求线性时不变系统的输出
学生实验报告开课学院及实验室： 电子楼317 2013年 3 月18 日学院机械与电气工程学院年级、专业、班姓名学号 实验课程名称 数字信号处理实验成绩 实验项目名称实验二 求线性时不变系统的输出指导老师一、实验目的1. 学习用递推法求解差分方程的方法。
2. 学习用线性卷积法求网络输出的方法。
二、实验原理已知()h n 和输入()x n ，求系统输出()y n ，也可以用线性卷积法进行。
线性卷积法的公式如下：()()()y n x n h n =* ()()m x m h n m ∞=-∞=-∑图2.1 用递推法求解差分方程的框图计算时，关键问题是根据()x n 和()h n 的特点，确定求和的上下限。
例如，()0.9()nh n u n =，10()()x n R n =，卷积公式为10()()0.9()n m m y n R m u n m ∞-=-∞=-∑根据上式中的10()R n ，限制非零区间为：09m ≤≤，由()u n 限制非零区间为：m n ≤。
由上面的不等式知道m 的取值和n 有关，可以分几种情况： 当n <0时，()0y n =当0n ≤<9时，()0.9nn m m y n -==∑当9n ≥时，9()0.9n m m y n -==∑最后得到：再用计算机计算()y n 。
如果给定的()x n 和()h n 是一些离散数据，更方便的是用MATLAB 语言的数字信号工具箱函数conv 计算两个n 的取值从零开始的有限长序列的线性卷积。
三、使用仪器、材料1、硬件：计算机2、软件：Matlab四、实验步骤1.已知系统的差分方程如式：()0.9(1)()y n y n x n =-+（1） 输入信号＝10()()x n R n =，初始条件(1)1y -=，求解输出()y n ； （2） 输入信号＝10()()x n R n =，初始条件(1)0y -=，求解输出()y n 。
数字信号处理实验实验二离散时间傅里叶变换.wps
实验二 离散时间傅里叶变换一、实验原理经由正、逆离散时间傅里叶变换表达的傅里叶表示式是信号分析的一个关键部分，下面是分析方程与综合方程。
∑∞-∞=-=n n j j e n x e X ωω][)( ωπππωωd e e X n X n j j ⎰-=)(21][由以上公式知，离散时间傅里叶变换是w 的周期复值函数，周期是π2， 并且周期常选为[-π, π].对离散时间傅里叶变换有两个问题：（1） DTFT 的定义对无限长信号是有效的。
（2） DTFT 是连续变量的ω函数。
第二个问题是频率抽样问题。
Matlab 擅长在有线网格点上计算DTFT 。
通常选择足够多的频率以使绘出的图平滑，逼近真实的DTFT 。
对计算有利的最好选择是在（-π，π）区间上一组均匀的隔开的频率，或者共轭对称变换选择【0，π】，采用上述抽样方法，DTFT 式变为1,...,1,0,][)()(10/2(/2-===∑-=-N k en x e X e X L n N k j N k j j n k ）ππω在对DTFT 进行抽样时，并不要求N=L ，尽管通常由DFT 进行计算时，如果N=L 计算很方便。
1. 实验内容（1） asic 的m 文件编写一个matlab 文件如asic （w ，L ），之间从式中计算在频格上的asinc （w ，L ），该函数有两个输入：L 和W ，函数必须检查被0除的情形。
直接计算混叠sinc 函数得到脉冲信号DTFT 绘出幅度，保存该图以便与dtft 得到的结果比较。
程序：函数文件是function a=asinc(w,L)if (w==0)a=L;else a=sin(1/2*w*L)/sin(1/2*w);end当输入asinc(0,2)ans =2混叠程序function a=L(m);w=-pi:pi/50:pi;b=sin(1/2*w*m)./sin(1/2*w);c=exp(-j*w*(m-1)/2);R=b.*c;plot(w,abs(R))Grid当输入L（8）有而dtft程序是function[H,W]=dtft(h,N)%DTFT calculate DTFT at N equally spaced frequencies %Usage%[H,W]=dtft(h,N)%h:finite-length input vector,whose length is L%N:nambei of frequencies for evaluation over [-pi,pi) %==>constraint:n>=L%H:DTFT values(complex)%W:(2nd output)vector of freqs where DTFT is computed %N=fix(N);L=length(h);h=h(:);if(N<L)error('DTFT:#data samples cannot exceed #freq samples')endW=(2*pi/N)*[0:(N-1)]';mid=ceil(N/2)+1;W(mid:N)=W(mid:N)-2*pi;W=fftshift(W);H=fftshift(fft(h,N));当输入n=0:7;r=ones(8,1);r(1)=1[X,W]=dtft(r,128);plot(W,abs(X))grid有结果分析：两种方法得到的图形基本一致，证明了公式的真确性。
数字信号处理实验二DFT和FFT
数字信号处理实验二DFT和FFT一．实验内容1.认真复习周期序列DFS、有限长序列DFT的概念、旋转因子的定义、以及DFS和DFT的性质等有关内容；复习基2-FFT的基本算法，混合基-FFT的基本算法、Chirp-Z变换的算法等快速傅立叶变换的方法。
2.掌握有限长序列的循环移位、循环卷积的方法，对序列共轭对称性的含义和相关内容加深理解和掌握，掌握利用DFT分析序列的频谱特性的基本方法。
3.掌握FFT算法的基本原理和方法、Chirp-Z变换的基本原理和方法，掌握利用FFT分析序列的频谱特性的方法。
4.熟悉利用MATLAB进行序列的DFT、FFT的分析方法。
二．实验内容a.设周期序列某(n)={…,0,1,2,3,0,1,2,3,0,1,2,3,….}，求该序列的离散傅立叶级数某",(k)=DFS[某(n)]，并画出DFS的幅度特性。
在matlab中新建函数df:function[某k]=df(某n,N)n=0:1:N-1;k=0:1:N-1;Wn=e某p(-j某2某pi/N);nk=n'某k;Wnk=Wn.^nk;某k=某n某Wnk;~~在matlab中输入以下代码：某n=[0,1,2,3];k=0:1:3;N=4;－1－某k=df(某n,N);y=ab(某k);tem(k,y);title('周期序列的离散傅立叶级数');生成图像如下：由定义可知，对于周期序列，根据离散傅里叶级数公式即可求出，实验中显示了一个周期的傅里叶级数。
b.设周期方波序列为某(n)=1(mNnmNL-1)(m=0,1,2,....)0(mNLn(m1)N-1)其中N为基波周期，L/N是占空比。
(1)用L和N求|某",(k)|的表达式；(2)当L和N分别为：L=5，N=20；L=5，N=40；L=5，N=60以及－2－L=7，N=60时画出DFS的幅度谱；(3)对以上结果进行讨论，总结其特点和规律。
实验二 模拟和数字信号光纤传输系统实验
三、实验仪器
1.光纤通信实验系统 1 台。 2.示波器 1 台。 3.光纤跳线 2 根。 4.电话 2 部
四、实验原理
1.模拟信号光纤传输系统 本实验中将模拟信号源输出的正弦波、三角波、方波信号通过光纤进行传输。模拟信号源的电路图如 下：
图中 P400 是输入的方波信号， 输入的方波信号有两种频率可选 1k、 2k。 P401 是三角波的输出端， P410 是正弦波的输出端。 模拟信号也可以通过 PCM 编码后变成数字信号。然后，再送入光发射模块数字信号端进行传输。接 收到信号后再送入 PCM 译码模块，得到模拟信号。这种传输方法将在后面的实验中进行。 2.电话语音光纤传输系统 本实验系统的电话系统采用了热线电话的模式，热线电话的工作模式：其中任意一路摘机后（假定是 甲路） ，另一路将振铃（假定是乙路）而电话甲将送回铃音。当乙路摘机后，双方进入通话状态。当其中
模拟信号源模块正弦波输出p410p500pcm编译码模块一pcm编码输入模拟信号源模块正弦波输出p410p512pcm编译码模块二pcm编码输入pcm编译码模块一pcm编码输出p503p643pcm编码复用解复用模块复用输入一pcm编译码模块二pcm编码输出p507p642pcm编码复用解复用模块复用输入二pcm编码复用解复用模块复用输出p641p2001550nm光发模块数字光发数据输入1550nm光收模块数字信号输出p206p640pcm编码复用解复用模块解复用输入pcm编码复用解复用模块解复用输出一p639p509pcm编译码模块二译码输入pcm编码复用解复用模块解复用输出二p638p506pcm编译码模块一译码输入pcm编码复用解复用模块位时钟输出p644p505pcm编译码模块一pcm译码位时钟pcm编码复用解复用模块位时钟输出p644p510pcm编译码模块二pcm译码位时钟pcm编码复用解复用模块帧同步信号输出p637p504pcm编译码模块一译码帧同步信号pcm编码复用解复用模块帧同步信号输出p637p511pcm编译码模块二译码帧同步信号pcm编译码模块一pcm译码输出p501p515电话甲音频输入pcm编译码模块二pcm译码输出p513p517电话乙音频输入5
数字信号处理 实验二用FFT做谱分析
数字信号处理实验二用F F T做谱分析《数字信号处理》实践报告题目：实验二用FFT做谱分析1. 实验目的(1) 进一步加深DFT算法原理和基本性质的理解（因为 FFT只是DFT的一种快速算法，所以FFT的运算结果必然满足DFT的基本性质）。
(2) 熟悉FFT算法原理和FFT子程序的应用。
(3) 学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析误差及其原因，以便在实际中正确应用FFT。
2. 实验步骤(1) 复习DFT 的定义、性质和用DFT 作谱分析的有关内容。
(2) 复习按时间抽选法FFT 算法原理及相应的运算流图(3) 编制信号产生子程序，产生以下典型信号供谱分析用：x1(n) = R4(n)x 2(n) = ⎪⎩⎪⎨⎧≤≤-≤≤+n n n n n 其他,074,830,1 x 3(n) = ⎪⎩⎪⎨⎧≤≤-≤≤-nn n n n 其他,074,330,4 x 4(n) = cos(πn /4)x 5(n) = sin(πn /8)x 6(t) = cos8πt + cos16πt + cos20πt应当注意，如果给出的是连续信号x a (t)，则首先要根据其最高频率确定抽样频率f s 以及由频率分辨率选择抽样点数N ，然后对其进行软件抽样(即计算 x(n)=x a (nT)，0≤n ≤N-1)，产生对应序列 x(n)。
对信x 6(t)，频率分辨率的选择要以能分辨开其中的三个频率对应的谱线为准则。
对周期序列，最好截取周期的整数倍进行谱分析，否则有可能产生较大的分析误差。
请实验者根据DFT 的隐含周期性思考这个问题。
(4) 编写主程序。
图2.1对2中所给出的信号逐个进行谱分析。
下面给出针对各信号的FFT 变换区间N 以及对连续信号x 6(t)的抽样频率f s ，供实验时参考。
x 1(n) , x 2(n) , x 3(n) , x 4(n) , x 5(n)：N = 8 , 16x 6(t)：f s = 64(Hz) , N = 16 , 32 , 643、实验内容（1）x1(n) = R4(n)程序代码：x1=ones(1,4);N=4;n=[0:1:3];stem(n,x1);X1=fft(x1,8);magX1=abs(X1);k=[0:7];stem(k,magX1);x12=ones(1,4);N=16;n1=[0:3];stem(n1,x12);X12=fft(x12,16);magX12=abs(X12);k1=[0:15];stem(k1,magX12);subplot(2,2,1):stem(n,x1); subplot(2,2,2):stem(n1,x12); subplot(2,2,3):stem(k,magX1); subplot(2,2,4):stem(k1,magX12);（2）x 2(n) = ⎪⎩⎪⎨⎧≤≤-≤≤+nn n n n 其他,074,830,1 程序代码：x2=[1:4,4:-1:1];n=[0:3,4:7];stem(n,x2);X2=fft(x2,8);magX2=abs(X2);k=[0:7];stem(k,magX2);x2=[1:4,4:-1:1];n1=[0:3,4:7];stem(n,x2);X21=fft(x2,16);magX21=abs(X21);k1=[0:15];stem(k1,magX21);subplot(2,2,1):stem(n,x2);subplot(2,2,2):stem(n1,x2);subplot(2,2,3):stem(k,magX2);subplot(2,2,4):stem(k1,magX21);（3）x 3(n) = ⎪⎩⎪⎨⎧≤≤-≤≤-nn n n n 其他,074,330,4程序代码：x3=[4:-1:1,1:4];n=[0:3,4:7];stem(n,x3);X3=fft(x3,8);magX3=abs(X3);k=[0:7];stem(k,magX3);x3=[4:-1:1,1:4];n1=[0:3,4:7];stem(n,x3);X31=fft(x3,16);magX31=abs(X31);k1=[0:15];stem(k1,magX31);subplot(2,2,1):stem(n,x3); subplot(2,2,2):stem(n1,x3); subplot(2,2,3):stem(k,magX3); subplot(2,2,4):stem(k1,magX31);（4）x4(n) = cos(πn /4)程序代码：n=[0:7];x4=cos(pi*n/4);stem(n,x4);X4=fft(x4,8);magX4=abs(X4);k=[0:7];stem(k,magX4);n1=[0:15];x41=cos(pi*n1/4);stem(n,x4);X41=fft(x4,16);magX41=abs(X41);k1=[0:15];stem(k1,magX41);subplot(2,2,1):stem(n,x4); subplot(2,2,2):stem(n1,x41); subplot(2,2,3):stem(k,magX4); subplot(2,2,4):stem(k1,magX41);（5）x5(n) = sin(πn /8)程序代码：n=0:1:7;x5=sin(pi*n/8);stem(n,x5);X5=fft(x5,8);magX5=abs(X5);k=[0:7];stem(k,magX5);n1=0:1:15;x51=sin(pi*n1/8);stem(n1,x51);X51=fft(x5,16);magX51=abs(X51);k1=[0:15];stem(k1,magX51);subplot(2,2,1):stem(n,x5); subplot(2,2,2):stem(n1,x51); subplot(2,2,3):stem(k,magX5); subplot(2,2,4):stem(k1,magX51);（6）x6(t) = cos8πt + cos16πt + cos20πtN=32程序代码：fs=64;T=1/fs;t=0:T:1/2-T;x6=cos(2*pi*4*t)+cos(2*pi*8*t)+cos(2*pi*10*t); stem(t,x6);X6=fft(x6,32);magX6=abs(X6);k=[0:31];stem(k,magX6);subplot(2,1,1):stem(t,x6);subplot(2,1,2):stem(k,magX6);N=16程序代码：fs=64;T=1/fs;t=0:T:1/4-T;x6=cos(2*pi*4*t)+cos(2*pi*8*t)+cos(2*pi*10*t); stem(t,x6);X6=fft(x6,16);magX6=abs(X6);k=[0:15];stem(k,magX6);subplot(2,1,1):stem(t,x6);subplot(2,1,2):stem(k,magX6);N=64程序代码：fs=64;T=1/fs;t=0:T:1-T;x6=cos(2*pi*4*t)+cos(2*pi*8*t)+cos(2*pi*10*t); stem(t,x6);X6=fft(x6,64);magX6=abs(X6);k=[0:63];stem(k,magX6);subplot(2,1,1):stem(t,x6);subplot(2,1,2):stem(k,magX6);思考题：1、在N=8时，x2（n）和x3（n）的幅频特性会相同吗？为什么？N=16呢？答：在N=8时，x2(n)和x3(n)的幅频特性相同，x3(n)= x2((n-4))8,0≤n≤7DFT(x3(n))= e-j(2π/8)k4X2[k]=e-jπkX2[k]，所以x2(n)和x3(n)的幅频特性相同。
数字信号处理实验二时域采样和频域采样
实验二-时域采样和频域采样一、实验目的时域采样理论与频域采样理论是数字信号处理中的重要理论。
要求掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样后的信号不丢失信息；要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。
二、实验原理及方法1、时域采样定理的要点：a)对模拟信号)(t x a 以间隔T 进行时域等间隔理想采样，形成的采样信号的频谱)(ˆΩj X 是原模拟信号频谱()aX j Ω以采样角频率s Ω（T s /2π=Ω）为周期进行周期延拓b)采样频率s Ω必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的频谱不产生频谱混叠。
利用计算机计算上式并不方便，下面我们导出另外一个公式，以便用计算机上进行实验。
2、频域采样定理的要点：a)对信号x(n)的频谱函数X(ej ω)在[0，2π]上等间隔采样N 点 则N 点IDFT[()N X k ]得到的序列就是原序列x(n)以N 为周期进行周期延拓后的主值区序列。
三、实验内容及步骤1、时域采样理论的验证程序：clear;clcA=444.128;a=50*sqrt(2)*pi;w0=50*sqrt(2)*pi;Tp=50/1000;F1=1000;F2=300;F3=200;T1=1/F1;T2=1/F2;T3=1/F3;n1=0:Tp*F1-1;n2=0:Tp*F2-1;n3=0:Tp*F3-1;x1=A*exp(-a*n1*T1).*sin(w0*n1*T1);x2=A*exp(-a*n2*T2).*sin(w0*n2*T2);x3=A*exp(-a*n3*T3).*sin(w0*n3*T3);f1=fft(x1,length(n1));f2=fft(x2,length(n2)); %f3=fft(x3,length(n3)); %k1=0:length(f1)-1;fk1=k1/Tp; %k2=0:length(f2)-1;fk2=k2/Tp; % k3=0:length(f3)-1;fk3=k3/Tp; % subplot(3,2,1)stem(n1,x1,'.')title('(a)Fs=1000Hz');xlabel('n');ylabel('x1(n)');subplot(3,2,3)stem(n2,x2,'.')title('(b)Fs=300Hz');xlabel('n');ylabel('x2(n)');subplot(3,2,5)stem(n3,x3,'.')title('(c)Fs=200Hz');xlabel('n');ylabel('x3(n)');subplot(3,2,2)plot(fk1,abs(f1))title('(a) FT[xa(nT)],Fs=1000Hz'); xlabel('f(Hz)');ylabel('·ù¶È')subplot(3,2,4)plot(fk2,abs(f2))title('(b) FT[xa(nT)],Fs=300Hz'); xlabel('f(Hz)');ylabel('·ù¶È')subplot(3,2,6)plot(fk3,abs(f3))title('(c) FT[xa(nT)],Fs=200Hz'); xlabel('f(Hz)');ylabel('·ù¶È')结果分析：由图2.2可见，采样序列的频谱的确是以采样频率为周期对模拟信号频谱的周期延拓。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
DSP实验二
信通8班7号
朱志强
实验4 离散信号的DTFT 和DFT
（1）
DFT代码：
n=0:15;
y=cos(5*pi*n/16);
Y=fft(y,16);
YY=fft(y,32);
subplot(2,1,1);
stem(abs(Y));title('16点幅度谱');
subplot(2,1,2);
stem(abs(YY));title('32点幅度谱');
图：
（2）
DTFT代码：
n=0:15;
y=cos(5*pi*n/16);
w=linspace(0,2*pi,1000);
X=y*exp(-j*n'*w);
plot(w,X);title('16点序列DTFT') 图：
结论：对DTFT进行满足一定关系的频域抽样即可得到DFT。
实验5 FFT 算法的应用
（1）
n=0:63;
Y1=cos(2*pi*7*2*n/64)+0.5*cos(2*pi*19*2*n/64);
Y2=cos(2*pi*7*(2*n+1)/64)+0.5*cos(2*pi*19*(2*n+1)/64);
Y3=Y1+j*Y2;
Y4=Y1-j*Y2;
X3=fft(Y3,64);
X4=fft(Y4,64);
X1=0.5*(X3+X4);
X2=-0.5*j*(X3-X4);
k=0:127;
X1=[X1,X1];
X2=[X2,X2];
X=X1+X2.*exp(-j*pi*k/64);
stem(abs(X));
结论：通过有效合理的利用共扼偶分量和奇分量，能显著减少运算量。
（2）
b=[1,0];a=[1,-0.8];
[H,w]=freqz(b,a,'whole');
subplot(2,2,1);plot(abs(H));
subplot(2,2,1);plot(abs(H));title('频率响应');
xn=ifft(H);
subplot(2,2,2);stem(xn); title('x(n)');
X1=freqz(b,a,64,'whole');
subplot(2,2,3);stem(abs(X1));title('64点频域响应');
x1=ifft(X1,64);
subplot(2,2,4);stem(x1);title('64点x(n)');
结论:在IDTFT中，时域抽样，频域连续；而在DTF中，是频域抽样，时域延拓。