短时傅里叶变换 matlab程序
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
clear all
%窗口函数%
n1=40;
window=boxcar(n1);
w1=window;
figure(1);
stem(w1);
% 非平稳信号产生%
fs=1000;
a=0:1/fs:1;
f0=0;
f1=150;
y1=chirp(a,f0,1,f1);
x=y1(1:510);
figure(2);
plot(x);
% 短时傅里叶变换%
t=1:length(x);
n=length(x);
[tfr,t,f]=tfrstft(x,t,n,w1,0); contour(t,f,abs(tfr))
clc
clear all
%窗口函数%
n1=40;
window=boxcar(n1);
w1=window;
figure(1);
stem(w1);
% 非平稳信号产生%
fs=1000;
a=0:1/fs:1;
f0=0;
f1=150;
y1=chirp(a,f0,1,f1);
x=y1(1:510);
figure(2);
plot(x);
% 短时傅里叶变换%
t=1:length(x);
n=length(x);
[S,F,T,P]=spectrogram(x,w1,32,256); surf(T,F,10*log10(abs(P)));
axis tight;
view(0,90);