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(幅度),xlabel