首页 > 学术百科

分数阶傅立叶变换中午matlab,怎么做短时分数阶傅里叶变换

分数阶傅⽴叶变换中午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;

rqsN = sqrt(N);

a = mod(a,4);

% do special cases

if (a==0), Faf = f; return; end;

if (a==2), Faf = flipud(f); return; end;济南817

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);

本文发布于:2024-09-23 14:33:14，感谢您对本站的认可！

本文链接：https://www.17tex.com/xueshu/436968.html

上一篇：傅里叶变换分析信号的缺点

下一篇：傅里叶变换本质及其公式解析

标签：变换分数阶傅燃烧

留言与评论（共有 0 条评论）