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

Załącznik 3 – kod programu - wyznaczanie efektywnego
współczynnika przewodzenia ciepła
Program napisano w języku Matlab/Comsol Script
function strAB=lapl3v2s(rIC,SP,plt)
%IRR=[1.8875,1.8262,0.26313,12,10,12];
IRR=[1.0628,6.1328,0.071875,9,9,31];
DZ_z=IRR(1);
DW_z=IRR(2);
DW_w=IRR(3);
PZ_z=IRR(4);
PZ_zw=IRR(5);
PZ_ww=IRR(6);
XAB=[];YAB=[];XBC=[];YBC=[];XCD=[];YCD=[];XDA=[];YDA=[];
XKW=[];YKW=[];XZR=[];YZR=[];errTW=[];errDTW=[];
XABZ=[];YABZ=[];XBCZ=[];YBCZ=[];XCDZ=[];YCDZ=[];XDAZ=[];YDAZ=[];
XKWZ=[];YKWZ=[];XKWZ2=[];YKWZ2=[];VNX=[];VNY=[];
xIC2=[];
yIC2=[];
varN=length(SP);
for j=2:2:varN
xIC(j/2)=SP(j-1);
yIC(j/2)=SP(j);
end
%{
xIC(1)=SP(j);
yIC(1)=SP(j);
xIC(2)=SP(3);
yIC(2)=SP(4);
xIC(3)=SP(5);
yIC(3)=SP(6);
%}
[v,KW_cnt]=size(xIC);
AP=[];
%M - ilośc punktów kollokacji na odc. jednostkowym
%MZ - ilośc punktów źródłowych na odc. jednostkowym
%NZ - Całkowita ilośc punktów źródłowych
%rIC - promien konturu wewnętrzengo
%xIC - współrzędna X śrdka konturu wew.
%YCkw - współrzędna Y śrdka konturu wew.
%L1 - wsp. przewodzenia osnowy
%L2 - wsp przewodzenia włókna
%T1 - Temp. na brzegu AB
%T2 - Temp. na brzegu CD
PK_z=80;
PK_w=50;
%PZ_z=9;
%PZ_zA=4*10;
%PZ_zw=12;
%PZ_ww=8;
%PZ_zw=15;
%PZ_ww=5;
%xIC=0.7;
%yIC=0.7;
L1=2.0;
L2=0.1;
T1=0;
118