matlab反正弦分布的代码,MATLAB里面怎么反正弦的值为虚数了呢？

clear

clc

time=cputime;

%% 初始化值

%%%%%%%%%%%%%%%%%%%%

dt=0.01*10^-9;%采样间隔

xt=0.06*10^-7;%子波长度

ht=100*10^-9; %时窗

fm=200*10^6; %天线频率

jpl=2*pi*fm;%天线角频率

c=0.3;%光速，单位ns/m

e0=8.85*10^-12;%真空介电常数

u0=4*pi*10^-7;%真空磁导率

e1=10;%湿粘土相对介电常数

v1=c/sqrt(e1);%shi粘土的电磁波速度

e2=1;%g管内填充物的介电常数(air)

v2=c;%g管内填充物的速度

zy=1e+2;%增益

Y=2;%管线埋深=Y-R1

ddl1=0.01;%湿粘土导电率，单位S/m

ddl2=1e-8;%air

O=0;%偏移距

ny=ht/dt;

k=0;

csd11=-0.9;

csd21=-0.4662;

csd22=0.2384;

csd23=0.0035;

%%%%%%%%%

gaga1=sqrt(1+(ddl1/(jpl*e0*e1))^2);

xsxs1=sqrt(0.5*u0*e0*e1*jpl^2*(gaga1-1));%第一层吸收系数

gaga2=sqrt(1+(ddl2/(jpl*e0*e2))^2);

xsxs2=sqrt(0.5*u0*e0*e2*jpl^2*(gaga2-1));%第二层吸收系数

fsxs1=(sqrt(e1)-sqrt(e2))/(sqrt(e1)+sqrt(e2));%电磁波在第一层中的反射系数

%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%for X0=-4:0.1:0

X0=-1;

X=X0-0.5*O;

absX=abs(0.5*O-X0);

ny=ht/dt;

for i=0:ny

h=0.0;

end %对振幅序列初始化</