function [xs,ys,zs] = bd_sat(t)
%*********************************
% ** 北斗星座 (地心地固坐標) 
% **********************************
%模擬衛星在地心地固坐標系中的運動，參數t是運行時間，單位秒（s），距離的單位（km）
we = 0;%7.2921151467e-5*360/(2*pi);        %0;%為了讓坐標不隨地球轉動，在此讓地球轉速we=0
trans=2*pi/360;      %用tran將度換算弧度
bomigaG0=0;      %參考時間的格林威治子午線赤經
%-------------GEO--------
a=42378137;      %a是軌道半徑
e=0.0;           %e是軌道離心率
i=0;             %i是軌道傾角
bomiga=0;%每個星的升交點赤經
w=0.0;           %w近地點角，在此取0
v0=[80,140];    %每個星的真近點角初始值
v=we*t+v0; 
bomigaG=(bomigaG0+we*t); %因為地球的自轉，坐標系也有一個旋轉
r=a.*(1-e.^2)./(1+e.*cos(trans*v));             %公式
xs(1:2)=r.*cos(trans*(v+w)).*cos(trans*(bomiga-bomigaG))-r.*sin(trans*(v+w)).*sin(trans*(bomiga-bomigaG)).*cos(trans*i);        %在地心地固坐標系下的軌道
ys(1:2)=r.*cos(trans*(v+w)).*sin(trans*(bomiga-bomigaG))+r.*sin(trans*(v+w)).*cos(trans*(bomiga-bomigaG)).*cos(trans*i);        %軌道方程
zs(1:2)=r.*sin(trans*(v+w)).*sin(trans*i);