%此函數用于產生一部分極化信號的B。2個角度分別是仰俯角、方位角
function a=v_a(angle_1,angle_2)
degrad=pi/180;
i=sqrt(-1);
j=i;
a=[cos(angle_2*degrad)*cos(angle_1*degrad) -sin(angle_2*degrad)
   sin(angle_2*degrad)*cos(angle_1*degrad)  cos(angle_2*degrad)
  -sin(angle_1*degrad)                      0
  -sin(angle_2*degrad)                     -cos(angle_2*degrad)*cos(angle_1*degrad)
   cos(angle_2*degrad)                     -sin(angle_2*degrad)*cos(angle_1*degrad)
   0                                        sin(angle_1*degrad)];