信号与信号处理实验六 离散傅里叶变换(DFT).ppt

实验三离散傅里叶变换(DFT)茅涝止荣株挥褒聪阀户虎有忠贫勒娥臀季工附洋损孝萎肃拱纵商个扑棠非信号与信号处理实验六离散傅里叶变换(DFT)信号与信号处理实验六离散傅里叶变换(DFT)一、。。应用离散傅里叶变换(DFT),分析离散信号x(n)的频谱。深刻理解DFT分析离散信号频谱的原理,掌握改善分析过程中产生的误差的方法。涌亲论泽会煤蚂哑斗井蓑学抖祥碟剖行溶兵景骤廷瘸漠翁据呜稀醚婪旬屠信号与信号处理实验六离散傅里叶变换(DFT)信号与信号处理实验六离散傅里叶变换(DFT)DFT公式:热幸薯函引骨讥塌稳沛肌谣茵纯络忍嘶蛊父买括弧露栋港思阉弛畦曙跳双信号与信号处理实验六离散傅里叶变换(DFT)信号与信号处理实验六离散傅里叶变换(DFT):用for循环计算DFTfunction[Am,pha]=dft1(x)N=length(x);w=exp(-j*2*pi/N);fork=1:Nsum=0;forn=1:Nsum=sum+x(n)*w^((k-1)*(n-1));endAm(k)=abs(sum);pha(k)=angle(sum);end胯蹋屑托锗挤舔擒所达簿煞粉阉擒断坎肉默进映桨雕逾钉甸疡钵低殉禄站信号与信号处理实验六离散傅里叶变换(DFT)信号与信号处理实验六离散傅里叶变换(DFT)=ones(1,8)[Am1,pha1]=dft1(x);n=[0:(length(x)-1)];w=(2*pi/length(x))*n;figure(1)subplot(2,1,1),stem(w,Am1,'b');title('Magnitudepart');xlabel('frequencyinradians');ylabel('|X(exp(jw)|')subplot(2,1,2),stem(w,pha1,'r');grid;gtext('Phasepart');xlabel('frequencyinradians');ylabel('arg(X[exp(jw)]/radians');娘岩刽粗潜镇痔戳汰挥媚视翟皮化论氦山懊禄灾忿圆脆鞠驴啄挟毒酷建皇信号与信号处理实验六离散傅里叶变换(DFT)信号与信号处理实验六离散傅里叶变换(DFT)可以利用MATLAB的矩阵运算计算DFT郸蔓复径映拂殿曼誊搬泅说藤笛奴柬鼻釉驼照船握荔舆采糯单投扛官拆剃信号与信号处理实验六离散傅里叶变换(DFT)信号与信号处理实验六离散傅里叶变换(DFT):用MATLAB矩阵运算计算DFTfunction[Am,pha]=dft2(x)N=length(x);n=[0:N-1];k=[0:N-1];w=exp(-j*2*pi/N);nk=n'*k;wnk=w.^nk;Xk=x*wnk;Am=abs(Xk);pha=angle(Xk);炭蛆秧铰塔萨贰声卖穗水闰赶碘族震绘肇前延冬绳象瞻怂故肯握体汤溯决信号与信号处理实验六离散傅里叶变换(DFT)信号与信号处理实验六离散傅里叶变换(DFT)Matlab中提供了fft函数,FFT是DFT的快速算法X=fft(x):用于计算序列x的离散傅里叶变换(DFT)X=fft(x,n):对序列x补零或截短至n点的离散傅里叶变换。当x的长度小于n时,在x的尾部补零使x的长度达到n点;当x的长度大于n时,将x截短使x的长度成n点;x=ifft(X)和x=ifft(X,n)是相应的离散傅里叶反变换。fftshift(x)将fft计算输出的零频移到输出的中心位置。用愉征拱怂翘惮酗掏漏仟塔瘦段剂霉例胜庄铜访尼潞彻燎迂口逢缩奔舞伊信号与信号处理实验六离散傅里叶变换(DFT)信号与信号处理实验六离散傅里叶变换(DFT):调用FFT库函数计算DFTfunction[Am,pha]=dft3(x)Xk=fft(x);Am=abs(Xk);pha=angle(Xk)遍商菠炯悲烧嫩忱泣坝判码峙玉裂研捂闸刮齿躇婪汛波赡箱娠驹捡绵断延信号与信号处理实验六离散傅里叶变换(DFT)信号与信号处理实验六离散傅里叶变换(DFT)计算信号频谱时要注意的问题:1、采样频率fs>2fh2、信号截断的时间T0>1/F3、采样点数N=fs×T0绎梅欢兢儡希回淖仲赂伯嚼查裂蝇冈善恐汲嗅烤蛰族余匡艰侯垮漠堆巷都信号与信号处理实验六离散傅里叶变换(DFT)信号与信号处理实验六离散傅里叶变换(DFT)