摘要
随着移动通信技术的发展，与日俱增的移动用户数量和日趋丰富的移动增值服务，使无线通信的业务量迅速增加，无限电波有限的带宽远远满足不了通信业务需求的增长。另一方面，由于移动通信系统中的同频干扰和多址干扰的影响严重，更影响了无线电波带宽的利用率。并且无线环境的多变性和复杂性，使信号在无线传输过程中产生多径衰落和损耗。这些因素严重地限制了移动通信系统的容量和性能。因此为了适应通信技术的发展，迫切需要新技术的出现来解决这些问题。这样智能天线技术就应运而生。智能天线是近年来移动通信领域中的研究热点之一，应用智能天线技术可以很好地解决频率资源匮乏问题，可以有效地提高移动通信系统容量和服务质量。开展智能天线技术以及其中的一些关键技术研究对于智能天线在移动通信中的应用有着重要的理论和实际意义。
论文的研究工作是在MATLAB 软件平台上实现的。首先介绍了智能天线技术的背景；其次介绍了智能天线的原理和相关概念，并对智能天线实现中的若干问题，包括:实现方式、性能度量准则、智能自适应算法等进行了分析和总结。着重探讨了基于MATLAB 的智能天线的波达方向以及波束形成，阐述了music 和capon 两种求来波方向估计的方法，并对这两种算法进行了计算机仿真和算法性能分析;
关键字:智能天线；移动通信；自适应算法；来波方向； MUSIC算法
Abstract
With development of mobile communication technology ，mobile users and communication，increment service are increasing ，this make wireless services increase so that bandwidth of wireless wave is unfit for development of communication ，On the other hand ，much serious Co-Channel Interruption and the Multiple Address interruption effect utilize rate of wireless wave’s bandwidth，so the transported signals are declined and wear down，All this has strong bad effect on the capacity and performance of question and be fit for the development of communication ，so smart
antenna arise Smart Antenna ，which is considered to be a solution to the problem of lacking frequency, becomes a hotspot in the Mobile Communication area ．With this technology, Capacity of Mobile Communication system can be increased effectively and the quality of service can be improved at the same time. To study Smart Antenna and its key technologies is important both in theory and in practice 。
All of the research work of this paper is based on the MA TLAB software environment．First ，in this paper，we make an introduction on the appearance background of the Smart Antenna technology and its relative theory and concept ．In addition, some problems about its realization such as modes of its realization, rules of its performance,
adaptive algorithms are analyzed. Focused on the smart antenna based on MA TLAB the DOA and beam forming， and capon on the music for two to wave the direction of the estimated method. In order to decrease calculating-time and complexity of the algorithm, a rule of maximum received signal is presented; Next several problems about realization of the algorithm are discussed。
Key Words : Smart Antenna；Mobile Communication；Adaptive algorithm Direction Of Arrival；Cyclic-MUSIC arithmetic
目录
第一章绪论 (6)
1.1 课题背景 (6)
1.2 主要概念 (7)
1.3 国内外研究现状 (8)
1.3.1 国外研究现状 (8)
1.3.2 国内研究现状 (10)
1.4 本文研究的意义 (11)
1.5 本论文研究的主要内容及章节安排 (11)
第二章智能天线 (12)
2.1 概述 (12)
2.2 移动通信环境与智能天线的信号模型 (12)
2.2.1 移动通信环境 (12)
2.2.2 移动通信中的空间信道模型 (14)
2.3 智能天线的基本原理及结构 (14)
2.3.1 智能天线的基本原理 (14)
2.3.2 智能天线系统的基本结构 (16)
2.4 智能天线的主要功能及应用 (18)
第三章基于MATLAB 的波达方向研究以及波束形成 (20)
3.1 MATLAB概述 (20)
3.1.1 MATLAB 语言简介 (20)
3.1.2 MATLAB 语言特点 (20)
3.2 天线阵的波达方向估计的MATLAB 仿真 (22)
3.2.1 原理 (22)
3.2.2波达方向估计的方法 . (24)
3.3 天线阵波束形成的MATLAB 仿真 . (27)
第四章总结.......................................... 31 参考文献 (32)
致谢 (32)
第一章绪论
1.1 课题背景
移动通信作为未来个人通信的主要手段，在全球通信业务中占据越来越重要的地位。随着移动通信用户数的迅速增长以及人们对通话质量要求的不断提高，要求移动通信网在大容量下仍保持较高的服务质量。而与此要求相对，目前移动通信中主要存在两大问题：第一，随着移动用户的增多，频谱资源日益匾乏；第二，由于信道传输条件较恶劣，所需信号在到达天线接收端前会经历衰减、衰落和时延扩展，另外还有来自其他用户的干扰，极大地限制了系统通信质量的提高。这两大问题是移动通信技术发展的主要矛盾，也是推动移动通信技术发展的原动力。必须采取有效方法对系统进行扩容并提高服务质量。为了解决系统容量问题，第二代数字蜂窝系统中主要采用时分多址(TDMA和码分多址(CDMA两种多址方式; 为了提高系统通信质量，在第二代系统中广泛采用了调制、信道编码、均衡(TDMA系统、RAKE 接收(CDMA系统等时、频域信号处理技术以及分集天线、扇形天线等简单空间处理技术。这些解决方法在发挥各自功效的同时，有着共同的不足，即无法对空域资源进行有效利用。
智能天线技术正是在这样的背景下被引入到移动通信中来的。理论研究和实测数据均表明:有用信号、其延时样本和干扰信号往往具有不同的到达角(DOA和空间信号结构，利用这种空域信息我们可以获得附加的信号处理自由度，从而能提高系统容量，并且能够更有效地对抗衰落和抑制干扰。应用于无线通信系统基站的智能天线技术正是充分利用了信号的空域信息，它能有效地扩充系统的容量, 大幅度提高系统的通信质量。
智能天线技术己经被公认为第三代移动通信系统的一项关键技术，并越来越受到人们的关注。在提交国际电联ITU 所有的3GRTT 标准中，几乎都附有一条:如果有可能，本建议将采用智能天线技术:在国际电联2000年3月份的会议上，更是提出要重视在CDMA 系统中使用智能天线技术，并在2000年8月份的会议上正式讨论了在CDMA 系统中使用智能天线的问题。可以预见，智能天线技术将在未来的移动通信体制中占据非常重要的地位。
1.2 主要概念
智能天线又称为自适应天线阵列，英文名为Smart Antenna或Intelligent Antenna 。智能大线技术的核心是阵列信号处理，早期应用集中于雷达和声纳信号处理领域，七十年代后期被引入到军事通信中，而应用于民用蜂窝移动通信则是近十儿年的事情。一般而言，智能天线是专指用于移动通信中的自适应天线阵列(这里的移动通信系统主要指数字蜂窝移动通信系统，它利用数字信号处理技术产生空间定向波束，使天线的主波束跟踪所需用户信号到达方向，旁瓣或零陷对准不希望的干扰信号到达方向，达到充分分离和有效利用用户信号并删除或抑制干扰信号的目的。在移动通信的基站中使用具有全向收发功能的智能天线，可以为每个用户提供一个窄的定向波束，使信号在有限的方向区域内发送和接收，这样就可以充分利用信号发射功率，降低信号全向发射带来的电磁干扰与相互干扰。智能天线是提高无线电数据通信，包括蜂窝通信、个人通信和第三代宽带CDMA 等系统容量的最佳选择，它超越了任何由信道复用和各种调制技术所达到的水平。
CDMA(Code Division Multiple Access是码分多址的英文缩写，它是在数字技术的分支—扩频通信技术上发展起来的一种崭新而成熟的无线通信技术。CDMA 技术的原理是基于扩频技术，即将需传送的具有一定信号带宽信息数据，用一个带宽远大于信号带宽的高速伪随机码进行调制，使原数据信号的带宽被扩展，再经载波调制并发送出去。接收端使用完全相同的伪随机码，与接收的带宽信号作相关处理，把宽带信号换成原信息数据的窄带信号即解扩，以实现信息通信。
SDMA(Space Division Multiple Access 是空分多址的英文缩写，移动通信中应用智能天线技术就产生了这种新的信道增容方式。它不同于传统的频分多址(FDMA,时分多址(TDMA或码分多址(CDMA，这种多址方式是利用用户空间位置的不同来区分不同用户，也就是说，在相同时隙、相同频率或相同地址码的情况下，仍然可以根据信号不同的空间传播路径而区分不同的信号。空分多址可以与其他多址方式完全兼容，从而可实现组合的多址方式，例如空分—码分多址(SD-CDMA、空分—时分多址(SD-TDMA等，这样可以成倍地增长系统容量。码间干扰(ISI Inter-Symbol Interference是数字通信系统中除噪声干扰
之外最主要的干扰，它与高斯分布的加性噪声干扰不同，是一种乘性干扰。信道的衰减和时延失真等都可能引起ISI ，实际上，只要传输信道的频带是有限的就不可避免地带来一定的ISI. 以一定速度传送的波形受到非理想信道的影响表现为各码元波形持续时间拖长，从而使相邻码元波形产生重叠，从而引起判决错误，当这种线性失真严重时，码间干扰显得尤为突出。
同信道干扰(CCI Common Channel Interference，又叫同频干扰，它是指使用相同频率的信道之间的干扰。在蜂窝移动通信中，同信道干扰主要指使用相同频率的小区间的干扰。
多址干扰(MAI Multiple Access Interference ，是在码分多址蜂窝移动通信中出现的一种干扰。由于在同一个小区内同时通信的用户是多个，多个用户均占同一时隙、同一频率，所不同的是选取的地址码不一样，而实际选用的地址码间的互相关函数不可能全为零，这样多个用户同时通信时必然会产生多址干扰。
天线增益——取定向天线主射方向上的某一点，在该点场强保持不变的情况下，此时用无方向性天线发射时天线所需的输入功率与采用定向天线时所需的输入功率之比称为天线增益，常用“G ”表示，天线增益可以用来描述天线往某一方向发射的能力。
1.3 国内外研究现状
1.3.1 国外研究现状
移动通信在经历了第一代模拟通信系统、第二代蜂窝数字通信系统和窄带CDMA 系统，正向第三代移动通信系统发展。目前正处于确立第三代移动通信技术标准之时，国外如欧美等发达国家非常重视智能天线技术在未来移动通信方案中的地位与作用，己经开展了大量的理论分析和研究，同时也建立了一些技术试验平台。
1. 欧洲
欧洲通信委员会在计划中实施了智能天线技术第一阶段研究，称之为TSUNAMI(The Technology in Smart antennas for Universal Advanced Mobile Infrastructure ，德国、英国、丹麦和西班牙合作完成。该项目组在DECT 基站
基础上构造智能天线试验模型，于1995年初开始现场试验。天线由8个阵元组成，射频工作频率为1.89GHz ，阵元间距可调，阵元分布分别有直线形、圆环形和平面形三种。模型用数字波束形成方法实现智能天线，采用ERA 技术有限公司的专用集成电路芯片DBF1108完成波束形成，系统评估了识别信号到达方向的多用户信号识别分类算法(MUSIC，采用的自适应算法有归一化最小均方算法(NLMS和递归最小平方算法。
实验验证了智能天线的功能，在2个用户4个空间信道(包括上行和下行链路情况下，试验系统比特差错率(BER优于310 。现场测试结果表明，圆形和平面形天线适用于室内通信环境，而市区环境则采用简单的直线阵更合适。
2. 日本
ATR光电通信研究所研制了基于波束空间处理方式的多波束智能天线。天线阵元布局为间距半波长的16阵元平面方阵，工作频率为1.545GHz. 阵元接收信号先经过模数变换，然后进行快速傅里叶变换处理，形成正交波束后，分别采用恒模算法或最大比值合并(MRC算法。数字信号处理部分由10片现场可编程门阵列完
成，整块电路板大小为23.3*34.0(CM。野外移动试验确认了采用恒模算法的多波束智能天线的功能。
根据用户所处环境不同，影响系统性能的主要因素(如噪声、同信道干扰或符号间干扰也不同，利用软件方法在不同环境应用不同算法。比如当噪声是主要因素时使用多波束最大比值合并算法; 当同信道干扰是主要因素时则使用多波束恒模算法:而要抵消符号间干扰时采用最小均方算法和最大似然连续估计算法(MLSE，以此提供算法分集。
3. 美国
美国对于智能天线技术的研究水平处于世界领先地位，并且许多电信设备生产商纷纷推出了自己的产品。最具代表性的是爱瑞通信公司(ArrayComm，它是一家拥有成熟的自适应智能天线技术的全球知名通信技术公司，在这一领域拥有多项专利技术，居世界领先水平。爱瑞通信公司拥有丰富的自适应智能天线的产品线，其中IntelliCell 技术在通信系统中的应用能有效地改善信号质量和频谱
利用率，使系统容量和覆盖范围增大，提高数据传输速率，从而获得最佳的语音质量。IntelliCell 处理器通过自适应处理算法，形成可以加权参数，在幅度，相位和信号空间到达角等多个指标_L进行每秒调整数百次的调整，从而完成上行处理和下行波束形成。该技术己经在全球超过7.5万个基站系统上得到应用，为450万名无线用户提供高质量的无线宽带(Wireless Bandwidth服务。这一技术支持第三代移动通(3G各种空中接口标准(IMT2000WCDMA,CDMA2000,TD-SCDMA，在容量、数据传输速率、覆盖范围及服务质量上都比传统移动通信系统具有明显的优势。据称，IntelliCell 可以使运营商的基站数量减少50%，由此可以减少大量的设备成本和营运成本。
另外，德州大学奥斯丁SDMA 小组建立了一套智能天线试验环境，并着手理论研究以及与实际系统相结合的研究。
4 .加拿大
加拿大McMaster 大学开发了四元I-4列天线，并进行了恒模(CMA Constant Module Algorithm算法的研究。
1.3.2 国内研究现状
国内对于智能天线的研究起步较晚，但也取得了一些成绩。如北京信威公司研制了应用于无线本地环路(WLI的智能天线系统，信威公司的智能天线采用8阵元环形自适应阵列，工作于1785-1805MHz, 采用时分双工方式，收发间隔10ms ，接收机灵敏度最大可提高9dB 。但该系统只能工作于无线本地环路中，并对用户位置、移动速率有一定的要求。
在国内一些大学和研究结构，如清华大学、西安交大、中国科技大学、西安电子科技大学、北方交通大学、北京邮电大学、电信科学技术研究院等相继开展了智能天线的理论研究一些大的电信设备生产企业如大唐电信、华为、中兴科技等也投入了很多的人力物力进行研发; 国家“八六三”、国家自然科学基金、博士点基金等也相应支持有关单位进行理论与技术平台的研究。
1.4 本文研究的意义
智能天线有着非常诱人的应用前景，许多国家都投入了大量的人力物力对该
技术进行研究，并取得了一些成就。概括地讲，目前研究主要解决了以下两方面
的问题:
1.研究论证了智能天线在不同移动通信系统中应用的可行性和有效性，建立
了一些技术试验平台，并且在一定的条件下（从目前情况来看，智能天线正逐步
应用在固定无线接入系统中，即用户固定无线传播环境不断变化的情况）实现了
智能天线技术，验证了智能天线在提高移动通信系统性能中表现出的强大优势。
2. 研究了智能天线基本结构以及功能模块，并提出了一些智能天线的性能度
量准则和自适应波束形成的算法.
智能天线应用于3G 系统时，同时也带来了相应的新问题，如：智能天线的校
准、智能天线和其它抗干扰技术的结合、波束赋形的速度问题、设备复杂性的考
虑、帧结构及有关物理技术等。
波束形成是智能天线的关键技术之一，本论文主要研究智能天线的波束形成
技术。
1.5 本论文研究的主要内容及章节安排
本论文主要研究智能天线波束形成技术，包括基本思想、技术原理，研究智
能天线波束形成的MA TLAB 仿真技术，包括天线阵波达方向估计MA TLAB 仿
真，以及天线阵波束形成的MATLAB 仿真。
本论文的章节安排如下：
第一章智能天线的研究现状以及本课题研究的意义；
第二章智能天线的原理以及结构和功能；
第三章利用MATLAB 仿真工具，重点对天线阵波达方向估计以及天线阵的
波束形成进行仿真。
第二章智能天线
2.1 概述
应用于具有复杂电磁波传播环境的移动通信系统中的智能天线是一种具有
上行接收和下行发射双向功能的天线系统，它能够有效提高天线系统的可靠性与
灵活性，增加通信系统容量和改善通信质量。本章将对这一技术展开研究，前半
部分介绍智能天线系统原理及应用，内容包括智能天线的应用环境、信道模型、
工作原理、结构、应用; 后半部分主要探讨智能天线技术实现中的若干问题，内
容包括智能天线的实现方式、性能度量准则等。
2.2 移动通信环境与智能天线的信号模型
移动通信系统都有着复杂的电磁波传播环境，从某种意义上说，智能天线的
工作原理可以理解为根据用户信号的不同空间传播方向调节用户信道参数，从而
减少干扰对系统的影响。因此，了解移动通信的环境以及复杂移动通信环境所引
起的干扰，建立可以有效而精确地预测无线系统的信道模型，对于智能天线技术
的研究和实现有着十分重要的意义。
2.2.1 移动通信环境
随着无线电技术的发展和广泛应用，空间传输的电磁波变得越来越复杂。应
用智能天线技术就是为了有效消除电波传输的不利影响。影响移动通信中信号传
输的主要因素有:通信信道的特点、电磁波传播方式、传输损耗及传输效应，下
面分别予以介绍。
1. 移动通信信道的主要特点
1传播的开放性，即一切无线信道都是基于电磁波在空间传播来实现信息传输的。
2接收地理环境的复杂性与多样性，有城市中心繁华区近郊小城镇农村及远郊区。
3通信用户的随机移动性，有慢速步行的通信和高速车载不间断的通信。
2. 电磁波是无线信号传输的介质，移动通信下的电磁波传播方式有:
1直射波 :它是指视距覆盖区无遮挡的传播，直射波信号最强。
2多径反射波：指不同建筑物或其它物反射后到达接收点的传播信号。
3绕射波：从较大的山丘或建筑物绕射后到达接收点的传播信号，其强度与放射波相当。
4散射波：由于空气中离子受激后二次反射后引起的漫反射后到达接收点的传播信号，其信号强度最弱。
3. 电磁波在传输过程中是有损的，移动通信电磁波传输可能会有下面三种不同的
损耗:
1路径传播损耗：又称为衰耗，它是指电波在空间传播所产生的损耗，它反映
了传播在宏观大范围的空间距离上的接收信号电平平均值的变化趋势。
2 慢衰落损耗：它是由于电磁波在传播路径上受到建筑物及山丘等的阻挡所
产生的阴影效应而产生的损耗。它反应了中等范围内数百波长量级接收电平的均
值变化而产生的损耗，一般遵从对数正态分布，其变化率较慢故又称为慢衰落。
3 快衰落损耗：它主要是由于多径传播而产生的衰落，它反应微观小范围内
数十波长量级接收电平的均值变化而产生的损耗，一般遵从瑞利分布或莱斯分
布，其变化率比慢衰落快，故称为快衰落。
4. 电磁波在传输过程中有不同的传播路径，且存在传输损耗; 另一方面移动通信
用户分布较为分散，用户与移动台的相对位置是变化的。因此，在移动通信中存
在三种不同效应:
1 阴影效应：由于大型建筑物和其他物体的阻挡而形成在传播接收区域上的
半盲区。
2 远近效应：由于接收用户的随机移动性，移动用户与基站间的距离也是随
机地变化，若各移动用户发射功率一样，那么到达基站的信号强弱不同，离基站
近信号强，离基站远信号弱。通信系统的非线性则进一步加重，出现前者更前，
弱者更弱和以强压弱的现象，通常称这类现象为远近效应。
3 多普勒效应：它是由于接收的移动用户高速运动而引起传播频率的扩散而
引起的，其扩散程度与用户运动的速度成正比。
智能天线技术正是应用在这样一个复杂的环境中，在智能天线的算法设计和
系统实现中我们必须充分考虑上面这些因素的影响。
2.2.2 移动通信中的空间信道模型
移动通信的电磁波传播环境复杂，那么究竟信号传输会经历那些过程，受到
那些因素的影响呢? 下面我们将研究移动通信中的空间信道模型。
无线移动通信的空间信道，指将信号分配单元输出端到信号组合单元输入端
之间的部分合并，简称信道。智能天线技术直接建立在信道参量的基础上，无线
移动信道的建模与估计是智能天线技术的基础，无论是算法描述还是算法的性能
分析以及仿真都必须依赖这一基石。而且，根据无线传播理论和对各种通信环境
的实际测量建立合理的无线移动信道模型，可以降低智能天线自适应算法的对实
时测量的要求。经典的模型仅考虑接收信号的功率水平分布和多普勒频移，而现
代空间信道模型建立在衰落和多普勒扩展的基础上，结合了时延扩展、到达方向
和自适应天线阵的几何特性。
2.3 智能天线的基本原理及结构
2.3.1 智能天线的基本原理
简单说，智能天线的基本原理就是根据一定的接收准则自动地调节天线阵元
的幅度和相位加权值，以实现最优接收和发射。从空间响应来看，智能天线是一
个空间滤波器，它在信号入射方向上增益最大，在干扰信号入射方向上形成零陷
或低陷。
下面以直线阵为例，说明智能天线的基本原理。假设满足天线传输窄带条件，
即某一入射信号在各天线单元的响应输出只有相位差异而没有幅度变化，这些相
位差异由入射信号到达各天线所走路线的长度差来决定。若入射信号为平面波
(只有一个入射方向，则这些相位差由载波波长、入射角度、天线位置分布唯一
确定。给定一组加权值，一定的入射信号强度，不同入射角度的信号由于在天线
间的相位差不同，合并后的输出信号强度也会不同。智能天线的工作原理可用下
图说明:
n
图2-1 智能天线基本原理图
对于N 元天线阵列，设信号自α方向入射，阵元间距为d ，接收信号功率为P ，
以第1个天线单元为参考，则第i 个天线单元的相位延迟为Φ=β(i-1dcos(α ，
其中β=2π/λ为电磁波波长。
智能天线的目的就是确定最佳权值矢量W ，己达到提取有用信号，抑制干扰
信号滤除噪声信号的目的。智能天线的方向图根据权值的变化而变化，它不同于
全向(omni-天线(理想时为一直线，而更接近定向(directional天线的方向
图，即有主瓣(main lobe.副瓣(side lobe等，但相比而言，智能天线通常有
较窄的主瓣，较灵活的主，副瓣大小、位置关系和较大的天线增益(无线术语，
天线的一项重要指标，是最强大力向的增益与各力向平均增益之比。它和固定
天线的最大区别是：不同的权值对应着不同的力向图，可以通过改变权值来调节
天线方向图，即天线模式(antenna Pattern。理想的智能天线就是要使天线方向图的主瓣对准目标用户方向，零瓣对准干扰信号方向，如下图所示
:
图2-2 理想智能天线工作的方向图
智能天线的方向图是随着权值的改变而动态变化的，智能天线正是通过自适应调整权值的来抑制干扰、提高信噪比，进而提高移动通信系统性能。
2.3.2 智能天线系统的基本结构
通常智能天线系统由3部分组成:实现信号空间采样的天线阵; 对各阵元输出进行加权合并的波束成型网络:更新合并权值的控制部分，其基本结构如图2-3所示。
天线阵列部分根据天线阵元之间的几何关系，阵列形状大致可划分为:线阵、面阵、圆阵等，甚至还可以组成三角阵、不规则阵和随机阵。天线阵的配置方式对智能天线性能有着直接的影响。在移动通信应用中天线阵多采用均匀线阵或均匀圆阵。因为直线阵己被证明更适合于市区移动通信环境，本文的讨论都是以直线阵为例。天线阵元数一般取4到16. 因为一方面天线阵元数越多，系统增益也就越高；但另一方面阵元数的增加会使射频通道相应增加，会导致基站成本上升过大，所以
智能天线的天线数不能过大。阵元间距一般为半个波长，因为如果阵元间距过大，接收信号的彼此相关程度会降低; 间距过小，会在天线的方向图上形
成不必要的栅瓣(有较大甚至和主瓣高度相同的旁瓣。
图2-3 典型的智能天线系统结构示意图
波束成型网络部分主要完成数模转换和天线方向图的自适应调整。每个天线阵阵元上都有ADC 和UAC ，将接收到的模拟信号转换为数字信号，将待发射的数字信号转换为模拟信号，完成模拟信号和数字信号的相互转换。所有收发数字信号都通过一组高速数字总线和基带数字信号处理器连接。天线方向图的调整是根据控制部分得到的权值调节天线输出来实现的。
控制部分(即算法部分是智能天线系统的核心部分，其功能是依据信号环境、按某种性能度量准则和自适应算法，选择或计算权值。
智能天线系统是由上面三部分组成的一个自适应控制系统，它根据一定的自适应算法自动调准天线阵方向图，使它在干扰方向形成零陷或低陷，在信号到达方向形成主瓣，从而达到加强有用信号，拟制干扰信号的目的。
智能天线系统的基本工作流程可以简单概括如下:
1、系统将首先对来自所有天线中的信号进行快照或取样，然后将它们转换成数字形式，并存储在内存中。
2、处理器将立即分析样本，对无线环境进行评估，确认有用用户所在的位置。
3、处理器对天线信号的组合方式进行计算，力争最佳地恢复用户的信号。
4、系统将进行模拟计算，使天线阵列可以有选择地向空间发送信号。
5、在上述处理的基础上，系统就能够在每条空间信道上发送和接收信号，从而使这些信道成为双向信道。
以上介绍了智能天线的基本结构、原理以及工作流程，在工程实际中智能天线的工作原理更复杂，并且每一部分的实现和结构往往根据所应用的系统不同而略有不同。
2.4 智能天线的主要功能及应用
上一节介绍了智能天线的原理，移动通信中引入智能天线技术究竟能实现那些功能，它又有那些优势；移动通信系统依据多址方式和实现方法的不同可以分为不同的类别，智能天线在不同的移动通信系统中的应用情况又是怎么样呢? 这一节里主要介绍这两个方面的内容。
智能天线的功能以及在移动通信中应用的强大优势主要表现在以下几个面:
1.抗衰落
基于M A T L A B的智能天 线及仿真 This model paper was revised by the Standardization Office on December 10, 2020
摘要 随着移动通信技术的发展，与日俱增的移动用户数量和日趋丰富的移动增值服务，使无线通信的业务量迅速增加，无限电波有限的带宽远远满足不了通信业务需求的增长。另一方面，由于移动通信系统中的同频干扰和多址干扰的影响严重，更影响了无线电波带宽的利用率。并且无线环境的多变性和复杂性，使信号在无线传输过程中产生多径衰落和损耗。这些因素严重地限制了移动通信系统的容量和性能。因此为了适应通信技术的发展，迫切需要新技术的出现来解决这些问题。这样智能天线技术就应运而生。智能天线是近年来移动通信领域中的研究热点之一，应用智能天线技术可以很好地解决频率资源匮乏问题，可以有效地提高移动通信系统容量和服务质量。开展智能天线技术以及其中的一些关键技术研究对于智能天线在移动通信中的应用有着重要的理论和实际意义。 论文的研究工作是在MATLAB软件平台上实现的。首先介绍了智能天线技术的背景；其次介绍了智能天线的原理和相关概念，并对智能天线实现中的若干问题，包括:实现方式、性能度量准则、智能自适应算法等进行了分析和总结。着重探讨了基于MATLAB的智能天线的波达方向以及波束形成，阐述了music和capon两种求来波方向估计的方法，并对这两种算法进行了计算机仿真和算法性能分析; 关键字:智能天线；移动通信；自适应算法；来波方向； MUSIC算法 Abstract With development of mobile communication technology，mobile users and communication，increment service are increasing，this make wireless services increase so that bandwidth of wireless wave is unfit for development of communication，On the other hand，much serious Co-Channel Interruption and the Multiple Address interruption effect utilize rate of wireless wave’s bandwidth，so the transported signals are declined and wear down，All this has strong bad effect on the capacity and
14元阵列天线方向图及其MATLAB仿真
阵列天线方向图及其MATLAB 仿真 1设计目的 1.了解阵列天线的波束形成原理写出方向图函数 2.运用MATLAB 仿真阵列天线的方向图曲线 3.变换各参量观察曲线变化并分析参量间的关系 2设计原理 阵列天线：阵列天线是一类由不少于两个天线单元规则或随机排列并通过适当激励获得预定辐射特性的特殊天线。 阵列天线的辐射电磁场是组成该天线阵各单元辐射场的总和—矢量和由于各单元的位置和馈电电流的振幅和相位均可以独立调整，这就使阵列天线具有各种不同的功能，这些功能是单个天线无法实现的。 在本次设计中，讨论的是均匀直线阵天线。均匀直线阵是等间距，各振源电流幅度相等，而相位依次递增或递减的直线阵。均匀直线阵的方向图函数依据方向图乘积定理，等于元因子和阵因子的乘积。 二元阵辐射场： 式中： 类似二元阵的分析，可以得到N 元均匀直线振的辐射场： 令 ，可得到H 平面的归一化方向图函数，即阵因子的方向函数： 式中：ζφθψ+=cos sin kd 均匀直线阵最大值发生在0=ψ 处。由此可以得出 ])[,(212121ζθθθ?θj jkr jkr m e r e r e F E E E E --+=+=12 cos ),(21jkr m e F r E E -=ψ?θθζ φθψ+=cos sin kd ∑-=+-=10)cos sin (),(N i kd ji jkr m e e r F E E ζ?θθ?θ2πθ=)2/sin()2/sin(1)(ψψψN N A =kd m ζ?-=cos
这里有两种情况最为重要。 1.边射阵，即最大辐射方向垂直于阵轴方向，此时 ，在垂直于阵轴的方向上，各元观察点没有波程差，所以各元电流不需要有相位差。 2.端射振，计最大辐射方向在阵轴方向上，此时0=m ?或π，也就是说阵的各元电流沿阵轴方向依次超前或滞后kd 。 3设计过程 本次设计的天线为14元均匀直线阵天线，天线的参数为：d=λ/2，N=14相位滞后的端射振天线。基于MATLAB 可实现天线阵二维方向图和三维方向图的图形分析。 14元端射振天线H 面方向图的源程序为： a=linspace(0,2*pi); b=linspace(0,pi); f=sin((cos(a).*sin(b)-1)*(14/2)*pi)./(sin((cos(a).*sin(b)-1)*pi/2)*14); polar(a,f.*sin(b)); title('14元端射振的H 面方向图 ，d=/2,相位=滞后'); 得到的仿真结果如图所示： 14元端射振天线三维方向图的源程序为： y1=(f.*sin(a))'*cos(b); z1=(f.*sin(a))'*sin(b); x1=(f.*cos(a))'*ones(size(b)); surf(x1,y1,z1); 2 π?±=m
2.3.1 阵列几何图 天线阵可以是各种排列，下图所示分别为圆阵(UCA)、线阵(ULA)、矩形阵(URA)排列方式与空间来波方向关系图，为简化整列分析，假设阵元间不考虑耦合，L 为天线数目，天线间距相等且均为d ，为入射在阵列上的水平波达角，为垂直波达角。 图2- 1 阵列排列方式与空间来波方向的关系 1） 圆阵排列方式的天线响应矢量为： 011cos() cos() cos() cos() (,)[,,...,,...,]l L j j j j T U C A a e e e e ξ?ψξ?ψξ?ψξ?ψ θ?-----= 公 式2- 1 其中2/,0,1,...,1l l L l L ψπ==-为第l 天线阵元的方位角，sin(),w w k r k ξθ=为波 数 2） 线阵排列方式的天线响应矢量为： cos sin (1)cos sin (,)[1,,...,]w w jk d jk d L T U LA a e e ?θ ?θ θ?-= 公式2- 2 3） 矩形阵列方式的天线响应矢量为： (1)()[(1)] (1)[(1)(1)](,)(()())[1,,...,,,,... ,...,,...,] T jv j p v ju j u v u URA N p j u p v j N u j N u p v T a vec a u a v e e e e e e e θ?-++---+-== 公式2- 3 ,N P 分别为x ，y 方向的天线数目，这里设x y d d =， (1)()[1,,...,]ju j N u T N a u e e -=； cos sin w x u k d ?θ=； (1)()[1,,...,]jv j p v T p a v e e -=；
天线阵代码 一、 clc clear all f=3e9; N1=4;N2=8;N3=12; a=pi/2; %馈电相位差 i=1; %天线电流值 lambda=(3e8)/f; %lambda=c/f 波长 d=lambda/2; beta=2.*pi/lambda; W=-2*pi:0.001:2*pi; y1=sin((N1.*W./2))./(N1.*(sin(W./2))); %归一化阵因子 y1=abs(y1); r1=max(y1); y2=sin((N2.*W./2))./(N2.*(sin(W./2))); %归一化阵因子 y2=abs(y2); r2=max(y2); y3=sin((N3.*W./2))./(N3.*(sin(W./2))); %归一化阵因子 y3=abs(y3); r3=max(y3); %归一化阵因子绘图程序, figure(1) subplot(311);plot(W,y1) ; grid on; %绘出N=4等幅等矩阵列的归一化阵因子 xlabel('f=3GHz,N=4,d=1/2波长,a=π/2') subplot(312);plot(W,y2) ; grid on; %绘出N=8等幅等矩阵列的归一化阵因子 xlabel('f=3GHz,N=8,d=1/2波长,a=π/2') subplot(313);plot(W,y3) ; grid on; %绘出N=12等幅等矩阵列的归一化阵因子 xlabel('f=3GHz,N=12,d=1/2波长,a=π/2') %--------------------- %只有参数N改变的天线方向图 t=0:0.01:2*pi; W=a+(beta.*d.*cos(t)); z1=(N1/2).*(W);
天线阵代码 .pudn./downloads164/sourcecode/math/detail750575.htm l 一、 clc clear all f=3e9; N1=4;N2=8;N3=12; a=pi/2; %馈电相位差 i=1; %天线电流值 lambda=(3e8)/f; %lambda=c/f 波长 d=lambda/2; beta=2.*pi/lambda; W=-2*pi:0.001:2*pi; y1=sin((N1.*W./2))./(N1.*(sin(W./2))); %归一化阵因子 y1=abs(y1); r1=max(y1); y2=sin((N2.*W./2))./(N2.*(sin(W./2))); %归一化阵因子 y2=abs(y2); r2=max(y2); y3=sin((N3.*W./2))./(N3.*(sin(W./2))); %归一化阵因子 y3=abs(y3);
r3=max(y3); %归一化阵因子绘图程序, figure(1) subplot(311);plot(W,y1) ; grid on; %绘出N=4等幅等矩阵列的归一化阵因子 xlabel('f=3GHz,N=4,d=1/2波长,a=π/2') subplot(312);plot(W,y2) ; grid on; %绘出N=8等幅等矩阵列的归一化阵因子 xlabel('f=3GHz,N=8,d=1/2波长,a=π/2') subplot(313);plot(W,y3) ; grid on; %绘出N=12等幅等矩阵列的归一化阵因子 xlabel('f=3GHz,N=12,d=1/2波长,a=π/2') %--------------------- %只有参数N改变的天线方向图 t=0:0.01:2*pi; W=a+(beta.*d.*cos(t)); z1=(N1/2).*(W); z2=(1/2).*(W); W1=sin(z1)./(N1.*sin(z2)); %非归一化的阵因子K1 K1=abs(W1); %---------------------- W=a+(beta.*d.*cos(t));
微波技术与天线作业 电工1001，lvypf（12） 1、二元阵天线辐射图matlab实现 1)matlab程序: theta = 0 : .01*pi : 2*pi; %确定θ的范围 phi = 0 : .01*pi : 2*pi; %确定φ的范围 f = input('Input f(Ghz)='); %输入频率f c = 3*10^8; %常量c lambda = c / (f*10^9); %求波长λ k = (2*pi) / lambda; %求系数k d = input('Input d(m)='); %输入距离d zeta = input('Input ζ='); %输入方向系数ζ E_theta=abs(cos((pi/2)*cos(theta))/sin(theta))*abs(cos((k*d*sin(theta)+zeta)/2)); %二元阵的E面方向图函数 H_phi=abs(cos((k*d*cos(phi)+zeta)/2)); %二元阵的H面方向图函数 subplot(2,2,1); polar(theta,E_theta); title('F_E_θ') subplot(2,2,2); polar(phi,H_phi); title('F_H_φ'); subplot(2,2,3); plot(theta,E_theta); title('F_E_θ'); grid xlim([0,2*pi]) subplot(2,2,4); plot(phi,H_phi); grid xlim([0,2*pi]) title('F_H_φ');
2)测试数据生成的图形： a)f=2.4Ghz,d=lambda/2,ζ=0 图1，f=2.4Ghz,d=lambda/2,ζ=0 b)f=2.4Ghz,d=lambda/2,ζ=pi 图2，f=2.4Ghz,d=lambda/2,ζ=pi
手把手教你天线设计—— 用MATLAB仿真天线方向图 吴正琳 天线是一种变换器，它把传输线上传播的导行波，变换成在无界媒介（通常是自由空间）中传播的电磁波，或者进行相反的变换。在无线电设备中用来发射或接收电磁波的部件。无线电通信、广播、电视、雷达、导航、电子对抗、遥感、射电天文等工程系统，凡是利用电磁波来传递信息的，都依靠天线来进行工作。此外，在用电磁波传送能量方面，非信号的能量辐射也需要天线。一般天线都具有可逆性，即同一副天线既可用作发射天线，也可用作接收天线。同一天线作为发射或接收的基本特性参数是相同的。这就是天线的互易定理。天线的基本单元就是单元天线。 1、单元天线 对称振子是一种经典的、迄今为止使用最广泛的天线，单个半波对称振子可简单地单独立地使用或用作为抛物面天线的馈源，也可采用多个半波对称振子组成天线阵。两臂长度相等的振子叫做对称振子。每臂长度为四分之一波长、全长为二分之一波长的振子，称半波对称振子。 对称振子是一种经典的、迄今为止使用最广泛的天线，单个半波对称振子可简单地单独立地使用或用作为抛物面天线的馈源，也可采用多个半波对称振子组成天线阵。两臂长度相等的振子叫做对称振子。每臂长度为四分之一波长、全长为二分之一波长的振子，称半波对称振子。
1.1用MATLAB画半波振子天线方向图 主要是说明一下以下几点： 1、在Matlab中的极坐标画图的方法： polar(theta,rho,LineSpec)； theta：极坐标坐标系0-2*pi rho：满足极坐标的方程 LineSpec：画出线的颜色 2、在方向图的过程中如果rho不用abs(f)，在polar中只能画出正值。也就是说这时的方向图只剩下一半。 3、半波振子天线方向图归一化方程： Matlab程序： clear all lam=1000;%波长 k=2*pi./lam;
1.均匀线阵方向图 %8阵元均匀线阵方向图，来波方向为0度 clc; clear all; close all; imag=sqrt(-1); element_num=8;%阵元数为8 d_lamda=1/2;%阵元间距d与波长lamda的关系 theta=linspace(-pi/2,pi/2,200); theta0=45/180*pi;%来波方向(我觉得应该是天线阵的指向) %theta0=0;%来波方向 w=exp(imag*2*pi*d_lamda*sin(theta0)*[0:element_num-1]'); for j=1:length(theta) %(我认为是入射角度，即来波方向,计算阵列流形矩阵A) a=exp(imag*2*pi*d_lamda*sin(theta(j))*[0:element_num-1]'); p(j)=w'*a; %(matlab中的'默认为共轭转置，如果要计算转置为w.'*a) end figure; plot(theta,abs(p)),grid on xlabel('theta/radian') ylabel('amplitude') title('8阵元均匀线阵方向图') 见张小飞的书《阵列信号处理的理论和应用2.3.4节阵列的方向图》
当来波方向为45度时，仿真图如下： 8阵元均匀线阵方向图如下，来波方向为0度，20log（dB）
随着阵元数的增加，波束宽度变窄，分辨力提高：仿真图如下：
2.波束宽度与波达方向及阵元数的关系clc clear all close all ima=sqrt(-1);
阵列天线方向图的MATLAB 实现课程名称：MATLAB程序设计与应用任课教师：周金柱 班级：04091202 姓名：黄文平 学号：04091158 成绩：
阵列天线方向图的MATLAB 实现 摘要：天线的方向性是指电磁场辐射在空间的分布规律，文章以阵列天线的方向性因子F（θ，φ）为主要研究对象来分析均匀和非均匀直线阵天线的方向性。讨论了阵列天线方向图中主射方向和主瓣宽度随各参数变化的特点，借助M ATLAB绘制出天线方向性因子的二维和三维方向图，展示天线辐射场在空间的分布规律，表现辐射方向图的特点。 关键词：阵列天线;;方向图；MATLAB 前言： 天线是发射和接收电磁波的重要的无线电设备，没有天线也就没有无线电通信。不同用途的天线要求其有不同的方向性，阵列天线以其较强的方向性和较高的增益在工程实际中被广泛应用。因此，对阵列天线方向性分析在天线理论研究中占有重要地位。阵列天线方向性主要由方向性因子F（θ，φ）表征，但F（θ，φ）在远区场是一组复杂的函数，如果对它的认识和分析仅停留在公式中各参数的讨论上，很难理解阵列天线辐射场的空间分布规律[ 1 ]。MATLAB以其卓越的数值计算能力和强大的绘图功能，近年来被广泛应用在天线的分析和设计中。借助MATLAB可以绘制出阵列天线的二维和三维方向图，直观地从方向图中看出主射方向和主瓣宽度随各参数的变化情况，加深对阵列天线辐射场分布规律的理解。 1 均匀直线阵方向图分析 若天线阵中各个单元天线的类型和取向均相同，且以相等的间隔d 排列在一条直线上。且各单元天线的电流振幅均为I，相位依次滞后同一数值琢，那么，这种天线阵称为均匀直线式天线阵，如图1 所示[ 2 ]: 均匀直线阵归一化阵因子为[ 3 ]： Fn（θ，φ）是一个周期函数，所以除§= 0 时是阵因子的主瓣最大值外，§= ±2 mπ
阵列天线方向图及其MATLAB 仿真 1设计目的 1.了解阵列天线的波束形成原理写出方向图函数 2.运用MATLAB 仿真阵列天线的方向图曲线 3.变换各参量观察曲线变化并分析参量间的关系 2设计原理 阵列天线：阵列天线是一类由不少于两个天线单元规则或随机排列并通过适当激励获得预定辐射特性的特殊天线。 阵列天线的辐射电磁场是组成该天线阵各单元辐射场的总和—矢量和由于各单元的位置和馈电电流的振幅和相位均可以独立调整，这就使阵列天线具有各种不同的功能，这些功能是单个天线无法实现的。 在本次设计中，讨论的是均匀直线阵天线。均匀直线阵是等间距，各振源电流幅度相等，而相位依次递增或递减的直线阵。均匀直线阵的方向图函数依据方向图乘积定理，等于元因子和阵因子的乘积。 二元阵辐射场： 式中： 类似二元阵的分析，可以得到N 元均匀直线振的辐射场： 令 ，可得到H 平面的归一化方向图函数，即阵因子的方向函数： ])[,(212121ζθθθ?θj jkr jkr m e r e r e F E E E E --+=+=12 cos ),(21jkr m e F r E E -=ψ?θθζ φθψ+=cos sin kd ∑-=+-=10)cos sin (),(N i kd ji jkr m e e r F E E ζ?θθ?θ2 πθ=) 2/sin() 2/sin(1)(ψψψN N A =
式中：ζφθψ+=cos sin kd 均匀直线阵最大值发生在0=ψ 处。由此可以得出 这里有两种情况最为重要。 1.边射阵，即最大辐射方向垂直于阵轴方向，此时 ，在垂直于阵轴的方向上，各元观察点没有波程差，所以各元电流不需要有相位差。 2.端射振，计最大辐射方向在阵轴方向上，此时 0=m ?或π，也就是说阵的 各元电流沿阵轴方向依次超前或滞后kd 。 3设计过程 本次设计的天线为14元均匀直线阵天线，天线的参数为：d=λ/2，N=14相位滞后的端射振天线。基于MATLAB 可实现天线阵二维方向图和三维方向图的图形分析。 14元端射振天线H 面方向图的源程序为： a=linspace(0,2*pi); b=linspace(0,pi); f=sin((cos(a).*sin(b)-1)*(14/2)*pi)./(sin((cos(a).*sin(b)-1)*pi/2)*14); polar(a,f.*sin(b)); title('14元端射振的H 面方向图 ，d=/2,相位=滞后'); 得到的仿真结果如图所示： kd m ζ?-=cos 2π ?±=m
天线阵代码 http://pudn/downloads164/sourcecode/math/detail750575.html 一、 clc clear all f=3e9; N1=4;N2=8;N3=12; a=pi/2; %馈电相位差 i=1; %天线电流值 lambda=(3e8)/f; %lambda=c/f 波长 d=lambda/2; beta=2.*pi/lambda; W=-2*pi:0.001:2*pi; y1=sin((N1.*W./2))./(N1.*(sin(W./2))); %归一化阵因子 y1=abs(y1); r1=max(y1); y2=sin((N2.*W./2))./(N2.*(sin(W./2))); %归一化阵因子 y2=abs(y2); r2=max(y2); y3=sin((N3.*W./2))./(N3.*(sin(W./2))); %归一化阵因子 y3=abs(y3); r3=max(y3);
%归一化阵因子绘图程序, figure(1) subplot(311);plot(W,y1) ; grid on; %绘出N=4等幅等矩阵列的归一化阵因子 xlabel('f=3GHz,N=4,d=1/2波长,a=π/2') subplot(312);plot(W,y2) ; grid on; %绘出N=8等幅等矩阵列的归一化阵因子 xlabel('f=3GHz,N=8,d=1/2波长,a=π/2') subplot(313);plot(W,y3) ; grid on; %绘出N=12等幅等矩阵列的归一化阵因子 xlabel('f=3GHz,N=12,d=1/2波长,a=π/2') %只有参数N改变的天线方向图 t=0:0.01:2*pi; W=a+(beta.*d.*cos(t)); z1=(N1/2).*(W); z2=(1/2).*(W); W1=sin(z1)./(N1.*sin(z2)); %非归一化的阵因子K1 K1=abs(W1); W=a+(beta.*d.*cos(t)); z3=(N2/2).*(W); z4=(1/2).*(W); W2=sin(z3)./(N2.*sin(z4)); %非归一化的阵因子K2
山东科技大学通信工程系 创新性实验设计报告 实验项目名称Matlab仿真天线阵的方向图及半功率宽度 组长姓名学号 手机 Email 成员姓名＿＿＿＿学号＿＿ 成员姓名＿＿＿＿学号＿＿ 专业通信工程班级 09—3 指导教师及职称 开课学期 2011 至 2012 学年＿ 2 ＿学期 提交时间 2012 年 7 月 13 日 一、实验摘要
熟悉课本第183页直线等距（LES）天线阵列的基本原理； 熟悉常用的算法，包括最小均方计算（LES），递归最小平方计算（RLS）和恒模算法（CMA）。 能够利用matlab软件仿真天线阵列的方向图，并且利用matlab计算其半功率宽度。 二、实验目的 1. 了解只能天线的组成，并能对天线矩阵进行分析 2. 运用MATLAB仿真阵列天线的方向图曲线 3. 变换各参量观察曲线变化并分析参量间的关系 4. 利用matlab计算LES阵列的半波宽度。 三、实验设计方案 1、实验原理 1. 阵列天线：阵列天线是一类由不少于两个天线单元规则或随机排列并 通过适当激励获得预定辐射特性的特殊天线。 阵列天线的辐射电磁场是组成该天线阵各单元辐射场的总和—矢量和由 于各单元的位置和馈电电流的振幅和相位均可以独立调整，这就使阵列天 线具有各种不同的功能，这些功能是单个天线无法实现的。
2.方向图原理：对于单元数很多的天线阵，用解析方法计算阵的总方向图相当繁杂。假如一个多元天线阵能分解为几个相同的子阵，则可利用方向图相乘原理比较简单地求出天线阵的总方向图。一个可分解的多元天线阵的方向图，等于子阵的方向图乘上以子阵为单元 阵列天线 天线阵的方向图。这就是方向图相乘原理。一个复杂的天线阵可考虑多次分解，即先分解成大的子阵，这些子阵再分解为较小的子阵，直至得到单元数很少的简单子阵为止，然后再利用方向图相乘原理求得阵的总方向图。这种情况适应于单元是无方向性的条件，当单元以相同的取向排列并自身具有非均匀辐射的方向图时，则天线阵的总方向图应等于单元的方向图乘以阵的方向图。 2、实验内容 假设阵元间距为波长的一半，天线元数目5个，导引向量a（0,0）； （1）仿真天线阵方向图实验主程序 thita=linspace(0,2*pi; fiai=0; lamda=0.4; d=lamda./2; beta=2.*pi./lamda; thita1=pi/6;
clc clear all f=3e9; N1=4;N2=8;N3=12; a=pi/2; % 馈电相位差 i=1; % 天线电流值 lambda=(3e8)/f; %lambda=c/f d=lambda/2; beta=2 、 *pi/lambda; W=-2*pi:0 、 001:2*pi; y1=sin((N1 、 *W 、 /2)) 、 /(N1 、 y1=abs(y1); r1=max(y1); y2=sin((N2 、 *W 、 /2)) 、 /(N2 、 y2=abs(y2); r2=max(y2); y3=sin((N3 、 *W 、 /2)) 、 /(N3 、 y3=abs(y3); r3=max(y3); %归一化阵因子绘图程序 figure(1) subplot(311);plot(W,y1) ; grid on; % 阵因子 xlabel('f=3GHz,N=4,d=1/2 波长,a= n /2') subplot(312);plot(W,y2) ; grid on; % 阵因子 xlabel('f=3GHz,N=8,d=1/2 波长,a= n /2') subplot(313);plot(W,y3) ; grid on; % 化阵因子 xlabel('f=3GHz,N=12,d=1/2 波长,a= n /2') % --------------- %只有参数N 改变的天线方向 图 t=0:0 、0 1 :2*pi; W=a+(beta 、*d 、*cos(t)); z1=(N1/2) 、*(W); z2=(1/2) 、*(W); W1=sin(z1) 、/(N1 、*sin(z2)); % 非归一化的阵因子 K1 K1=abs(W1); % --------------- 天线阵代码 波长 *(sin(W 、/2))); % *(sin(W 、/2))); % *(sin(W 、/2))); % 归一化阵因子 归一化阵因子 归一化阵因子 绘出N=4等幅等矩阵列的归一化 绘出N=8等幅等矩阵列的归一化 绘出N=12等幅等矩阵列的归一
matlab仿真天线阵代码 天线阵代码 tail750575.html 一、 clc clear all f=3e9; N1=4;N2=8;N3=12; a=pi/2; %馈电相位差 i=1; %天线电流值 lambda=(3e8)/f; %lambda=c/f 波长 d=lambda/2; beta=2.*pi/lambda; W=-2*pi:0.001:2*pi; y1=sin((N1.*W./2))./(N1.*(sin(W./2))); %归一化阵因子y1=abs(y1); r1=max(y1); y2=sin((N2.*W./2))./(N2.*(sin(W./2))); %归一化阵因子y2=abs(y2); r2=max(y2); y3=sin((N3.*W./2))./(N3.*(sin(W./2))); %归一化阵因子y3=abs(y3); r3=max(y3);
%归一化阵因子绘图程序, figure(1) subplot(311);plot(W,y1) ; grid on; %绘出N=4等幅等矩阵列的归一化阵因子 xlabel('f=3GHz,N=4,d=1/2波长,a=π/2') subplot(312);plot(W,y2) ; grid on; %绘出N=8等幅等矩阵列的归一化阵因子 xlabel('f=3GHz,N=8,d=1/2波长,a=π/2') subplot(313);plot(W,y3) ; grid on; %绘出N=12等幅等矩阵列的归一 化阵因子 xlabel('f=3GHz,N=12,d=1/2波长,a=π/2') %--------------------- %只有参数N改变的天线方向图 t=0:0.01:2*pi; W=a+(beta.*d.*cos(t)); z1=(N1/2).*(W); z2=(1/2).*(W); W1=sin(z1)./(N1.*sin(z2)); %非归一化的阵因子K1 K1=abs(W1); %---------------------- W=a+(beta.*d.*cos(t)); z3=(N2/2).*(W); z4=(1/2).*(W); W2=sin(z3)./(N2.*sin(z4)); %非归一化的阵因子K2 K2=abs(W2); %------------------------- W=a+(beta.*d.*cos(t)); z5=(N3/2).*(W);
方向图 clear lamda=10;%自由空间的波长 n=1.125; l=n*lamda; k=2*pi/lamda;%自由空间的相移常数 theta0=[0.0001:0.1:360]; theta=theta0*pi/180; for i=1:length(theta0) fe(i)=abs((cos(k*l*cos(theta(i)))-cos(k*l))/sin(theta(i))); % fh(i)=1-cos(k*l); end % figure polar(theta,fe/max(fe)); xlabel('l=1.125λ'); % figure % polar(theta,fh/max(fh)); 辐射阻抗 clear lamda=10;%自由空间的波长 n=0.1:0.01:0.9; l=n*lamda; k=2*pi/lamda;%自由空间的相移常数 gama=0.5772156; for i=1:length(n) Rr1(i)=sin(2*k*l(i))*(Si(4*k*l(i))-2*Si(2*k*l(i))); Rr2(i)=cos(2*k*l(i))*(gama+log(k*l(i))+Ci(4*k*l(i))-2*Ci(2*k*l(i))); Rr3(i)=2*(gama+log(2*k*l(i))-Ci(2*k)); Rr(i)=(Rr1(i)+Rr2(i)+Rr3(i))*30; end plot(n,Rr) 阻抗特性 clc clear lamda=100;%自由空间的波长 n=0.1:0.01:0.9; l=n*lamda; a=[30,300,3000];%l和a的比值 k=2*pi/lamda;%自由空间的相移常数 gama=0.5772156; for j=1:length(a) for i=1:length(n) Rr1(i)=sin(2*k*l(i))*(Si(4*k*l(i))-2*Si(2*k*l(i))); Rr2(i)=cos(2*k*l(i))*(gama+log(k*l(i))+Ci(4*k*l(i))-2*Ci(2*k*l(i)));
北京邮电大学基于Matlab的MIMO通信系统仿真 专业：信息工程 班级：2011211126 姓名： 学号：
目录 一、概述 (1) 1、课题的研究背景 (1) 2、课程设计的研究目的 (1) 3、MIMO系统 (1) 【1】MIMO的三种主要技术 (1) 【2】MIMO系统的概述 (2) 【3】MIMO系统的信道模型 (2) 二、基本原理 (3) 1、基本流程 (3) 2、MIMO原理 (3) 3、空时块码 (4) 三、仿真设计 (5) 1、流程图 (5) 2、主要模块及参数 (5) 3、信源产生 (5) 4、信道编码 (6) 5、调制 (6) 6、AWGN信道 (6) 7、输出统计 (7) 四、程序块设计 (7) 1、代码 (7) 五、仿真结果分析 (11) 1、仿真图 (11) 2、结果分析 (12) 六、重点研究的问题 (12) 七、心得与体会 (12) 八、参考文献 (12)
一、概述 1、背景 MIMO 表示多输入多输出。在第四代移动通信技术标准中被广泛采用，例如IEEE 802.16e (Wimax)，长期演进(LTE)。在新一代无线局域网(WLAN)标准中，通常用于IEEE 802.11n，但也可以用于其他 802.11 技术。MIMO 有时被称作空间分集，因为它使用多空间通道传送和接收数据。只有站点（移动设备）或接入点（AP）支持 MIMO 时才能部署MIMO。 MIMO 技术可以显著克服信道的衰落，降低误码率。该技术的应用，使空间成为一种可以用于提高性能的资源，并能够增加无线系统的覆盖范围。 通常，多径要引起衰落，因而被视为有害因素。然而研究结果表明，对于MIMO系统来说，多径可以作为一个有利因素加以利用。MIMO系统在发射端和接收端均采用多天线(或阵列天线)和多通道，MIMO的多入多出是针对多径无线信道来说的。传输信息流s(k)经过空时编码形成N个信息子流ci(k)，I=1，……，N。这N个子流由N个天线发射出去，经空间信道后由M个接收天线接收。多天线接收机利用先进的空时编码处理能够分开并解码这些数据子流，从而实现最佳的处理。 2、课程设计的研究目的 设计一个基于Matlab的MIMO通信系统仿真。 3、MIMO系统 【1】 MIMO的三种主要技术 当前，MIMO技术主要通过3种方式来提升无线传输速率及品质： （1）空间复用（Spatial Multiplexing）: 系统将数据分割成多份，分别在发射端的多根天线上发射出去，接收端接收到多个数据的混合信号后，利用不同空间信道间独立的衰落特性，区分出这些并行的数据流。从而达到在相同的频率资源内获取更高数据速率的目的。 （2）传输分集技术，以空时编码（Space Time Coding）为代表： 在发射端对数据流进行联合编码以减小由于信道衰落和噪声所导致的符号错误率。空时编码通过在发射端增加信号的冗余度，使信号在接收端获得分集增益。 （3）波束成型（Beam Forming）： 系统通过多根天线产生一个具有指向性的波束，将信号能量集中在欲传输的方向，从而提升信号质量，并减少对其他用户的干扰。
天线阵代码 https://www.doczj.com/doc/ea16837287.html,/downloads164/sourcecode/math/de tail750575.html 一、 clc clear all f=3e9; N1=4;N2=8;N3=12; a=pi/2; %馈电相位差 i=1; %天线电流值 lambda=(3e8)/f; %lambda=c/f 波长 d=lambda/2; beta=2.*pi/lambda; W=-2*pi:0.001:2*pi; y1=sin((N1.*W./2))./(N1.*(sin(W./2))); %归一化阵因子 y1=abs(y1); r1=max(y1); y2=sin((N2.*W./2))./(N2.*(sin(W./2))); %归一化阵因子 y2=abs(y2); r2=max(y2); y3=sin((N3.*W./2))./(N3.*(sin(W./2))); %归一化阵因子 y3=abs(y3); r3=max(y3); %归一化阵因子绘图程序, figure(1) subplot(311);plot(W,y1) ; grid on; %绘出N=4等幅等矩阵列的归一化阵因子 xlabel('f=3GHz,N=4,d=1/2波长,a=π/2') subplot(312);plot(W,y2) ; grid on; %绘出N=8等幅等矩阵列的归一化阵因子 xlabel('f=3GHz,N=8,d=1/2波长,a=π/2') subplot(313);plot(W,y3) ; grid on; %绘出N=12等幅等矩阵列的归一化阵因子 xlabel('f=3GHz,N=12,d=1/2波长,a=π/2') %--------------------- %只有参数N改变的天线方向图 t=0:0.01:2*pi; W=a+(beta.*d.*cos(t)); z1=(N1/2).*(W); z2=(1/2).*(W); W1=sin(z1)./(N1.*sin(z2)); %非归一化的阵因子K1
阵列天线方向图及其MATLAB仿真一．实验目的 1.了解阵列天线的波束形成原理写出方向图函数 2.运用MATLAB仿真阵列天线的方向图曲线 3.变换各参量观察曲线变化并分析参量间的关系 二．实验原理 1.阵列天线：阵列天线是一类由不少于两个天线单元规则或随机排列并 通过适当激励获得预定辐射特性的特殊天线。 阵列天线的辐射电磁场是组成该天线阵各单元辐射场的总和—矢量和由于各单元的位置和馈电电流的振幅和相位均可以独立调整，这就使阵列天线具有各种不同的功能，这些功能是单个天线无法实现的。 ^ 2.方向图原理：对于单元数很多的天线阵，用解析方法计算阵的总方向图相当繁杂。假如一个多元天线阵能分解为几个相同的子阵，则可利用方向图相乘原理比较简单地求出天线阵的总方向图。一个可分解的多元天线阵的方向图，等于子阵的方向图乘上以子阵为单元 阵列天线 天线阵的方向图。这就是方向图相乘原理。一个复杂的天线阵可考虑多次分解，即先分解成大的子阵，这些子阵再分解为较小的子阵，直至得到单元数很少的简单子阵为止，然后再利用方向图相乘原理求得阵的总方向图。这种情况适应于单元是无方向性的条件，当单元以相同的取向排列并自身具有非均匀辐射的方向图时，则天线阵的总方向图应等于单元的方向图乘以阵的方向图。
三．源程序及相应的仿真图1.方向图随n变化的源程序 clear; sita=-pi/2::pi/2; lamda=; ] d=lamda/4; n1=20; beta=2*pi*d*sin(sita)/lamda; z11=(n1/2)*beta; z21=(1/2)*beta; f1=sin(z11)./(n1*sin(z21)); F1=abs(f1); figure(1); plot(sita,F1,'b'); hold on; n2=25; ： beta=2*pi*d*sin(sita)/lamda; z12=(n2/2)*beta; z22=(1/2)*beta; f2=sin(z12)./(n2*sin(z22)); F2=abs(f2); plot(sita,F2,'r'); hold on; n3=30; beta=2*pi*d*sin(sita)/lamda; z13=(n3/2)*beta; z23=(1/2)*beta; > f3=sin(z13)./(n3*sin(z23)); F3=abs(f3); plot(sita,F3,'k') hold off; grid on; xlabel('theta/radian'); ylabel('amplitude'); title('方向图与阵列个数的关系'); legend('n=20','n=25','n=30');
摘要 随着移动通信技术的发展，与日俱增的移动用户数量和日趋丰富的移动增值服务，使无线通信的业务量迅速增加，无限电波有限的带宽远远满足不了通信业务需求的增长。另一方面，由于移动通信系统中的同频干扰和多址干扰的影响严重，更影响了无线电波带宽的利用率。并且无线环境的多变性和复杂性，使信号在无线传输过程中产生多径衰落和损耗。这些因素严重地限制了移动通信系统的容量和性能。因此为了适应通信技术的发展，迫切需要新技术的出现来解决这些问题。这样智能天线技术就应运而生。智能天线是近年来移动通信领域中的研究热点之一，应用智能天线技术可以很好地解决频率资源匮乏问题，可以有效地提高移动通信系统容量和服务质量。开展智能天线技术以及其中的一些关键技术研究对于智能天线在移动通信中的应用有着重要的理论和实际意义。 论文的研究工作是在MATLAB 软件平台上实现的。首先介绍了智能天线技术的背景；其次介绍了智能天线的原理和相关概念，并对智能天线实现中的若干问题，包括:实现方式、性能度量准则、智能自适应算法等进行了分析和总结。着重探讨了基于MATLAB 的智能天线的波达方向以及波束形成，阐述了music 和capon 两种求来波方向估计的方法，并对这两种算法进行了计算机仿真和算法性能分析; 关键字:智能天线；移动通信；自适应算法；来波方向； MUSIC算法 Abstract With development of mobile communication technology ，mobile users and communication，increment service are increasing ，this make wireless services increase so that bandwidth of wireless wave is unfit for development of communication ，On the other hand ，much serious Co-Channel Interruption and the Multiple Address interruption effect utilize rate of wireless wave’s bandwidth，so the transported signals are declined and wear down，All this has strong bad effect on the capacity and performance of question and be fit for the development of communication ，so smart