不同矩形窗长度的修正短时自相关函数MATLAB仿真 语音信号处理

用MATLAB绘制不同矩形窗长度的修正短时自相关函数，包含M文件和语音文件。

fid=fopen（‘test.txt‘‘rt‘）
b=fscanf（fid‘%f‘）;
b1=b（1:640）;
N=320;                                     %窗长设定为320
A=[];
for k=1:320;
sum=0;
for m=1:N;
sum=sum+b1（m）*b1（m+k-1）;
end
A（k）=sum;
end
for k=1:320
A1（k）=A（k）/A（1）;                           %归一化
end
figure（1）
subplot（311）
plot（A1）;
xlabel（‘延时 k‘）
ylabel（‘R（k）‘）
legend（‘N=320‘）
axis（[0320-11]）;

b2=b（1:320）;
N=160;                                     %窗长设定为160
B=[];
for k=1:160;
sum=0;
for m=1:N;
sum=sum+b2（m）*b2（m+k-1）;
end
B（k）=sum;
end
for k=1:160
B1（k）=B（k）/B（1）;%归一化B（k）
end
figure（1）
subplot（312）
plot（B1）;
xlabel（‘延时 k‘）
ylabel（‘R（k）‘）
legend（‘N=160‘）
axis（[0320-11]）;

b3=b（1:140）;                           
N=70;                                      %窗长设定为70
C=[];
for k=1

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2014-06-30 16:25  No.1\
     文件        1080  2014-06-30 16:42  No.1\pro_1.m
     文件        9033  2007-11-06 16:35  No.1\test.txt