function ges = gr12yagi7(az,el,z)
% ges = gr12yagi7(az,el,z)
% Gruppe von yagi7 mit minimalen Nebenzipfeln und 22.6 dBi
% az, el :  Azimut und Elevation in Grad
% z      :  vertikale Positionen der Yagis:
% z(1)   :  Abstand der oberen und unteren inneren Yagis von der Elevationsachse
% z(2)   :  Abstand der unteren inneren Yagis von der Elevationsachse
% z(3)   :  Abstand der oberen  inneren Yagis von der Elevationsachse
% z(4)   :  Abstand der unteren ueren Yagis von der Elevationsachse
% z(5)   :  Abstand der mittleren ueren Yagis von der Elevationsachse
% z(6)   :  Abstand des unteren Querrohrs von der Elevationsachse
%           Default:  z = [3.75 0.92 1.64 2.05 0.50 1.40]
% !!! Kreuzschellen-Versatz in x nicht modelliert (alles auf x = 0) !!!
% !!! Seil-Abspannung nicht modelliert !!!

% Klaus von der Heide DJ5HG,  v.d.heide@on-line.de

if nargin<1,  az = 0; end
if nargin<2,  el = 1; end
if nargin<3,   z = [3.75 0.92 1.64 2.05 0.50 1.40]; end

hg  = 10.00;  % Masthhe
b   = 1.05;   % Abstand der inneren Trger vom Mast
c   = 3.00;   % Abstand der ueren Trger vom Mast
d   = 0.045;  % Abstand Yagi-Boom vom Trger (Kreuzschelle)
uz  = 0.050;  % Trger-berstand ber Yagi-Boom
uy  = 0.100;  % Rohr-berstand (horizontal)
b1  = b + d;  % untere Yagis auen am Trger fr Freiheit zum Gittermast
b2  = b - d;  % obere Yagis innen am Trger fr besseres Richtdiagramm
c1  = c + d;  % Abstand der ueren Yagis vom Mast

yag = group(yagi7,[...
                   -b1 -z(1) ;  +b1 -z(1);  % ganz unten innen
                   -b2 +z(1) ;  +b2 +z(1);  % ganz oben innen
                   -c1 -z(4) ;  +c1 -z(4);  % unten auen
                   -c1 +z(5) ;  +c1 +z(5);  % oben auen 
                   -b1 -z(2) ;  +b1 -z(2);  % unten innen 
                   -b2 +z(3) ;  +b2 +z(3)   % oben innen
                  ]);

tub = wires([0 -(c+uy) 0; 0 -c -(z(4)+uz); 0 +c -(z(4)+uz); 0 -b -(z(1)+uz); 0 +b -(z(1)+uz); 0 +b1 -z(6)],...
            [0 +(c+uy) 0; 0 -c +(z(5)+uz); 0 +c +(z(5)+uz); 0 -b +(z(1)+uz); 0 +b +(z(1)+uz); 0 -b1 -z(6)],...
            [    0.05            0.045            0.045            0.05             0.05        0.04],...
            [      31              19               19               31               31          19]);
kon = construct(yag,tub);    % vertikal drehbares System
dre = rotation(translation(rotation(kon,cos(pi/180*[el 90 90+el;90 0 90; 90-el 90 el])),[0 0 hg]),[0 0 az]);   % Drehung und Anhebung
mst = wire([0.30 0 0],[0.30 0 hg],0.15,81);
ges = comment(construct(dre,mst),{'84-Element-Gruppe';'144MHz'}); 
