DSP实验四离散傅里叶变换.docx

DSP实验四离散傅里叶变换
装 订 线
实验报告
课程名称： 数字信号处理 指导老师： 成绩：
实验名称： 实验四离散傅里叶变换 实验类型： 程序实验
DSP实验四离散傅里叶变换
装 订 线
实验报告
课程名称： 数字信号处理 指导老师： 成绩：
实验名称： 实验四离散傅里叶变换 实验类型： 程序实验 同组学生姓名：
一、实验目的和要求（必填）
三、主要仪器设备（必填）
五、实验数据记录和处理
七、讨论、心得
二、实验内容和原理（必填）
四、操作方法和实验步骤
六、实验结果与分析（必填）
一、实验目的和要求


二、实验内容和原理
1. 编写函数DFT_new(x,N,win)
参数x表示输入序列
N表示DFT长度
win控制窗函数类型（条件语句）

求X3(n)的DFT，验证其线性
3. 序列的循环移位
求出的离散傅立叶变换X(k)，再对x(n)进行循环移位，移位量为m，并对其结果进行DFT变换得到X1(k)，看是否与成对应关系，m=6,N=11

L=P=6
xk=xn*A;%矩阵相乘得到xn的离散傅里叶变换
end
验证：N=10时，即N>序列长度
x=[1,1,1,1,1];
N=10;
win=1;%boxcar窗
xk=dft_new(x,N,win);
stem(0:1:N-1,abs(xk));
得到图形：
验证：N=8,N<序列长度时，需要用到数
x=[1,2,3,4,5,6,6,5,4,3,2,1];
N=8;
win=3；%hanning窗
xk=dft_new(x,N,win);
stem(0:1:N-1,abs(xk));

求X3(n)的DFT，验证其线性
n=0:1:10;
for N=0:10
x1(N+1)=10*^N;%x1序列
end
x2=[1 1 1 1 1 1 1 1 1 0 0];%x2序列
y11=dft_new(x1,length(x1),1);
y12=dft_new(x2,length(x1),1);%x1,x2分别DFT
y1=y11+y12; %DFT后求和
figure
subplot(2,1,1)
stem(n,abs(y1));
title('DFT(X1)+DFT(X2)');
x=x1+x2;
y2=dft_new(x,length(x1),1);%x1,x2的和进行DFT
subplot(2,1,2)
stem(n,abs(y2));
title('DFT(X1+X2)');
比较两幅图，是一样的，所以说明了DFT的线性性质
3. 序列的循环移位
求出的离散傅立叶变换X(k)，再对x(n)进行循环移位，移位量为m，并对其结果进行DFT变换得到X