DC带你玩Matlab——我不是教你诈

上午看到@清华南都 童鞋发了这么一个微博：

http://www.vjianke.com/board/b5e51221f40e42489627af17814aae89.clip#ZQ0YL 

下午DC准备用Matlab小试牛刀:

主料：360vs百度的视频

辅料：IDM，格式工厂，Wavosaur，Matlab

做法：先IDM抓取视频（flv格式），再用格式工厂转成wav格式，然后用Wavosaur剪辑关键帧，最后使用Matlab进行音频信号处理（画时域和频域图即可）。

以下是本人的拙劣代码：

clc,clear,close all
%读取波形
[x1, sampFreq, nBits] = wavread('361.wav');
%取双声道的第一个声道
x1=x1(:,1);                  
%作时域波形图
figure
plot(x1)
grid on
x=x1(256000:266000);%以最后一段（号码为8）为例
 
figure
plot(x)
grid on
 
N =sampFreq; w1 = 2*pi/N; k = 0:N-1; 
X = fft(x,N);%对信号进行FFT 
magX=abs(X(1:N/2+1)); 
figure                             
plot(magX); %作幅度谱
title('FFT');
grid on
xlabel('k');
ylabel('abs(X)');
axis();

最后附上几张图，周总的号码我就不公布了。

参考资料：DTMF