以前随手写的,不过有错误(波的原理性错误),需要修改才对
a=1;
b=0;
s=a*exp(j*b/180*pi);
f=2e9;
lemada=4e8/f;
t=0:0.01e-9:4e-9;
z=-1+0.005:0.005:0;
wave_in=zeros(200,401);
wave_re=zeros(200,401);
for n=1:100
temp=2*n;
wave_in(1:temp,n)=exp(j*(2*pi*f*t(n)-2*pi/lemada*z(1:temp)'));
end;
for n=101:401
wave_in(:,n)=exp(j*(2*pi*f*t(n)-2*pi/lemada*z'));
end;
for n=1:100
temp=200-2*n+2;
wave_re(temp:200,100+n)=exp(j*(2*pi*f*t(n)+2*pi/lemada*z(temp:200)'))*s;
end;
for n=201:401
wave_re(:,n)=exp(j*(2*pi*f*t(n)+2*pi/lemada*z'))*s;
end;
wave=wave_in+wave_re;
wave_RE=real(wave);
for k = 1:401
L{k} =real(wave(: