分数阶傅立叶变换中午matlab,怎么做短时分数阶傅里叶变换
我的分数阶傅里叶变换如下,怎么将它和短时傅里叶变换结合起来,变成短时分数阶傅里叶变换,有没有大佬指点一下下!!!
frft:
function Faf = frft(f, a)
% The fast Fractional Fourier Transform
% input: f = samples of the signal
% a = fractional power
% output: Faf = fast Fractional Fourier transform
error(nargchk(2, 2, nargin));
f = f(:);
N = length(f);
shft = rem((0:N-1)+fix(N/2),N)+1;
sN = sqrt(N);
a = mod(a,4);
% do special cases
if (a==0), Faf = f; return; end;
if (a==2), Faf = flipud(f); return; end;
if (a==1), Faf(shft,1) = fft(f(shft))/sN; return; end
if (a==3), Faf(shft,1) = ifft(f(shft))*sN; return; end
% reduce to interval 0.5 < a < 1.5
if (a>2.0), a = a-2; f = flipud(f); end
if (a>1.5), a = a-1; f(shft,1) = fft(f(shft))/sN; end
if (a<0.5), a = a+1; f(shft,1) = ifft(f(shft))*sN; end
% the general case for 0.5 < a < 1.5
alpha = a*pi/2;
tana2 = tan(alpha/2);
sina = sin(alpha);
f = [zeros(N-1,1) ; interp(f) ; zeros(N-1,1)];
% chirp premultiplication
chrp = exp(-i*pi/N*tana2/4*(-2*N+2:2*N-2)'.^2);
f = chrp.*f;
% chirp convolution
c = pi/N/sina/4;
Faf = fconv(exp(i*c*(-(4*N-4):4*N-4)'.^2),f);
Faf = Faf(4*N-3:8*N-7)*sqrt(c/pi);
% chirp post multiplication
Faf = chrp.*Faf;
% normalizing constant
Faf = exp(-i*(1-a)*pi/4)*Faf(N:2:end-N+1);
%%%%%%%%%%%%%%%%%%%%%%%%%
function xint=interp(x)
% sinc interpolation
N = length(x);
y = zeros(2*N-1,1);
y(1:2:2*N-1) = x;
xint = fconv(y(1:2*N-1), sinc([-(2*N-3):(2*N-3)]'/2));
xint = xint(2*N-2:end-2*N+3);
%%%%%%%%%%%%%%%%%%%%%%%%%
function z = fconv(x,y)
% convolution by fft
N = length([x(:);y(:)])-1;
P = 2^nextpow2(N);
z = ifft( fft(x,P) .* fft(y,P));
z = z(1:N);
分数阶傅立叶变换中午matlab,怎么做短时分数阶傅里叶变换相关推荐
- 分数阶傅立叶变换 matlab,【综述】分数阶傅里叶变换(FRFT)
原标题:[综述]分数阶傅里叶变换(FRFT) 作者:WTT整理 傅里叶级数(傅里叶变换)几乎在所有科学和工程领域发挥着重要作用.黎曼积分和勒贝格积分均起源于对傅里叶级数的研究,傅里叶级数(傅里叶变换) ...
- 分数阶傅立叶变换程序汇总(转载)
转载链接http://forum.vibunion.com/forum.php?mod=viewthread&action=printable&tid=41561 声振论坛 标题: 分 ...
- 时域脉冲傅立叶变换(MATLAB)
采用MATLAB编程,对太赫兹脉冲波的时域波形进行傅立叶变换,画出幅度谱线和相位谱线.在网上没能找到相关代码的教程,只能自己根据师兄的指导来做一遍了. 1. MATLAB代码 clear variab ...
- matlab实现傅里叶变换_傅立叶变换求解偏微分方程和积分方程
本编文章探讨一下数学物理方法中一个常见的问题,即如何用傅立叶变换得到PDE或者积分方程的解.在文章的最后,会附上MATLAB的实现代码.学疏才浅,欢迎大家指点! 1. 什么是傅立叶变换? 从数学形式上 ...
- 全面解析傅立叶变换(非常详细)
前言 第一部分. DFT 第一章.傅立叶变换的由来 第二章.实数形式离散傅立叶变换(Real DFT) 从头到尾彻底理解傅里叶变换算法.下 第三章.复数 第四章.复数形式离散傅立叶变换 前言: &q ...
- # 傅立叶变换 画曲线_让你永远忘不了的傅立叶变换解析
使用联想链条和几何直观,辅以从实际需求衍生概念的思考模式,详解什么是傅立叶变换,为什么要做傅立叶变换等,帮助记忆和理解,目的当然是标题所说:让你永远忘不了傅立叶变换这个公式.另,这篇博客还从侧面一定程 ...
- 傅立叶变换,时域,频域一
参考文献: 信号完整性分析 "信息传输调制和噪声"P31, "傅立叶变换的数学再认识"及若干网上博客. 目录 信号分析方法概述 时域 频域 ...
- 二维小波变换_让你永远忘不了的傅立叶变换解析
使用联想链条和几何直观,辅以从实际需求衍生概念的思考模式,详解什么是傅立叶变换,为什么要做傅立叶变换等,帮助记忆和理解,目的当然是标题所说:让你永远忘不了傅立叶变换这个公式.另,这篇博客还从侧面一定程 ...
- 线性代数28——复矩阵和快速傅立叶变换
原文 | https://mp.weixin.qq.com/s/YzPoPnRb-gEm_EiV9et0TA 实矩阵也可能碰到复特征值,因此无可避免地在矩阵运算中碰到复数. 矩阵当然也有可能包含复数, ...
最新文章
- “三高”程序员谈:“三高”集群架构 哈哈哈!!!
- Weak Pair HDU - 5877 树状数组+离散化+DFS遍历
- 1-10:学习shell之进程
- 嗖嗖移动业务大厅代码_移动云应用性能监控,掌控云时代的多变应用
- iOS UISegmentedControl
- 程序员入门:三本必看的书
- 用nmap扫描内网conficker
- H5+CSS3移动商城界面.七天从零实战课程-会员_登录
- CSS3 页面布局
- 从思维走向实践,数字化转型 IT 经营的成功路径
- win10开发环境搭建之wsl2(Ubuntu)+Terminal+docker
- [转贴]馒头的N种吃法
- 一站式数字藏品收款系统开发解决方案
- opencv-python给图片添加半透明mask
- 炒股:如何安全地追热点
- 线性插值函数的基函数构造
- css-背景渐变色动画
- 关于.SFV文件的知识
- C语言数据结构之树超全详解
- iOS开发中,PNG图片的各种压缩软件,以及压缩率对比
热门文章
- ei会议论文录用但不参加会议_科学网—推荐一个录用论文全部发表在EI期刊上的会议 - 徐庆征的博文...
- scuba 报表_是否想了解JavaScript的for循环? 这个动画的SCUBA潜水员可以提供帮助!...
- 最简单DIY基于ESP32CAM的物联网相机系统⑥(用上位机VS2013 MFC实现WIFI图传)
- 箭头函数与普通函数,以及使用场景
- 计算机图标右击管理打不开,电脑图标打不开怎么办6
- python命名规则数字开头的成语_浅谈Python中带_的变量或函数命名
- Linux 中 TCL 和 Expect语法
- oracle12免安装版配置,配置免安装版Oracle客户端
- 大火的扩散模型终于杀到了NLP领域
- 使用sunshine+moonlight 实现电脑串流到电视(Android 设备)低延迟投屏