%本程序意在实现四项单参数WFRFT的基本功能,仅适用于单参数WFRFT
%参考公式为王晓鲁硕士论文公式(2-3)和(2-5)
function y=spwfrft(data,a)
%function y=spwfrft(sr,si,a)
%%
%信号读取部分
% L=length(si); %测试信号的长度
% sr_in=sr;
% si_in=si; %读入实部和虚部
%%
%四项分数阶傅立叶变换
% x0=sr_in+1i*si_in; %组成原始信号
x0=data;
x1= 1/sqrt(length(x0))*fft(x0);%1阶DFT
x2= 1/sqrt(length(x1))*fft(x1);%2阶DFT
x3= 1/sqrt(length(x2))*fft(x2);%3阶DFT
%%
%
wl=zeros(1,4); %预留加权系数存储空间
for j=1:1:4
wl(j)=cos(pi*(a-j+1)/4)*cos(2*pi*(a-j+1)/4)*exp(-1i*3*pi*(a-j+1)/4);%加权系数的组成
end;
y=wl(1)*x0+wl(2)*x1+wl(3)*x2+wl(4)*x3; %组成最终的输出信号