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

发布时间:2025-03-09 15:38

学习如何进行简单的财务报表分析,理解自己的钱去了哪里。 #生活常识# #金融理财常识#

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

短时傅里叶变换进行时频分析
在这里插入图片描述

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)

1234567891011121314151617181920

在这里插入图片描述
在这里插入图片描述

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)

123456789101112131415161718192021 扫频信号谱阵分析
在这里插入图片描述

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)

123456789101112131415161718192021 扫频信号生成 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) 1234567 频率混叠
在这里插入图片描述

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

网址:matlab简单分析短时傅里叶变换STFT https://www.yuejiaxmz.com/news/view/815857

相关内容

python与matlab代码等价切换
matlab中for循环的简单使用
MATLAB图像处理(包括图像类型转换)
相量与正弦量的数学变换原理
PCA (主成分分析)详解 (写给初学者) 结合matlab
频谱函数X(jω)=cos(ω)的傅里叶逆变换x(t)为( )。
Matlab 透视变换原理 代码实现
智能客服的语音识别技术实践
在智能微电网中,如何应用贝叶斯博弈理论来优化需求侧管理,并通过MATLAB进行模拟验证?
在便携式电子设备中,电荷泵如何优化DC/DC电压转换以提高效率?它与传统DC/DC变换器相比有哪些优势?

随便看看