13.04.2014 Views

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

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

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

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!