极大似然估计程序代码

% x=A+noise其中A为直流信号大小未知，noise 为AWGN均值为0，noise~N（1sigma^2）;
% 其中sigma=A;
% 估计出A来；
%x[n]=A+w[n]A is unknownw[n]is a AWGN NOISE
%w为高斯白噪声，方差为未知参数A；
%估计A
A=1;M=1000;
s=A*ones（1M）;
w=randn（1M）.*A;
x=w+s;
x_square=x.^2;

for N=1:M
I=ones（1N）;
I（M）=0;
A_estim（N）=-0.5+sqrt（（I*x_square‘）/N+0.25）;
end
A_estim_aver（1）=A_estim