八、DFT的原理及實現(xiàn)VIP

實驗主題DFT的原理與實現(xiàn)實驗時間實驗臺編號實驗的目的DFT算法原理序列的DFT運算實驗儀器PC臺式機，軟件MATLAB預(yù)習實驗實驗預(yù)習：1、DTFT和IDTFT式；2.DFT轉(zhuǎn)換公式3、DFT的性質(zhì)實驗原理：1、DTFT :改寫如下：2 .反變換IDTFT :3、DFT :實驗步驟函數(shù)for ()；首先，使用定義法定義DFT運算函數(shù)1、二次循環(huán)法：function X=DFTfor(xn )N=length(xn )X=zeros(1，n )for k=0:N-1；for n=0:N-1；x (k1)=x (k1) xn (n1) * exp (-j *2* pi * n * k/n )；結(jié)束結(jié)束2、矩陣乘法運算DFT :function Xk=DFTmat(xn )N=length(xn )n=0:N-1；k=n；nk=n*k；WN=exp(-j*2*pi/N )； Wnk=WN.nk；Xk=xn*Wnk；主題8-1 :你計算11點的DFT嗎？n=-5:5； xn=(-0.9).n；k1=03336508； w=(pi/500)*k1；X=xn*(exp(-j*pi/500).(n*k1) %求連續(xù)譜magX=abs(X) angX=angle(X )Xk=DFTmat(xn) %調(diào)用矩陣乘法N=length(xn) k=0:N-1；Xk1=Xk.*exp(j*2*pi*5*k/N) %xn在實際取值范圍內(nèi)求解magXk=abs(Xk1) angXk=angle(Xk1 )subplot (2，1，1 ) plot (w/pi、magX、-)hold on； stem(2*k/N，magXk )hold off； axis ( 0，2，0，15 ) grid on；xlabel(omega(xpi ) ) ylabel (寬度)subplot (2，1，2 ) plot (w/pi、angX/pi、-)hold on； stem(2*k/N，angXk/pi )hold off； axis ( 0，2，- 1，1 ) grid on；xlabel(omega(xpi ) ) ylabel (相位(弧度/pi )； %結(jié)果8-1主題8-2 :二、用直接法計算其四點DFT？n=0:3； xn= 1，1，1，1 ；k1=03336508； w=(pi/500)*k1；X=xn*(exp(-j*pi/500).(n*k1)magX=abs(X) angX=angle(X )N=length(xn )n1=0:N-1； k=n1； nk=n1*k；WN=exp(-j*2*pi/N )； Wnk=WN.nk；Xk=xn*Wnk；magXk=abs(Xk) angXk=angle(Xk )subplot (2，1，1 ) plot (w/pi，magX，k- )；hold on； stem(2*k/N，magXk )hold off； axis ( 0，2，0，5 ) grid on；xlabel(omega(xpi ) ) ylabel (寬度|X(k)|)subplot (2，1，2 ) plot (w/pi，angX/pi，k- )；hold on； stem(2*k/N，angXk/pi )hold off； axis ( 0，2，- 1，1 ) grid on；xlabel(omega(xpi ) ) ylab