%%%?????Rx（0）?=?mean（x（1:n）^2）;%%%%%%%%%%%%%
%%%?????Ry（0）?=?mean（y（1:n）^2）;%%%%%%%%%%%%%
%%%?????Rxy（0）?=?mean（x（1:n）?.*?y（1:n））;%%%%%%%%%%%%%

f?=?50;
w?=?2*pi*f;
N?=?100;
delt_t?=?1/f/100;
t?=?0?:?delt_t?:?（N-1）*delt_t;


noise1?=?0.5*rand（