GNURadio探密.docx

1、GNURadio探密GNU Radio 探密2010-03-04 20:30作者：Eric Blossom简介：软件无线电是一门让软件代码尽可能靠近接收天线的技术。它把无线电的硬件问题放到软件中来解决。软件无线电最基本的特征是软件定义（调制）无线电传输的波形，软件解调接收到的波形。这个和目前大多数的由模拟电路或者由模拟数字组合电路组成的无线电有鲜明的区别。GNU Radio是一个开源的可以构建软件无线电平台的软件包。软件无线电是无线电设计领域的一次革命，使得建立无线电设备更加灵活，给用户带来更多的机遇。软件无线电能够很好的实现传统无线电所能实现的功能，另外，软 件无线电最精彩的特性是由软件提供

2、给我们的灵活性，一些年以后，我们将会看到传统的由很多固定元件组成的无线电设备将会被通用的通信设备来取代。可以想 像，一个软件无线电设备能够变体成为蜂窝电话，能够灵活的使用GPRS,802.11 Wi-Fi,802.16 WiMax,卫星链路，或者将出现的其他通讯标准进行通讯，而硬件设备无需改动，只需要软件设置。你能够决定使用GPS，GLONASS或者两个一起使用来给自己定位。也许，最令人激动的是使用软件无线电来建立一个分布式个人通讯系统。当今的通讯系统，大多数是采用自上而下的结构，广播和电视系统提供一个单向链路，内容被 严格控制在一小部分人的手里。蜂窝电话系统给人们带来极大的便利，但是你的手机

3、提供的功能是由运营商来控制的，而不是你自己。这种集中控制的系统限制了人们的创新，代替做蜂窝电话系统的二等公民，我们可以建造一个智能设备，这些设备能够自组织的，在使用者之间建成一个网络。处理流程图图一图一，表示一个典型的软件无线电处理流程图。为了理解无线电的软件模块，首先需要理解和其关联的硬件。在这个图中的接收路径上，能够看到一个天线，一个神奇的RF前端，一个模拟数字转换器ADC和一堆代码。ADC是一个连接连续模拟的自然世界和离散的数字世界的桥梁。ADC有两个主要特性，抽样率和动态范围。抽样率是ADC测量模拟信号的速度，动态范围是ADC区别最低信号值和最大信号值的精度，这决定ADC数字信号输出的

4、比特数（位数）。例如，8位的AD转换器最多能代表256个信号层次，而一个16位的转换器能够代表65536个层次信号。总的来说，ADC的物理特性和价格决定了抽样率和动态范围。在我们深入研究软件之前，先来了解一些理论知识，在1927年，出生于瑞典的物理和电子学家Harry Nyquist提出了如果AD转换想没有混叠现象发生，那么抽样率至少是目标信号带宽的2倍。混叠现象就像是车子重复的碾在过去的车轮印迹上一样，让你分不清楚是本次的还是以前碾过的车轮印迹。假设我们要处理一个低通信号，我们感兴趣的信号带宽是 0到fmax，按照Nyquist理论，抽样率必须至少是2*fmax。如果我们的ADC工作在20M

5、HZ，但是我们想收听92.1MHZ的FM电台，我 们该怎么办呢？答案是使用RF前端，接收机的RF前端能够把它接收到的高频段信号下变频到一个低频段信号后输出。例如，我们能让RF前端把90- 100MHZ频段内的信号下变频到0-10MHZ的低频范围内，那么我们的20MHZ的ADC就能够派上用场了。大多数情况下，我们把RF前端当作一个信号控制的黑盒子，负责处理输入信号的中心频率。举一个具体例子，一个调制解调器的调制模块能够把50M到800M 之间的6M带宽的中心频率下变频到一个中心频率是5.75MHZ的输出。这个输出的中心频率通常称作中频(IF)。按照越简单越易用的原则，RF前端最好也能够被去掉，有

6、一个GNU Radio用户已经成功的使用一个100英尺（译注：等于12英寸，合0.305米）的天线直接连接到一个20M抽样率的ADC上收听到了AM和短波广播。深入软件内部 GNU Radio提供一个信号处理模块的库，并且有个机制可以把单个的处理模块连接在一起形成一个系统。编程者通过建立一个流向图（flow graph）就能搭建成一个无线电系统。信号处理模块是使用C+来实现的，理论上说，信号数据流不停的从输入端口流入从输出端口流出。信号处理块 (blocks)的属性包括输入和输出的端口数，流过它们的数据的类型，经常使用的数据流的类型是短整型(short),浮点型(float),和复数 （comp

7、lex）类型。一些处理模块仅仅有输入端口或者输出端口，它们分别成为信号源(data source)和信号接收器(sink)。有的信号源从文件或者ADC读入数据，信号接收器写入文件或者DAC或者PC的多媒体接口。GNU Radio提供了超过100个信号处理块，并且扩展新的处理模块也是非常容易的。软件图形化接口和信号处理模块的链接机制是通过python脚本语言来进行的，例1，是一个GNUadio “Hello World”的例子。它产生两个sine波形并且把他们输出到声卡，一个输出到声卡的左声道，一个输出到右声道。例子1.输出拨号音#!/usr/bin/env pythonfrom gnuradi

8、o import grfrom gnuradio import audiodef build_graph ():sampling_freq = 48000ampl = 0.1fg = gr.flow_graph ()src0 = gr.sig_source_f (sampling_freq, gr.GR_SIN_WAVE, 350, ampl)src1 = gr.sig_source_f (sampling_freq, gr.GR_SIN_WAVE, 440, ampl)dst = audio.sink (sampling_freq)fg.connect (src0, 0), (dst, 0)

9、fg.connect (src1, 0), (dst, 1)return fgif _name_ = _main_:fg = build_graph ()fg.start ()raw_input (Press Enter to quit: )fg.stop ()我们开始建立一个流向图(flow graph)把所有的信号处理模块连接到一起，调用gr_sig_source_f产生两个sine波形，后缀f表面这个信号源的数据类型是浮点型的，一 个波形是350HZ，另外一个是440HZ，合在一起他们听起来像一个美国电话拨号音。audio_sink是一个接收器，它把接收到的信号输入到声卡 中，我们把3个

10、信号处理模块用流向图的connect方法连接到一起。connect方法有两个参数，源端点和目的端点，用来建立一个从源到目的处理模块 的链路。每个端点(endpoint)有两个成员：一个信号处理模块和一个端口号。端口号表示哪个输入或者输出端口应该被连接的。通常端点使用 python语言的tuple来表示，像：（block,port_number）。当端口号是0时，portt_number可以被省略的。下面的两个表示方法是一样的:fg.connect (src1, 0), (dst, 1)fg.connect (src1, (dst, 1)一旦流向图被建立了，我们调用start生成一个或者多个线程

11、去运行它，按下任意键程序控制权就会返回给调用者。一个完整的FM接收机例2 是一个简单化的但是是一个完整的FM广播接收机，它包括对RF前端的控制和所有FM信号的处理。这个例子的RF前端是使用一个cable调制解调器和一个20M抽样率的ADC来实现的。例2.FM广播接收机#!/usr/bin/env pythonfrom gnuradio import grfrom gnuradio import audiofrom gnuradio import mc4020import sysdef high_speed_adc (fg, input_rate):# return gr.file_source

12、 (gr.sizeof_short, dummy.dat, False)return mc4020.source (input_rate, mc4020.MCC_CH3_EN | mc4020.MCC_ALL_1V)# return a gr.flow_graph#def build_graph (freq1, freq2):input_rate = 20e6cfir_decimation = 125audio_decimation = 5quad_rate = input_rate / cfir_decimationaudio_rate = quad_rate / audio_decimat

13、ionfg = gr.flow_graph ()# use high speed ADC as input sourcesrc = high_speed_adc (fg, input_rate)# compute FIR filter taps for channel selectionchannel_coeffs = gr.firdes.low_pass (1.0, # gaininput_rate, # sampling rate250e3, # low pass cutoff freq8*100e3, # width of trans. bandgr.firdes.WIN_HAMMING

14、)# input: short; output: complexchan_filter1 = gr.freq_xlating_fir_filter_scf (cfir_decimation,channel_coeffs,freq1, # 1st station freqinput_rate)(head1, tail1) = build_pipeline (fg, quad_rate, audio_decimation)# sound card as final sinkaudio_sink = audio.sink (int (audio_rate)# now wire it all toge

15、therfg.connect (src, chan_filter1)fg.connect (chan_filter1, head1)fg.connect (tail1, (audio_sink, 0)return fgdef build_pipeline (fg, quad_rate, audio_decimation):Given a flow_graph, fg, construct a pipelinefor demodulating a broadcast FM signal. Theinput is the downconverted complex basebandsignal.

16、The output is the demodulated audio.build_pipeline returns a two element tuplecontaining the input and output endpoints.fm_demod_gain = 2200.0/32768.0audio_rate = quad_rate / audio_decimationvolume = 1.0# input: complex; output: floatfm_demod = gr.quadrature_demod_cf (volume*fm_demod_gain)# compute

17、FIR filter taps for audio filterwidth_of_transition_band = audio_rate / 32audio_coeffs = gr.firdes.low_pass (1.0, # gainquad_rate, # sampling rateaudio_rate/2 - width_of_transition_band,width_of_transition_band,gr.firdes.WIN_HAMMING)# input: float; output: floataudio_filter = gr.fir_filter_fff (audi

18、o_decimation, audio_coeffs)fg.connect (fm_demod, audio_filter)return (fm_demod, 0), (audio_filter, 0)def main (args):nargs = len (args)if nargs = 1:# get station frequency from command linefreq1 = float (args0) * 1e6else:sys.stderr.write (usage: fm_demod freqn)sys.exit (1)# connect to RF front endrf

19、_front_end = gr.microtune_4937_eval_board ()if not rf_front_end.board_present_p ():raise IOError, RF front end not found# set front end gainrf_front_end.set_AGC (300)# determine the front ends Intermediate FrequencyIF_freq = rf_front_end.get_output_freq () # 5.75e6# Tell the front end to tune to fre

20、q1.# I.e., freq1 is translated down to the IF frequencyrf_front_end.set_RF_freq (freq1)# build the flow graphfg = build_graph (IF_freq, None)fg.start () # fork thread(s) and returnraw_input (Press Enter to quit: )fg.stop ()if _name_ = _main_:main (sys.argv1:)和Hello World的例子一样，这个例子建立了一个流向图，把所有的信号处理模块

21、连接到了一起然后开始运行。在这个例子中，mc4020.source是一个和 Measurement Computing PCI-DAS 4020/12 高速 ADC的连接接口。在它后面连接一个gr.freq_xlating_fir_filter_scf，这是一个有限脉冲响应（FIR）滤波器，负责选择 FM广播的电台频段并且把它变频到基带（0HZ，直流）。然后用20M的converter和解调器，我们就能够接收到邻近的6MHZ带宽的频率了，在这 个频带中大概能包含10个或者更多频道的FM电台。接着使用gr.freq_xlating_fir_filter_scf来选择一个我们想收听的频道。 在这个

22、例子中，我们选择了RF前端精确的中频频率5.75MHZ，gr.freq_xlating_fir_filter_scf的输出是一个经过 160,000次/秒的复采样数据流，然后把复基带信号输出给gr.quadrature_demod_cf进行积分取模，这个模块负责真正的FM解调。 gr.quadrature_demod_cf通过比较每个邻近的复采样信号的角度（by subtracting the angle of each adjacent complex sample）来有效的区分频率，gr.quadrature_demod_cf的输出包含了FM立体声的左加右声道信号（left-plus- r

23、ight FM mono audio signal），19kHZ的立体声的控制音（pilot tone）,中心频率是38KHZ的左减右（left-minus-right）信号和其他的在这个上面的次载波（sub-carriers）信号。在这 个简单的接收机中，我们使用低通滤波和抽取滤波，提取出了left-plus-right信号信息， 把这个信号送给声卡，经过声卡的32，000HZ的DAC，就能够把数字信号还原成模拟声音，到此我们可以听到FM广播信号了。如果想进一步了解FM接收 机的原理，请参考：Listening to FM, Step by Step.图形化用户接口GNU Radio应用程序

24、的图形化接口是使用python来实现的，接口能使用pyhone的任何toolkit来实现，我们推荐使用wxPython，它能最大化的实现跨平台应用。GNU Radio提供由C+到Python的连接机制。硬件要求GNU Radio是硬件独立的，今天带有单独浮点运算单元的上GHZ的CPU已经很常见了，这给台式机实现数字信号处理功能带来了可能。一个3G的 Pentium 或者 Athlon处理器能够每秒处理30亿次浮点FIR运算。我们现在能够在PC上建立软件通讯系统，这是几年前所不敢想象的。 你的硬件要求依赖于你想做什么。总的来说，一个1G或者2G带有256M内存的机器应该是足够了。你也需要一些其他

25、的模拟外设连接在你的PC上，包括内置 的声卡或者96 kHz, 24-bit高保真声卡。使用这些模拟设备，你只能处理有限的窄带信号。另外的方案是使用高速的PCI模拟到数字的外设，这些设备可以达到20M的抽样率，但是他们是比较昂贵的，差不多是买一台PC的价钱。使用这些高速板，调制解调器可以作为RF前端来使用。为了有一个更满意的替代产品，现在我们已经设计了一个通用软件无线电外设（USRP）。通用软件无线电外设我们首选的硬件解决方案是USRP,USRP是Matt Ettus的杰作，USRP是一个非常灵活的USB设备，它把你的PC连接到RF世界。USRP包含一个小的母板，母板包含4个12bit /64

26、M抽样率的ADC，4个14bit/128M DAC,一个百万门的FPGA芯片和一个可编程的USB2.0控制器。每个USRP母板支持4个子板，2个接收，2个发射。RF前端是实现在子板上的，不 同的子板处理不同的频率带宽。作为业余无线电使用，可以选择低能量的能够接收和发送440MHZ和1.24GHZ的子板。一个仅有接收功能的基于 cable modem 调谐的子板覆盖50MHZ到800MHZ频率范围被设计成容易手工来做原型的，目的是为了方便实验。图二 通用无线电外设USRP 的灵活性得益于2个可编程的元件，通过他们和PC上的host端交互。为了对USRP有一个初步映像，让我们看一下它的启动过程。U

27、SRP本身不含有 ROM，仅仅有一个存储VendorID和productID和版本号信息的EEPROM。当USRP插到PC的USB口上以后，主机上的程序通过 VID，PID和版本号识别这是一个未配置的USRP，主机上的程序第一步要做的是下载一个8051固件到USB控制芯片上，这个固件控制USB的行为。 当USB固件下载好后，USRP模拟一个USB设备的重枚举过程，此后主机识别到一个不同的设备，VID,PID和版本号都不同了。现在这个USB固件定 义USB端口，接口和用户自定义命令。其中一个命令是load FPGA，收到这个命令后USB设备就能够把FPGA配置bitstream下载到FPGA芯片

28、中开始工作。FPGA是一个通用硬件，它的行为完全由配置的bitstream来决定，你可以把bitstream看作是一个目标码。这个bitsteam是由一个高 级硬件描述语言编译得到的，在USRP里面这是由verilog硬件描述语言来实现的。这些代码是开源的，和其他的GNU Radio代码一样，是基于GNU的GPL协议的。FPGA职责FPGA 像一个小的，高性能的并行计算机一样，可以完成你设计的任务，设计FPGA需要一些技能，并且如果不慎还会烧坏你的板子。还好我们已经提供一个标准的适应 性很广的FPGA配置。通过使用一个好的USB控制器，USRP能够有32M/s的处理能力，USB是半双工的，基于

29、你的需要，USB能够在传输和接收之 间转换。在接收方式下，标准的配置能够允许你选择你感兴趣的频率，同时完成基带化和抽取滤波的工作。和RF前端处理方法一样，但现在做的是在数字域的采样。执行这些功能的代码叫做数字下变频转换器DDC（Digital Down Converter Block Diagram），利用这些处理我们能够在数字域快速的改变中心频率。图三 数字下变频转换器在传输路径上，执行一个完全反过来的过程。FPGA包含多个数字上下变频实例，根据需要这些实例连接到相同或者不同的ADC上面。鉴于篇幅，这里不能介绍所有的理论，在GNU Radio的wiki上有更加多的信息可以参考。GNU Rad

30、io应用除了上面提到的例子，GNU Radio还实现了一个完全的HDTV传输和接收的例子，一个频谱分析仪，一个示波器，以及并发多通道接收器和很多调制解调的例子。目前利用GNU Radio在研究和进行中的工程有： 一个和TiVo类似的无线电设备，但是它能够同时记录多个电台的信号. TDMA 波形 利用广播TV信号作为信号源的无源雷达系统。利用它你也许可以侦测到从你上空飞过的飞机。 天文无线电（译注: 射电望远镜？） TETRA 收发机 Digital Radio Mundial (DRM). 软件GPS 分布式传感器网络 分布式频谱利用测量系统 业余无线电收发机 Ad Hoc自组织网络 RFID

31、 系统 多输入多输出系统(MIMO)政策问题每项革新都会遇到政治上的阻碍，构建软件无线电也会遇到麻烦，在美国，我们不小心站到了全美音像学会（Motion Picture Association of America）的对立面，他们对TV数字信号进行加密以限制不同的接收器。另外，美国联邦通讯委员会也已经发布一个感知无线电和软件无线电的制造规则提议（Proposed Rule Making (NPRM)）。在NPRM中有几个麻烦的问题被提出，包括限制高速DAC的销售，软件无线电发出的信号需要被授权，和限制未授权的软件无线电在业余无线电市场上的销售。总结软件无线电是一个很有意义的领域，GNU Radio提供了一个探索软件无线电领域的工具，想对软件无线电有一个深入的了解需要很多学科的知识积累。我们所做的就是尽可能的降低这个入门的门槛。