moon_scan

Moon Scan

0001 %Moon Scan
0002 close all
0003 
0004 %d = read_arc('26-jan-2010:05:57:00','26-jan-2010:06:12:00'); %moon_scan
0005 %d = read_arc('26-jan-2010:06:25:00','26-jan-2010:06:34:00'); %taua_scan
0006 d = read_arc('26-jan-2010:19:05:00','26-jan-2010:19:12:00'); %casa_scan
0007 
0008 
0009 d = reshapeRx(d);
0010 d = reshapeAz(d);
0011 d = interpAz(d);
0012 
0013 %sourceScanMap(d, 0.25)
0014 
0015 
0016 
0017 figure;
0018 plot(d.antenna0.receiver.utc, d.antenna0.receiver.dataNew(:,1))
0019 % now we calculate which azimuth and elevation the source should be at at
0020 % every time stamp in d.
0021 time = d.antenna0.receiver.utc;
0022 az_scan = d.antenna0.servo.interpAz;
0023 el_scan = d.antenna0.servo.interpEl;
0024 [az_src, el_src] = calcAzEl(d);
0025 az_poff = -21.75;
0026 el_poff = -1.65;
0027 
0028 az_scan_mended=az_scan-az_poff;
0029 el_scan_mended=el_scan-el_poff;
0030 rx=d.antenna0.receiver.dataNew;
0031 %plot(az_scan_mended,rx(:,1));
0032 %title('Azimuth Scan and Rx Power');
0033 %figure
0034 %plot(el_scan_mended,rx(:,1));
0035 %title('Elevation Scan and Rx Power');
0036 
0037 
0038 %plot3(az_scan_mended,el_scan_mended,d.antenna0.receiver.dataNew(:,1))
0039 
0040 deltaAz = az_scan - az_src-az_poff;
0041 deltaEl = el_scan - el_src-el_poff;
0042 plot3(deltaAz,deltaEl,d.antenna0.receiver.dataNew(:,1))