clear;
clc;

fp=24;%構(gòu)造函數(shù)
a=-10;
dt=0.002;
TT=0.5;
t=0:dt:TT;
s=exp（a*t）.*sin（2*pi*fp*t）;


nsmp=length（t）;
N=512;
for?k=nsmp+1:N
????s（k）=0.0;
end

xf=fft（sN）;

for?i=1:N
????if?i==1
????????xf（i）=xf（i）;
????elseif?i>1&&i????????xf（i）=2*xf（i）;
????else
????????xf（i）=0;
????end
end
i=1:N;

%?
%?x=ifft（xfN）;


%?xr=real（x）;
%?xi=imag（x）;
%?aa=abs（x）;

%?ts=i*dt;

%?p=-511:1:512;
%?[PI]=meshgrid（pi）;
%?ans=P.*（I.^2）;
%?XF=meshgrid（xf）;
%?G=XF.*exp（0.5j.*ans）;
%?g=ifft（GN）;
%?D=sqrt（1+P.^2）.*（abs（g））.^2;
%?
PP=500;
p=-PP:1:PP;
np=length（p）;
G=zeros（Nnp）;
Z=xf;
[PI]=meshgrid（pi）;
for?w=1:N
?????for?p=1:PP
?????????G（wPP+p）=Z（w）.*exp（0.5j.*（p）.*w^2）;
?????end
end
?for?w=1:N
?????for?p=1:PP+1
?????????G（wp）=Z（w）.*exp（0.5j.*（p-PP-1）.*w^2）;
?????end
?end
g=ifft（GN）;
M=zeros（12*PP+1）;
D=zeros（12*PP+1）;
for?p=1:PP
????M（p+PP）=sqrt（1+p.^2）;
end
for?p=1:PP+1
????M（p）=sqrt（1+（p-PP-1）.^2）;
end
for?w=1:N
?????for?p=1:PP
?????????D（wPP+p）=M（p）.*（abs（g（wPP+p）））.^2;
?????end
end
?for?w=1:N
?????for?p=1:PP+1
?????????D（wp）=M（p）.*（abs（g（wp）））.^2;
?????end
?end


%?D=M.*（abs（g））.^2;

%?angH=angle（x）;
%?dxdt_angH=diff（angH）/dt;
%?dxdt_angH（512）=0;
%?
%?figure
%?plot（ixf）;axis（[-100600-3030]）;

figure
plot（iabs（xf））;
axis（[0300050]）;
%?
%?title（‘頻域分布X（k）‘）;grid?on;
%?figure
%?plot（tsxr‘-‘tsxi‘:‘tsaa‘--‘）;axis（[00.5-11]）;title（‘解析函數(shù)及其包絡(luò)‘）;grid?on;
%?figure
%?plot（tsangH）;axis（[00.5-55]）;title（‘瞬時(shí)相位‘）;grid?on;
%?figure
%?plot（tsdxdt_angH）;axis（[00.5-4000800]）;title（‘瞬時(shí)頻率‘）;grid?on;

figure
surf（IPD）;colormap（jet）;shading?interp;



%?figure
%?subplot（311）;plot（tsxr‘-‘tsxi‘:‘tsaa‘--‘）;
%?axis（[01.1-11]）;
%?subplot（312）;plot（tsangH）;
%?%?axis（[01.1-55]）;
%?subplot（313）;plot（tsdxdt_angH）;
%?%?axis（[01.1-40001000]）;

%?for?k=1:nsmp
%?????bb（k）=b（k）;
%?????yr（k）=real（x（k））;
%?????yi（k）=imag（x（k））;
%?????tts（k）=k*dt;
%?????aaa（k）=abs（x（k））;
%?????aangH（k）=angH（k）;
%?????adxdt_angH（k）=dxdt_angH（k）;
%?end
%?figure
%?
%?subplot（311）;plot（ttsyr‘-‘ttsyi‘:‘ttsaaa‘--‘）;?axis（[00.6-11]）;
%?subplot（312）;plot（ttsaangH）;axis（[00.6-55]）;
%?subplot（313）;plot（ttsadxdt_angH）;axis（[00.6-4000500]）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????2355??2013-05-17?20:09??rwtransform.m