武汉科技大学《信号与系统》实验四 实验报告

信息科学与工程学院 《信号与系统》 实验报告四

专业班级 电信(DB)1102 班 姓 名 X X 学 号 20111213505X

实验时间 2013 年 6 月 18 日 指导教师 陈华丽 成 绩

注：“实验记录及个人小结”部分可另附页或在背面续写 第 页

注：“实验记录及个人小结”部分可另附页或在背面续写 第 页

二： function y=x(n) %定义函数文件 y=cos(0.48*pi*n)+cos(0.52*pi*n); endx(n)(n=0:10) 8

(1)n=0:10; f1=x(n); f2=(fft(f1)); stem(n,f2) title('x(n)(n=0:10)') axis ([0,10,-2,8])

7 6 5 4 3 2 1 0 -1 -2

0

1

2

3

4

5

6

7

8

9

10

x(n)(n=0:20) 8

(2)f3=zeros(1,20); for n=0:10 a=2*n+1; f3(a)=n; end f4=x(f3); f5=fft(f4) stem(f5)

6 4 2 0 -2 -4 -6 -8

0

5

10

15

20

25

实验记录 及个人小 结(包括： 实 验 源 程 (3)n=0:100; f1=x(n); 序、注释、 f2=(fft(f1)); 结果分析 stem(n,f2) 与讨论等)

x(n)(n=0:100) 40

30

20

10

0

-10

-20

0

10

20

30

40

50

60

70

80

90

100

x(n) 的幅度曲线是否一致三： (1) N=8;n=0:N-1; x1=(n+1).*(u(n-4)-u(n))+(8-n).*(u(n-8)-u(n-4)); x2=cos(pi*n/4); x3=sin(pi*n/8); %x4=cos(8*pi.*t)+cos(16*pi.*t)+cos(20*pi.*t); y1=fft(x1,N);subplot(2,2,1); stem(n,abs(y1)); y2=fft(x2,N);subplot(2,2,2);

注： “实验记录及个人小结”部分可另附页或在背面续写

第

页

注：“实验记录及个人小结”部分可另附页或在背面续写 第 页

n=0:63;f=64; x4=cos(8*pi*n/f)+cos(16*pi*n/f)+cos(20*pi*n/f); x4k64=fft(x4,64); x4FFT变 换 N=16 subplot(3,1,3) 15 stem(n,abs(x4k64)) title('x4FFT变换 N=64') axis([0,63,0,40])10

5

0

0

2

4

6

8

10

12

14

x4FFT变 换 N=32 20

N 不同,进行的周期延拓序 列也不相同,所以傅里叶级 数也不相同,所以变换后的 结果也不相同.

15 10 5 0

0

5

10

15 x4FFT变 换 N=64

20

25

30

40 30 20 10

实验记录 0 0 10 20 30 40 50 60 及个人小 结(包括

: 实验源程 序、注释、 个人小结： 结果分析 与讨论等) 通过本次实验我首先掌握离散信号谱分析的方法：序列的傅里叶变换、离散傅

里叶级数、 离散傅里叶变换、 快速傅里叶变换， 进一步理解了这些变换之间的关系。 通过 Matlab 作图，观察到了频谱混叠的现象。其次掌握序列的傅里叶变换、离散傅 里叶级数、离散傅里叶变换、快速傅里叶变换的 Matlab 实现。最后还了解到 FFT 对连续信号和离散信号进行谱分析的方法，了解到了可能出现的分析误差及其原 因，以便在实际中正确应用 FFT。

注： “实验记录及个人小结”部分可另附页或在背面续写

第

页

注：“实验记录及个人小结”部分可另附页或在背面续写第 页