%?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