窗口傅里叶正反变换matlab程序

窗口傅里叶正反变换matlab程序，同时包含使用的例子程序

%Program for plotting spectrum and error of windowed Fourier transform
clear;clc;

%Add path
addpath（[cd‘\lib‘]）;

%Define signal 
L=2^10;
dt=1/L;
t=（0:L-1）*dt;
f0=100; %Constrat sine signal
s=sin（2*pi*f0*t.^3）; 
% f1=100;f2=200;f3=300;f4=400;%Construnct four frequencies signal
% jump1=fix（L/4）;jump2=fix（2*L/4）;jump3=fix（3*L/4）;
% s=[sin（2*pi*f1*t（1:jump1））sin（2*pi*f2*t（jump1+1:jump2））sin（2*pi*f3*t（jump2+1:jump3））sin（2*pi*f4*t（jump3+1:end））];
% f0=100; %Construct Ricker signal
% t0=fix（L/2）*dt;
% s=（1-2*（pi*f0*（t-t0））.^2）.*exp（-（pi*f0*（t-t0））.^2）;
% s=zeros（1L）; %Construct Delta signal
% s（fix（L/2））=1; 

%Define window function
epsilon=1e-3;
alpha=2*log（1/epsilon）/（（fix（sqrt（L））*dt）^2）; %Assume a point radius of fix（sqrt（L））
w=@（x）exp（-alpha/2*x.^2）;  %构造高斯窗
S=wft（swdt）;   %进行窗口傅里叶变换

%Plot reconstruction error
subplot（212）;
s0=iwft（Swdt）;
xIndex=（0:L-1）*dt;
plot（xIndexs-s0）;
axis tight;
xlabel（‘t（s）‘）;ylabel（‘error‘）;
title（‘重构误差图‘）;

%Plot spectrum
subplot（211）;
S=S（1:fix（L/2）+1:）;
[mn]=size（S）;
x=（1:n）*dt;
df=1/（L*dt）;
y=（1:m）*df;
%mesh（xyabs（S）.^2）;
pcolor（xyabs（S）.^2）;
shading interp;
axis tight;
xlabel（‘t（s）‘）;
ylabel（‘f（Hz）‘）;
title（‘谱图‘）;
colormap（‘jet‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2018-12-10 17:24  窗口傅里叶变换\
     目录           0  2018-12-10 17:24  窗口傅里叶变换\lib\
     文件         580  2018-12-11 16:58  窗口傅里叶变换\lib\iwft.m
     文件         610  2018-12-11 16:58  窗口傅里叶变换\lib\wft.m
     文件        1268  2018-12-14 11:27  窗口傅里叶变换\plot_spectrum_and_error.m