skyplot
구문
설명
skyplot(는 행렬로 지정된 방위각과 고도 데이터(단위: 도)를 사용하여 스카이 플롯을 만듭니다. 방위각은 북쪽 방향에서 시계 방향으로 양의 각도로 측정됩니다. 고도각은 수평선으로부터 90도 수직선 위에서 측정됩니다. 스카이 플롯 Figure 요소의 세부 정보는 스카이 플롯 기본 요소 항목을 참조하십시오.azdata,eldata)
skyplot(는 방위각 데이터와 고도 데이터를 필드 status)SatelliteAzimuth와 SatelliteElevation을 갖는 구조체로 지정합니다.
skyplot(___,는 위에 열거된 구문의 입력 인수 외에 하나 이상의 이름-값 인수를 사용하여 옵션을 지정합니다. 이름-값 인수는 Name,Value)SkyPlotChart 객체의 속성입니다. 속성의 목록을 보려면 SkyPlotChart Properties 항목을 참조하십시오.
skyplot(는 parent,___)parent에 의해 지정된 Figure, 패널 또는 탭에 스카이 플롯을 만듭니다.
h = skyplot(___)SkyPlotChart 객체인 h로 반환합니다. 차트를 만든 후 속성을 수정하려면 h를 사용하십시오. 속성의 목록을 보려면 SkyPlotChart Properties 항목을 참조하십시오.
예제
GNSS 센서 모델을 gnssSensor System object™로 만듭니다.
gnss = gnssSensor;
센서의 위치와 속도를 지정합니다. 센서 측정값을 시뮬레이션하고 가시 위성에서 상태를 가져옵니다. 방위각과 고도각을 벡터로 저장합니다.
pos = [0 0 0]; vel = [0 0 0]; [~, ~, status] = gnss(pos, vel); satAz = status.SatelliteAzimuth; satEl = status.SatelliteElevation;
로컬 환경을 추가하기 위해 최대 고도가 30도인 무작위 로컬 고도 마스크를 만듭니다.
rng(8) terrainMaskElevations = 30*rand(1,12); % elevations (degrees) terrainMaskEdges = [0 24 48 100 132 180 204 240 276 288 300 312 360]; % azimuth edges (degrees)
위성 위치를 고도 마스크와 함께 플로팅합니다.
skyplot(satAz,satEl,MaskElevation=terrainMaskElevations,MaskAzimuthEdges=terrainMaskEdges);
GNSS 센서로부터 시간 경과에 따른 위성 위치의 궤적을 애니메이션으로 표시합니다.
스카이 플롯 Figure를 초기화합니다. 관련 시간 스텝 정보를 지정합니다.
skyplotHandle = skyplot(0,0);
numHours = 12; dt = 100; numSeconds = numHours * 60 * 60; numSimSteps = numSeconds/dt;
GNSS 센서 모델을 gnssSensor System object™로 만듭니다.
gnss = gnssSensor('SampleRate', 1/dt); 시간 스텝을 반복하고 다음을 수행합니다.
센서 측정값을 시뮬레이션합니다. 정지 상태인 센서의 0 위치와 속도를 지정합니다.
방위각과 고도각을 벡터로 저장합니다.
SkyPlotChart핸들의AzimuthData속성과ElevationData속성을 직접 설정합니다.
for i = 1:numSimSteps [~, ~, status] = gnss([0 0 0],[0 0 0]); satAz = status.SatelliteAzimuth; satEl = status.SatelliteElevation; set(skyplotHandle,'AzimuthData',satAz,'ElevationData',satEl); drawnow end
Adafruit® GPS 위성 센서에서 생성된 로그 파일로부터 방위각 데이터와 고도 데이터를 불러옵니다. 이 예제에서 제공하는 데이터는 각 위성의 방위각과 고도, 의사 랜덤 잡음(PRN) 코드를 포함합니다. 이러한 값을 벡터로 저장합니다.
load('gpsHWInfo','hwInfo') satAz = hwInfo.SatelliteAzimuths; satEl = hwInfo.SatelliteElevations; prn = hwInfo.SatellitePRNs;
PRN 코드를 기준으로 위성을 구분합니다. 각 위치마다 그룹을 연결하기 위해 categorical형 배열을 만듭니다. 이 위성 집합의 경우 PRN이 32 미만인 위성만 측위 결정에 사용됩니다.
isUnused = (prn > 32); group = categorical(isUnused,[false true],["Used in Positioning Solution" "Unused"]);
위성을 시각화하고 GroupData 이름-값 인수에 이 범주 그룹을 지정합니다. PRN을 각 점에 대한 레이블로 지정합니다. 범례를 표시합니다.
skyplot(satAz,satEl,prn,GroupData=group) legend('Used','Unused')
수신기 위치, RINEX 내비게이션 파일, 마스크 각도, 시간 스텝 크기, RINEX 파일에서 샘플링할 데이터의 시간을 지정합니다.
recPos = [42 -71 50]; navfile = "GODS00USA_R_20211750000_01D_GN.rnx"; maskAngle = 25; dt = 60; % seconds numHours = 4;
내비게이션 파일을 읽어서 이 파일에 캡처된 모든 위성의 GPS 데이터를 가져옵니다.
data = rinexread(navfile); [~,satIdx] = unique(data.GPS.SatelliteID); navmsg = data.GPS(satIdx,:);
시작 시간을 내비게이션 메시지의 초기 시간으로 설정합니다. 그런 다음 시간 벡터 t를 만듭니다.
startTime = navmsg.Time(1); secondsPerHour = 3600; timeElapsed = 0:dt:(secondsPerHour*numHours); t = startTime + seconds(timeElapsed);
방위각과 고도에 대한 벡터를 초기화합니다. 그런 다음 모든 위성에 대해 시간 t에서 방위각 데이터와 고도 데이터를 수집합니다.
numSats = numel(navmsg.SatelliteID); allAz = NaN(numel(t),numSats); allEl = allAz; for idx = 1:numel(t) [satPos,~,satID] = gnssconstellation(t(idx),RINEXData=navmsg); [az,el,vis] = lookangles(recPos,satPos,maskAngle); allAz(idx,:) = az; allEl(idx,:) = el; end
고도가 0 미만인 수평선 아래의 모든 위성을 누락값으로 표시합니다.
allEl(allEl < 0) = missing;
스카이 플롯을 만들고 AzimuthData 속성과 ElevationData 속성을 업데이트하여 위성 궤적을 애니메이션으로 표시합니다.
figure sp = skyplot(allAz(1,:),allEl(1,:),satID,MaskElevation=maskAngle); for idx = 1:size(allAz, 1) set(sp,AzimuthData=allAz(1:idx,:),ElevationData=allEl(1:idx,:)); drawnow limitrate end
