【DSP】数字信号FFT分析

一、实验内容及要求

(1)离散信号的频谱分析:

设信号

        此信号的和0.302pi两根谱线相距很近,谱线 0.45pi 的幅度很小,请选择合适的序列长度 N 和窗函数,用 DFT 分析其频谱,要求得到清楚的三根谱线。

(2) DTMF 信号频谱分析

用计算机声卡采用一段通信系统中电话双音多频(DTMF)拨号数字 0~9的数据,采用快速傅立叶变换(FFT)分析这10个号码DTMF拨号时的频谱。

二、实验结果 & 代码注释

2.1 离散信号的频谱分析:

2.1.1 实验结果

【DSP】数字信号FFT分析_第1张图片

 

        如图所示,上图为输入信号的时域波形图。下图为要求的时域波形。0.3附件可以看到两个相距很近的两根分离的谱线。在0.45附近可以看到一根幅值较小的谱线。

        符合题目的要求。可以看到三根分立的谱线。

2.1.2 代码及注释

close all;%关闭所有的图形窗口
N=1000;   %采样点数
n=1:1:N;
x=0.001*cos(0.45*n*pi)+sin(0.3*n*pi)-cos(0.302*n*pi-pi/4);     %输入的信号
y=fft(x,N); %进行FFT变换
a=abs(y);%FFT求模
k=0:1:N-1;
w=2*pi/N*k; %数字角频率
subplot(2,1,1)
stem(n,x,'.');
title('时域波形');
xlabel('n');ylabel('x(n)');
subplot(2,1,2);
stem(w/pi,a); 
axis([0.25 0.46 0 2]); %控制坐标范围使谱线幅度合适 
title('1000DFT频谱分析'); xlabel('数字频率'); ylabel('X(k)'); 

2.2 DTMF 信号频谱分析

电话按键进行下图编码:

低频./Hz

高频./Hz

1209

1336

1477

1633

697

1

2

3

A

770

4

5

6

B

852

7

8

9

C

941

*

0

#

D

2.2.1 运行结果

【DSP】数字信号FFT分析_第2张图片

 

2.2.2 代码及分析

clear;close all;  %关闭所有的窗口
column=[1209,1336,1477,1633];  %设置高频数组
line=[697,770,852,941];         %设置低频数组
fs=8000;                           %采样率为8000
N=1024;                            %变换点数为1024
ts=1/fs;
n=0:N-1;
f=0:fs/N:fs/N*(N-1);
 
key=zeros(10,N);  %建立全部为0的矩阵
key(1,:)=cos(2*pi*column(1)*ts*n)+cos(2*pi*line(1)*ts*n); %1
key(2,:)=cos(2*pi*column(2)*ts*n)+cos(2*pi*line(1)*ts*n); %2
key(3,:)=cos(2*pi*column(3)*ts*n)+cos(2*pi*line(1)*ts*n); %3
key(4,:)=cos(2*pi*column(1)*ts*n)+cos(2*pi*line(2)*ts*n); %4
key(5,:)=cos(2*pi*column(2)*ts*n)+cos(2*pi*line(2)*ts*n); %5
key(6,:)=cos(2*pi*column(3)*ts*n)+cos(2*pi*line(2)*ts*n); %6
key(7,:)=cos(2*pi*column(1)*ts*n)+cos(2*pi*line(3)*ts*n); %7
key(8,:)=cos(2*pi*column(2)*ts*n)+cos(2*pi*line(3)*ts*n); %8
key(9,:)=cos(2*pi*column(3)*ts*n)+cos(2*pi*line(3)*ts*n); %9
key(10,:)=cos(2*pi*column(2)*ts*n)+cos(2*pi*line(4)*ts*n); %0
figure;
for i=1:9
    subplot(3,4,i) %把整个页面分为12部分,按序进行排列频谱图
    plot(f,abs(fft(key(i,:))));
    xlabel('频率(Hz)');ylabel('幅值');
    title(i);
    grid;
    subplot(3,4,10) %绘制 “0”的频谱图
    plot(f,abs(fft(key(10,:))));
    xlabel('频率(Hz)');ylabel('幅值');
    title(0);
    grid;
end

你可能感兴趣的:(资料,matlab)