调试SQL、igrfgridData编程辅导、辅导Java编程、调试SQL语言

load igrfgridData.txt % 1
[a,b,c]=size(igrfgridData)
load igrfdata.txt % 1
[d,e,f]=size(igrfdata)
for i=1:a
X(i,1)=igrfgridData(i,1);
X(i,2)=igrfgridData(i,2);
X(i,3)=igrfgridData(i,3);
xw(i)=X(i,1);
yw(i)=X(i,2);
end
for i=1:d
U(i,1)=igrfdata(i,1);
U(i,2)=igrfdata(i,2);
U(i,3)=igrfdata(i,3);
U(i,4)=igrfdata(i,4);
U(i,5)=igrfdata(i,5);
U(i,6)=igrfdata(i,6);
U(i,7)=igrfdata(i,7);
U(i,8)=igrfdata(i,8);
U(i,9)=igrfdata(i,9);
U(i,10)=igrfdata(i,10);
U(i,11)=igrfdata(i,11);
U(i,12)=igrfdata(i,12);
U(i,13)=igrfdata(i,13);
U(i,14)=igrfdata(i,14);
U(i,15)=igrfdata(i,15);
U(i,16)=igrfdata(i,16);
U(i,17)=igrfdata(i,17);
end
xwmax=max(xw);
xwmin=min(xw);
ywmax=max(yw);
ywmin=min(yw);
nx=round((xwmax-xwmin).*111000./a);
ny=round((ywmax-ywmin).*75000./a);
m1=0:16;
n1=0:16;
m=zeros(d,1);
n=zeros(d,1);
K1=(mod(d,2)==0);
K2=(mod(d,2)==0);
switch K1
case 0
for i=1:d
if m1(i)<=fix(d/2)
m(i)=m1(i);
else
m(i)=m1(i)-d;
end
end
case 1
for i=1:a
if m1(i)<=a/2-1
m(i)=m1(i);
else
m(i)=m1(i)-a;
end
end
end
switch K2
case 0
for i=1:d
if n1(i)<=fix(d/2)
n(i)=n1(i);
else
n(i)=n1(i)-d;
end
end
case 1
for i=1:b
if n1(i)<=b/2-1
n(i)=n1(i);
else
n(i)=n1(i)-b;
end
end
end
Lx=d*nx;
Ly=d*ny;
kx=m/(Lx);
ky=n/(Ly);
A=zeros(d,d);
H=input('')
for ix=1:d
for iy=1:d
A(ix,iy)=exp(-2*pi*H*sqrt(kx(ix)^2+ky(iy)^2));
end
end
UF=fft2(U);
UF=UF.*A;
UA=ifft2(UF);