1、% % wavrecord。ey 垃 e? o 送卩？% Fs = 11025;% y =wavrecord(5*Fs,Fs,int16;%?5? mic euo?% wavplay(y,Fs; %2?3? dfs = 11025 ; % e?2ee? u? ? ey =wavrecord(1*fs ,fs ,uint8 ; %?ae?2 e? 8?6?dlD?o?, e? ?5ssubplot(411;plot (y ;wavplay(y ,fs ; % ? ? u2e?-戰(zhàn) 6?0?Y1 =fft (y ;% ?D探o? /?fft ?Y2 =fftshift(Y1;f=(0:le ngt

2、h(Y2-1/le ngth(Y2*fs-fs/2 ;subplot(412plot(f,abs(Y1;subplot(413plot(abs(Y2 subplot(414 N=le ngth(Y2; plot(-N/2:N/2-1*fs/N,abs(Y2*2/N; figure(2 fs=1000;f=100;t=0:1/fs:1;y=s in (2*pi*f*t;subplot(211plot(t,y;Y1=fft(y;Y2 =fftshift (Y1;f=(0:le ngth(Y1-1/le ngth(Y1*fs-fs/2; subplot(212plot(f,abs(Y2;axis(0

3、 fs 0 700subplot(313plot(abs(Y1% load sounddat y; %?ae?o?卩 ey?YFs=32000;Bits=16; channels=1; recorder = audiorecorder(Fs,Bits,cha nn els;%3e?o?sound(y,Fs;%2? 0?卩recordblock in g(recorder,4.0;x = getaudiodata(recorder;% i=o ? ?2?卩？ e ?(io_eybinopd% ?e? bin opdf K n, p ? p ?a? ?e? e e?t A -辰up? e 垃K?a

4、e?t A-getK，?n?ae? e Xu ?eyclear;px=binopdf(45,100,0.5 % ? ?x=45 p? e%x+?x=1:100;p=bi nopdf(x,100,0.5;plot(x,p;title(?e ?2? i?%?i e? i? 0? D?o?x1 X63 x22Z? oe? eu垃 e?ooo ? a ? ? ? p? p?x1o k2 戰(zhàn)%?dl?oey?：e?6 ? eu?x1 ?a?X?d?o eyclear alln=99;T=6;t=0:2*pi/( n-1:2*T*pi;N=le ngth(t;fs=1/( n-1;% x1=4*sin(t

5、; %? D?o?%y1=4*sin(t+pi/6; %? D?o? 6 3 u pi/6?x=awgn(x1,20; %? 6?1X? euD? eSNR=20dBy=awgn(y1,0.001; % ? x2=awg n( x1./4,-3;subplot(311,plot(t,xsubplot(312,plot(t,ytic%tico toc o?e 6 ?6u ?3Dd? ?e?Cc,r=xcorr(x,y; %xcorr 6 ?6U?6?6oey 垃Cc?0?d? 6 ,r? 6 3 ue?a?Cc1,r1=xcorr(x1,y1; % ? e? isubplot(313,plot(r*2*pi/( n-1,Cc,r,hold on; plot(Cc1,-;ticS仁fft(x,2*N-1; %? x p?2*N-1 QfftS2=fft(y,2*N-1;X=S1.* conj(S2; %? xGDy p? (?e?XR仁real(fftshift(ifft(X;% o ?ifft? ($? p? ?c1?o ey 垃 o ?fftshift Xa toc%? o tic? 0? ?e?figure(2subplot(211,plot(Cc,title(