bd503b752b2ccf3d8e0912b7522098d7.m

clear;
clf reset;
close all;

q=2;
dt=0.1;
w=[3 8];
z=exp（1i*w*dt）;
N=100;
C=[4*exp（1i*pi/8） exp（1i*pi/2.5）]‘;
sigma2=2;

for k=1:N
    H（k:）=z.^k;
end

n=sqrt（sigma2）*randn（N1）+1i*sqrt（sigma2）*randn（N1）;
r=H*C+n;

C_estim=（[H（1:）;H（2:）]‘*[H（1:）;H（2:）]）^（-1）*[H（1:）;H（2:）]‘*[r（1:）;r（2:）]; % C_estim=theta chapeau
J=（1/sigma2）*H（1:）‘*H（1:）; % J（1）
J=J（1）+（1/sigma2）*H（2:）‘*H（2:）; % J（2）
C_ALL=[];
for i=3:N
    J=J+（1/sigma2）*H（i:）‘*H（i:）;
    C_estim=C_estim+（sigma2+H（i:）*J^（-1）*H（i:）‘）^（-1）*J^（-1）*H（i:）‘*（r（i:）-H（i:）*C_estim）;
    C_ALL（:i-2） = C_estim;
end

figure（2）
clf; subplot（211）; hold on
plot（[3 N] [real（C（1