Politechnika Poznańska Zastosowanie algorytmów genetycznych do ...
Politechnika Poznańska Zastosowanie algorytmów genetycznych do ...
Politechnika Poznańska Zastosowanie algorytmów genetycznych do ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
colormap(jet(128))<br />
% surf(iX,iY,iU);<br />
shading interp;<br />
colorbar();<br />
hold<br />
scatter(XZR,YZR,3,[0 0 0],'filled')<br />
scatter(XAB,YAB,3,[0 0 1],'filled')<br />
scatter(XBC,YBC,3,[0 0 1],'filled')<br />
scatter(XCD,YCD,3,[0 0 1],'filled')<br />
scatter(XDA,YDA,3,[0 0 1],'filled')<br />
scatter(XKWA,YKWA,3,[0 0 1],'filled')<br />
view(2);<br />
end<br />
%<br />
disp({'errAB=',errAB,'errBC=',errBC,'errCD=',errCD,'errDA=',errDA,'errT=',errT,'errDT=',<br />
errDT})<br />
%}<br />
% rozwiązanie podstawowe<br />
function v=f1(x,y,xs,ys)<br />
r=realsqrt((x-xs)^2+(y-ys)^2);<br />
v=reallog(r);<br />
% rozwiązanie podstawowe - symetryczne<br />
function v=f2(x,y,xs,ys)<br />
r1=realsqrt((x-xs)^2+(y-ys)^2);<br />
r2=realsqrt((x+xs)^2+(y-ys)^2);<br />
v=(reallog(r1)-reallog(r2));<br />
function v=df1dx(x,y,xs,ys)<br />
r2=(x-xs)^2+(y-ys)^2;<br />
v=(x-xs)/r2;<br />
function v=df1dy(x,y,xs,ys)<br />
r2=(x-xs)^2+(y-ys)^2;<br />
v=(y-ys)/r2;<br />
function v=df1dx2(x,y,xs,ys)<br />
r2=(x-xs)^2+(y-ys)^2;<br />
v=-((x - xs)^2 - (y - ys)^2)/r2^2;<br />
function v=df1dy2(x,y,xs,ys)<br />
r2=(x-xs)^2+(y-ys)^2;<br />
v=((x - xs)^2 - (y - ys)^2)/r2^2;<br />
function v=df1dn(x,y,xs,ys,xIC,yIC)<br />
r=realsqrt((x-xIC)^2+(y-yIC)^2);<br />
vnx=(x-xIC)/r;<br />
vny=(y-yIC)/r;<br />
fx=df1dx(x,y,xs,ys);<br />
fy=df1dy(x,y,xs,ys);<br />
w=vnx*fx+vny*fy;<br />
%disp({x,y,vnx,vny,vnx^2,vny^2});<br />
v=w;<br />
function v=roz(x,y,KW_cnt,PZ_z,PZ_ww,PZ_zw,XR,XZR,YZR,XKW,YKW,XKWZ,YKWZ)<br />
vT=0;<br />
p=0;<br />
for k=1:KW_cnt<br />
p=k;<br />
xrow=XKW(k,:);<br />
yrow=YKW(k,:);<br />
a=inpolygon(x,y,xrow,yrow);<br />
if(a),break,end<br />
end<br />
if(a)<br />
for i=1:PZ_zw<br />
vT=vT+XR(4*PZ_z+KW_cnt*PZ_ww+(p-1)*PZ_zw+i)*f1(x,y,XKWZ(p,i),YKWZ(p,i));<br />
end<br />
else<br />
122