文档介绍：姓名:谢英明 班级:学号:
实验题目:用FFT对信号作频谱分析
实验目的
学****用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析 误差及其原因，以便正确应用FFT。
实验原理与方法
 用FFT对信号作频谱为x3(n）=x2((n+3））8R8（n），所以，x3(n)与x2(n)的8点DFT的模相等,如图（2a）和（2b），但是当N=16时，x2(n）与 x3(n）不满足循环移位关系，模值不相等。下图所示：
对以下周期序列进行谱分析。
选择FFT的变换区间N为8和16 两种情况分别对以上序列进行频谱分析。分别打印其幅频特性曲线。并进行对比、分析和讨论。
解：>> syms n a N
for i=1:4
disp('请输入变换区间N:');
b=input('N=');
N=b;
n=0:N-1;
disp('请输入序列xn:');
a=input('xn=','s');%输入字符串x1n
xn=eval(a); %字符串转数字
xk=fft(xn,N);
subplot(3,2,i);
if i==1
mstem(xk);xlabel('w/pi');ylabel('幅度');title('x4n的8点DFT');grid on;
end
if i==2
mstem(xk);xlabel('w/pi');ylabel('幅度');title('x4n的16点DFT');grid on;
end
if i==3
mstem(xk);xlabel('w/pi');ylabel('幅度');title('x5n的8点DFT');grid on;
end
if i==4
mstem(xk);xlabel('w/pi');ylabel('幅度');title('x5n的16点DFT');grid on;
end
i=i+1;
end
请输入变换区间N:
N=8
请输入序列xn:
xn=cos(pi/4*n)
请输入变换区间N:
N=16
请输入序列xn:
xn=cos(pi/4*n)
请输入变换区间N:
N=8
请输入序列xn:
xn=cos(pi/4*n)+cos(pi/8*n)
请输入变换区间N:
N=16
请输入序列xn:
xn=cos(pi/4*n)+cos(pi/8*n)
分析：对于x4(n）=cos(pi*n/4),周期T=8，N=8和16是周期的整数倍，所以得到正确的单一频率正弦波的频谱，
π处有1根单一谱线。如图（4a）和（4b）所示 
对于x5n=cos(pi*n/4)+cos(pi*n/8)，周期T=16，但是N=8不是该周期的整数倍，所以8点的dft和16点的dft不一样。
（3）对模拟周期信号进行谱分析： 选择采样频率Fs=64Hz变换区间N=16,32,64 三种情况进行谱分析。分别打印其幅频特性，并进行分析和讨论。
解：>> syms n a N k
disp('请选择采样频率:');
b=input('Fs=');
Fs=b;
T=1/Fs;
for i=1:3
disp('请输入变换区间N