FFT对连续信号进行谱分析.docVIP

二○一二～二○一三学年第二学期

电子信息工程系

课程设计方案书

班级：电信1103班

学号：201113136068

姓名：陈凯

课程名称：信号与系统课程设计

学时学分：1周1学分

指导教师：盛玉霞

二○一三年六月十七日

1、课程设计目的：

综合运用本课程的理论知识对信号进行采样，频谱分析，掌握线性卷积和循环卷积的区别和联系，利用MATLAB编程实现，并通过理论推导得出相应结论，从实践上初步实现对信号的处理

2、课程设计时间安排：

第一天：教师布置设计任务，讲解设计要求，提示设计要点。

第二天：查阅资料〔在图书馆或上网〕，弄清题目要求，提出解决方案。

第三～五天：根据题目要求，将理论推导与编程实现相结合，写出设计报告。

第六天：验收。

3、课程设计内容及要求：

题目一、FFT对连续信号进行谱分析

假设信号由20Hz，220Hz和750Hz三个正弦频率成分构成，试用FFT计算其频谱并画出其幅度谱〔要求横坐标以Hz为单位〕，验证信号包含的频率成分。

实验内容包括：

Matlab实现fft并画出幅度谱。

c语言编程实现FFT算法。

改变采样点数分析对结果的影响。

分析：

设

由于是周期信号，首先要确定该信号的周期，如果截取信号时没有至少一个完整周期的信号，那么会发生泄漏现象。

∵

∴信号的周期为

因此，截取信号的持续时间最好为0.1s的整数倍。

信号的最高频率为，因此，采样频率或采样间隔需满足

，

取采样频率和采样间隔分别为：,

DFT的点长N从持续时间中算得

令t=nT,n=0,1,2,…,255,代入原信号，得

利用FFT计算200点的DFT。

matlab实现fft以及频谱图

源程序：n=0:255;

x=0.25*sin(pi*0.02*n)+0.5*sin(pi*0.2*n)+sin(pi*0.75*n);

c=fft(x,200);

m=abs(c);

subplot(2,1,1);

stem(x);

holdon;

plot(n,x,g--);

axis([040-22])

title(x=0.25sin(pi*0.02n)+0.5sin(pi*0.2n)+sin(pi*0.75n));

gridon

subplot(2,1,2);

stem(m);

gridon

运行结果

〔2〕C语言实现FFT算法

源程序如下：

#includeiom128.h

#includeintrinsics.h

#includemath.h

#defineFFT_N128

structcompx{floatreal,imag;};

structcompxs[FFT_N];

floatSIN_TAB[FFT_N/4+1];

structcompxEE(structcompxa,structcompxb)

{

structcompxc;

c.real=a.real*b.real-a.imag*b.imag;

c.imag=a.real*b.imag+a.imag*b.real;

return(c);

}

voidcreate_sin_tab(float*sin_t)

{

inti;

for(i=0;i=FFT_N/4;i++)

sin_t[i]=sin(2*PI*i/FFT_N);

}

floatsin_tab(floatpi)

{

intn;

floata;

n=(int)(pi*FFT_N/2/PI);

if(n=0n=FFT_N/4)

a=SIN_TAB[n];

elseif(nFFT_N/4nFFT_N/2)

{

n-=FFT_N/4;

a=SIN_TAB[FFT_N/4-n];

}

elseif(n=FFT_N/2n3*FFT_N/4)

{

n-=FFT_N/2;

a=-SIN_TAB[n];

}

elseif(n=3*FFT_N/4n3*FFT_N)

{

n=FFT_N-n;

a=-SIN_TAB[n];