请老师指点
 

clear all;
clc;
close all;
r=36;
lamd=1.5;k=2*pi/lamd;
th=pi/5; 角度可以任意设定
omiga=2000*pi;
p0=1;
t=0:1/5000:(19/5000);
p=p0*exp(j*k*r-j*omiga*t); 这里的高斯白噪声 加的对吗？
S=awgn(p,20);
vx=cos(th)*S;
vy=sin(th)*S;
v1=0;T=4/1000;
dt=1/5000;
for i=1:length(t);
c=vx(i)*dt;
v1=v1+c;
end
gangvx=v1/T
v2=0;
for i=1:length(t);
d=vy(i)*dt;
v2=v2+d;
end
gangvy=v2/T
a=gangvy/gangvx;
TH=abs(atan(a))
JDWC=TH-th;
XDWC=(th-TH)/th;


这个程序的问题是 高斯白噪声跟本没有影响结果 我十分困惑