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

AP( PK_z+i, 4*PZ_z+KW_cnt*PZ_ww+ j ) = 0;
AP( 2*PK_z+i, 4*PZ_z+KW_cnt*PZ_ww+ j ) = 0;
AP( 3*PK_z+i, 4*PZ_z+KW_cnt*PZ_ww+ j ) = 0;
end
end
for k=1:KW_cnt % pętla przez KONTURY WEWNĘTRZNE
for i=1:PK_w
% pętla przez PKT. KOL. na KONTURZE WEWNĘTRZNYM
iter=4*PZ_z+KW_cnt*PZ_ww;
for j=1:iter %pętla przez PKT. ZR. dla KONTURU ZEWNĘTRZNEGO
AP( 4*PK_z+(k-1)*PK_w+i, j ) = f1(XKW(k,i), YKW(k,i), XZR(j), YZR(j)); % TEMP.
CONT.
AP( 4*PK_z+(KW_cnt+k-1)*PK_w+i, j ) = L1*df1dn(XKW(k,i), YKW(k,i), XZR(j),
YZR(j),xIC(k),yIC(k)); % DERIV. CONT.
end
for j=1:PZ_zw
%pętla przez PKT. ZR. dla KONTURU WEWNĘTRZNEGO
AP( 4*PK_z+(k-1)*PK_w+i, 4*PZ_z+KW_cnt*PZ_ww+(k-1)*PZ_zw+j ) = -f1(XKW(k,i),
YKW(k,i), XZR(4*PZ_z+KW_cnt*PZ_ww+(k-1)*PZ_zw+j), YZR(4*PZ_z+KW_cnt*PZ_ww+(k-
1)*PZ_zw+j)); % TEMP. CONT.
AP( 4*PK_z+(KW_cnt+k-1)*PK_w+i, 4*PZ_z+KW_cnt*PZ_ww+(k-1)*PZ_zw+j ) = -
L2*df1dn(XKW(k,i), YKW(k,i), XZR(4*PZ_z+KW_cnt*PZ_ww+(k-1)*PZ_zw+j),
YZR(4*PZ_z+KW_cnt*PZ_ww+(k-1)*PZ_zw+j),xIC(k),yIC(k)); % DERIV. CONT.
end
end
end
%xlswrite('c:\APv3.xls',AP);
%5. Warunki brzegowe
BP=zeros(4*PK_z+2*KW_cnt*PK_w,1); %wektor pełen zer, więc nie jest zmieniany jeśli na
danym brzegu jest 0.
for i=1:PK_z
BP(0*PK_z+i) = T1;
BP(2*PK_z+i) = T2;
BP(1*PK_z+i) = 0;
BP(3*PK_z+i) = 0;
end
X=linsolve(AP,BP);
iX=[];iY=[];iU=[];
aa=0;
for i=-1:0.02:1
aa=aa+1;
bb=0;
for j=-1:0.02:1
bb=bb+1;
dfu=roz(i,j,KW_cnt,PZ_z,PZ_ww,PZ_zw,X,XZR,YZR,XKW,YKW,XKWZ2,YKWZ2);
iX(aa,bb)=i;
iY(aa,bb)=j;
iU(aa,bb)=dfu;
end
end
% obliczenie wartości strumienia ciepła
step2=0.01;
str=0;
for i=-1:step2:1
str=str+stream(i,-1,KW_cnt,PZ_z,PZ_ww,PZ_zw,X,XZR,YZR,XKW,YKW,XKWZ2,YKWZ2,L1,L2);
end
strAB=(str*(yC-yA))/((2/step2)*(T2-T1));
%sprawdzenie warunków
step=0.01;
%brzeg AB
errAB=0;
for j=-1:step:1
vT=0;
iter=4*PZ_z+KW_cnt*PZ_ww;
for i=1:iter
vT=vT+X(i)*f1(j,-1,XZR(i),YZR(i));
end
errAB=errAB+(T1-vT)^2;
end
errAB=(errAB^0.5)/(2/step);
%brzeg BC
errBC=0;
for j=-1:step:1
120