应用FFT实现信号频谱分析.doc

文档介绍：应用FFT实现信号频谱分析
应用FFT实现信号频谱分析
应用FFT实现信号频谱分析
信号、系统与信号处理实验Ⅱ
实验报告
实验名称:应用FFT实现信号频谱分析
实验目的
1、能够熟练掌握快速离散傅里叶变换FFT的原理及应用FFgXK1)-1;
stem(k,magXK1);
xlabel('k');
ylabel('|X(k)|');
title('X(k) N = 20');
subplot(4,1,3);
k=0:length(magXK2)-1;
stem(k,magXK2);
xlabel('k');
ylabel('|X(k)|');
title('X(k) N = 40');
subplot(4,1,4);
k=0:length(magXK3)-1;
stem(k,magXK3);
xlabel('k');
ylabel('|X(k)|');
title('X(k) N = 22');
复合函数序列:v3(n) = 0、9sin(2*pi*n/N) + 0、6sin(2*pi*n*3/N)。
clear all
应用FFT实现信号频谱分析
应用FFT实现信号频谱分析
应用FFT实现信号频谱分析
N = 40;
n = 0:N-1;
xn = 0、9*sin(2*pi、*n/N)+0、6*sin(2*pi*n*3/N);
XK = fft(xn,N);
magXK = abs(XK);
phaXK = angle(XK);
subplot(1,2,1);
stem(n,xn);
xlabel('n');
ylabel('x(n)');
title('x(n) N =40');
subplot(1,2,2);
k = 0:length(magXK)-1;
stem(k,magXK);
xlabel('k');
ylabel('|X(k)|');
title('X(k) N =40');
已知信号x(t) = 0、15sin(2*pi*f1*t)+sin(2*pi*f2*t)-0、1sin(2*pi*f3*t),其中f1=1Hz,f2=2Hz,f3=3Hz,取样频率为32Hz。① 32 点FFT,画出其幅度谱。② 64 点FFT,画出其幅度谱,比较两者间的差异,思考实际频率与离散频谱图中横坐标k的对应关系。
clear all
N=32;
N1=32;
N2=64;
n=0:N-1;
t=n/N;
xn=0、15*sin(2*pi*t)+sin(4*pi*t)-0、1*sin(6*pi*t);
XK1=fft(xn,N1);
应用FFT实现信号频谱分析
应用FFT实现信号频谱分析
应用FFT实现信号频谱分析
magXK1=abs(XK1);
phaXK1=angle(XK1);
XK2=fft(xn,N2);
magXK2=abs(XK2);
phaXK2=angle(XK2);
subplot(3,1,1);
stem(n,xn);
xlabel('n');
ylabel('x(n)');
title('x(n)');
subplot(3,1,2);
k=0:length(magXK1)-1;