matlab的dc程序,DC带你玩Matlab——我不是教你诈

上午看到&＃64;清华南都 童鞋发了这么一个微博&＃xff1a;

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

主料&＃xff1a;360vs百度的视频

辅料&＃xff1a;IDM&＃xff0c;格式工厂&＃xff0c;Wavosaur&＃xff0c;Matlab

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

以下是本人的拙劣代码&＃xff1a;

clc,clear,close all

%读取波形

[x1, sampFreq, nBits] &＃61; wavread(&＃39;361.wav&＃39;);

%取双声道的第一个声道

x1&＃61;x1(:,1);

%作时域波形图

figure

plot(x1)

grid on

x&＃61;x1(256000:266000);%以最后一段(号码为8)为例

figure

plot(x)

grid on

N &＃61;sampFreq; w1 &＃61; 2*pi/N; k &＃61; 0:N-1;

X &＃61; fft(x,N);%对信号进行FFT

magX&＃61;abs(X(1:N/2&＃43;1));

figure

plot(magX); %作幅度谱

title(&＃39;FFT&＃39;);

grid on

xlabel(&＃39;k&＃39;);

ylabel(&＃39;abs(X)&＃39;);

axis();

最后附上几张图&＃xff0c;周总的号码我就不公布了。

参考资料&＃xff1a;DTMF