1、通信实验报告通信系统综合设计实训实验报告姓名:王福海班级:01011406学号:2014210509实验时间:2016-09-29指导老师:郑丹玲实验 用 FFT 进行谱分析一、实验原理及目的。 实验内容:1.编写 matlab M 文件对信号做8点和16 点的 FFT.clcclear allclose allN=8; %定义8点FFT长度 x=1 1 1 1 0 0 0 0; %表示离散序列x1xk=fft(x,N); %调用FFT函数figure; %显示窗口 subplot(211); %将窗口分成两部分显示第一部分stem(0:length(x)-1,x,.); %表示x图形titl
2、e(x1的波形); %标题subplot(212); %显示第二个部分stem(0:N-1,abs(xk),.); %表示x的8点fft图形title(x1的8点FFT); %标题clcclear allclose allN=16; %定义16点FFT长度 x=1 1 1 1 0 0 0 0;%表示离散序列x1xk=fft(x,N); %调用FFT函数figure; %显示窗口subplot(211); %将窗口分成两部分显示第一部分stem(0:length(x)-1,x,.);%表示x图形title(x1的波形); %标题subplot(212); %显示第二个部分stem(0:N-1,a
3、bs(xk),.); %表示x的8点fft图形title(x1的16点FFT); %标题2.编写matlab M 文件对信号做 8 点和 16 点的 FFT.clcclear allclose allN=8; %定义8点FFT长度 x=1:4 4:-1:1;%表示离散序列x2xk=fft(x,N); %调用FFT函数figure(1); %显示第一个窗口subplot(211); %将窗口分成两部分显示第一部分stem(0:length(x)-1,x,.); %表示x2图形title(x2的波形); subplot(212);stem(0:N-1,abs(xk),.);title(x2的8点F
4、FT); N=16;x=1:4 4:-1:1;xk=fft(x,N);figure(2);subplot(211);stem(0:length(x)-1,x,.);title(x2的波形);subplot(212);stem(0:N-1,abs(xk),.);title(x2的16点FFT);3.编写matlab M 文件对信号做 8 点和 16 点的FFT.clcclear allclose allN=8;n=0:7;x4=cos(pi/4*n); %cos离散函数表示xk=fft(x4,N); %调用FFT函数 figure(1); subplot(211);stem(0:length(x
5、4)-1,x4,.); %长度定义title(x4的波形); subplot(212); stem(0:N-1,abs(xk),.);title(x4的8点FFT); N=16;n=0:15;x4=cos(pi/4)*n); %cos离散函数表示xk=fft(x4,N); %调用FFT函数figure(2); subplot(211); stem(0:length(x4)-1,x4,.); title(x4的波形);subplot(212);stem(0:length(x4)-1,abs(xk),.);title(x4的16点FFT);4.、编写 matlab M 文件对信号以fs=64(Hz
6、)采样后做 N=16、32、64点的FFT.clcclear allclose allf=64;N=32; %定义fft的长度可以为16、32或64n=0:N-1;t=n/f; %离散化,确定采样点x=cos(8*pi*t)+cos(16*pi*t)+cos(20*pi*t);%表示出采样后函数xk=fft(x,N);figure;subplot(211);stem(0:length(x)-1,x,.);title(x6的波形);subplot(212);stem(0:N-1,abs(xk),.);title(x6的32点FFT);二、结合实验中所得给定典型序列幅频特性曲线,与理论结合比较,并
7、分析说明误差产生的原因以及用 FFT作谱分析时有关参数的选择方法。 三、总结实验所得主要结论。四、简要回答思考题。 思考题1、在 N=8 和 N=16 两种情况下x2(n),x3(n)的幅频特性会相同吗?为什么? 2、如果周期信号的周期预先不知道,如何用 FFT 进行分析? 3、试使用函数fft(x)近似画出x(n)=R10(n)在(-4,4)上的幅频响应曲线|FT(X(n)|。clcclear allclose allN=512; %定义点FFT长度 x=1 1 1 1 1 1 1 1 1 1; %表示x序列y1=fft(x,N); %fftN点变换y2=y1,y1,y1,y1; %表示y2 plot(linspace(-4*pi,4*pi,4*N),abs(y2); %显示y2在-4到4上的axis(-4*pi 4*pi 0 10); %框定表示范围title(幅频响应曲线|FT(X(n)|); %标题
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1