برنامه درونیابی محلی به روش idw در متلب

برنامه را در mfile کپی کنید

%daroonyabi mahalli-Grid Based-idw%
%zp=?%
clc
clear all
format long g
a=input('xyz 4 nogte abcd[xa ya za;xb yb zb;xc yc zc;xd yd zd=');
p=input('xy p [x y]=');
for i=1:4;
d(i)=sqrt((p(1,1)-a(i,1))^2+(p(1,2)-a(i,2))^2);
end;
w1=(1/d(1).^2)/sum(1./(d.^2));
w2=(1/d(2).^2)/sum(1./(d.^2));
w3=(1/d(3).^2)/sum(1./(d.^2));
w4=(1/d(4).^2)/sum(1./(d.^2));
zp=(a(1,3)*w1)+(a(2,3)*w2)+(a(3,3)*w3)+(a(4,3)*w4)