myffttotestpcp.m

来自「自己定义了一种有关声音的特征」· M 代码 · 共 23 行

function myfftToTestpcp()
N=512*16;
fs=8000*2;
fc1=261;
s=30;
figure(1);
t = 1/fs:1/fs:s;
x=zeros(size(t));
for i=0:12
  x=x+sin(2*pi*fc1*2^(i/12)*t);   
end
y = x + 2*randn(size(t));
plot(fs*t(1:100),y(1:100));
title('Signal Corrupted with Zero-Mean Random Noise');
xlabel('time (milliseconds)');
%y=y(1:N*floor(s*fs/N));
%y=reshape(y,N,floor(s*fs/N));
Y=enframe(y,N,N/2).';
flag=1;
[pcpdata,lodd,leven]=pcp(Y,flag,fs);
fprintf('lodd=%d  leven=%d\n',lodd,leven);
%Pyy = abs(Y).^2 / 512;