%?Main?Function
%%
clear;
clc;
%%
%?畫出原始信號;
r?=?1;?L=512;n?=1:r:L;
x?=?5*exp（j*（0.15*n.^2））+6*exp（j*（300*n-0.15*n.^2））;
figure（1）;plot（nx）;
xlabel（‘n‘）;?ylabel（‘x‘）;xlim（[0512]）;
title（‘原始信號X‘）;
%%
%?創(chuàng)建四個不同寬度的高斯窗口默認singma=2.5
N1=7;?w1=gausswin（N1）;?????????
N2=15;w2=gausswin（N2）;????????
N3=21;w3=gausswin（N3）;????????
N4=31;w4=gausswin（N4）;?????????
%%
%?畫出不同寬度的高斯窗口函數(shù).
figure（2）
subplot（221）;plot（1:N1w1）;title（[‘寬度N1=‘?num2str（N1）]）?;xlim（[032]）;
subplot（222）;plot（1:N2w2）;title（[‘寬度N2=‘?num2str（N2）]）?;xlim（[032]）;
subplot（223）;plot（1:N3w3）;title（[‘寬度N3=‘?num2str（N3）]）?;xlim（[032]）;
subplot（224）;plot（1:N4w4）;title（[‘寬度N4=‘?num2str（N4）]）?;xlim（[032]）;
%%
%?用不同的窗口函數(shù)對x進行短時傅里葉變換
z=x‘;
[tfr1t1f1]=tfrstft（z1:LLw1）;?
[tfr2t2f2]=tfrstft（z1:LLw2）;?
[tfr3t3f3]=tfrstft（z1:LLw3）;?
[tfr4t4f4]=tfrstft（z1:LLw4）;?
%%
%?%?畫出不同寬度窗口函數(shù)的時頻分析2維與3維圖.
figure（3）
subplot（221）;imagesc（abs（tfr1））;?xlabel（‘t‘）;ylabel（‘?f‘）;xlim（[0512]）;title（[‘寬度N1=‘?num2str（N1）]）;?
subplot（222）;imagesc（abs（tfr2））;?xlabel（‘t‘）;ylabel（‘?f‘）;xlim（[0512]）;title（[‘寬度N2=‘?num2str（N2）]）;?
subplot（223）;imagesc（abs（tfr3））;?xlabel（‘t‘）;ylabel（‘?f‘）;xlim（[0512]）;title（[‘寬度N3=‘?num2str（N3）]）;?
subplot（224）;imagesc（abs（tfr4））;?xlabel（‘t‘）;ylabel（‘?f‘）;xlim（[0512]）;title（[‘寬度N4=‘?num2str（N4）]）;
figure（4）
subplot（221）;mesh（t1f1abs（tfr1））;xlabel（‘t‘）;ylabel（‘?f‘）;zlabel（‘?A‘）;xlim（[0512]）;title（[‘寬度N1=‘?num2str（N1）]）;?
subplot（222）;mesh（t2f2abs（tfr2））;xlabel（‘t‘）;ylabel（‘?f‘）;zlabel（‘?A‘）;xlim（[0512]）;title（[‘寬度N2=‘?num2str（N2）]）;?
subplot（223）;mesh（t3f3abs（tfr3））;xlabel（‘t‘）;ylabel（‘?f‘）;zlabel（‘?A‘）;xlim（[0512]）;title（[‘寬度N3=‘?num2str（N3）]）;?
subplot（224）;mesh（t4f4abs（tfr4））;xlabel（‘t‘）;ylabel（‘?f‘）;zlabel（‘?A‘）;xlim（[0512]）;title（[‘寬度N4=‘?num2str（N4）]）;
%%
%?進行短時傅里葉反變換
s1=tfristft（tfr11:Lw1）;
s2=tfristft（tfr21:Lw2）;
s3=tfristft（tfr31:Lw3）;
s4=tfristft（tfr41:Lw4）;
%%
%?畫出重構(gòu)信號與原始信號的差圖
figure（5）
subplot（221）;plot（nabs（s1-x‘））;?xlabel（‘n‘）;ylabel（‘△‘）;xlim（[0512]）;title（[‘寬度N1=‘?num2str（N1）]）;
subplot（222）;plot（nabs（s2-x‘））;?xlabel（‘n‘）;ylabel（‘△‘）;xlim（[0512]）;title（[‘寬度N2=‘?num2str（N2）]）;
subplot（223）;plot（nabs（s3-x‘））;?xlabel（‘n‘）;ylabel（‘△‘）;xlim（[0512]）;title（[‘寬度N3=‘?num2str（N3）]）;
subplot（224）;plot（nabs（s4-x‘））;?xlabel（‘n‘）;ylabel（‘△‘）;xlim（[0512]）;title（[‘寬度N4=‘?num2str（N4）]）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

????.CA....??????2399??2012-12-13?19:28??tfristft.m

????.CA....??????2572??2012-12-13?19:28??tfrstft.m

????.CA....??????2556??2012-12-13?19:28??my_tfrstft2.m

-----------?---------??----------?-----??----

?????????????????7527????????????????????3