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.
end<br />
Temp=I1/I2;<br />
newline=NaN;<br />
plik='C:\wyniki\kolo4v2.mat'; %definicja ścieżki <strong>do</strong> pliku wynikowego<br />
save(plik,'cnt','newline','newline','Temp','newline','newline','x','newline','newline','<br />
hN','newline','newline','hN_tmp','newline','newline','x','newline','newline','newline','<br />
newline','-ascii','-append','-tabs');<br />
t=I1/I2;<br />
%funkcja testowa<br />
function tN=yvalN(x,varLim,bitsNV,popN,varN)<br />
kN=zeros(1,popN);<br />
for i=1:popN<br />
for j=1:varN<br />
hN_tmp(i,j)=(x(i,j)/2^bitsNV(j))*(varLim(j,2)-varLim(j,1))+varLim(j,1);<br />
end<br />
for j=1:2:7<br />
jjk=(j-1)/2;<br />
hN(i,j)=hN_tmp(i,1)*cos(pi/4+jjk*pi/2);<br />
hN(i,j+1)=hN_tmp(i,1)*sin(pi/4+jjk*pi/2);<br />
end<br />
%disp(hN);<br />
flclear fem<br />
% Geometry<br />
g1=ellip2(1,1,'base','center','pos',[0,0]);<br />
g2=ellip2(0.05,0.05,'base','center','pos',[hN(i,1),hN(i,2)]);<br />
g3=ellip2(0.05,0.05,'base','center','pos',[hN(i,3),hN(i,4)]);<br />
g4=ellip2(0.05,0.05,'base','center','pos',[hN(i,5),hN(i,6)]);<br />
g5=ellip2(0.05,0.05,'base','center','pos',[hN(i,7),hN(i,8)]);<br />
[g10,ctx]=geomcomp({g1,g2,g3,g4,g5},'ns',{'g1','g2','g3','g4','g5'},'sf','g1-g2-g3-g4-<br />
g5','edge','none','out',{'ctx'});<br />
[allmp,nbs] = geominfo(g10,'out',{'mp','nbs'},'od',[0;1]);<br />
dd=ctx{1};<br />
for ii=1:nbs<br />
for jj=1:4<br />
wb=dd(ii,jj);<br />
if(~(wb==0))<br />
tc(jj)=ii;<br />
end<br />
end<br />
end<br />
clear s<br />
s.objs={g10};<br />
s.name={'CO1'};<br />
s.tags={'g10'};<br />
fem.draw=struct('s',s);<br />
fem.geom=geomcsg(fem);<br />
% Initialize mesh<br />
fem.mesh=meshinit(fem);<br />
% Application mode 1<br />
clear appl<br />
appl.mode.class = 'HeatTransfer';<br />
appl.assignsuffix = '_ht';<br />
clear prop<br />
prop.analysis='static';<br />
appl.prop = prop;<br />
clear bnd<br />
bnd.type = 'T';<br />
bnd.T0 = {343,303};<br />
bnd.ind = [1,1,1,1, 1,1,1,1, 1,1,1,1, 1,1,1,1, 1,1,1,1];<br />
for jj=1:4<br />
bnd.ind(tc(jj))=2;<br />
end<br />
appl.bnd = bnd;<br />
clear equ<br />
equ.init = 303;<br />
equ.ind = [1];<br />
appl.equ = equ;<br />
fem.appl{1} = appl;<br />
fem.border = 1;<br />
fem.outform = 'general';<br />
fem.units = 'SI';<br />
112