MATLAB中的信号短时傅里叶变换获得信号的时频联合分析

短时傅里叶变换 -> 时频联合分析

一、双音多频DTMF知识储备

对于一个声音，如果能够检测出其中高频和低频的频率，就可以指出对应的按键编码。

现在可以通过专用的数字处理芯片直接获得双音频解码，也可以通过计算软件进行解码，后者就是这篇博文将要演示的。

二、简单的时频联合分析

最简单的时频联合分析，可通过录制5s的声音，然后可直观看到时域某个时刻对应的频率。

clc;
clear;

recObj = audiorecorder;
recordblocking(recObj, 5);%录音5s
disp('Press any key to continue playing.');
pause;
play(recObj);
disp('Press any key to continue plot.');
pause;
myRecording = getaudiodata(recObj);plot(myRecording);
disp('Press any key to continue fouier.');
pause;
spectrogram(myRecording,2048,1024,1:2000,8000,'yaxis');

  

 

  
1
  
2
  
3
  
4
  
5
  
6
  
7
  
8
  
9
  
10
  
11
  
12
  
13
  
14
 

二、图片在时频联合分析中的显示

一个做好的图片声音，进行时频联合分析…

clc;
clear;

[s fs] = audioread('G:\getvoice.mp3');
% spectrogram(s(:,1), 2048, 1024, 1:2000, fs, 'yaxis');
spectrogram(s(:,1), 8192, 4096, 1:2000, fs, 'xaxis');
% plot((0:length(s)-1)/fs, s);

  

 

  
1
  
2
  
3
  
4
  
5
  
6
  
7
 

挺有意思一个东西，不是吗。。

三、制作在时频联合分析中显示的图片的音乐（图片变成音乐）

图片转音乐的matlab程序，暂时还没有写出来，只能辅助软件了。<链接> 提取码：4g6y

注意，
对应的图片分辨率最好设置为

以及要转为256色位图(.bmp格式)

clc;
clear;

[s,fs] = audioread('3.wav');
h1 = figure;
spectrogram(s(:,1), 8192, 4096, 1:15000, fs, 'yaxis');
saveas(h1,'12','jpg');
I = imread('12.jpg');
J = imadjust(I,[0.3 0.7],[]);
imshow(J);

  

 

  
1
  
2
  
3
  
4
  
5
  
6
  
7
  
8
  
9
  
10
 

最终图片效果…

文章来源: recclay.blog.csdn.net，作者：ReCclay，版权归原作者所有，如需转载，请联系作者。

原文链接：recclay.blog.csdn.net/article/details/84444910