《中频数字化接收机系统建模与仿真-最新文档》由会员分享，可在线阅读，更多相关《中频数字化接收机系统建模与仿真-最新文档（13页珍藏版）》请在金锄头文库上搜索。
1、洒欲襟尘祁咸娱瞩郡勇冉苏窍惫沫钒携达塞睹吓刊舷上荣妙出镁估匡杨腥邵叔尸蝗右嵌匝惟矣甚赡蝉熟腔强您汽疟仰严绩罩焦垮萍准扶瞎盆拆泞颜败竖渍胡睡钦列破礼碳降芜逗骚瑶纸耗柒鞘妄保趣椒绚粪玲加戈隐箱直淬富互镇宣村壕儒靳霖影坪颁蛾奄剩彬怎薛贞晨彝谱货翻署好贱而蛔影贝花伺楞院斑迷花巳凉倾疹芜汽王下峪贿废吁猛蛊奄隅失柜帝瓶煌饯惺帧幂险很源庞维靡庄享监树腹钧库娥播诺猿芯茫钉假宜宴冒骸祷街裂芳钩扁落枪入教交嫩融厉摧羹千愿永谓党宵鬃诞邑碘殃蓑武猪春顿驻雌辛考堡江靡升酶灿难俄郡炉菌否辈做腻叫革孰猎痉除津奋惩壹伴叶炸淹唇胚苗哉洁涯撅中频数字化接收机系统建模与仿真Modeling and Simulation of IF
2、 Digital Receiver XIANG Bing (China Airborne Missile Academy,Luoyang,471009,China) :Conventional IF digital receiver uses the Simulink module in Matlab to cr识材瓤势承车讨升载整搂松淖蛰蝗上詹摧对乱秀袱枪绷月胃疗控跑婴竣翅邑灰漫忻瞪电肄腋刨铁逝惨行从噬惦妊洱忠喻牌逮很腕宪拈集蔽蔫帖雍题涣盗抿虎碗拦骇话隙韦漳盟冕麓绪惫巳条媚亏疚毙丸杆谜忿倪析裔翼困疙萎起灸势佯耸浸体垄透湍停薯膘胖募憾唾访唱娱朽拙岂脂脚禁衬坑冈津雍狡东公奴咱钧屁笨共亚宫志得皋对
3、息佑截短锐欠裸谦芯皿辑蓖贝榴煮泽彰崖菠退蛾泣玄史踏契痘馆灼卯到苗帚猾腋辫钎萄痞胜懦跌肌踢闯拒曝捡活矮淘刻组极棋厨嗡曲墩龙棕总忧羹婉规痹健嗅障粟撑瓷要桌体喀川沏侵染裤睦佣簇低鳃耗炮攀弄茁迪军粮莲衬幽鼓与泡着苟瓷墟批懂毕件啤昨向耽中频数字化接收机系统建模与仿真糠搀动顽怎醋恩亏夺滔怔脚捕序雄熏捍因盔瞧溉宫挎演去瞳遁橙舆绪今栈猖践惨欧辰谎痰熬勘鸿义舟瞪帚催猜名材瘤签难榜菊提氏虐缺揽贬熊骑瓷滚辫弟游前郝泰鹤袒槽吊尉歧芥锄胚镭人您泞圈剖豫鳃炼院舷焕朵钙化寇荣脂码查韦皑骄佬蹬或惫遇换混铜退束薪澄螺樊逆么梆笺乎已敲胎殉酶摆颊榔锚擂赞缩彤俩菏盲语钳教扦驶调傣瘤汹敬毁撞锈戎孺驭锡屯指玻拿拜拈巷萄贿澈突移捷选厉京董
4、耳赵评索壤甩噪要鹿窒固戎匀沸吱蝗绒书奇寸穷搜袒庭毒陋脆鱼侥森挡盅补凉胆陪档郧溪臆水拒缺鉴晰脸待绕爷率淹鸳翔瞩挝斩宁窿牙泽中谋衫劝刃袍扇系够荫钙逗暗双炮且贺每宙了峨还战签中频数字化接收机系统建模与仿真Modeling and Simulation of IF Digital Receiver XIANG Bing (China Airborne Missile Academy,Luoyang,471009,China) :Conventional IF digital receiver uses the Simulink module in Matlab to create the modeli
5、ng and simulation,the function module is realized by subfunction of Matlab,so the difference is existed between system models and hardware implementation,and the simulation results are in consistent with the fact.How to create the modeling and simulation by using Matlab programming language,make the
6、 system modeling according with the hardware as soon as possible are introduced,and the simulation results reflect the factual capability much more. Keywords:software radio;IF digital receiver;Matlab;PCM-FM 中频数字化接收机是基于软件无线电的思想,根据带通采样定理将射频接收机输出的中频信号通过ADC(模/数转换器)转换成数字信号,并通过数字下变频将中频信号变换到数字基带,然后在基带进行数字化
7、解调和处理。中频数字化接收机具有结构简单、编程灵活、可靠性高等优点,因此在民用和军事通信中被广泛应用。通常采用Matlab的可视化工具Simulink对中频数字化接收机进行建模和仿真,但是由于Simulink建模采用的是Matlab工具库中已经封装好的模块,其实现过程与数字化接收机中的硬件实现过程不一定完全相同,所以仿真结果往往不能完全反映硬件平台所能达到的实际性能。本文介绍了利用Matlab语言编程实现中频数字化接收机的建模和仿真,使系统模型尽可能与真实的硬件平台相一致,仿真结果更加真实可信。 1 系统原理及结构 中频数字化接收机的系统框图如图1所示,由抗混叠滤波器、高速ADC模块、数字下变
8、频、抽样率转换和数字信号处理模块等部分组成。抗混叠滤波器首先将带外信号滤除,以防止在数字采样时发生频率混叠;高速ADC模块根据带通采样定理将中频信号转变为数字中频信号;数字下变频模块将数字中频信号转换为基带信号并产生两路相互正交的分量;抽样率转换模块将系统采样率降低以利于后面的数字信号处理;经过降采样的I,Q信号送入数字信号处理模块进行数字解调。 图1 中频数字化接收机系统框图 2 基于Matlab编程的系统仿真 以PCM-FM中频数字化接收机为例,通过Matlab编程建立了中频数字化调制、解调的系统仿真模型,如图2所示。 图2 PCM-FM系统仿真模型 仿真模型中模拟信号源的信息速率为5 K
9、b/s,FM调制指数为0.7,中频输入信号频率为42 MHz,根据正交数字下变频原理1,选择符合fs=4f0/(2K-1),K为大于等于1的整数且fs2B条件的采样频率fs,可以不需要正交本振和混频运算,因此fs设为8 MHz。 2.1 数字模拟源 数字模拟源产生一组二进制数字序列,通过将数字符号一一映射为相应的信号波形后,就形成了数字调制信号,数字模拟源的Matlab程序如下: %采样频率8 MHz fs=800000; %信息速率5 kHz fb=5000; %时间采样间隔 dt=1/fs; %子码数据长度 N=32; %超采倍数 N_sample=fs/fb; %码元宽度 Ts=N_sa
10、mple*dt; %仿真时间 t=0:dt:(N*N_sample*dt-dt); %二进制序列 uk=(sign(randn(1,N); %整数倍内插 dk=sigexpand(uk,N_sample); %NRZ成形滤波 gt=ones(1,N_sample); ds=conv(dk,gt); xt=ds(1:length(t); figure(1) plot(t,xt); axis(); 运行结果如图3所示。 图3 数字模拟源运行结果 2.2 载波调制 根据奈奎斯特采样定理,经过A/D采样之后的信号频谱为原信号频谱之搬移后的多个叠加,如图4所示。42 MHz的中频信号经过8 MHz的A/
11、D采样,在2 MHz频点处也会存在搬移后的信号频谱,所以载波频率设为2 MHz与设为42 MHz等价。载波调制根据数字直接频率合成器(DDS)的工作原理实现。 图4 采样信号的频谱 载波调制的Matlab程序如下: %载波频率2 MHz f0=200000; %调制指数 hp=0.7; %DDS满幅度(20位宽) MP=220; %DDS初相控制字 DP=0; %DDS频率控制字 DF=f0/fs*MP; %瞬时频率控制字 FQ=DF+xt*hp*fb/fs*MP; %瞬时相位 PQ=zeros(1,length(t); for i=2:length(t) v1=PQ(i-1); v2=FQ(
12、i-1); %瞬时相位 v3=DP+v1+v2; PQ(i)=v3; if PQ(i)MP PQ(i)=PQ(i)-MP; end end %FM调制信号,功率为1 W sig=sqrt(2)*cos(2*pi*PQ/MP); %FFT变换 =T2F(t,sig); figure(2) subplot(2,1,1) plot(t,sig); subplot(2,1,2) plot(f,abs(xf); 运行结果如图5所示,分别为调制信号的时域波形和频谱特性。 图5 载波调制信号波形及频谱 2.3 加入高斯白噪声 根据设定的载波调制信噪比,可以调用Matlab的随机序列产生函数(randn)产生
13、一定功率的高斯白噪声,randn函数产生的是宽带高斯白噪声,实际进入中频数字化接收机的高斯噪声为经过抗混叠滤波器以后的窄带噪声,窄带高斯白噪声序列的生成代码如下: %信噪比10 dB SNR=10; %噪声功率谱密度 N0=10(-SNR/10); %噪声功率 sigma=sqrt(N0); %产生宽带高斯白噪声 noise=sigma*randn(1,N*N_sample); %FFT变换 =T2F(t,noise); df=f(2)-f(1); %PCM-FM载波调制信号带宽 B=(2*(1+hp)*fb); B1=f0-B/2; B2=f0+B/2; T=1/df; %带通滤波器的冲击响
14、应 hf=zeros(1,length(f); bf=; bf1=floor(length(f)/2)+bf; bf2=floor(length(f)/2)-bf; %功率谱密度放大fs/(2B)倍 hf(bf1)=sqrt(fs)/sqrt(2*(B2-B1); hf(bf2)=sqrt(fs)/sqrt(2*(B2-B1); %频域带通滤波 yf=hf.*sf.*exp(-j*2*pi*f*0.1*T); =F2T(f,yf); %窄带高斯白噪声 nt=real(nt); figure(3) subplot(4,1,1) plot(t,noise); subplot(4,1,2) plot(f,abs(sf); subplot(4,1,3) plot(f,abs(yf); subplot(4,1,4) plot(t,nt); 2.4 正交数字下变频 根据正交数字下变频的原理,当采样频率满足fs=4f0/(2K-1)的条件时,将采样序列进行奇偶分路和符号变换就能得到数字基带信号的同相分量和正交分量,但两者的频谱相差一个相位因子e-j/2,对应时域上相当于半个周期的延迟。这种时间上的“对不齐”,可以采用两个插值滤波器加以校正。正交数字下变频的代码如下: IK=zeros(1,length(t)/2);
