Un.rar_Matlab 录音

fs = 16000 ; % \ ?豕?????迄???? y =wavrecord(3*fs ,fs ,'uint8'); %????????8??????????,?㊣????5s subplot(411); plot (y) ; wavplay(y ,fs) ; % ?????迄?????????? Y1 =fft (y) ; %℅???????fft ㊣??? Y2 =fftshift(Y1); f=((0:length(Y2)-1)/length(Y2))*fs-fs/2; subplot(412) plot(f,abs(Y1)); subplot(413) plot(abs(Y2)) subplot(414) N=length(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=sin(2*pi*f*t); subplot(211) plot(t,y); Y1=fft(y); Y2 =fftshift (Y1); f=(0:length(Y1)-1)/length(Y1)*fs-fs/2; subplot(212) plot(f,abs(Y2)); axis([0 fs 0 700]) subplot(313) plot(abs(Y1)) Fs=32000;Bits=16; channels=1; recorder = audiorecorder(Fs,Bits,channels);%???????? sound(y,Fs);%???????? recordblocking(recorder,4.0); x = getaudiodata(recorder); clear; px=binopdf(45,100,0.5) % ????x=45?????? %℅‾?? x=1:100; p=binopdf(x,100,0.5); plot(x,p);title('??????????') clear all n=99; T=6; t=0:2*pi/(n-1):2*T*pi; N=length(t); fs=1/(n-1); x=awgn(x1,20); %??????∼℅???迄??????㊣?SNR==20dB y=awgn(y1,0.001); %???? x2=awgn(x1./4,-3); subplot(311),plot(t,x) subplot(312),plot(t,y) tic %tic??toc???????????????辰?????㊣?? [Cc,r]=xcorr(x,y); %xcorr?????車???角?????車??????Cc?????角?????車,r?????????? [Cc1,r1]=xcorr(x1,y1); %???? subplot(313),plot(r*2*pi/(n-1),Cc,'r'),hold on; plot(Cc1,'-'); tic S1=fft(x,2*N-1); %?車x??2*N-1??fft S2=fft(y,2*N-1); X=S1.* conj(S2); %?車x??y?????????℅ R1=real(fftshift(ifft(X)));%????ifft?車??㊣??????????角??????????fftshift??℅? toc %????tic?車?????㊣?? figure(2) subplot(211),plot(Cc),title('xcorr') subplot(212),plot(R1); title('ifft') figure(3) plot(xcorr(x2,'unbiased'));