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

6.12

clc;

clear;

M=15;

Lb=10;

%hb=[0.4070.8150.407];

hb=[0.04-0.050.07-0.21-0.50.720.3600.210.030.07];

Hb=zeros(2*M+1,2*M+Lb+1);

fork=1:2*M+1

Hb(k,k:1:k+Lb)=hb;

end

EA1=zeros(2000,1);

EA2=zeros(2000,1);

fork=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);

fork=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;

fork=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;

fork=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.4070.8150.407];

Hb=zeros(2*M+1,2*M+Lb+1);

fork=1:2*M+1

Hb(k,k:1:k+Lb)=hb;

end

EA3=zeros(2000,1);

EA4=zeros(2000,1);

fork=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);

fork=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;

fork=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;

fork=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;

EA4=EA4+MSEb_RLS;

end

%figure

plot(EA1/100);

holdon

plot(EA2/100);

holdon

plot(EA3/100);

holdon

plot(EA4/100);

6.15

clc;

clear;

EA1=zeros(999,1);

A1=zeros(2,1000);

fori=1:100

a1=0.99;

sigma=0.995;

N=1000;

vn=sqrt(sigma)*randn(N,1);

nume=1;

deno=[1a1];

u0=zeros(length(deno)-1,1);

xic=filtic(nume,deno,u0);

un=filter(nume,deno,vn,xic);

n0=1;

M=2;

b=un(n0+1:N);

L=length(b);

un1=[zeros(M-1,1).',un.'];

A=zeros(M,L);

fork=1:L

A(:,k)=un1(M-1+k:-1:k);

end

delta=0.004;

lambda=0.98;

w=zeros(M,L+1);

epsilon=zeros(L,1);

P1=eye(M)/delta;

fork=1:L

PIn=P1*A(:,k);

denok=lambda+A(:,k)'*PIn;

kn=PIn/denok;

epsilon(k)=b(k)-w(:,k)'*A(:,k);

w(:,k+1)=w(:,k)+kn*conj