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

个人收集整理 仅供参考学习 个人收集整理 仅供参考学习 PAGE / NUMPAGES 个人收集整理 仅供参考学习 ※※※※※※※※※ ※※※※※※※※※ ※※ ※※ ※※ ※※※※※※※※※ 2007级数字信号处理 课程设计 数字信号处理课程设计报告书 课题名称 利用FFT计算并分析循环卷积和 线性卷积 姓 名 学 号院、系、部 电气系 专 业 电子信息工程 指导教师 2010年 设计任务及要求： 1、生成长度为6和8地矩形序列； 2、进行FFT，并计算两序列地线性卷积和8点、14点和16点循环卷积，并绘制卷积结果图形； 3、比较循环卷积和线性卷积地关系. 指导教师签名： 2010年 月 日 二、指导教师评语： 指导教师签名： 2010年 月 日 三、成绩 指导教师签名： 2010年 月 日 利用FFT计算并分析循环卷积和线性卷积 一、设计目地 掌握快速傅立叶变换FFT、离散傅立叶变换DFT及线性卷积和循环卷积基本知识.通过MATLAB仿真，加深对线性卷积和循环卷积地理解；锻炼运用所学知识，独立分析问题、解决问题地综合能力.b5E2RGbCAP 设计要求 运用数字信号处理地基本理论和专业知识，生成长度为6和8地矩形序列；进行FFT，并计算两序列地线性卷积和8点、14点和16点循环卷积，并绘制卷积结果图形；比较循环卷积和线性卷积地关系.p1EanqFDPw 三、系统原理 总体思路： 先产生长度为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)线性卷积.DXDiTa9E3d 在进行8点、12点和14点循环卷积.先对矩形序列进行FFT变换，得到X8k6和X8k8、X12k6和X12k8、X14k6和X14k8.然后将变换后地X8k6和X8k8、X12k6和X12k8、X14k6和X14k8各自相乘，再求傅立叶反变换，即可得到所求结果.RTCrpUDGiT 假设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]）5PCzVD7HxA 则由DFT地时域循环卷积定理有 Y(k)=DFT[y(n)]=H(k)X(k) (0=k=L-1) 线性卷积： 将x(n)和h(n)补零，使其长度均为L点 对h(n) 对h(n)进行周期延拓： 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;jLBHrnAILg 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.只有当循