matlab 三维时频幅值图

1、内容简介

略

269
2、内容说明

时频分析（JTFA）即时频联合域分析，作为分析时变非平稳信号的有力工具，清楚地描述了信号频率随时间的变化关系。

时频图清楚地描述了信号频率随时间的变化关系，并且还通过颜色表达了对应的振幅大小。

（属于三维空间 H（w，t））

答1疑2咨3询4企q鹅q号：1762016542

3、仿真分析

clc
close all
clear
data1 = load(‘1.txt’);t1 = data1(:,1);x1 = data1(:,2);y1 = data1(:,3);z1 = data1(:,4);
data2 = load(‘2.txt’);t2 = data2(:,1);x2 = data2(:,2);y2 = data2(:,3);z2 = data2(:,4);

Fs = 8000;
figure
subplot(131)
plot(t1,x1)
title x
subplot(132)
plot(t1,y1)
title y
subplot(133)
plot(t1,z1)
title z
figure
subplot(131)
plot(t2,x2)
title x
subplot(132)
plot(t2,y2)
title y
subplot(133)
plot(t2,z2)
title z

[t_x_1,frequency_x_1,f_spectrum_x_1]=fft_s(x1,0.8,Fs);
figure,surf(t_x_1,frequency_x_1,f_spectrum_x_1); % 绘图
shading interp;
xlabel(‘时间/s’);ylabel(‘频率/Hz’);zlabel(‘幅值’);title 1信号x信号

[t_y_1,frequency_y_1,f_spectrum_y_1]=fft_s(y1,0.8,Fs);
figure,surf(t_y_1,frequency_y_1,f_spectrum_y_1); % 绘图
shading interp;
xlabel(‘时间/s’);ylabel(‘频率/Hz’);zlabel(‘幅值’);title 1信号y信号

[t_z_1,frequency_z_1,f_spectrum_z_1]=fft_s(z1,0.8,Fs);
figure,surf(t_z_1,frequency_z_1,f_spectrum_z_1); % 绘图
shading interp;
xlabel(‘时间/s’);ylabel(‘频率/Hz’);zlabel(‘幅值’);title 1信号z信号

[t_x_2,frequency_x_2,f_spectrum_x_2]=fft_s(x2,0.5,Fs);
figure,surf(t_x_2,frequency_x_2,f_spectrum_x_2); % 绘图
shading interp;
xlabel(‘时间/s’);ylabel(‘频率/Hz’);zlabel(‘幅值’);title 2信号x信号

[t_y_2,frequency_y_2,f_spectrum_y_2]=fft_s(y2,0.5,Fs);
figure,surf(t_y_2,frequency_y_2,f_spectrum_y_2); % 绘图
shading interp;
xlabel(‘时间/s’);ylabel(‘频率/Hz’);zlabel(‘幅值’);title 2信号y信号

[t_z_2,frequency_z_2,f_spectrum_z_2]=fft_s(z2,0.5,Fs);
figure,surf(t_z_2,frequency_z_2,f_spectrum_z_2); % 绘图
shading interp;
xlabel(‘时间/s’);ylabel(‘频率/Hz’);zlabel(‘幅值’);title 2信号z信号
 

4、参考论文

略