参数化重采样时频变换（PRTF变换）附matlab代码

​✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，matlab项目合作可私信。

🍎个人主页：Matlab科研工作室

🍊个人信条：格物致知。

更多Matlab仿真内容点击👇

智能优化算法       神经网络预测       雷达通信       无线传感器        电力系统

信号处理              图像处理               路径规划       元胞自动机        无人机 

⛄ 内容介绍

许多信号包含多个具有时变瞬时频率 (IF) 的分量，这些分量具有共同的轨迹趋势，例如机械振动信号、语音信号和生物医学信号。为了分析此类信号并实现高时频分辨率，本文提出了一种称为参数化重采样时频变换（PRTF变换）的方法。采用一般参数化时频变换（GPTF变换）的思想，我们使用参数化核来表示重采样函数，并进一步构造时变和时不变重采样算子来消除中频变化并重新定位中频位置。这些算子可以同时提高时频表示（TFR）中多个分量的能量集中度。

⛄ 核心代码

<span style="color:#333333"><span style="background-color:#fafafa"><code>%   Numerical example <span style="color:#ca7d37">for</span> W-PRTF transform</code><code>​</code><code>clc</code><code>clear</code><code><span style="color:#ca7d37">close</span> all</code><code>​</code><code>%% Initialize</code><code>SampFreq = <span style="color:#0e9ce5">5000</span>;</code><code>t = <span style="color:#0e9ce5">0</span>:<span style="color:#0e9ce5">1</span>/SampFreq:<span style="color:#0e9ce5">1</span>;</code><code>​</code><code>RTF = <span style="color:#0e9ce5">200</span>*(<span style="color:#0e9ce5">1.5</span>+<span style="color:#ca7d37">sin</span>(<span style="color:#0e9ce5">10</span>*t)+<span style="color:#ca7d37">sin</span>(<span style="color:#0e9ce5">30</span>*t)/<span style="color:#0e9ce5">3</span>+<span style="color:#ca7d37">sin</span>(<span style="color:#0e9ce5">50</span>*t)/<span style="color:#0e9ce5">5</span>+<span style="color:#ca7d37">sin</span>(<span style="color:#0e9ce5">70</span>*t)/<span style="color:#0e9ce5">7</span>);</code><code>phase = <span style="color:#0e9ce5">200</span>*(<span style="color:#0e9ce5">1.5</span>*t-<span style="color:#ca7d37">cos</span>(<span style="color:#0e9ce5">10</span>*t)/<span style="color:#0e9ce5">10</span>-<span style="color:#ca7d37">cos</span>(<span style="color:#0e9ce5">30</span>*t)/<span style="color:#0e9ce5">9</span>/<span style="color:#0e9ce5">10</span>-<span style="color:#ca7d37">cos</span>(<span style="color:#0e9ce5">50</span>*t)/<span style="color:#0e9ce5">25</span>/<span style="color:#0e9ce5">10</span>-<span style="color:#ca7d37">cos</span>(<span style="color:#0e9ce5">70</span>*t)/<span style="color:#0e9ce5">49</span>/<span style="color:#0e9ce5">10</span>);</code><code>​</code><code>R = [(<span style="color:#0e9ce5">1</span>:-<span style="color:#0e9ce5">0</span>.<span style="color:#0e9ce5">1</span>:<span style="color:#0e9ce5">0</span>)<span style="color:#dd1144">';zeros(11,1)];</span></code><code>G = [(1:-0.1:0)';(<span style="color:#0e9ce5">0</span>:<span style="color:#0e9ce5">0</span>.<span style="color:#0e9ce5">1</span>:<span style="color:#0e9ce5">1</span>)<span style="color:#dd1144">'];</span></code><code>B = [ones(11,1);(1:-0.05:0.5)'];</code><code><span style="color:#ca7d37">map</span> = [R G B];</code><code>​</code><code>%% Generate the signal</code><code>Sig = <span style="color:#0e9ce5">0</span>;</code><code><span style="color:#ca7d37">for</span> i = <span style="color:#0e9ce5">1</span>:<span style="color:#0e9ce5">6</span></code><code>    Sig = Sig + <span style="color:#ca7d37">cos</span>(<span style="color:#0e9ce5">2</span>*pi*(<span style="color:#0e9ce5">0</span>.<span style="color:#0e9ce5">6</span>+i*<span style="color:#0e9ce5">0</span>.<span style="color:#0e9ce5">4</span>*phase));</code><code>end</code><code>​</code><code>%% Waveform</code><code>figure</code><code>set(gcf,<span style="color:#dd1144">'Position'</span>,[<span style="color:#0e9ce5">20</span> <span style="color:#0e9ce5">100</span> <span style="color:#0e9ce5">320</span> <span style="color:#0e9ce5">250</span>]);      </code><code>set(gcf,<span style="color:#dd1144">'Color'</span>,<span style="color:#dd1144">'w'</span>); </code><code>plot(t,Sig,<span style="color:#dd1144">'linewidth'</span>,<span style="color:#0e9ce5">0</span>.<span style="color:#0e9ce5">5</span>,<span style="color:#dd1144">'color'</span>,<span style="color:#dd1144">'b'</span>);</code><code>xlabel(<span style="color:#dd1144">'Time (Sec)'</span>,<span style="color:#dd1144">'FontSize'</span>,<span style="color:#0e9ce5">12</span>,<span style="color:#dd1144">'FontName'</span>,<span style="color:#dd1144">'Times New Roman'</span>);</code><code>ylabel(<span style="color:#dd1144">'Amplitude'</span>,<span style="color:#dd1144">'FontSize'</span>,<span style="color:#0e9ce5">12</span>,<span style="color:#dd1144">'FontName'</span>,<span style="color:#dd1144">'Times New Roman'</span>);</code><code>set(gca,<span style="color:#dd1144">'FontSize'</span>,<span style="color:#0e9ce5">12</span>)</code><code>set(gca,<span style="color:#dd1144">'linewidth'</span>,<span style="color:#0e9ce5">1</span>);</code><code>axis([<span style="color:#0e9ce5">0</span> <span style="color:#0e9ce5">1</span> -<span style="color:#0e9ce5">8</span> <span style="color:#0e9ce5">8</span>])</code><code>​</code><code>%% Short-<span style="color:#ca7d37">time</span> frequency transform</code><code>figure</code><code>[Spec,f] = STFT(Sig<span style="color:#dd1144">',SampFreq,1024,1024);</span></code><code>imagesc(t,f,Spec); </code><code>colormap(map)</code><code>axis([0 t(end) 0 1500]);</code><code>set(gcf,'Position<span style="color:#dd1144">',[20 100 320 250]);   </span></code><code>xlabel('Time (Sec)<span style="color:#dd1144">','</span>FontSize<span style="color:#dd1144">',12,'</span>FontName<span style="color:#dd1144">','</span>Times New Roman<span style="color:#dd1144">');</span></code><code>ylabel('Frequency (Hz)<span style="color:#dd1144">','</span>FontSize<span style="color:#dd1144">',12,'</span>FontName<span style="color:#dd1144">','</span>Times New Roman<span style="color:#dd1144">');</span></code><code>set(gca,'YDir<span style="color:#dd1144">','</span>normal<span style="color:#dd1144">')</span></code><code>set(gca,'FontSize<span style="color:#dd1144">',12);</span></code><code>set(gcf,'Color<span style="color:#dd1144">','</span>w<span style="color:#dd1144">');</span></code><code>​</code><code>%% Parameterized resampling time-frequency transform with Fourier Series kernel</code><code>figure</code><code>coef = [10 1.5 0 0 0 0 0 0 0 1 0 1/3 0 1/5 0 1/7];</code><code>[Spec,f] = W_PRTFtransform(Sig',SampFreq,coef,<span style="color:#0e9ce5">1024</span>);</code><code>%[Spec,f] = W_PRTFtransform2(Sig<span style="color:#dd1144">',SampFreq,coef,1024,0.99);</span></code><code>imagesc(t,f,abs(Spec)); </code><code>colormap(map)</code><code>axis([0 1 0 1500]);</code><code>set(gcf,'Position<span style="color:#dd1144">',[20 100 320 250]);   </span></code><code>xlabel('Time (Sec)<span style="color:#dd1144">','</span>FontSize<span style="color:#dd1144">',12,'</span>FontName<span style="color:#dd1144">','</span>Times New Roman<span style="color:#dd1144">');</span></code><code>ylabel('Frequency (Hz)<span style="color:#dd1144">','</span>FontSize<span style="color:#dd1144">',12,'</span>FontName<span style="color:#dd1144">','</span>Times New Roman<span style="color:#dd1144">');</span></code><code>set(gca,'YDir<span style="color:#dd1144">','</span>normal<span style="color:#dd1144">')</span></code><code>set(gca,'FontSize<span style="color:#dd1144">',12);</span></code><code>set(gcf,'Color<span style="color:#dd1144">','</span>w<span style="color:#dd1144">');</span></code><code>​</code></span></span>

⛄ 运行结果

⛄ 参考文献

T. Li, Z. Peng, H. Xu and Q. He, "Parameterized domain mapping for order tracking of rotating machinery," in IEEE Transactions on Industrial Electronics, 2022, doi: 10.1109/TIE.2022.3201311

⛳️ 代码获取关注我

❤️部分理论引用网络文献，若有侵权联系博主删除

❤️ 关注我领取海量matlab电子书和数学建模资料

🍅 仿真咨询

1 各类智能优化算法改进及应用

生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化

2 机器学习和深度学习方面

卷积神经网络（CNN）、LSTM、支持向量机（SVM）、最小二乘支持向量机（LSSVM）、极限学习机（ELM）、核极限学习机（KELM）、BP、RBF、宽度学习、DBN、RF、RBF、DELM、XGBOOST、TCN实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断

2.图像处理方面

图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知

3 路径规划方面

旅行商问题（TSP）、车辆路径问题（VRP、MVRP、CVRP、VRPTW等）、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划、天线线性阵列分布优化、车间布局优化

4 无人机应用方面

无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配

、无人机安全通信轨迹在线优化

5 无线传感器定位及布局方面

传感器部署优化、通信协议优化、路由优化、目标定位优化、Dv-Hop定位优化、Leach协议优化、WSN覆盖优化、组播优化、RSSI定位优化

6 信号处理方面

信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号、信号配时优化

7 电力系统方面

微电网优化、无功优化、配电网重构、储能配置

8 元胞自动机方面

交通流 人群疏散 病毒扩散 晶体生长

9 雷达方面

卡尔曼滤波跟踪、航迹关联、航迹融合