clear;
P=0.1;
mu=0.01;
N=50000;
x=rand(1,N);
temp=var(x);
c=sqrt((P-mu^2)/var(x));
d=mu-c*mean(x);
u=c*x+d;
mean(u)
plot(u(1:200));
grid on;
ylabel('u(n)')
xlabel('n')
power_u=dot(u,u)
最后一句求的是平均功率；

dot()做内积，要求两个向量参数大小相同，先将对应元素相乘，然后求和