10 VOR信号生成

1、VOR信号生成 2015.2 VOR的一般特性的一般特性 1. 工作频段工作频段 f =108.00117.95MHz， f =0.05MHz，共，共200个频道，个频道，VOR只用其只用其 中的中的160个频道。个频道。 2. 频道分配（频道分配（1） 200个频道个频道 (108.oo117.95MHz) VOR：160个频道个频道 (108.00117.95MHz) LOC：40个频道个频道 (108.10111.95MHz) 航路航路VOR：120个频道个频道 (112.00117.95MHz) 终端终端VOR：40个频道个频道 (108.00111.85MHz) ICAO建议，对于航

2、路建议，对于航路VOR，应优先选用，应优先选用112112117.95117.95MHz范围内范围内 小数点后第一位为奇数的频率，然后考虑该频率范围内小数点后第一位小数点后第一位为奇数的频率，然后考虑该频率范围内小数点后第一位 为偶数的频率。为偶数的频率。 中国民航大学 CAUC VOR在导航中在导航中的应用的应用 N N VOR-A VOR-B 60 90 30 中国民航大学 CAUC VOR系统的基本原理系统的基本原理 VOR信标辐射信号信标辐射信号 可变相位信号可变相位信号 基准相位信号基准相位信号 识别信号识别信号 话音信号话音信号 中国民航大学 CAUC 基准相位信号的发射 先用30

3、Hz对9960Hz副载波调频，然后调频副载波再对载波调幅。 而30Hz调频信号的相位在 VOR台周围 360方位上是相同的。 CVOR 基准相位信号波形 可变相位信号的发射 用30Hz对载波调幅，相位随VOR台的径向方位而变化 合成辐射场为旋转的心形方向性图 中国民航大学 CAUC V OR空间合成信号波形图 中国民航大学 CAUC 基准和可变 30HZ的相位关系 设当心形方向性图 的最大值在磁北时， 9960HZ副载波的 频率为10440HZ。 北 10440HZ 东 9960HZ 南 9480HZ 西 9960HZ 3.5 3.5 DVOR信标的工作原理信标的工作原理 图片说明图片说明：D

4、VOR常与测距机（常与测距机（DME）信标合装在同一台站。图中在圆周上布）信标合装在同一台站。图中在圆周上布 置的天线为置的天线为DVOR天线，旁边直立的有一个是天线，旁边直立的有一个是DME信标的天线。信标的天线。 仿真波形生成思路 利用DDS（查表法）生成正弦波 先调频 再调幅 DDS：数字直接频率合成 具体仿真参数 采样速率fs=48000S/s 查表法实现 余弦表cos_tab：01023 相位累加器位宽：20bit(整数部分10bit，小数 部分10bit） 20bit时对应频率分辨率=fs/220=1/22Hz 查表前右移位数：20-10=10bit DDS工作原理 计算数字频率（

5、相位增量）：Phase_inc 累加相位增量Phase_inc得到当前相位：Phase_acc 根据波形表格的大小对当前相位进行移位：右移20-10=10bit 查表得到波形 举例：fs=48000Hz，acc_width=20，则： fc=1000Hz,则Phase_inc=220*1000/48000=21845 Matlab程序说明(参数设置) %生成1024点余弦表 cos_tab=round(32767*cos(2*pi*(0:1023)/1024); %主要参数设置 fs=48000; %采样频率 Sample/s phase_inc_1Hz=(220/fs); %1Hz对应的数字

6、频率（相位累加器宽度=20bit） shift_r=10; %查表前右移位数 fc=1000; %载波频率（Hz） f_30=30; %调制信号频率（Hz） fc_phase_dlt=round(fc*phase_inc_1Hz); %1kHz载波对应的数字频率 phase_dlt_30Hz=round(f_30*phase_inc_1Hz); %30Hz调制信号对应的数字频率 phase_FM=0; %FM调频信号(载波)瞬时频率 phase_30Hz_am=0; phase_30Hz_fm=0; %FM parameter MaxFre=9960; MinFre=9480; FREQDEV

7、=(MaxFre-MinFre)/2; phase_fm_delta=round(FREQDEV*phase_inc_1Hz); %FM频偏对应的数字频率 Matlab程序说明-生成FM/AM信号 for sampNo=1:round(fs*1/30) %1/30Hz %gen 30Hz_fm phase_30Hz_fm=phase_30Hz_fm+phase_dlt_30Hz; tmp1=round(phase_30Hz_fm/2shift_r); tmp=mod(tmp1,1024); y_30Hz_fm=cos_tab(tmp+1); y_fm(sampNo)=y_30Hz_fm; %g

8、en FM signal tmp1=(y_30Hz_fm/32767)*phase_fm_delta; phase_FM=phase_FM+fc_phase_dlt+tmp1; tmp1=round(phase_FM/2shift_r); tmp=mod(tmp1,1024); y_carry=cos_tab(tmp+1); %gen 30H_am phase_30Hz_am=phase_30Hz_am+phase_dlt_30Hz; tmp1=round(phase_30Hz_am/2shift_r); tmp=mod(tmp1,1024); y_30Hz_am=cos_tab(tmp+1)

9、; %gen am y_fc(sampNo)=round(y_30Hz_am*0.8/2+16384)*y_carry)/215); end Matlab程序说明(运行结果) 02004006008001000120014001600 -4 -3 -2 -1 0 1 2 3 4 x 10 4 plot(y_fc,b+:); hold on; plot(y_fm,r*:); TMS320C5509算法说明-主程序 int BlockCnt=0; long DMARec_p; void main(void) initial_sys(); while(1) if(DataGeted=1) DataG

10、eted=0; *IODATA=(*IODATA)/B7=0, GPIO7=0, 23J0=1, for test CPU use ratio BlinkLed(BlockCnt/100)%2); BlockCnt+; Deal_data(); *IODATA=(*IODATA)|(116; /通道目的起始地址Upper 数据地址为0 x20000 *DMA_CSSA_L2 =(unsigned int)(long)BufferB*2); /通道目的起始地址Low DMARec_p = BufferA; else *DMA_CSSA_U2 =(long)BufferA*2)16; /通道目的起

11、始地址Upper 数据地址为0 x20000 *DMA_CSSA_L2 =(unsigned int)(long)BufferA*2); /通道目的起始地址Low DMARec_p = BufferB; *DMA_CCR2 =0 x1986;/通道使能，不重复，利用McBSP1发送事件进行同步 LBDS1=uLED(numTransfers_Tx5) 问题： 1、该函数的中断周期=？ 2、DataGeted变量的作用？ 波形生成1 for(i=0;i15); /generate cos(2*pi*fc*t) phase_fc+=fc_phase_dlt+dat3; dat1=Get_cos(p

12、hase_fc);/FM /generate AM dat2=(int)(long)dat2*102)8)+16384;/m=0.8, 0.4*28=102.4 dat1=(int)(long)dat1*(long)dat2)16);/AM */ Buffer1rec_blk_cnt*240+i=dat1; *(int *)(DMARec_p+i*2+0)=dat1; *(int *)(DMARec_p+i*2+1)=dat1; 波形生成2 int Get_cos(unsigned long phase_curren) int dat1,dat2; unsigned long tmp1; tmp1=(phase_currenshift_r);/ shift_r=10 dat1=(int)tmp1)/mod 1024 dat2=cos_tabdat1;/look up table, get cos r