MATLAB数字式自己振荡器及其应用.doc

%提高题 数字式自激振荡器及其应用clc,clear,close all,N1=6;% - 差分方程y(n)-y(n-1)*2cos(2pi/N)+y(n-2)=x(n)-cos(2pi/N)*x(n-1)a1 =1 -2*cos(2*pi/N1) 1; %系统函数H(z)分母多项式系数b1 =1 -cos(2*pi/N1) 0; %系统函数H(z)分子多项式系数uN1=1,zeros(1,79999); x1=uN1; ycos=filter(b1,a1,x1)N2=6;% - 差分方程y(n)-y(n-1)*2cos(2pi/N)+y(n-2)=x(n)+sin(2pi/N)*x(n-1)a2 =1 -2*cos(2*pi/N2) 1; %系统函数H(z)分母多项式系数b2 =0 sin(2*pi/N2) 0; %系统函数H(z)分子多项式系数uN2=1,zeros(1,79999); x2=uN2; ysin=filter(b2,a2,x2)figure(1);subplot(221)plot(ycos)xlabel(x),ylabel(cosx),title(cosx的时域);subplot(222)plot(abs(fft(ycos)ylabel(幅度),xlabel(频率/Hz),title(cosx的频域谱);subplot(223)plot(ysin)xlabel(x),ylabel(sinx),title(sinx的时域);subplot(224)plot(abs(fft(ysin)ylabel(幅度),xlabel(频率/Hz),title(sinx的频域谱); figure(2);N=80000;Fs=22000;Ts=Fs/1;n=0:N-1;ss1,fs1,nbit1 = wavread(aaa.wav);%语音信号调用s1=ss1(1:N);subplot(221);plot(n/fs1,s1);title(原信号1);ylabel(幅度),xlabel(t/s);f1=fft(s1);%傅里叶变换f=n/N*fs1;subplot(222);F1=abs(f1);plot(f,F1);title(原信号1的幅度频谱);ylabel(幅度),xlabel(频率/Hz);sound(s1,fs1,nbit1);ss2,fs2,nbit2 = wavread(bbb.wav);%语音信号调用s2=ss2(1:N);t2=1/fs2;subplot(223)plot(n/fs2,s2)title(原信号2);ylabel(幅度),xlabel(t/s);f2=fft(s2);%傅里叶变换f=n/N*fs2;subplot(224);F2=abs(f2);plot(f,F2);title(原信号2的幅度频谱);ylabel(幅度),xlabel(频率/Hz);sound(s2,fs2,nbit2); figure(3);y=s1.*ycos+s2.*ysin;%调制Y=fft(y);%傅里叶变换subplot(211)plot(n/N*Ts,abs(Y);title(调制后幅度谱)ylabel(幅度),xlabel(频率/Hz);subplot(212)plot(n/Ts,y)title(调制后信号);ylabel(幅度),xlabel(t/s);sound(y) figure(4)a1=ones(1,15000);%低通滤波器aa=zeros(1,N-15000*2);a2=a1 aa a1;z1=y.*ycos;%解调，时域相乘z2=y.*ysin;Z1=fft(z1);%傅里叶变换v1=Z1.*a2;%低通滤波V1=real(ifft(v1);%傅里叶反变换E1=fft(V1);E1=fftshift(E1);subplot(221);plot(n/fs1,V1);title(还原后信号1),ylabel(幅度),xlabel(t/s);subplot(222)plot(n/N*Ts,abs(E1);title(还原后信号1幅度频谱);xlabel(频率/Hz),ylabel(幅度);sound(V1,fs1)Z2=fft(z2);v2=Z2.*a2;V2=real(ifft(v2);%取实数部分E2=fft(V2);subplot(223);plot(n/fs2,V2);title(还原后信号2),ylabel(幅度),xla