WYKŁAD 14

INTERPOLACJA FUNKCJI DWÓCH ZMIENNYCH

....................................................
N:=50;
h:=2*Pi/N;
for i:=0 to N do begin
x:=-Pi+i*h;
for j:=0 to N do begin
y:=-Pi+j*h;
fun[i,j]:=Sin(x-1.5*y)*Sqrt(Abs(x*y)/3);
end;
end;
N1:=N-1;
{.......................................}
ZX:=GetMaxX div 2; ZY:=GetMaxY;
a:=ZY-50; b:=ZY/2;
xw[1]:=ZX-a/2; xw[2]:=ZX+a/2;
xw[3]:=xw[2]; xw[4]:=xw[1];
yw[1]:=b+a/2; yw[2]:=yw[1];
yw[3]:=b-a/2; yw[4]:=yw[3];
h:=a/N; M:=Round(h);
SetBkColor(White); SetColor(Blue);
min:=fun[0,0]; max:=fun[0,0];
for i:=0 to N do
for j:=0 to N do begin
pp:=fun[i,j];
if pp<min then min:=pp;
if pp>max then max:=pp;
end;
del:=(max-min)/25;

for j:=0 to N1 do begin
yp:=yw[1]-j*h;
for i:=0 to N1 do begin
xp:=xw[1]+i*h;
for jj:=0 to M do begin
s:=jj/h;
s1:=1-s;
for ii:=0 to M do begin
r:=ii/h;
r1:=1-r;
pp:=s1*(r1*fun[i,j]+r*fun[i+1,j])+
s*(r1*fun[i,j+1]+r*fun[i+1,j+1]);
K:=Trunc((pp-min)/del);
if Odd(K) then

PutPixel(Round(xp+ii),Round(ypjj),

LightGray)
else
PutPixel(Round(xp+ii),Round(yp-jj),White);
end;
end;
end;
end;
..................................................................



 



,

5

.

1

,

5

.

1

5

.

1

,

5

.

1

,

2

)

,

(

)

(

2

2

2













D

y

x

e

y

x

f





,

]

,

[

]

,

[

,

2

sin

)

,

(

2

2





















D

y

x

y

x

f

KONIEC