现代数字信号处理及其应用 第六章仿真题.docxVIP

6.12 clc; clear; M=15; Lb=10; % hb=[0.407 0.815 0.407]; hb=[0.04 -0.05 0.07 -0.21 -0.5 0.72 0.36 0 0.21 0.03 0.07]; Hb=zeros(2*M+1,2*M+Lb+1); for k=1:2*M+1 Hb(k,k:1:k+Lb)=hb; end EA1 = zeros(2000, 1); EA2 = zeros(2000, 1); for k=1:100 sigma=1e-3; N=2000; s=randsrc(2*M+Lb+N,1); vn=sqrt(sigma)*randn(2*M+N,1); S=zeros(2*M+Lb+1,N); V=zeros(2*M+1,N); for k=1:N S(:,k)=s(2*M+Lb+k:-1:k); V(:,k)=vn(2*M+k:-1:k); end Ub=Hb*S+V; errb_LMS=zeros(N,1); wb_LMS=zeros(2*M+1,N); wb_LMS(M+1,1)=1; dn=S(M+Lb+1,:); errb_LMS(1)=dn(1)-wb_LMS(:,1)*Ub(:,1); mu=0.025; for k=1:N-1 wb_LMS(:,k+1)=wb_LMS(:,k)+mu*Ub(:,k)*conj(errb_LMS(k)); errb_LMS(k+1)=dn(k+1)-wb_LMS(:,k+1)*Ub(:,k+1); end MSEb_LMS=abs(errb_LMS).^2; EA1=EA1+MSEb_LMS; lambda=0.99; delta=0.004; wb_RLS=zeros(2*M+1,N+1); wb_RLS(M+1,1)=1; epsilon=zeros(N,1); P1=eye(2*M+1)/delta; for k=1:N PIn=P1*Ub(:,k); deno=lambda+Ub(:,k)*PIn; kn=PIn/deno; epsilon(k)=dn(k)-wb_RLS(:,k)*Ub(:,k); wb_RLS(:,k+1)=wb_RLS(:,k)+kn*conj(epsilon(k)); P1=P1/lambda-kn*Ub(:,k)*P1/lambda; end MSEb_RLS=abs(epsilon).^2; EA2=EA2+MSEb_RLS; end M=15; Lb=2; hb=[0.407 0.815 0.407]; Hb=zeros(2*M+1,2*M+Lb+1); for k=1:2*M+1 Hb(k,k:1:k+Lb)=hb; end EA3 = zeros(2000, 1); EA4 = zeros(2000, 1); for k=1:100 sigma=1e-3; N=2000; s=randsrc(2*M+Lb+N,1); vn=sqrt(sigma)*randn(2*M+N,1); S=zeros(2*M+Lb+1,N); V=zeros(2*M+1,N); for k=1:N S(:,k)=s(2*M+Lb+k:-1:k); V(:,k)=vn(2*M+k:-1:k); end Ub=Hb*S+V; errb_LMS=zeros(N,1); wb_LMS=zeros(2*M+1,N); wb_LMS(M+1,1)=1; dn=S(M+Lb+1,:); errb_LMS(1)=dn(1)-wb_LMS(:,1)*Ub(:,1); mu=0.025; for k=1:N-1 wb_LMS(:,k+1)=wb_LMS(:,k)+mu*Ub(:,k)*conj(errb_LMS(k)); errb_LMS(k+1)=dn(k+1)-wb_LMS(:,k+1)*Ub(:,k+1); end MSEb_LMS=abs(errb_LMS).^2; EA3=EA3+MSEb_LMS; lambda=0.99; delta=0.004; wb_RLS=zeros(2*M+1,N+1); wb_RLS(M+1,1)=1; epsilon=zeros(N,1); P1=eye(2*M+1)/delta; for k=1:N PIn=P1*Ub(:,k); deno=lambda+Ub(:,k)*PIn; kn=PIn/deno; epsilon(k)=dn(k)-wb_RLS