文章目录
原信号: x ( t ) = 2.5 + s i n ( 2 π ∗ 100 t ) + 3 s i n ( 2 π ∗ 900 t ) x(t)=2.5+sin(2\pi *100t)+3sin(2\pi *900t) x(t)=2.5+sin(2π∗100t)+3sin(2π∗900t)
1. 确定最小采样频率
- 根据奈奎斯特采样定理保证不混叠且能还原出原信号信息,确定最小采样频率 f s m i n = 1800 H z f_{smin}=1800 Hz fsmin=1800Hz。
2. 确定最小采样总时间
- 为了能够区分出原信号中的三个信号分量:0Hz、100Hz以及9000Hz。首先确定频率分辨率最大为 100Hz,则最小采样总时间 L m i n = 1 / 100 秒 L_{min} =1/100 秒 Lmin=1/100秒。
- 频率分辨率的选取举例:如果有三个信号分量,分别是 201Hz、202Hz、500Hz,那么频率分辨率最大为1Hz,这是为了能够区分出 201Hz 和 202 Hz的信号分量。
3. 确定最小采样点数
- 根据 L = N / f s L=N/f_s L=N/fs,得到最小采样点数 N m i n = L m i n × f s m i n = 1800 / 100 = 18 N_{min}=L_{min}\times f_{smin}=1800/100=18 Nmin=Lmin×fsmin=1800/100=18。
- 基2-FFT要求点数N为2的指数倍(不考虑补零或截断,采样点数直接等于FFT运算点数),因此可得到大于 N m i n N_{min} Nmin的 N 的可取值序列 N ( n ) = 32 , 64 , 128 , 256.... N(n)={32,64,128,256....} N(n)=32,64,128,256....。
4. 确定采样频率
- 为了保证频率和幅度的准确性(博主理解为在该段采样时间内获取的都是整数倍原周期信号,频谱尽可能不泄露),需满足
m 1 × f s / N = f 1 = 100 , m 1 ∈ 正整数 m_1\times f_s/N=f_{1}=100,m_1∈正整数 m1×fs/N=f1=100,m1∈正整数
m 2 × f s / N = f 2 = 900 , m 2 ∈ 正整数 m_2\times f_s/N=f_{2}=900,m_2∈正整数 m2×fs/N=f2=900,m2∈正整数 - 联立上两式得到 m 2 m 1 = 9 \frac{m_2}{m_1}=9 m1m2=9
- 对 m 1 、 m 2 和 N m_1、m_2和N m1、m2和N依次取值,凑数。
令 m 1 = 1 , m 2 = 9 , N = 32 m_1=1,m_2=9,N=32 m1=1,m2=9,N=32,得到