通信实验报告.docx
《通信实验报告.docx》由会员分享，可在线阅读，更多相关《通信实验报告.docx（9页珍藏版）》请在冰豆网上搜索。
通信实验报告
通信系统综合设计实训实验报告
姓名：
王福海
班级：
01011406
学号：
2014210509
实验时间：
2016-09-29
指导老师：
郑丹玲
实验用FFT进行谱分析
一、实验原理及目的。
实验内容：
1.编写matlabM文件对信号做8点和16点的FFT.
clc
clearall
closeall
N=8;%定义8点FFT长度
x=[11110000];%表示离散序列x1
xk=fft（x,N）;%调用FFT函数
figure;%显示窗口
subplot（211）;%将窗口分成两部分显示第一部分
stem（0:
length（x）-1,x,'.'）;%表示x图形
title（'x1的波形'）;%标题
subplot（212）;%显示第二个部分
stem（0:
N-1,abs（xk）,'.'）;%表示x的8点fft图形
title（'x1的8点FFT'）;%标题
clc
clearall
closeall
N=16;%定义16点FFT长度
x=[11110000];%表示离散序列x1
xk=fft（x,N）;%调用FFT函数
figure;%显示窗口
subplot（211）;%将窗口分成两部分显示第一部分
stem（0:
length（x）-1,x,'.'）;%表示x图形
title（'x1的波形'）;%标题
subplot（212）;%显示第二个部分
stem（0:
N-1,abs（xk）,'.'）;%表示x的8点fft图形
title（'x1的16点FFT'）;%标题
2.编写matlabM文件对信号
做8点和16点的FFT.
clc
clearall
closeall
N=8;%定义8点FFT长度
x=[1:
44:
-1:
1];%表示离散序列x2
xk=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点FFT'）;
N=16;
x=[1:
44:
-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.编写matlabM文件对信号做8点和16点的FFT.
clc
clearall
closeall
N=8;
n=0:
7;
x4=cos（pi/4*n）;%cos离散函数表示
xk=fft（x4,N）;%调用FFT函数
figure
（1）;
subplot（211）;
stem（0:
length（x4）-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.、编写matlabM文件对信号
以fs=64（Hz）采样后做N=16、32、64点的FFT.
clc
clearall
closeall
f=64;
N=32;%定义fft的长度可以为16、32或64
n=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'）;
二、结合实验中所得给定典型序列幅频特性曲线，与理论结合比较，并分析说明误差产生的原因以及用FFT作谱分析时有关参数的选择方法。
三、总结实验所得主要结论。
四、简要回答思考题。
思考题
1、在N=8和N=16两种情况下x2（n）,x3（n）的幅频特性会相同吗？
为什么？
2、如果周期信号的周期预先不知道，如何用FFT进行分析？
3、试使用函数fft（x）近似画出x（n）=R10（n）在（-4π,4π）上的幅频响应曲线|FT[（X（n）]|。
clc
clearall
closeall
N=512;%定义点FFT长度
x=[1111111111];%表示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*pi4*pi010]）;%框定表示范围
title（'幅频响应曲线|FT[（X（n）]|'）;%标题
升级会员 