v10=17.0;
n=0.01:0.01:10;
xn=1200*n./v10;
k=0.03;
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矩阵,分别求出R0,R1,R2,R3,R4
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(i,j)=quadl(H0,0.001,10,0.001,0,k,dx,dz,v10);
R0(j,i)=R0(i,j);
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