MATLAB之频谱分析（七）

时域信号频谱分析

如果需要对时域内信号的频谱进行分析，需要了解几个基本的概念：

时域：时域（Time domain）是描述数学函数或物理信号对时间的关系，比如说一个数学表达式：$y=t^2+t$,这个表达式y是随着t的改变而改变，也就是说因变量y随着自变量t的变化而变化。这个式子也就是在时域上的。再例如一个信号的时域波形可以表达信号随着时间的变化。时域是真实世界，是惟一实际存在的域。因为我们的经历都是在时域中发展和验证的，已经习惯于事件按时间的先后顺序地发生。而评估数字产品的性能时，通常在时域中进行分析，因为产品的性能最终就是在时域中测量的。

频域：频域（frequency domain）是描述信号在频率方面特性时用到的一种坐标系。在电子学，控制系统工程和统计学中，频域图显示了在一个频率范围内每个给定频带内的信号量。频域，尤其在射频和通信系统中运用较多，在高速数字应用中也会遇到频域。频域最重要的性质是：它不是真实的，而是一个数学构造。时域是惟一客观存在的域，而频域是一个遵循特定规则的数学范畴。

正弦波是频域中唯一存在的波形，这是频域中最重要的规则，即正弦波是对频域的描述，因为频域中的任何波形都可用正弦波合成。这是正弦波的一个非常重要的性质。然而，它并不是正弦波的独有特性，还有许多其他的波形也有这样的性质。正弦波有四个性质使它可以有效地描述其他任一波形：

（1）频域中的任何波形都可以由正弦波的组合完全且惟一地描述。

（2）任何两个频率不同的正弦波都是正交的。如果将两个正弦波相乘并在整个时间轴上求积分，则积分值为零。这说明可以将不同的频率分量相互分离开。

（3）正弦波有精确的数学定义。

（4）正弦波及其微分值处处存在，没有上下边界。

使用正弦波作为频域中的函数形式有它特别的地方。若使用正弦波，则与互连线的电气效应相关的一些问题将变得更容易理解和解决。如果变换到频域并使用正弦波描述，有时会比仅仅在时域中能更快地得到答案。

关系：时域分析是以时间轴为坐标表示动态信号的关系；频域分析是把信号变为以频率轴为坐标表示出来。

转换：动态信号从时间域变换到频率域主要通过傅立叶级数和傅立叶变换实现。周期信号靠傅立叶级数，非周期信号靠傅立叶变换。时域越宽，频域越短。

采样频率：采样频率，也称为采样速度或者采样率，定义了单位时间内从连续信号中提取并组成离散信号的采样个数，它用赫兹（Hz）来表示。采样频率的倒数是采样周期或者叫作采样时间，它是采样之间的时间间隔。通俗的讲采样频率是指计算机单位时间内能够采集多少个信号样本。

奈奎斯特采样定理（香农定理）：采样是将一个信号（即时间或空间上的连续函数）转换成一个数值序列（即时间或空间上的离散函数）。采样定理指出，如果信号是无限的，并且采样频率高于信号带宽的两倍，那么，原来的连续信号可以从采样样本中完全重建出来。

带限信号变换的快慢受到它的最高频率分量的限制，也就是说它的离散时刻采样表现信号细节的能力是有限的。

采样定理是指，如果信号带宽不到采样频率的一半（即奈奎斯特频率），那么此时这些离散的采样点能够完全表示原信号。高于或处于奈奎斯特频率的频率分量会导致混叠现象。大多数应用都要求避免混叠，混叠问题的严重程度与这些混叠频率分量的相对强度有关。

采样频率必须大于被采样信号带宽的两倍，另外一种等同的说法是奈奎斯特定律必须大于被采样信号的带宽。如果信号的带宽是100Hz，那么为了避免混叠现象采样频率必须大于200Hz。换句话说就是采样频率必须至少是信号中最大频率分量频率的两倍，否则就不能从信号采样中恢复原始信号。

这里只对香农定理概念作一定的描述，具体的内容可以参考《信号与系统》，里面有详细的描述。

实例：假设数据采样频率为1400Hz,一个信号包含频率为65Hz,振幅为0.75的正弦波和一个频率为145Hz，振幅为0.8的正弦波，噪声为零平均值的随机噪声，利用傅里叶变换分析其频谱。

Fs =1400; %采样频率
T= 1/Fs; %采样时间
L= 1400; %信号长度
t = (0:L-1)*T %时间向量
x = 0.75*sin(2*pi*65*t)+0.8*sin(2*pi*145*t); %有用信号
y = x+2*randn(size(t)); %加噪声正弦信号
plot(Fs*t(1:50),y(1:50)) %以时间为x轴，零平均值噪声信号为y轴画出信号
title("零平均噪声信号");
xlabel("time");
NFFT = 2^nextpow2(L);
Y = fft(y,NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2);
plot(f,2*abs(Y(1:NFFT/2)))
title("单边振幅频谱")
xlabel('frequency')
ylabel('|Y(f)|')

运行结果：