Diplomska naloga (.pdf)

51
M=repmat(eye(2), [1, 1, int]);
alpha=zeros(int,st_plasti);
calpha=zeros(int,st_plasti);
salpha=zeros(int,st_plasti);
for j=1:st_plasti
alpha(:,j)=(2*pi*interval.*(mi(:,P(1,j,k))*E(P(1,j,k))).^0.5.*(P(2,j,k)))/3e8;
calpha(:,j)=cos(alpha(:,j));
salpha(:,j)=sin(alpha(:,j));
for n=1:int
N=M(:,:,n);
M(:,:,n)=N*[calpha(n,j) -i*(mi(n,P(1,j,k))/E(P(1,j,k))).^0.5.*salpha(n,j);
-i*salpha(n,j)./(mi(n,P(1,j,k))/E(P(1,j,k))).^0.5 calpha(n,j)];
end
end
q=0;
for n=1:int
R(n)=(-M(1,2,n)-M(2,2,n))/(-M(1,2,n)+M(2,2,n));
R(n)=20*log10(abs(R(n)));
%Preverimo, ali za vsako vrednost iz intervala ustreza dopustnemu pogoju.
if R(n)>pogoj
q=1;
end
end
%Če optimiziramo po debelini absorberja.
if opt==0
if q==1
%Označimo nedopustno rešitev.
f(k)=1;
else
%Vrnemo debelino absorberja.
f(k)=sum(P(2,:,k));
end
%Če optimiziramo po širini frekvenčnega pasu absorberja.
elseif opt==1
clear y
clear z
if R(int1)>=pogoj
f(k)=0; %Označimo nedopustni osebek.
else
y=int1+1;
while y