p3_8.m

来自「基于MMSE算法的自适应均衡LMS算法仿真」· M 代码 · 共 38 行

%%%%%%%-----Problem3.8-----%%%%%%%
%%Date  2007.4.17

%--------Settings--------%
clear all;
close all;
A=[-0.195 0.95;-1.5955 0.95;-1.9114 0.95];
segma_v2=1;                 %variance of noise
mu=0.001;                   %stepsize
N=1000;

Rxx=zeros(2,2);
Rxd=zeros(2,1);
I=eye(2);
w=zeros(2,N);

%-----Main process-----%
for n_c=1:3
    r0=(1+A(n_c,2))*segma_v2/((1-A(n_c,2))*((1+A(n_c,2))^2-A(n_c,1)^2));
    r1=(-1)*A(n_c,1)*segma_v2/((1-A(n_c,2))*((1+A(n_c,2))^2-A(n_c,1)^2));
    Rxx=[r0 r1;r1 r0];
    Rxd=[r0; r1];
    
    for n_i=1:N-1
        w(:,n_i+1)=(I-2*mu.*Rxx)*w(:,n_i)+2*mu.*Rxd;   %MSD
    end
       
    figure;   %plotting curves
    plot(1:N,w(1,:),'r-');hold on;
    plot(1:N,w(2,:));
    legend('w1','w2')
    xlabel('Iteration');
    ylabel('Weight');
    t_c=num2str(A(n_c,:));
    title({['Curve of a1,a2=[',t_c,']']}); 
    grid on;
end