test.asv

很多matlab的源代码

[y,fs,bits]=wavread('Laser');
x=y(:,1);
subplot(4,2,1);
plot(x);
N1=2000;
N2=2*N1;
d=[1 zeros(1,N1-1) 0.7 zeros(1,N1-1) 0.5];
n=[1];
y1=conv(x,d);
subplot(4,2,2);
plot(y1);
y2=filter(n,d,y1);
sound(y2,fs);
subplot(4,2,3);
plot(y2);
subplot(4,2,5);
nn=0:length(d)-1;
h=stem(nn,d);
set(h,'Linewidth',2);
set(h,'MarkerSize',1);
xx=[1 zeros(1,20000)];
h2=filter(n,d,xx);
nn=0:length(h2)-1;
subplot(4,2,6);
h_2=stem(nn,h2);
set(h_2,'Linewidth',2);
set(h_2,'Linewidth',1);

PW1=abs(fft(d,2048));
N=length(PW1);
HN=fix(N/2);
PW1=PW1(1:HN);
nn=0:HN-1;
w=nn*pi/HN;
subplot(4,2,7);
plot(w,PW1);

PW2=abs(freqz(n,d,1024));
HN=length(PW2);
PW2=PW2(1:HN);
nn=0:HN-1;
w=nn*pi/HN;
subplot(4,2,8);
plot(w,PW2);
figure
P1=PW1(1:200);
P2=PW2(1:200);
Y=P1.*P2';
plot(Y);