sa = spiralArchimedean('Turns',9, 'InnerRadius',1e-3, ...
    'OuterRadius',11e-3);
CBSA = cavityCircular('Exciter',sa,'Radius',12e-3,'Height',...
    25e-3,'Spacing',23e-3);
    N2D=7;
    p(:,1) = [  0.0000;   0.0000; 0]; 
    p(:,2) = [ 25;   7.5; 0]; 
    p(:,3) = [  7.5;  25; 0]; 
    p(:,4) = [-19;  17.5; 0]; 
    p(:,5) = [-26;  -8.3; 0]; 
    p(:,6) = [ -8.3; -26; 0]; 
    p(:,7) = [ 17.5; -19; 0]; 
        p = p*1e-3;
        a=zeros(1,N2D);e=90*ones(1,N2D);
        DA03 = phased.ConformalArray('Element',CBSA,'ElementPosition',p,...
                                                        'ElementNormal',[a;e]);                
    
        fc = 10*1e9;
        Tc=1/fc;
        df = 10e6;    
        WL = physconst('LightSpeed')/fc; 
        t = linspace(0,2*Tc,400)';
        dt=t(2);      
            SigAz = 10;
            SigEl = 20;            
    
            S1 = 10*cos(2*pi*fc*t);
            S1=awgn(S1,20,'measured');
    
            collector = phased.Collector('Sensor',DA03,'OperatingFrequency',fc,...
                                    'Wavefront','Plane','Polarization','none');
      yDa = collector(S1,[SigAz;SigEl]);
            svDA = phased.SteeringVector('SensorArray',DA03,...
                                                'IncludeElementResponse',true);   
            A=svDA(fc,[SigAz;SigEl]);