数字信号处理课程设计报告书.docVIP

数字信号处理课程设计报告书 课题名称 利用FFT计算并分析循环卷积和 线性卷积 姓 名 学 号 院、系、部 电气系 专 业 电子信息工程 指导教师 2010年 6 月30日 设计任务及要求： 1、生成长度为6和8的矩形序列； 2、进行FFT，并计算两序列的线性卷积和8点、14点和16点循环卷积，并绘制卷积结果图形； 3、比较循环卷积和线性卷积的关系。 指导教师签名： 2010年 月 日 二、指导教师评语： 指导教师签名： 2010年 月 日 三、成绩 指导教师签名： 2010年 月 日 利用FFT计算并分析循环卷积和线性卷积 一、设计目的 掌握快速傅立叶变换FFT、离散傅立叶变换DFT及线性卷积和循环卷积基本知识。通过MATLAB仿真，加深对线性卷积和循环卷积的理解；锻炼运用所学知识，独立分析问题、解决问题的综合能力。 设计要求 运用数字信号处理的基本理论和专业知识，生成长度为6和8的矩形序列；进行FFT，并计算两序列的线性卷积和8点、14点和16点循环卷积，并绘制卷积结果图形；比较循环卷积和线性卷积的关系。 三、系统原理 总体思路： 先产生长度为6和8的矩形序列x6n和x8n。 x8n=[1 1 1 1 1 1 1 1]; x6n=[1 1 1 1 1 1 ]; 接着进行线性卷积。MATLAB提供了conv函数，用于计算两个有限长序列的卷积。利用conv函数即可实现yn=conv(x6n, x8n)线性卷积。 在进行8点、12点和14点循环卷积。先对矩形序列进行FFT变换，得到X8k6和X8k8、X12k6和X12k8、X14k6和X14k8。然后将变换后的X8k6和X8k8、X12k6和X12k8、X14k6和X14k8各自相乘，再求傅立叶反变换，即可得到所求结果。 假设h(n)和x(n)都是有限长序列，长度分别为N,M，其中 L点循环卷积： 且H(k)=DFT[h(n)] 与X(k)=DFT[x(n)] （0=k=L-1,L=max[N,M]） 则由DFT的时域循环卷积定理有 Y(k)=DFT[y(n)]=H(k)X(k) (0=k=L-1) 线性卷积： 将x(n)和h(n)补零，使其长度均为L点 L点循环卷积是线性卷积以L为周期的周期延拓序列的主值序列。 当 L=M+N-1 时 ，线性卷积以L为周期进行周期延拓无混叠现象；L点的循环卷积能代表线性卷积。 四、程序设计 clc; clear; x8n=[1 1 1 1 1 1 1 1]; x6n=[1 1 1 1 1 1 ]; yn=conv(x6n,x8n); subplot(2,2,1); n=0:length(yn)-1; stem(n,yn,.); title(x8(n)与x6(n)线性卷积结果); xlabel(n);ylabel(y(n)); %当L=8时候的循环卷积 X8k6=fft(x6n,8); X8k8=fft(x8n,8); Y8k=X8k6.*X8k8; y8n=ifft(Y8k,8); subplot(2,2,2); n=0:length(y8n)-1; stem(n,y8n,.); title(当L=8时候的循环卷积); xlabel(n);ylabel(y8(n)); %当L=14时候的循环卷积 X14k6=fft(x6n,14); X14k8=fft(x8n,14); Y14k=X14k6.*X14k8; y14n=ifft(Y14k,14); subplot(2,2,3); n=0:length(y14n)-1; stem(n,y14n,.); title(当L=14时候的循环卷积); xlabel(n);ylabel(y14(n)); %当L=16时候的循环卷积 X16k6=fft(x6n,16); X16k8=fft(x8n,16); Y16k=X16k6.*X16k8; y16n=ifft(Y16k,16); subplot(2,2,4); n=0:length(y16n)-1; stem(n,y16n,.); title(当L=16时候的循环卷积); xlabel(n);ylabel(y16(n)); 五、仿真结果及分析线性卷积的长度为L=N+M-1只有当循环卷积长度L=N+M-1时，线性卷积以L为周期进行周期延拓时才无时域混叠现象。也就是说循环卷积等于线性卷积的条件是L=N+M-1。 此设计中分别进行了8点、14和16点的循环卷积。8点的循环卷积因为长度小于13，出现时域混