資源描述：
《實驗二離散傅立葉變換及譜分析》由會員上傳分享，免費在線閱讀，更多相關內(nèi)容在行業(yè)資料-天天文庫。
1、實驗二離散傅立葉變換及譜分析三、實驗步驟例1、本例檢驗實序列的性質(zhì)DFT[xec(n)]=Re[X(k)]DFT[xoc(n)]=Im[X(k)]設x(n)=10*(0.8).^n0<=n<=10將x(n)分解為共扼對稱及共扼反對稱部分實驗結(jié)果：例2、本例為計算序列的圓周卷積程序，運行之前應在命令窗口輸入x1,x2,N的值。程序：N=8;x1=[1,5,-2,1];x2=[2,2,1,-1];iflength(x1)>Nerror('Nmustbe>=thelengthofx1')endiflength(x2)>Nerror('Nmust
2、be>=thelengthofx2')endx1=[x1zeros(1,N-length(x1))];x2=[x2zeros(1,N-length(x2))];m=[0:1:N-1];x2=x2(mod(-m,N)+1);H=zeros(N,N);forn=1:1:NH(n,:)=cirshftt(x2,n-1,N);endy=x1*H'實驗結(jié)果：y=21272-53-1例3、本例驗證采樣定理令，繪制其傅立葉變換。用不同頻率對其進行采樣，分別畫出離散時間傅立葉變換。（1）f=5k時，實驗結(jié)果：（2）f=1k時，程序：Dt=0.00005;
3、%步長為0.00005st=-0.005:Dt:0.005;xa=exp(-1000*abs(t));%取時間從-0.005s到0.005s這段模擬信號Wmax=2*pi*2000;%信號最高頻率為2*2000K=500;%頻域正半軸取500個點進行計算k=0:1:K;W=k*Wmax/K;%求模擬角頻率Xa=xa*exp(-j*t'*W)*Dt;%計算連續(xù)時間傅立葉變換（利用矩陣運算實現(xiàn)）Xa=real(Xa);%取實部W=[-fliplr(W),W(2:501)];%將角頻率范圍擴展為從-到+Xa=[fliplr(Xa),Xa(2:5
4、01)];subplot(2,2,1);plot(t*1000,xa);%畫出模擬信號，橫坐標為時間（毫秒），縱坐標為幅度xlabel('time(millisecond)');ylabel('xa(t)');title('anologsignal');subplot(2,2,2);plot(W/(2*pi*1000),Xa*1000);%畫出連續(xù)時間傅立葉變換xlabel('frequency(kHZ)');%橫坐標為頻率（kHz）ylabel('xa(jw)');%縱坐標為幅度title('FT');%下面為采樣頻率5kHz時的程序T
5、s=0.001;%采樣間隔為n=-25:1:25;x=exp(-1000*abs(n*Ts));%離散時間信號K=500;k=0:1:K;w=pi*k/K;%w為數(shù)字頻率X=x*exp(-j*n'*w);%計算離散時間傅立葉變換（序列的傅立葉變換）X=real(X);w=[-fliplr(w),w(2:K+1)];X=[fliplr(X),X(2:K+1)];subplot(2,2,3);stem(n*Ts*1000,x);%畫出采樣信號（離散時間信號）xlabel('time(millisecond)');gtext('Ts=0.2ms
6、');%該語句可以將引號中的內(nèi)容放置在figure中的任何地方，只需%將十字的中心放在想放置內(nèi)容的地方，然后按鼠標即可。ylabel('x1(n)');title('discretesignal');subplot(2,2,4);plot(w/pi,X);%畫出離散時間傅立葉變換xlabel('frequency(radian)');%橫坐標為弧度ylabel('x1(jw)');title('DTFT');實驗結(jié)果：例4、本例說明補零序列的離散傅立葉變換序列，已給出序列的傅立葉變換程序和將原序列補零到10長序列的DFT。（1）補零到10
7、長，實驗結(jié)果：（2）補零到20長，程序：n=0:4;x=[ones(1,5)];%產(chǎn)生矩形序列k=0:999;w=(pi/500)*k;X=x*(exp(-j*pi/500)).^(n'*k);%計算離散時間傅立葉變換Xe=abs(X);%取模subplot(3,2,1);stem(n,x);ylabel('x(n)');%畫出矩形序列subplot(3,2,2);plot(w/pi,Xe);ylabel('
8、X(ejw)
9、');%畫出離散時間傅立葉變換N=20;x=[ones(1,5),zeros(1,N-5)];%將原序列補零為20長
10、序列n=0:1:N-1;X=dft(x,N);%進行DFTmagX=abs(X);k=(0:length(magX)'-1)*N/length(magX);subplot(3,2,3);ste
