%
%
%10.4   Multibandantennen %<'FontWeight',MSVTITLWEIGHT,'FontSize',MSVTITLSIZE,'Color',MSVTITLCOLOR>
%
%
%
% Bei Multibandantennen spielt der Gewinn eine untergeordnete Rolle gegenber der 
% Anpassbarkeit auf vielen Bndern.  Die Kostenfunktion knnte z.B. die mittlere 
% quadratische Abweichung der Impedanz des Antennenspeisepunktes von einem 
% vorgegebenen Wellenwiderstand  {\fontname{Courier}z_0}  bei Frequnenzen  {\fontname{Courier}f}  sein.  Die Formulierung 
% in Matlab ist dann
%
% {\fontname{Courier}costfcn = 'sum((imp(antenne,f) - z0).^2)'}
%
% Wir versuchen die Lngen der beiden Enden der FD4 dahingehend zu optimieren, dass 
% das Stehwellenverhltnis im Sinne des quadratischen Mittels an den Bandgrenzen  
% {\fontname{Courier}[3.5  3.8  7.0  7.1] MHz} fr eine Speisung mit  {\fontname{Courier}150 \Omega}  am kleinsten ist:
%
costfcn = 'sum(swr(imp(fd4(0,x(1),x(2)),[3.5  3.8  7.0  7.1]),x(3)).^2)';
x0 = [13  27; 15  29; 0.2  0.2];
[cm,ppo,c] = gridopt(costfcn,[0 0 150],[1 2],x0,0,'0',1,1);
imp(fd4(0,ppo(1),ppo(2)),3:0.05:8,150,[1.5 2 3 5]);
%
% Offenbar ist die FD4 schon gut optimiert, so dass dieser Versuch keine nennenswerte 
% Verbesserung bringt.  Dennoch machen wir einen weiteren Versuch, um zu zeigen, dass 
% die Kostenfunktion auch anders formuliert werden kann.  Wir wnschen nmlich, dass 
% das maximale Stehwellenverhltnis an den Grenzen smtlicher mit der FD4 nutzbaren 
% Bnder am kleinsten ist:
%
costfcn = 'min(swr(imp(fd4(0,x(1),x(2)),[3.5 3.8 7.0 7.1 14.0 14.35 18.068 18.168 24.89 24.99 28.0 29.7]),x(3)))';
x0 = [13.9  27.8; 0.2  0.2];
[cm,ppo] = randopt(costfcn,[0 0 150],[1 2],x0,150);
imp(fd4,3:0.1:30,150,[1.5 2 3 5]);
imp(fd4(0,ppo(1),ppo(2)),3:0.1:30,150,[1.5 2 3 5]);
%
% Das gut gemeinte Vorhaben hat allerdings zur Folge, dass versucht wird, die sowieso 
% schlecht passenden Bandgrenzen des 10m-Bandes zu verbessern ohne Rcksicht auf die 
% anderen Bnder.  Die Bandgrenzen ersetzen wir deshalb durch die Bandmitten:
%
costfcn = 'min(swr(imp(fd4(0,x(1),x(2)),[3.65 7.05 14.175 18.118 24.94 28.85]),x(3)))';
x0 = [13.9  27.8; 0.2  0.2];
[cm,ppo] = randopt(costfcn,[0 0 150],[1 2],x0,150);
imp(fd4,3:0.1:30,150,[1.5 2 3 5]);
imp(fd4(0,ppo(1),ppo(2)),3:0.1:30,150,[1.5 2 3 5]);
%
% Man erkennt auch hier, dass der Betrieb auf  {\fontname{Courier}18 MHz}  verbessert wird, dies aber auf 
% Kosten der anderen Bnder.
%