function y=swerling2pu(n)
%由數據量為n的高斯白噪聲產生向量為n，功率譜為高斯型的斯維凌II型隨機向量
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
r=3;
[z1,z2]=gaussian(n);
x=gaussianpu(z1);
y=gaussianpu(z2);
y=x.^2+y.^2;
y=r*y;
b=canshu(y);
y=y-b(1);     %去掉直流分量
y=conv(y,y);
y=fft(y);
y=abs(y);
i=1:2*n-1;
plot(i,y)