【PRTF 变换】参数化重采样时频变换(Matlab代码实现）

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

2.1 P-PRTF变换的数值示例

2.2 W-PRTF 变换的数值示例

🎉3 参考文献

🌈4 Matlab代码、文章

---

💥1 概述

摘要：
许多信号包含多个具有时变瞬时频率（IFs）的组成部分，它们共同具有一个轨迹趋势，例如机械振动信号、语音信号和生物医学信号。为了分析这种信号并获得高时频分辨率，本文提出了一种称为参数化重采样时频变换（PRTF 变换）的方法。借鉴通用参数化时频变换（GPTF 变换）的思想，我们使用参数化核函数来表示重采样函数，并进一步构建时变和时不变的重采样算子来消除 IF 的变化并重新定位 IF 的位置。这些算子可以同时提高时频表示（TFR）中多个组成部分的能量集中度。针对不同类型的信号，提供了包含多项式函数和傅里叶级数的典型核函数。提出了相应的核估计方法，通过利用多个组成部分检测 IF 的共享趋势，并递归逼近核参数。数值模拟和实际实验证明了该方法在提高非平稳多组分信号的时频分辨率方面的有效性。

详细文章见第四部分。

📚2 运行结果

2.1 P-PRTF变换的数值示例

2.2 W-PRTF 变换的数值示例

部分代码：

%%%%%%%%%%%%%%%%%%%%%%%  input %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   Sig:                The signal to be analyzed
%   SampFreq:     Sampling frequency 
%   coef_f:            Coefficients of polynomial
%   WinLen:          The length of window used to locate the signal in time.
%   overlap:          Overlap rate of the window
%%%%%%%%%%%%%%%%%%%%%%% output %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   Spec:              The absolute value of the spectrogram result of the transform
%   f:                    Frequency axis

%%  Initializing
f_wait = waitbar(0,'1',"Name","P-PRTF transform");

Sig = Sig(:);
SigLen = length(Sig);
Sig = real(Sig);
t = (0:SigLen-1)/SampFreq;

%   Calculate parameters
coef_r = coef_f;
shift = ceil((1-overlap)*WinLen);
for i = 1:length(coef_f)
    coef_r(i) = coef_f(i)/i;
end
coef_r = [0 coef_r];
coef_r = flip(coef_r);
coef_f = flip(coef_f);

%   Zero padding
r_t = polyval(coef_r,t);
f_t = polyval(coef_f,t);

🎉3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

🌈4 Matlab代码、文章