fengsumoxing.m

v10=20;
n=0.01:0.01:10;
xn=1200*n./v10;
k=0.001;
ti=0.1;
s1=4*k*v10^2*xn.^2./n./（1+xn.^2）.^（4/3）;
%定义空间点坐标
for i=1:20
    x（i）=5+i;
    z（i）=5+i;
end
%AR模型取4阶，p=4
%求R矩阵，分别求出R0R1R2R3R4
syms f
R0=zeros（20）;
for i=1:20
      for j=i:20
H0=inline（‘（4*k*v10^2*（1200*f/v10）.^2）./f./（1+（1200*f/v10）.^2）.^（4/3）.*exp（-sqrt（dx^2/50^2+dz^2/60^2））‘‘f‘‘k‘‘dx‘‘dz‘‘v10‘）;
       dx=x（i）-x（j）;
       dz=z（i）-z（j）;
       R0（ij）=quadl（H00.001100.0010kdxdzv10）;
       R0（ji）=R0（ij）;
      end
end
R1=zeros（20）;
for i=1:20
    for j=i:20
H1=inline（‘（4*k*v10^2*（1200*f/v10）.^2）./f./（1+（1200*f/v10）.^2）.^（4/3）.*exp（-sqrt（dx^2/50^2+dz^2/60^2））.*cos（2*pi*f*1*ti）‘‘f‘‘k‘‘dx‘‘dz‘‘ti‘‘v10‘）;
      dx=x（i）-x（j）;
      dz=z（i）-z（j）;
      R1（ij）=quadl（H10.001100.0010kdxdztiv10）;
      R1（ji）=R1（ij）;
     end
end
R2=zeros（20）;
for i=1:20
    for j=i:20
H2=inline（‘（4*k*v10^2*（1200*f/v10）.^2）./f./（1+（1200*f/v10）.^2）.^（4/3）.*exp（-sqrt（dx^2/50^2+dz^2/60^2））.*cos（2*pi*f*2*ti）‘‘f‘‘k‘‘dx‘‘dz‘‘ti‘‘v10‘）;
      dx=x（i）-x（j）;
      dz=z（i）-z（j）;
      R2（ij）=quadl（H20.001100.0010kdxdztiv10）;
      R2（ji）=R2（ij）;
     end
end
R3=zeros（20）;
for i=1:20
    for j=i:20
       H3=inline（‘（4*k*v10^2*（1200*f/v10）.^2）./f./（1+（1200*f/v10）.^2）.^（4/3）.*exp（-sqrt（dx^2/50^2+dz^2/60^2））.*cos（2*pi*f*3*ti）‘‘f‘‘k‘‘dx‘‘dz‘‘ti‘‘v10‘）;
      dx=x（i）-x（j）;
      dz=z（i）-z（j）;      R3（ij）=quadl（H30.001100.0010k