随机实验报告2CJ.doc

程序 2 N=1024; Ts=1/50000; n=1:1:N; w=randn(1,N); a=5; f1=a*1000; f2=a*1000*2; f3=a*1000*3; x=sin(2*pi*f1*n*Ts)+sin(2*pi*f2*n*Ts)+sin(2*pi*f3*n*Ts)+w(n); R=xcorr(x); [S,W]=periodogram(x); m=-(N-1):1:(N-1); %f=n/(N*Ts); figure(1) subplot(3,1,1);plot(n,x); subplot(3,1,2);plot(m,R); subplot(3,1,3);plot(W,S); fp=f1;fs=f2;ap=1;as=40; %数字滤波器指标 T=Ts,fsa=1/T; %采样频率与间隔 wp=2*pi*fp/fsa;ws=2*pi*fs/fsa; %转换为数字角频率 Wp=2/T*tan(wp/2);Ws=2/T*tan(ws/2); %由数字角频率转换为模拟角频率 [N1,Wc]=buttord(Wp,Ws,ap,as,s); %获取模拟滤波器的阶数和3dB截止频率 [Z,P,K]=buttap(N1); %归一化模拟滤波器模型的零极点形式参数 [B,A]=zp2tf(Z,P,K); %归一化模拟滤波器传递函数的系数 [Bl,Al]=lp2lp(B,A,Wc); %把模拟滤波器原型转换成截至频率为Wc的低通滤波器 [b,a]=bilinear(Bl,Al,fsa);%用双线性变换法实现模拟滤波器到数字滤波器的转换 [H,w]=freqz(b,a); %获取频率响应 figure(2) plot(w*fsa/(2*pi),abs(H));grid; %绘制频率响应曲线 xlabel(频率(Hz));ylabel(频率响应幅度); y=filter(b,a,x); Ry=xcorr(y); [Sy,Wy]=periodogram(y); m=-(N-1):1:(N-1); %f=n/(N*Ts); figure(3) subplot(3,1,1);plot(n,y); subplot(3,1,2);plot(m,Ry); subplot(3,1,3);plot(Wy,Sy); 结果 程序2 clc; R_x=zeros(1,81);R_x(41)=1; % S_x=fftshift(abs(fft(R_x))); % No=39; % r=1-1/(No+1); h0=zeros(1,40); i=1:41; h1=r.^i; h=[h0,h1]; % H=fftshift(abs(fft(h)));% m_x=0; % sigma_x=R_x(41); P_x=R_x(41); figure(1), subplot(221),stem(R_x),title(Rx); subplot(222),stem(S_x),title(Sx); subplot(223),stem(h),title(h); subplot(224),stem(H),title(H); %%%%%%%%%%%%%%%%%% %时域法求解 R_xy=conv(R_x,h);R_xy=R_xy(41:121); R_yx=conv(R_x,fliplr(h));R_yx=R_yx(41:121); R_y=conv(R_yx,h);R_y=R_y(41:121); figure(2), subplot(221),stem(R_x);title(Rx); subplot(222),stem(R_xy);title(Rxy); % subplot(223),stem(R_yx);title(Ryx); subplot(224),stem(R_y);title(Ry); % S_xy=abs(fft(R_xy));S_xy=fftshift(S_xy); S_yx=fftshift(abs(fft(R_yx)));S_y=fftshift(abs(fft(R_y))); figure(3), subplot(221),stem(S_x);title(Sx); subplot(222),stem(S_xy);title(Sxy); % subplot(223),stem(S_yx);title(Syx); subplot(224),stem(S_y);title(Sy);