%?Code?function:?the?performance?analysis?of?diversely?polarization?array
%?Code?writer:?Borne?Zhang
%?Date:?2012-04-23
%?scalar-sensor?array‘s?steering?vector?matrix:?As
%?vector-sensor?array‘s?steering?vector?matrix:?A

%-----clc?&?clear?all-----%
clc;
clear?all;

%-----the?partial?differentiation?of?As-----%
N=6;
M=2;
lamda=1;
d=0.25*lamda;
syms?theta1?theta2?gama1?gama2;
phi=90*（pi/180）;
eta=0*（pi/180）;
as_1=exp（1i*2*pi*d*（0:（N-1））.‘*sin（theta1）*sin（phi）/lamda）;
as_2=exp（1i*2*pi*d*（0:（N-1））.‘*sin（theta2）*sin（phi）/lamda）;
As=[as_1as_2];
As_theta_1=diff（Astheta1）;
As_theta_2=diff（Astheta2）;
As_theta=As_theta_1+As_theta_2;

%-----the?partial?differentiation?of?A------%
Lh=diag（[cos（gama1）cos（gama2）]）;
Lv=diag（[sin（gama1）sin（gama2）]）;
Qp=1i*diag（[1-11-11-1]）;
A=As*Lh+Qp*As*Lv;
A_theta=diff（Atheta1）+diff（Atheta2）;
A_gama=diff（Agama1）+diff（Agama2）;

%-----CRB?for?scalar-sensor?array-----%
K=1;
Pn=1;
SNR=0:5:20;
CRB_scalar=zeros（1length（SNR））;
CRB_vector=zeros（1length（SNR））;
for?s=1:length（SNR）
Ps=Pn*10^（0.1*SNR（s））;
Rss=Ps*eye（MM）;
As=subs（As{theta1theta2}{015*（pi/180）}）;
Rxx=As*Rss*As‘+Pn*eye（NN）;
As_theta=subs（As_the