資源簡介
用MATLAB繪制不同矩形窗長度的修正短時(shí)自相關(guān)函數(shù),包含M文件和語音文件。
代碼片段和文件信息
fid=fopen(‘test.txt‘‘rt‘)
b=fscanf(fid‘%f‘);
b1=b(1:640);
N=320;?????????????????????????????????????%窗長設(shè)定為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(‘延時(shí)?k‘)
ylabel(‘R(k)‘)
legend(‘N=320‘)
axis([0320-11]);
b2=b(1:320);
N=160;?????????????????????????????????????%窗長設(shè)定為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(‘延時(shí)?k‘)
ylabel(‘R(k)‘)
legend(‘N=160‘)
axis([0320-11]);
b3=b(1:140);???????????????????????????
N=70;??????????????????????????????????????%窗長設(shè)定為70
C=[];
for?k=1
?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????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
評(píng)論
共有 條評(píng)論