matlab 注解一下这个程序里的几个地方

来源：百度知道 编辑：UC知道 时间：2024/06/26 01:27:31
x=wavread('F:\a.wav',1024);%读取该音频文件前1024个采样点
plot(x);

y=wavread('F:\a.wav',1000);
a1=hamming(1000);设计长度1000的窗函数
a2=fft(y.*a1);％对语音信号滤波
a3=real(a2);％取滤波后信号的实部
plot(a3);

z=wavread('F:\a.wav',[1,500]);
han=hamming(100);
En=conv(han,abs(z));％通过卷积实现滤波
plot(En);

w=wavread('F:\a.wav',1024);
h=hamming(1000);
En=conv(h,w.*w);
plot(En);

k=wavread('F:\a.wav',[1000,2000]);
han=hamming(1000);
for i=1:1000
if k(i)>=0
b(i)=1;
else
b(i)=-1;
end
end
for i=2:1000
x(i)=b(i)-b(i-1);
end
z=conv(abs(x),han);
Zn=z./2000;
plot(Zn);

x=wavread('F:\a.wav',1000);
a1=xcorr(x);
plot(a1);

yc(1) = 0;
step = 0.0001;
x=wavread('F:\a.wav',1000);
a1=xcorr(x);
plot(a1);
for n=1:1000
d(n

close all,clear all,clc
figure
x=wavread('a.wav');%读取该音频文件前1024个采样点
plot(x);

figure
y=wavread('a.wav',1000);
a1=hamming(1000);%设计长度1000的窗函数
a2=fft(y.*a1);%对语音信号加窗后求频谱
a3=real(a2);%取频谱的实部
plot(a3);
figure

z=wavread('a.wav',[1,500]);
han=hamming(100);
En=conv(han,abs(z));%通过卷积实现低通滤波
plot(En);

w=wavread('a.wav',1024);
h=hamming(1000);
En=conv(h,w.*w); %通过卷积实现w.^2的低通滤波
plot(En);

k=wavread('a.wav',[1000,2000]);
han=hamming(1000);
for i=1:1000
if k(i)>=0
b(i)=1;
else
b(i)=-1;
end
end
for i=2:1000
x(i)=b(i)-b(i-1); %x相当于对k做了一个简单的差分编码
end
z=conv(abs(x),han); %z是abs(x)通过低通滤波器后的样子
Zn=z./2000;
plot(Zn);

x=wavread('a.wav',1000);
a1=xcorr(x); %对信号x作自相关，能用来求语音的频率
plot(a1);

yc(1) = 0; %自适应增量调制编码
s