Politechnika Poznańska Zastosowanie algorytmów genetycznych do ...

T2=100;
%DZ_z=1.4;
%DW_z=0.2;
%DW_w=0.19;
%1. współrzędne narożników
xA=-1.0; yA=-1.0;
xB=1.0; yB=-1.0;
xC=1.0; yC=1.0;
xD=-1.0; yD=1.0;
del=2/PK_z;
% 2. kontur zewnętrzny
% 2.1 punkty kollokacji na konturze zewnętrznym
for i=1:PK_z
XAB(i)=xA+0.5*del+(i-1)*del; YAB(i)=yA;
XBC(i)=xB;
YBC(i)=yB+0.5*del+(i-1)*del;
XCD(i)=xC-0.5*del-(i-1)*del; YCD(i)=yC;
XDA(i)=xD;
YDA(i)=yD-0.5*del-(i-1)*del;
end;
% 2.2 punkty źródłowe wokół konturu zewnętrznego
del = (2 + 2*DZ_z)/PZ_z;
for i=1:PZ_z
XABZ(i)= xA-DZ_z+0.5*del+(i-1)*del; YABZ(i)= yA-DZ_z;
XBCZ(i)= xB+DZ_z ;
YBCZ(i)= yB-DZ_z+0.5*del+(i-1)*del;
XCDZ(i)= xA-DZ_z+0.5*del+(i-1)*del; YCDZ(i)= yC+DZ_z;
XDAZ(i)= xD-DZ_z;
YDAZ(i)= yB-DZ_z+0.5*del+(i-1)*del;
XZR(i) = XABZ(i); YZR(i) = YABZ(i);
XZR(PZ_z+i) = XBCZ(i);
YZR(PZ_z+i) = YBCZ(i);
XZR(2*PZ_z+i) = XCDZ(i);
YZR(2*PZ_z+i) = YCDZ(i);
XZR(3*PZ_z+i) = XDAZ(i);
YZR(3*PZ_z+i) = YDAZ(i);
end;
% 3. kontury wewnętrzne
% 3.1 punkty kollokacji na konturach zewnętrznych
for k=1:KW_cnt
for i=1:PK_w
rd=i*2*pi/PK_w;
XKW(k,i)=xIC(k)+rIC*sin(rd);
YKW(k,i)=yIC(k)+rIC*cos(rd);
XKWA((k-1)*PK_w+i)=xIC(k)+rIC*sin(rd);
end
YKWA((k-1)*PK_w+i)=yIC(k)+rIC*cos(rd);
VNY(k,i)=cos(rd); %składowe wekt. norm.
VNX(k,i)=sin(rd);
end;
% 3.2.1 punkty źródłowe wokół konturu wewnętrznego (wew kont. wew.)
for i=1:PZ_ww
rd=i*2*pi/PZ_ww;
XKWZ(k,i)=xIC(k)+(rIC-1*DW_w)*sin(rd);
YKWZ(k,i)=yIC(k)+(rIC-1*DW_w)*cos(rd);
XZR(4*PZ_z+(k-1)*PZ_ww+i) = XKWZ(k,i);
YZR(4*PZ_z+(k-1)*PZ_ww+i) = YKWZ(k,i);
end;
% 3.2.1 punkty źródłowe wokół konturu wewnętrznego (na zew. kont. wew.)
for i=1:PZ_zw
rd=i*2*pi/PZ_zw;
XKWZ2(k,i)=xIC(k)+(rIC+1*DW_z)*sin(rd);
YKWZ2(k,i)=yIC(k)+(rIC+1*DW_z)*cos(rd);
XZR(4*PZ_z+KW_cnt*PZ_ww+(k-1)*PZ_zw+i) = XKWZ2(k,i);
YZR(4*PZ_z+KW_cnt*PZ_ww+(k-1)*PZ_zw+i) = YKWZ2(k,i);
end;
%4. Tworzenie macierzy układu
for i=1:PK_z
% pętla przez PKT. KOL. na KONTURZE ZEWNĘTRZNYM
iter=4*PZ_z+KW_cnt*PZ_ww;
for j=1:iter %pętla przez PKT. ZR. dla KONT. ZEW.
AP( i, j ) = f1(XAB(i), YAB(i), XZR(j), YZR(j)); % T1
AP( PK_z+i, j ) = df1dx(XBC(i), YBC(i), XZR(j), YZR(j)); % DT/DX=0
AP( 2*PK_z+i, j ) = f1(XCD(i), YCD(i), XZR(j), YZR(j)); % DT/DX=0
AP( 3*PK_z+i, j ) = df1dx(XDA(i), YDA(i), XZR(j), YZR(j)); % SM
end
for j=1:KW_cnt*PZ_zw
%pętla przez PKT. ZR. dla KONTURU WEWNĘTRZNEGO
AP( i, 4*PZ_z+KW_cnt*PZ_ww+ j ) = 0; % Te punkty źródłowe nie dotyczą KONT.
ZEW.
119