beamgsc_forgraduate.m

在此提供波束成型與方向性的指標天線供研究者參考

%% copy from SPCOM
%% method by GSC
clear all;
close all;
for  run_times=1:1
    run_times
    
    
    D_A=[0];									    % DOA of desired signal
    I_A=[20 50 -35 -20 -45 60 30 80 -75 -55];			     	% DOA of interference
  
    
    N_D=length(D_A);								% number of desired signal
    N_I=length(I_A);								% number of interference signal
    P_ant=12;									    % number of antenna
    pow_I=20;									    % power of interference (db)
    pow_N=0;								       	% power of noise (db)
    SNR=0;                                          % SNR (db)
    R_A=1;								         	% the number of reference antenna
    order=1;
    
    symbols=500*3;                % length of one symbol for each frame
    frame=1;                      % length of frame for each simulation, don't to large
    ini_ad=0;%0.01;               % initinal value for adaptive
    step_gsc=0.00001;
    
    uncor_gsc=zeros(1,1);
    uncor_lcmv=zeros(1,1);
    uncor_mmse=zeros(1,1);
    
    conv_gsc=zeros(symbols,frame);
    conv_gsc_nar=zeros(symbols,frame);
    
    w_gsc=zeros(P_ant,1);
    w_nar_lcmv=zeros(P_ant,1);
    w_nar_mmse=zeros(P_ant,1);
    
    %%% caculated the amplitude
    amp_I=sqrt(10^(pow_I/10));
    amp_N=sqrt(10^(pow_N/10));
    amp_D=sqrt(10.^(SNR/10))*amp_N;
    
    for ii=1:P_ant
        steer_D(ii,:)=exp(-j*(ii-R_A)*pi*sin(D_A*pi/180));
        steer_I(ii,:)=exp(-j*(ii-R_A)*pi*sin(I_A*pi/180));
    end
    
    k0=fix(P_ant/2)+1;
    k=(1:P_ant).'-k0*ones(P_ant,1);
    
    C=[];
    for ii=1:(order+1)
        C=[C (k.^(ii-1))];
    end % "for ii=1:(order+1)"
    lon_B=length(C(1,:));
    fnar_const=[1 zeros(1,lon_B-1)]';
    [u,d,v]=svd(C);
    B=u(:,lon_B+1:P_ant);                 % generate the blocking matrix
    wq=C*inv(C'*C)*fnar_const;       % wq