Exploring-gnuradio(中文).doc

Exploring-gnuradio简介： 软件无线电是一门让软件代码尽可能靠近接收天线的技术。它把无线电的硬件问题放到软件中来解决。软件无线电最基本的特征是软件定义（调制）无线电传输的波形，软件解调接收到的波形。这个和目前大多数的由模拟电路或者由模拟数字组合电路组成的无线电有鲜明的区别。GNU Radio是一个开源的可以构建软件无线电平台的软件包。软 件无线电是无线电设计领域的一次革命，使得建立无线电设备更加灵活，给用户带来更多的机遇。软件无线电能够很好的实现传统无线电所能实现的功能，另外，软 件无线电最精彩的特性是由软件提供给我们的灵活性，一些年以后，我们将会看到传统的由很多固定元件组成的无线电设备将会被通用的通信设备来取代。可以想 像，一个软件无线电设备能够变体成为蜂窝电话，能够灵活的使用GPRS,802.11 Wi-Fi,802.16 WiMax,卫星链路，或者将出现的其他通讯标准进行通讯，而硬件设备无需改动，只需要软件设置。你能够决定使用GPS，GLONASS或者两个一起使用来给自己定位。 也 许，最令人激动的是使用软件无线电来建立一个分布式个人通讯系统。当今的通讯系统，大多数是采用自上而下的结构，广播和电视系统提供一个单向链路，内容被 严格控制在一小部分人的手里。蜂窝电话系统给人们带来极大的便利，但是你的手机提供的功能是由运营商来控制的，而不是你自己。这种集中控制的系统限制了人们的创新，代替做蜂窝电话系统的二等公民，我们可以建造一个智能设备，这些设备能够自组织的，在使用者之间建成一个网络。处理流程图图一，表示一个典型的软件无线电处理流程图。为了理解无线电的软件模块，首先需要理解和其关联的硬件。在这个图中的接收路径上，能够看到一个天线，一个神奇的RF前端，一个模拟数字转换器ADC和一堆代码。ADC是一个连接连续模拟的自然世界和离散的数字世界的桥梁。ADC有两个主要特性，抽样率和动态范围。抽样率是ADC测量模拟信号的速度，动态范围是ADC区别最低信号值和最大信号值的精度，这决定ADC数字信号输出的比特数（位数）。例如，8位的AD转换器最多能代表256个信号层次，而一个16位的转换器能够代表65536个层次信号。总的来说，ADC的物理特性和价格决定了抽样率和动态范围。 在我们深入研究软件之前，先来了解一些理论知识，在1927年，出生于瑞典的物理和电子学家 Harry Nyquist提出了如果AD转换想没有混叠现象发生，那么抽样率至少是目标信号带宽的2倍。 混叠现象就像是车子重复的碾在过去的车轮印迹上一样，让你分不清楚是本次的还是以前碾过的车轮印迹。假设我们要处理一个低通信号，我们感兴趣的信号带宽是 0到fmax，按照Nyquist理论，抽样率必须至少是2*fmax。如果我们的ADC工作在20MHZ，但是我们想收听92.1MHZ的FM电台，我 们该怎么办呢？答案是使用RF前端，接收机的RF前端能够把它接收到的高频段信号下变频到一个低频段信号后输出。例如，我们能让RF前端把90- 100MHZ频段内的信号下变频到0-10MHZ的低频范围内，那么我们的20MHZ的ADC就能够派上用场了。 大多数情况下，我们把RF前端当作一个信号控制的黑盒子，负责处理输入信号的中心频率。举一个具体例子，一个调制解调器的调制模块能够把50M到800M 之间的6M带宽的中心频率下变频到一个中心频率是5.75MHZ的输出。这个输出的中心频率通常称作中频(IF)。 按照越简单越易用的原则，RF前端最好也能够被去掉，有一个GNU Radio用户已经成功的使用一个100英尺（译注：等于12英寸，合0.305米）的天线直接连接到一个20M抽样率的ADC上收听到了AM和短波广播。 深入软件内部 GNU Radio提供一个信号处理模块的库，并且有个机制可以把单个的处理模块连接在一起形成一个系统。编程者通过建立一个流向图（flow graph）就能搭建成一个无线电系统。信号处理模块是使用C+来实现的，理论上说，信号数据流不停的从输入端口流入从输出端口流出。信号处理块 (blocks)的属性包括输入和输出的端口数，流过它们的数据的类型，经常使用的数据流的类型是短整型(short),浮点型(float),和复数 （complex）类型。一些处理模块仅仅有输入端口或者输出端口，它们分别成为信号源(data source)和信号接收器(sink)。有的信号源从文件或者ADC读入数据，信号接收器写入文件或者DAC或者PC的多媒体接口。GNU Radio提供了超过100个信号处理块，并且扩展新的处理模块也是非常容易的。软件图形化接口和信号处理模块的链接机制是通过python脚本语言来进行的，例1，是一个GNUadio “Hello World”的例子。它产生两个sine波形并且把他们输出到声卡，一个输出到声卡的左声道，一个输出到右声道。 例子1.输出拨号音#!/usr/bin/env pythonfrom gnuradio import grfrom gnuradio import audiodef build_graph (): sampling_freq = 48000 ampl = 0.1 fg = 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) 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个信号处理模块用流向图的connect方法连接到一起。connect方法有两个参数，源端点和目的端点，用来建立一个从源到目的处理模块 的链路。每个端点(endpoint)有两个成员：一个信号处理模块和一个端口号。端口号表示哪个输入或者输出端口应该被连接的。通常端点使用 python语言的tuple来表示，像：（block,port_number）。当端口号是0时，portt_number可以被省略的。 下面的两个表示方法是一样的:fg.connect (src1, 0), (dst, 1)fg.connect (src1, (dst, 1) 一旦流向图被建立了，我们调用start生成一个或者多个线程去运行它，按下任意键程序控制权就会返回给调用者。一个完整的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 (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 = 20e6 cfir_decimation = 125 audio_decimation = 5 quad_rate = input_rate / cfir_decimation audio_rate = quad_rate / audio_decimation fg = gr.flow_graph () # use high speed ADC as input source src = high_speed_adc (fg, input_rate) # compute FIR filter taps for channel selection channel_coeffs = gr.firdes.low_pass (1.0, # gain input_rate, # sampling rate 250e3, # low pass cutoff freq 8*100e3, # width of trans. band gr.firdes.WIN_HAMMING) # input: short; output: complex chan_filter1 = gr.freq_xlating_fir_filter_scf (cfir_decimation, channel_coeffs, freq1, # 1st station freq input_rate) (head1, tail1) = build_pipeline (fg, quad_rate, audio_decimation) # sound card as final sink audio_sink = audio.sink (int (audio_rate) # now wire it all together fg.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 pipeline for demodulating a broadcast FM signal. The input is the downconverted complex baseband signal. The output is the demodulated audio. build_pipeline returns a two element tuple containing the input and output endpoints. fm_demod_gain = 2200.0/32768.0 audio_rate = quad_rate / audio_decimation volume = 1.0 # input: complex; output: float fm_demod = gr.quadrature_demod_cf (volume*fm_demod_gain) # compute FIR filter taps for audio filter width_of_transition_band = audio_rate / 32 audio_coeffs = gr.firdes.low_pass (1.0, # gain quad_rate, # sampling rate audio_rate/2 - width_of_transition_band, width_of_transition_band, gr.firdes.WIN_HAMMING) # input: float; output: float audio_filter = gr.fir_filter_fff (audio_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 line freq1 = float (args0) * 1e6 else: sys.stderr.write (usage: fm_demod freqn) sys.exit (1) # connect to RF front end rf_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 gain rf_front_end.set_AGC (300) # determine the front ends Intermediate Frequency IF_freq = rf_front_end.get_output_freq () # 5.75e6 # Tell the front end to tune to freq1. # I.e., freq1 is translated down to the IF frequency rf_front_end.set_RF_freq (freq1) # build the flow graph fg = build_graph (IF_freq, None) fg.start () # fork thread(s) and return raw_input (Press Enter to quit: ) fg.stop ()if _name_ = _main_: main (sys.argv1:) 和Hello World的例子一样，这个例子建立了一个流向图，把所有的信号处理模块连接到了一起然后开始运行。在这个例子中，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来选择一个我们想收听的频道。 在这个例子中，我们选择了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- right 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应用程序的图形化接口是使用python来实现的，接口能使用pyhone的任何toolkit来实现，我们推荐使用wxPython，它能最大化的实现跨平台应用。GNU Radio提供由C+到Python的连接机制。 硬件要求GNU Radio是硬件独立的，今天带有单独浮点运算单元的上GHZ的CPU已经很常见了，这给台式机实现数字信号处理功能带来了可能。一个3G的 Pentium 或者 Athlon处理器能够每秒处理30亿次浮点FIR运算。我们现在能够在PC上建立软件通讯系统，这是几年前所不敢想象的。 你的硬件要求依赖于你想做什么。总的来说，一个1G或者2G带有256M内存的机器应该是足够了。你也需要一些其他的模拟外设连接在你的PC上，包括内置 的声卡或者96 kHz, 24-bit高保真声卡。使用这些模拟设备，你只能处理有限的窄带信号。另外的方案是使用高速的PCI模拟到数字的外设，这些设备可以达到20M的抽样率，但是他们是比较昂贵的，差不多是买一台PC的价钱。使用这些高速板，调制解调器可以作为RF前端来使用。 为了有一个更满意的替代产品，现在我们已经设计了一个通用软件无线电外设（USRP）。 通用软件无线电外设 我们首选的硬件解决方案是USRP,USRP是Matt Ettus的杰作，USRP是一个非常灵活的USB设备，它把你的PC连接到RF世界。USRP包含一个小的母板，母板包含4个12bit /64M抽样率的ADC，4个14bit/128M DAC,一个百万门的FPGA芯片和一个可编程的USB2.0控制器。每个USRP母板支持4个子板，2个接收，2个发射。RF前端是实现在子板上的，不 同的子板处理不同的频率带宽。作为业余无线电使用，可以选择低能量的能够接收和发送440MHZ和1.24GHZ的子板。一个仅有接收功能的基于 cable modem 调谐的子板覆盖50MHZ到800MHZ频率范围被设计成容易手工来做原型的，目的是为了方便实验。 Figure 2. Universal Software Radio Peripheral USRP 的灵活性得益于2个可编程的元件，通过他们和PC上的host端交互。为了对USRP有一个初步映像，让我们看一下它的启动过程。USRP本身不含有 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芯片中开始工作。 FPGA是一个通用硬件，它的行为完全由配置的bitstream来决定，你可以把bitstream看作是一个目标码。这个bitsteam是由一个高 级硬件描述语言编译得到的，在USRP里面这是由verilog硬件描述语言来实现的。这些代码是开源的，和其他的GNU Radio代码一样，是基于GNU的GPL协议的。 FPGA职责 FPGA 像一个小的，高性能的并行计算机一样，可以完成你设计的任务，设计FPGA需要一些技能，并且如果不慎还会烧坏你的板子。还好我们已经提供一个标准的适应 性很广的FPGA配置。通过使用一个好的USB控制器，USRP能够有32M/s的处理能力，USB是半双工的，基于你的需要，USB能够在传输和接收之 间转换。在接收方式下，标准