matlab 三维时频图

1、内容简介

略
384-可以交流、咨询、答疑

2、内容说明

略

3、仿真分析

clc
close all
clear
data = load('data.TXT');
t1 = data(:,1);x1 = data(:,2);
Fs = 500e3;% 采样频率500KHZ
figure
plot(t1,x1)
title 冲击波压力
%% 
close all
[t_x_1,frequency_x_1,f_spectrum_x_1]=fft_s(x1,0.01,Fs);
figure,surf(t_x_1,frequency_x_1,f_spectrum_x_1); % 绘图
shading interp;
xlabel('时间/s');ylabel('频率/Hz');zlabel('幅值');title 冲击波压力
[X,Y]=meshgrid(t_x_1,frequency_x_1);
figure
pcolor(X,Y,f_spectrum_x_1)
colorbar;
shading interp;
ylim([0 5000]) % 显示的频率范围
xlabel('时间/s')
ylabel('频率/Hz')
%% 
[f,y0]=change_fft(Fs,x1);
figure
plot(f,y0)
xlabel('频率/Hz')
ylabel('幅值')
fre = [20 100 250 500 1000 2000 5000 10000 20000 40000]; % 频率范围
power = zeros(1,length(fre)+1);
for i = 1:length(fre)
    power(i) = sum( y0 (f<=fre(i)))/sum(y0);
end
power(end) = 1;
power_diff =[power(1) diff(power(1:end))];
figure
bar(power_diff)
xlabel 频率段
ylabel 归一化能量占比
figure
bar(power)
xlabel 频率段
ylabel 归一化能量占比

4、参考论文

略