正弦频率调制傅里叶变换

afgc223

于 2023-10-18 20:53:34 发布

阅读量369

点赞数

文章标签：学习傅立叶分析算法

本文链接：https://blog.csdn.net/afgc223/article/details/133914441

版权

文章目录

1、原理
2、仿真
3、总结
4、参考

1、原理

相比于傅里叶变换，正弦频率调制傅里叶变换的原理本质上是对信号的相位项进行傅里叶变换，变换的目的是将原始信号的相位项分解为傅里叶级数的展开，假设离散的原始信号为x(n)，可以表示为
在这里插入图片描述

系数ak可以通过傅里叶级数的正交性得到，即
在这里插入图片描述

2、仿真

仿真信号的采样率设置为2500Hz，采样时间长度0.1s，所采样的信号为

s_t=exp(1j*2*pi*cos(2*pi*100*t+pi/3)).*exp(1j*4*pi*sin(2*pi*50*t));

通过相位解缠后所提取的信号相位与原始信号的理论相位对比结果如下
在这里插入图片描述

两条曲线能够很好的吻合，计算信号的正弦频率调制傅里叶变换的幅度谱为
在这里插入图片描述

代码如下：

%% 用于离散信号正弦调频傅里叶变换
close all; clear; clc;
%% 参数设置
T=0.1;  %设置信号长度
fs=2500;%设置采样率
%t=(0.5*T:1/fs:1.5*T);
t=0:1/fs:T;
f=linspace(-fs/2,fs/2,length(t));
%% 信号采样
s_t=exp(1j*2*pi*cos(2*pi*100*t+pi/3)).*exp(1j*4*pi*sin(2*pi*50*t));
%% 正弦调频傅里叶变换
ln_s_t=unwrap(real(-1j*log(s_t)));
ln_s_t=ln_s_t-mean(ln_s_t);
figure;plot(t,ln_s_t,t,2*pi*cos(2*pi*100*t+pi/3)+4*pi*sin(2*pi*50*t),'--','linewidth',2);
a_k=fftshift(fft(ln_s_t))/length(f);
figure;plot(f,abs(a_k));

3、总结

该变换比较关键的一点是需要解决信号的相位解缠问题，这需要信号的采样率需要满足一定的条件。

4、参考

B. Peng, X. Wei, B. Deng, H. Chen, Z. Liu and X. Li, “A Sinusoidal Frequency Modulation Fourier Transform for Radar-Based Vehicle Vibration Estimation,” in IEEE Transactions on Instrumentation and Measurement, vol. 63, no. 9, pp. 2188-2199, Sept. 2014, doi: 10.1109/TIM.2014.2308031.