《数字信号处理上机指导》电信09级(实验1)
数字信号处理上机一
电子科技大学生命科学与技术学院标准实验报告（实验）课程名称数字信号处理2012-2013-第2学期电子科技大学教务处制表电子科技大学实验报告实验地点：清水河校区实验时间：2013年4月10 日一、实验室名称：科二 504 机房二、实验名称：三、实验学时：2学时四、实验原理：五、实验目的：（详细填写）1、熟悉MATLAB编程；2、了DFT的性质及其应用。
六、实验内容：（详细填写）1、利用DFT实现两个序列的线形卷积运算，并研究DFT点数与计算结果的关系；2、利用FFT/IFFT实现简单的频率滤波。
3、验证DFT的线形性质。
七、实验器材（设备、元器件）：八、实验步骤：九、实验数据及结果分析：（详细填写）（包括程序、图、结果等）实验一：“程序：x=[2 0 1 1 0 9 1 0 2 0 0 0 4];%长度13y=[1 9 9 4 0 1 1 0];%长度8k=0;for i=1:3N=17+i;xk=fft(x,N);yk=fft(y,N);f=ifft(xk.*yk);k=k+1;figure(k);stem(f);title(['用DFT计算线性卷积: N=' num2str(N)]);ycov=conv(x,y);k=k+1;figure(k);stem(ycov,'r');title('调用CONV计算线性卷积');k=k+1;figure(k);stem(f);hold on;stem(ycov,'r');title('两种方法比较')；endN=180102030405060708090100用DFT 计算线性卷积: N=180102030405060708090100调用CONV 计算线性卷积0102030405060708090100两总方法比较N=190102030405060708090100用DFT 计算线性卷积: N=1901020304050607080901000102030405060708090100两总方法比较N=2001020304050607080901000102030405060708090100调用CONV 计算线性卷积0102030405060708090100两总方法比较实验二 程序：n=101; x1=[1:n]; y1=sin(x1.*0.3); noise=rand(1,n)-0.5; y2=y1+noise; figure(1);plot(x1,y1,'r',x1,y2,'b:'); xxx=fft(y2); figure(2); plot(abs(xxx)); title('信号的幅度谱'); for i = 10 : n-20 xxx(i) = 0; endf1 = ifft(xxx); figure(3);plot(x1,y1,'r',x1,f1,'b:'); title('滤波的效果');020406080100120-1.5-1-0.50.511.52040608010012005101520253035404550信号的幅度谱020406080100120-1.5-1-0.50.511.5滤波的效果实验三： 程序：x=[1 2 3 4];y=[5 6 7 8]; xk=fft(x); yk=fft(y); zk1=3*xk+4*yk; stem(zk1);title('3*X(k)+4*Y(K)'); figure(2); z=3*x+4*y; zk2=fft(z); stem(zk2);title('DFT(3*x(k)+4*y(k))');3*X(k)+4*Y(K)1 1.52 2.53 3.54DFT(3*x(k)+4*y(k))1 1.52 2.53 3.54十、实验结论：（详细填写）实验一：当N为两个序列长度和减一是，循环卷积等于线形卷积。
数字信号处理实验报告
《数字信号处理》实验报告课程名称：《数字信号处理》学院：信息科学与工程学院专业班级：通信1502班学生姓名：侯子强学号：02指导教师：李宏2017年5月28日实验一离散时间信号和系统响应一. 实验目的1. 熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解2. 掌握时域离散系统的时域特性3. 利用卷积方法观察分析系统的时域特性4. 掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对离散信号及系统响应进行频域分析二、实验原理1. 采样是连续信号数字化处理的第一个关键环节。
对采样过程的研究不仅可以了解采样前后信号时域和频域特性的变化以及信号信息不丢失的条件，而且可以加深对离散傅里叶变换、Z 变换和序列傅里叶变换之间关系式的理解。
对连续信号()a x t 以T 为采样间隔进行时域等间隔理想采样，形成采样信号：ˆ()()()a a xt x t p 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 为信号的最高频率，才可能由采样信号无失真地恢复出原模拟信号计算机实现时，利用计算机计算上式并不方便，因此我们利用采样序列的傅里叶变换来实现，即而()()j j n n X e x n e ωω∞-=-∞=∑为采样序列的傅里叶变换()()n P t t nT δ∞=-∞=-∑µ1()()*()21()n a a a s X j X j P j X j jn T π∞=-∞Ω=ΩΩ=Ω-Ω∑µ()()|j a TX j X e ωω=ΩΩ=2. 时域中，描述系统特性的方法是差分方程和单位脉冲响应，频域中可用系统函数描述系统特性。
数字信号处理实验指导书(需印刷)
实验一 常见离散信号的MATLAB 产生和图形显示1、实验目的：（１）熟悉MATLAB 应用环境，常用窗口的功能和使用方法。
（２）加深对常用离散时间信号的理解。
（３）掌握简单的绘图命令。
（４）掌握线性卷积的计算机编程方法。
2、实验原理：（１）单位抽样序列⎩⎨⎧=01)(n δ0≠=n n如果)(n δ在时间轴上延迟了k 个单位，得到)(k n -δ即：⎩⎨⎧=-01)(k n δ≠=n k n（２）单位阶跃序列⎩⎨⎧=01)(n u 00<≥n n（３）矩形序列 ⎩⎨⎧=01)(n R N 其他10-≤≤N n（４）正弦序列)sin()(ϕ+=wn A n x（５）复正弦序列jwnen x =)(（６）指数序列na n x =)(（７）线性时不变系统的响应为如下的卷积计算式：∑∞-∞=-=*=m m n h m x n h n x n y )()()()()(3、实验内容及步骤：（１）复习常用离散时间信号的有关内容。
（２）编制程序产生上述6种序列（长度可输入确定，对(4) (5) (6)中的参数可自行选择），并绘出其图形。
（３）已知系统的单位脉冲响应),(9.0)(n u n h n=输入信号)()(10n R n x =，试用卷积法求解系统的输出)(n y ，并绘出n n x ~)(、n n h ~)(及n n y ~)(图形。
4、实验用MATLAB 函数介绍（１）数字信号处理中常用到的绘图指令(只给出函数名，具体调用格式参看help)figure(); plot(); stem(); axis(); grid on; title(); xlabel(); ylabel(); text(); hold on; subplot()（２）离散时间信号产生可能涉及的函数zeros(); ones(); exp(); sin(); cos(); abs(); angle(); real(); imag(); （３）卷积计算可能涉及的函数conv(); length()注：实验过程中也可以使用自己编制的自定义函数，如impseq()、stepseq()等。
数字信号处理上机实验报告
数字信号处理上机实验报告实验一 熟悉MATLAB 环境一、实验目的1、 熟悉MATLAB 的主要操作命令。
2、 学会简单的矩阵输入和数据读写。
3、 掌握简单的绘图命令。
4、 用MATLAB 编程并学会创建函数。
5、 观察离散系统的频率响应。
二、实验内容认真阅读本章附录，在MATLAB 环境下重新做一遍附录中的例子，体会各条命令的含义。
在熟悉MATLAB 基本命令的基础上，完成以下实验。
上机实验内容：1、 数组的加减乘除和乘方运算，输入[]4 3 2 1=A ，[]6 5 4 3=B ，求B A C +=，B A D -=，B A E *=.，B A F /.=，B A G .^=，并用stem 语句画出A 、B 、C 、D 、E 、F 、G 。
程序：>> 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;subplot(2,4,1);stem(A,'.'); subplot(2,4,2);stem(B,'.'); subplot(2,4,3);stem(C,'.'); subplot(2,4,4);stem(D,'.'); subplot(2,4,5);stem(E,'.'); subplot(2,4,6);stem(F,'.');subplot(2,4,7);stem(G,'.')2、 用MATLAB 实现下列序列。
a) 150 8.0)(≤≤=n n x nb) 150 )()32.0(≤≤=+n en x n jc) 150 )1.025.0sin(2)2.0125.0cos(3)(≤≤+++=n n n n x ππππ 程序： A) clear;clc; n=[0:15]; x1=0.8.^n;subplot(3,1,1),stem(x1) title('x1=0.8^n')xlabel('n'); ylabel('x1');B)clear;clc;n=[0:15];x2=exp((0.2+3j)*n);subplot(3,1,1),stem(x2)title('x2=exp((0.2+3j)*n)')xlabel('n'); ylabel('x2');C)clear;clc;n=[0:15];x3=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi);subplot(3,1,1),stem(x3)title('x3=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi)') xlabel('n'); ylabel('x3');3、 绘出下列时间常数的图形，对x 轴，y 轴以及图形上方均须加上适当的标注： a) s t t t x 100 )2sin()(≤≤=π b) 4s t 0 )sin()100cos()(≤≤=t t t x ππ >> m=0:0.01:10; n=0:0.01:4; x1t=sin(2*pi*m);x2t=cos(100*pi*n).*sin(pi*n); subplot(2,1,1);plot(m,x1t); subplot(2,1,2);plot(n,x2t);4、 给定一因果系统H(z)=(1+2-1-z z 2+)/(2167.0-1--+z z )，求出并绘制H(z)的幅频响应与相频响应。
数字信号处理上机实验dft
本科实验报告实验名称：数字信号处理上机实验作业1：用DFT 分析周期序列的频谱任务：设周期序列()cos(0.48)cos(0.52)xn n n ππ=+ 截取 N 点长得到 ()()()N x n x n R n = （1）N=10，做10点DFT ，得到 X1(k);（2）N=10，做100点补零DFT ，得到 X2(k); （3）N=100，做100点DFT ，得到 X3(k)。
要求：针对以上三种情况，分别输出|X1(k)|、|X2(k)|、|X3(k)|的图形，并进行比较、分析和讨论。
程序：clear all ; n=0:1000;xn=cos(pi*0.48*n)+cos(pi*0.52*n); Xk1=fft(xn(1:10),10); X1=abs(Xk1); subplot(3,1,1); stem(X1,'.'); xlabel('k'); ylabel('|X1(k)|'); title('N=10,10点DFT'); Xk2=fft(xn(1:10),100); X2=abs(Xk2); subplot(3,1,2); stem(X2,'.'); xlabel('k'); ylabel('|X2(k)|');title('N=10,100点补零DFT'); Xk3=fft(xn(1:100),100); X3=abs(Xk3); subplot(3,1,3); stem(X3,'.'); xlabel('k'); ylabel('|X3(k)|'); title('N=100,100点DFT');运行结果：k|X 1(k )|N=10,10点DFTk|X 2(k )|N=10,100点补零DFTk|X 3(k )|N=100,100点DFT分析：从幅度谱中我们可以明显看出，X1(k)的相邻谱线间隔大，栅栏效应明显，频率分辨率低。
数字信号处理实验报告一二
数字信号处理课程实验报告实验一 离散时间信号和系统响应一. 实验目的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. 时域中，描述系统特性的方法是差分方程和单位脉冲响应，频域中可用系统函数描述系统特性。
已知输入信号，可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应。
实验一数字信号分析与处理实验指导书
数字信号分析与处理实验指导书实验一用FFT做频谱分析一、实验目的：1.进一步加深对DFT算法原理和基本性质的了解。
2.学习用FFT对连续时域信号进行频谱分析的方法。
了解可能出现的分析误差及其原因，以便在实际中正确应用FFT。
二、实验内容：1.用FFT对连续信号进行频谱分析，先对Xa(t)进行时域采样得到X(n)=Xa(nt)。
再对X(n)进行FFT变换得到Xk(k)。
2.观察和分析泄露原因；改变程序参数观察混叠现象。
3.如何将频谱的横坐标定为频率刻度？三、实验设备与材料：PC计算机；MATLAB软件。
四、实验程序：T1=0.01；N1=40；n1=0:(N1-1);t1=n1*T1;X1=2*s i n(4*p i*t1)+5*c o s(8*p i*t1);X k1=f f t(x1,N1);M1=a b s(X k1):K1=(0:l o n g t h(M1)-1)*N1/l e n g t h(M1)’s u b p l o t(2,4,1);p l o t(t1，x1);a x i s([0,0.4,-7.5.,7]);Ti t l e(‘T1=0.01s,t1=0.4s’);y l a b e l(‘x1(t)’);s u b p l o t(2,4,5);s t e m(k1,M1);t i t l e(‘T1=0.01s,N1=40’);y l a b e l(‘x1(k));T2=0.01；N2=50；n2=0:(N2-1);t2=n2*T2;X2=2*s i n(4*p i*t2)+5*c o s(8*p i*t2);X k2=f f t(x2,N2);M2=a b s(X k2):K2=(0:l o n g t h(M2)-1)*N2/l e n g t h(M2)’s u b p l o t(2,4,2);p l o t(t2，x2);a x i s([0,0.4,-7.5.,7]);Ti t l e(‘T2=0.01s,t2=0.5s’);y l a b e l(‘x2(t)’);s u b p l o t(2,4,6);s t e m(k2,M2);t i t l e(‘T2=0.01s,N2=40’);y l a b e l(‘x2(k));T3=0.005；N3=40；n3=0:(N3-1);t3=n3*T3;X3=2*s i n(4*p i*t3)+5*c o s(8*p i*t3);X k3=f f t(x3,N3);M3=a b s(X k3):K3=(0:l o n g t h(M3)-1)*N3/l e n g t h(M3)’s u b p l o t(2,4,3);p l o t(t3，x3);a x i s([0,0.4,-7.5.,7]);Ti t l e(‘T3=0.005s,t3=0.2s’);y l a b e l(‘x3(t)’);s u b p l o t(2,4,7);s t e m(k3,M3);t i t l e(‘T3=0.005s,N3=40’);y l a b e l(‘x3(k));T4=0.005；N4=50；n4=0:(N4-1);t4=n4*T4;X4=2*s i n(4*p i*t4)+5*c o s(8*p i*t4);X k4=f f t(x4,N4);M4=a b s(X k4):K4=(0:l o n g t h(M4)-1)*N4/l e n g t h(M4)’s u b p l o t(2,4,4);p l o t(t4，x4);a x i s([0,0.4,-7.5.,7]);Ti t l e(‘T4=0.005s,t4=0.25s’);y l a b e l(‘x4(t)’);s u b p l o t(2,4,8);s t e m(k4,M4);t i t l e(‘T4=0.005s,N4=50’);y l a b e l(‘x4(k));。
《数字信号处理》实验指导书(完整)
《数字信号处理》实验指导书通信教研室安阳工学院二零零九年三月第1章 系统响应及系统稳定性1.1 实验目的● 学会运用MATLAB 求解离散时间系统的零状态响应；● 学会运用MATLAB 求解离散时间系统的单位取样响应；● 学会运用MATLAB 求解离散时间系统的卷积和。
1.2 实验原理及实例分析1.2.1 离散时间系统的响应离散时间LTI 系统可用线性常系数差分方程来描述，即∑∑==-=-Mj jN i i j n x b i n y a 00)()( （1-1） 其中，i a （0=i ，1，…，N ）和j b （0=j ，1，…，M ）为实常数。
MATLAB 中函数filter 可对式（13-1）的差分方程在指定时间范围内的输入序列所产生的响应进行求解。
函数filter 的语句格式为y=filter(b,a,x)其中，x 为输入的离散序列；y 为输出的离散序列；y 的长度与x 的长度一样；b 与a 分别为差分方程右端与左端的系数向量。
【实例1-1】 已知某LTI 系统的差分方程为)1(2)()2(2)1(4)(3-+=-+--n x n x n y n y n y试用MATLAB 命令绘出当激励信号为)()2/1()(n u n x n=时，该系统的零状态响应。
解：MATLAB 源程序为>>a=[3 -4 2];>>b=[1 2];>>n=0:30;>>x=(1/2).^n;>>y=filter(b,a,x);>>stem(n,y,'fill'),grid on>>xlabel('n'),title('系统响应y(n)')程序运行结果如图1-1所示。
1.2.2 离散时间系统的单位取样响应系统的单位取样响应定义为系统在)(n 激励下系统的零状态响应，用)(n h 表示。
DB09级数字信号处理课程设计
二○一一～二○一二学年第一学期电子信息工程系课程设计计划书班级：电信（DB）2009级班课程名称：数字信号处理课程设计学时学分：1周1学分指导教师：二○一一年十二月十六日1、课程设计目的：通过本课程设计使学生进一步巩固数字信号处理的基本概念、理论、分析方法和实现方法；使学生掌握的基本理论和分析方法方面的知识得到进一步扩展；使学生能有效地将理论和实际紧密结合；增强学生软件编程实现能力和解决实际问题的能力。
要求学生能够熟练地用Matlab/Simulink 等工具编写数字信号处理的应用程序；熟悉用FFT 对连续信号和离散信号进行谱分析的理论和方法；编程实现IIR 数字滤波器和FIR 数字滤波器；了解各种窗函数对滤波器特性的影响等，进一步明确数字信号处理理论的工程应用。
2、课程设计时间安排：第一天：教师布置设计任务，讲解设计要求，提示设计要点。
第二天：查阅资料（在图书馆或上网），弄清题目要求，提出解决方案。
第三～五天：根据题目要求，将理论推导与编程实现相结合，写出设计报告。
第六天：验收。
3、课程设计内容及要求：(1) 已知某测试系统观测到的信号12()()()()x t s t s t n t =++，其中 1()20sin(70)s t t π=， 2()40sin(200)s t t π=， n(t) 为高斯白噪声，均值为0，方差为0.1，试设计滤波系统，输入为 ()x t ，输出为1()s t 和2()s t 。
要求：① 确定采样频率s f 、采样长度N ，对()x t 进行采样； ② 设计IIR 数字滤波器；③ 通过计算机仿真对滤波器的性能进行分析。
(2) 针对一个含有5Hz 、15Hz 和30Hz 的混合正弦波信号()x t ，试设计滤波系统，滤除5Hz 和30Hz 的正弦分量，阻带的最小衰减不小于50dB 。
要求：① 确定采样频率s f 、采样长度N ，对()x t 进行采样； ② 选择合适的窗函数，设计FIR 带通数字滤波器； ③ 通过计算机仿真对滤波器的性能进行分析。
数字信号处理实验指导讲解
数字信号处理实验指导编写马令坤郭文强电气与电子信息学院实验一离散时间的信号和系统一、实验目的1、复习离散时间的信号和系统，复习离散时间重要类型的信号和它们的运算的实现。
2、复习离散时间信号理论中一些重要的结果，它们在数字信号处理中很有用。
二、实验原理1、典型序列单位采样序列；单位阶跃序列；实数指数序列；复数指数序列；正余弦序列；随机序列：MATLAB可用rand(1,N)和randn(1,N)来生成；周期序列。
2、序列的运算信号加；信号乘；改变比例；移位；折叠：fliplr(x)；取样和：sum(x(n1:n2)) 取样积：prod(x(n1:n2))；信号能量：sum(abs(x)^2)；信号功率：sum(abs(x)^2)/length(x)3、一些有用的结果单位采样合成：奇偶合成：几何级数：序列相关：卷积运算：∑∞-∞=-=kknkxnx)()()(δ)()()(nxnxn xoe+=1||,11<-→∑∞=aaann对∑∞-∞=-=nyxllnynxlr称为移位),()()(,),(yxconv∑∑==---=MmNkkmknyamnxbny01) ()()(差分方程： 在Matlab 中： 三、实验内容 1.典型序列的实现单位阶跃序列；实数指数序列；复数指数序列；正余弦序列；随机序列用rand(1,N)和randn(1,N)来生成； 2.序列的运算给定序列x1=[1 2 3 4 5 6 7 8 9], ns1=-4; x2=[9 8 7 6 5 4 3 2 1], ns2=4 求：1) x1+x2; 2) y3=x1×x2；3) y1=0.5×x1+0.8×x2；4) y2=0.3×x1(n)×δ(n-6)+0.8×δ(n-5)×x2(n); 5) x1和x2的反折序列; 6) x1(n)和x2(n)的功率； 7) y3=x1*x2 (线性卷积)；四、本实验用到的一些MATLAB 函数Stem(x,y),Plot(x,y): x 轴和y 轴均为线性刻度（Linear scale ） Xlable (‘x’), Ylable(‘y’) Title(‘x’)axis([xmin,xmax,ymin,ymax]):调整图轴的范围 Subplot(x,y,z):同时画出数个小图形在同一视窗中。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1
《数字信号处理》
上机实验指导书
光电信息与通信工程学院
电子信息工程专业
北京信息科技大学
2011年11月
2
目 录
实验一 离散傅里叶变换的性质及应用............................................................. 3
3
实验一 离散傅里叶变换的性质及应用
一、实验目的
1．了解DFT的性质及其应用
2．熟悉MATLAB编程
二、实验仪器及材料
计算机，MATLAB软件
三、实验内容及要求
1．验证 DFT性质。
（1）已知序列)3()2(2)1(3)(4)(nnnnnx＝，计算序列的6点
DFT的幅度、相位、实部及虚部，并用图形表示相应的)],([,)(),(kXanglekXnx
)(Im,)(RekXkX
。
（2）编程求解)(nxep和)(nxop ，验证
)](Re[)(DFTkXnxep
)](Im[)(DFTkXjnxop
（3）（选作）对序列)3()2(2)1(3)(4)(nnnnnx＝，自行设计
程序验证DFT时移及频移性质。
2．利用DFT实现两序列的卷积运算，并研究圆周卷积与线性卷积的关系。
（1）已知两序列
41103)(n
n
nx
3,2,1,0,2)(nnh
n
用直接法（即用线性卷积的定义计算，见下式）计算线性卷积y(n)=x(n)*h(n)的结果，
并以图形方式表示结果；
20),()()(2110NNnmnhmxnyNm
（2）利用DFT计算圆周卷积；分别令圆周卷积的点数为L=6，8，10，以图形方式表
示结果。
3．利用DFT进行谱分析。
（1）对模拟信号进行谱分析
4
)20cos(5.0)16cos(2)8cos()(ttttx
请问信号)(tx的最高频率是多少？
（2）假定抽样频率Hz64sf，编程实现记录数据点数为N=16,32,64 三种情况下的谱
分析。请问三种情况下谱分析的频率分辨率是多少？分别画出其幅度谱，讨论三种情况幅
度谱的区别、幅度谱与频率分辨率的对应关系。
四、思考题
1． 从在实验内容1中都可以验证DFT什么性质？
2． 说明为什么在实验内容2中L=6，8，10时采用DFT方法对两序列计算圆周卷积
的三种结果与线性卷积的关系如何？
3．
说明为什么在实验内容2中为N=16,32,64 三种情况下的谱分析结果不同？周期信
号的谱分析的点数应该如何选取？
五、实验报告要求
1．列出本实验编写的所有文件及各项实验结果曲线，加注必要的说明；
2．对实验结果作理论计算，解释实验结果；
3．回答思考题；
4．总结实验体会及实验中存在的问题。
六、本实验所用的部分MATLAB函数
abs(x) : 对复数取模（对实数取绝对值）。
angle(x) : 对复数取相角。
real(x) : 对复数取实部。
imag(x) : 对复数取虚部。
y=conv(x1,x2): 将序列x1(n)与序列x2(n)做线性卷积，结果赋给序列y(n)。
y=fft(x,L) : 将序列x(n)做L点快速傅立叶变换，结果赋给序列y(n)。
y=ifft(x,L) : 将序列x(n)做L点快速傅立叶反变换，结果赋给序列y(n)。
stem(x,y) : 以x为横坐标，y为纵坐标画离散曲线。