matlab简单分析短时傅里叶变换STFT

华中科技大学《数字信号分析理论实践》学习总结记录

短时傅里叶变换进行时频分析

Fs = 5120.0;
dt = 1/Fs;
N = 32768;
df = 500.0/N;
T = dt * N;
t = linspace(0,T,N);
for i = 1:8192x(i) = sin(2*pi*200*i*dt);x(i+8192) = sin(2*pi*400*i*dt);x(i+16384) = sin(2*pi*600*i*dt);x(i+24576) = sin(2*pi*800*i*dt);
end
Z = spectrogram(x,1024,512);
P = 20*log10(sqrt(Z.*conj(Z)));
X = linspace(0,Fs/2,513);
Y = linspace(0,T,63);
mesh(X,Y,P');
xlabel('频率')
ylabel('时间')
view(15,70)

调幅信号谱阵分析

Fs = 5120.0;
dt = 1/Fs;
N = 32768;
df = 500.0/N;
T = dt * N;
t = linspace(0,T,N);
for i = 1:Nx(i) = sin(2*pi*0.3*i*dt)*sin(2*pi*(500+df*i)*i*dt);
end
plot(x)
figure
Z = spectrogram(x,1024,512);
P = 20*log10(sqrt(Z.*conj(Z)));
S = size(P);
% 513 = 1024/2+1
X = linspace(0,Fs/2,513);
Y = linspace(0,T,S(2));
mesh(X,Y,P');
xlabel('频率')
ylabel('时间')
view(15,70)

扫频信号谱阵分析

Fs = 5120.0;
dt = 1/Fs;
N = 32768;
df = 500.0/N;
T = dt * N;
t = linspace(0,T,N);
for i = 1:Nx(i) = sin(2*pi*(500+df*i)*i*dt);
end
plot(x)
figure
Z = spectrogram(x,1024,512);
P = 20*log10(sqrt(Z.*conj(Z)));
S = size(P);
% 513 = 1024/2+1
X = linspace(0,Fs/2,513);
Y = linspace(0,T,S(2));
mesh(X,Y,P');
xlabel('频率')
ylabel('时间')
view(15,70)

扫频信号生成 chirp()

Fs = 1000;
T = 0:0.001:12;
X = chirp(T,1,12,400,'q');
subplot(211)
plot(T,X);
subplot(212)
spectrogram(X,512,256,512,Fs)

频率混叠

Fs = 1000;
T = 0:0.001:12;
X = chirp(T,1,12,900,'q');
spectrogram(X,512,256,512,Fs)

matlab简单分析短时傅里叶变换STFT相关推荐

短时傅里叶变换(STFT)及matlab
笔记~自用版~ 短时傅里叶变换的基础理论短时傅里叶变换(Short-Time Fourier Transform, STFT)是一种时频分析方法,它将信号在时间域上分成若干个短时段,对每个短时段进行 ...
matlab简单分析数字滤波器FIR
时域FIR滤波器华中科技大学<数字信号分析理论实践>第六单元 FIR滤波器学习总结记录 FIR 滤波器脉冲响应函数设计法 MATLAB 中 FIR 滤波器设计函数 fir2() f = ...
matlab简单分析频域滤波和时域滤波
华中科技大学<数字信号分析理论实践>第六单元信号数字滤波的概念学习总结记录频域滤波对信号进行低通滤波的一种快速方法,低频信号幅度设为1,负频率镜像过去,把虚部设为0 Fs = 20 ...
matlab简单分析矩阵乘法——一阶行向量乘一阶列向量
A=[1,2,3],B=[321]A=\begin{bmatrix} 1,&2,&3 \end{bmatrix},B=\begin{bmatrix} 3\\ 2\\ 1 \end{bm ...
matlab简单分析离散小波变换
华中科技大学<数字信号分析理论实践> 学习总结记录离散小波变换进行多尺度分解离散二进小波变化进行小波分解 h1 = [-0.014434 0.050359 0.017710 -0.09 ...
matlab简单分析其他类型滤波器（陷波尖峰梳状半带希尔伯特）
其他类型滤波器华中科技大学<数字信号分析理论实践>滤波器学习总结记录陷波滤波器 iirnotch Fs = 200; dt = 1.0/Fs; T = 2; N = T/dt; t ...
matlab简单分析模拟滤波器 IIR
IIR 滤波器华中科技大学<数字信号分析理论实践>第六单元 IIR滤波器学习总结记录 IIR滤波器一般先借助模拟滤波器进行设计,然后再将模拟滤波器转为数字滤波器脉冲响应不变法 &am ...
2021-05-10 Matlab短时傅里叶变换和小波变换的时频分析
Matlab短时傅里叶变换和小波变换的时频分析简介本文主要给定一小段音频,通过短时傅里叶变换和小波变换制作时频图.音频的采样率为44100, 短时傅里叶变换在matlab中,短时傅里叶变换的分析 ...
librosa.stft() 短时傅里叶变换
librosa 短时傅里叶变换 import numpy as np # pip install numpy import librosa # pip install librosa y, sr = ...

matlab简单分析短时傅里叶变换STFT

matlab简单分析短时傅里叶变换STFT相关推荐

最新文章

热门文章