三相异步电动机Matlab仿真

————————————————————————————————作者：————————————————————————————————日期：

中国石油大学胜利学院综合课程设计总结报告

题目：三相异步电机直接启动特性实验模型

学生姓名：潘伟鹏

系别：机械与电气工程系

专业年级： 2012级电气工程专业专升本2班

指导教师：王铭

2013年 6 月 27日

一、设计任务与要求

普通异步电动机直接起动电流达到额定电流的6--7倍，起动转矩能达到额定转矩的1.25倍以上。过高的温度、过快的加热速度、过大的温度梯度和电磁力，产生了极大的破坏力，缩短了定子线圈和转子铜条的使用寿命。但在电网条件和工艺条件允许的情况下，异步电动机也可以直接启动。本次课程设计通过MATLAB软件建模模拟三相异步电动机直接启动时的各个元器件上的电量变化。

参考：

电力系统matlab仿真类书籍

电机类教材

二、方案设计与论证

三相异步电动机直接起动就是利用开关或接触器将电动机的定子绕组直接接到具有额定电压的电网上。

由《电机学》知三相异步电动机的电磁转矩M与直流电动机的电磁转矩有相似的表达形式。它们都与电机结构（表现为转矩常数）和每级下磁通有关，只不过在三相异步电动机中不再是通过电枢的全部电流，而是点数电流的有功分量。三相异步电机电磁转矩的表达式为：

（1-1）式中——转矩常数

——每级下磁通

——转子功率因数

式(1-1)表明，转子通入电流后，与气隙磁场相互作用产生电磁力，因此，反映了电机中电流、磁场和作用力之间符合左手定则的物理关系，故称为机械特性的物理表达式。该表达式在分析电磁转矩与磁通、电流之间的关系时非常方便。

从三相异步电动机的转子等值电路可知，

（1-2）

（1-3）将式(1-2)、(1-3)代入(1-1)得：

（1-4）

一：我们做如下分析：

1.当s=0时，，

M=0，说明电动机的理想空载

转速为同步转速。

2.当s很小时，有

，，说明

电磁转矩T近似与s呈线性关系，即随着M的增加，略有下降。

因而，类似直流电动机的机械特性，是一条下倾的直线。

3.当s很大时，有，，说明电磁转矩M近

似与s成反比，即M增加时n反而升高。

4.当s=1时，n=0，=常数，此即三相异步电机的启动转

矩。

从上述可见，三相异步电动机的机械特性由两段组成：当s较小（n较高）时，n与M近似呈线性关系；当s较大（n较低）时，n随M增大而升高。将两部分机械特性圆滑连接，既得三相异步电动机机械特性，如图1-1所示。

由上述公式可得三相异步电动机直接起动时，起动电流较大，可达额定电流的 5～7 倍，同时会引起电网电压波动，该方法只适合于小容量轻负载的异步电

动机的起动。

二：异步电动机直接启动的电路图如图1-2：

图1-1 三相异步电动机机械特性

图1-2

三、电路设计与参数计算

MATLAB简介

八十年代以来，计算机仿真成为交流电机及其调速系统分析，研究和设计的有利工具。应用计算机的仿真技术，我们可以用软件建立起电机及其传动、控制的仿真模型，再以这个模型在计算机内人为模拟的环境或条件下的运行研究，替代真实电机在实际场合下的运行实验，既可得到可靠的数据，又节约了研究的时间及费用。MATLAB 语言具有以下特点：1功能强大、2界面友好、编程效率高3扩展性强

启动SIMULINK只需在MATLAB的命令窗口键人“SIMULNIK”命令，此时出现一个SIMULINK 窗口。这个窗口包含7个模块库，它们分别是信号源模块库(sources)、输出模块库(Sinks)、离散模块库(Diserete)、线性模块库(Linear)、非线性模块库(Nonlinear)、连接与接口模块库(Connections)和扩展模块库(Extrax)。

1、仿真电路设计和参数

应用MatLab/simulink 建立三相异步电动机直接起动仿真模型如图1 所示。主要包括三相异步电动机模块（Asynchronous Machine）、电源模块（Power Source）、选择器模块（Selector）和测量模块（Machinemeasurement）等。其中，异步电机模块参数设置如下：Pn =110 KW，Un =400 V，fn =50 Hz, Un =1 487 r/min, 定子电阻Rs=0.02155Ω，定子电感Ls=0.000 226H，转子电阻Rr=0.012 31Ω，转子电感Ls=0.000 226H，漏抗X1+X’2=0.254Ω,励磁互感L m=0.010

38H，转动惯量J=2.3kg·m，摩擦系数F=0.05421N·m·s，极对数p=2。

2、原理图（含仿真模块名称与型号）

图2三相异步电动机直接起动模型

3、元件清单

元件序号元件

符号

元件

型号

主要参数数量备注

1 U 3-Phase

Programmable

Voltage

Source

400v 1

2 M Asynchronous

Machine

SI Units

110kw 400v 1

3 V Voltage

Measurement

1 4 Constant 500 1

Machines Measurement Demux

Selector

ia 2

Scope

Gain

四、仿真过程与仿真结果

根据原件清单使用matalab6.5连线仿真结果如下

五、调试优化

（1）静态调试图如下，负载为零时。

转子电流

转子电流定子电流

转速

（2）动态调试图，带负载运行时

转子电

流

定子电

流

转速

六、仿真正确性验证

(1)转子电流验证

带入上述已知参数，利用M语言程序画出如下的转速与转子电流特性图；可得电流与转速的曲线关系Y/X(转子电流（A）/转速（r/min）)：

由上图曲线变化可得知电流随转速的增加而下降，异步电动机直接起动电流达到额定电流的6--7倍，当转速达到1500r/min时电流大约为300A,经过验证所得结果与理论值相同。

（2）机械转矩验证

将电机参数和数据代入电磁转矩公式，利用M语言程序画出如下的转速与机械转矩特性图；

Y/X(转速/机械转矩)

由图可知电机刚启动时起动转矩能达到额定转矩的1.25倍以上，与仿真图相符。

七、结论与心得

经过两个个星期的学习与仿真，使我在整个电机设计过程中，学到了不少东西，更深一步掌

握了MatLab在电机电路设计的原理，深入了解了我们专业的基础课程。

本次的实验目的是通过解决比较简单的实际问题，巩固和加深在《电机学》课程中所学的理论知识和实验技能。训练学生综合运用学过的电机及拖动

基础知识，在教师指导下完成查找资料，选择、论证方案，科学实验，分析结果等工作，通过理论联系实际提高和培养学生分析、解决实际问题的能力和创新能力，为后续课程的学习、毕业设计和毕业后的工作打下一定的基础。加强了理论知识与实践统一的能力，加强了自己动手操作的能力。与此同时，也培养了我们的团队精神，有助于我们更快的适应社会，适应工作。我相信通过这次课程设计会给我以后的学习生活中更大的帮助。最后感谢在这次训练中给予我帮助的的老师和同学。

八、参考文献

[1] 顾绳谷.电机及拖动基础[M].4版.北京：机械工业出版社，2007

[2] 李维，matlab在电气中的应用.北京：中国电力出版社，2007

[3] 汤蕴缪.电机学[M].北京：机械工业出版社，2000.

[4] 刘启新，电机与拖动基础北京中国电力出版社，2007.2

[5] 唐孝镐.实心砖子异步电机及其应用.机械工业出版社，1991.3

基于MATLAB的智能天线及仿真

基于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

认知无线电频谱切换源码matlab仿真

clear clc %rand('twister',1); blockpu=[]; blocksu=[]; for N=3:2:7 block=[]; for lambdap =0.01:0.05:0.5 %***************************************** %假设 1. CR网络和主网络(授权网络)共同存在于同一区域，并且使用同一频段。假设该频段共有N个信道，每个主用户或CR用户每次接入只占用一个信道。 % 若所有信道均被主用户占用，此时CR用户到达就被阻塞。若CR用户正在使用的信道有主用户出现，此时CR用户被迫中断，并进入缓存区排队等待 % 空闲可用信道以继续刚被中断的通信，若等待超过一定时限，则判定CR用户强制中断退离缓存区。 % 故共有三个队列，分别表示如下： % X队列——主用户队列，抢占优先，优先级最高 % Y队列——次用户队列，优先级最低 % Z队列——次用户切换队列，优先级次高，若在时延Tao内，则较次用户队列优先接入可用信道 % 2. 主用户和次用户的到达服从泊松分布，参数分别为lambdap和lambdas，平均服务时间服从参数为mup和mus的负指数分布 % 3. 对次用户而言，主用户抢占优先。总共有N个信道，也就是最多可以有N个主用户抢占所有信道， % 故Z队列的长度不会超过N，这里给定Z队列长度为N。 % 4. 假设初始状态所有N个信道均空闲，次用户理想感知，感知延时为0.005 %***************************************** % 2009年10月12日10月25日 %***************************************** %初始化 %***************************************** a = 100; %主用户数量 b = 100; %次用户数量 %N =3 %Z队列最大长度/总的信道数 %Tao=5 %切换时延门限Tao A = [ ]; %某主用户到达时刻占用信道序号的集合 B = [ ]; %某次用户到达时刻占用信道序号的集合 C = [ ]; %切换用户占用的当前所有信道序号集合 D = [ ]; %某次用户到达时刻主用户占用信道集合 member = [ ]; member_CR = [ ]; j1=1; %主用户参数*****************************************

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

用matlab 仿真不同天线阵列个天线的相关系数

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 -=；

泊松分布及其应用研究

泊松分布及其应用研究 Prepared on 22 November 2020

湖南科技大学信息与电气工程学院《课程论文》题目：泊松分布及其应用研究专业：通信工程班级： 13级3班姓名：黄夏妮学号：目录一、摘要 (1) 二、泊松分布的概念 (2) 三、计数过程为广义的泊松过程 (4) 四、泊松分布及泊松分布增量 (5) 五、泊松分布的特征 (5) 六、泊松分布的应用 (6) 七、基于MATLAB的泊松过程仿真 (8) 八、参考文献 (12)

摘要作为一种常见的离散型随机变量的分布,泊松分布日益显示其重要性,成为概率论中最重要的几个分布之一。服从泊松分布的随机变量是常见的,它常与时间单位的计数过程相联系。在现实生活中应用更为广泛,如数学建模、管理科学、运筹学及自然科学、概率论等等。并且在某些函数关系起着一种重要作用。例如线性的、指数的、三角函数的等等。同样, 在为观察现象构造确定性模型时, 某些概率分布也经常出现。泊松分布作为大量试验中稀有事件出现的频数的概率分布的数学模型, 它具有很多性质。为此本文讲述了泊松分布的一些性质, 并讨论了这些性质在实际生活中的重要作用。

二、泊松分布的概念：定义1 设随机变量X 的可能取值为,,2,1,0 且 {}0,,2,1,0,! >===-λλ k e k x k X P k 为常数。则称X 服从参数为λ的泊松分布,记作X ～ D(λ) 。定义2 设ε是任意一个随机变量,称 )t (- e t)(it +∞<<∞=Φε是ε的特征函数。主要结论：定理1 如果X 是一个具有以λ为参数的泊松分布,则E( X) = λ且D ( X) =λ。证明设X 是一随机变量,若 ] X) E( - X [ E{2}存在,则称它为X 的方差,记作D( X) ,即 ] X) E( - X [ E{ X) D(2}=。设X 服从泊松分布D ( X) ,即有: 则()()λλλλλλλλ λ=?=-==- ∞ =--∞ =-∑∑ e e k e k e k X E k k k k 11 0!1! 从而()() () λλλλλλλ λ +=-+-==-∞ =-∞ =--∞ =∑ ∑ ∑2122 2 2 !1!2! e k e k e k k X E k k k k k k 故λλλλ - X) E( - ) X E( X) D(2222=+== 定理2 设随机变量) , ,2 1 n ( x n =服从二项分布,其分布律为 {}n k p p C k x P k n n k n k n n ,,2,1,0,)1( =-==-。又设0>=λn np 是常数,则{}λλ-∞ →==e k k x P k n n ! lim 。证明由λ=n np 得: 显然,当k = 0 时,故λ-n e k} x P{→=。当k ≥1 且k → ∞时,有

MATLAB仿真天线阵代码

天线阵代码一、 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);

MATLAB仿真天线阵代码

天线阵代码 .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));

matlab仿真天线辐射图

微波技术与天线作业电工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的泊松分布的仿真

泊松过程样本轨道的MATLAB 仿真一、 Poisson Process 定义若有一个随机过程{:0}t N N t =≥是参数为λ>0的Poisson 过程，它满足下列条件： 1、0N = 0； 2、对任意的时间指标0s t ≤<，增量()()t s N N t s ω-ωλ(-)服从参数为泊松分布。 3、对任意的自然数n ≥2和任意的时间指标0120n t t t t =<<

手把手教你天线设计——用MATLAB仿真天线方向图

手把手教你天线设计—— 用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;

基于MATLAB的数字模拟仿真..

基于MATLAB的数字模拟仿真摘要：本文阐述了计算机模拟仿真在解决实际问题时的重要性，并较为系统的介绍了使用计算机仿真的原理及方法。对于计算机模拟仿真的三大类方法：蒙特卡罗法、连续系统模拟和离散事件系统模拟，在本文中均给出了与之对应的实例及基于MATLAB模拟仿真的相关程序，并通过实例深入的分析了计算机模拟解决实际问题的优势及不足。关键词：计算机模拟；仿真原理；数学模型；蒙特卡罗法；连续系统模拟；离散事件系统模拟在实际问题中，我们通常会面对一些带随机因素的复杂系统，用分析方法建模常常需要作许多简化假设，这样进行处理过后的模型与我们面临的实际问题可能相差很远，以致求解得到答案根本无法应用，这时，计算机模拟几乎成为唯一的选择。本文通过对计算机模拟仿真进行系统地介绍，寻求利用模拟仿真来解决问题的一般方法，并深入探讨了这些方法的长处和不足。我们定义一些具有特定的功能、相互之间以一定的规律联系的对象所组成的总体为一个系统，模拟就是利用物理的、数学的模型以系统为问题解决对象，来类比、模仿现实系统及其演变过程，以寻求过程规律的一种方法。模拟的基本思想是建立一个实验的模型，这个模型包含所研究系统的主要特点，这样做的目的就是通过对这个实验模型的运行，获得所要研究系统的必要信息。另外，系统的运行离不开算法，仿真算法是将系统模型转换成仿真模型的一类算法，在数字仿真模型中起核心和关键作用。 1、所谓计算机仿真计算机仿真是利用计算机对一个实际系统的结构和行为进行动态演示，以评价或预测该系统的行为效果。它是解决较复杂的实际问题的一条有效途径。针对一个确定的系统，根据运行的相似原理，利用计算机来逼真模仿研究对象（研究对象可以是真实的系统，也可以是设想中的系统），计算机仿真是将研究对象进行数学描述，建模编程，且在计算机中运行实现。对比于物理模拟通常花费较大、周期较长，且在物理模型上改变系统结构和系数都较困难的诸多缺陷，计算机模拟不怕破坏、易修改、可重用，有更强的系统适应能力。但是计算机模拟也有缺陷，比如受限于系统建模技术，即系统数学模型不易建立、程序调试复杂等。计算机仿真可以用于研制产品或设计系统的全过程中，包括方案论证、技术指标确定、设计分析、生产制造、试验测试、维护训练、故障处理等各个阶段。 2、计算机仿真的目的对于一个系统，是否选择进行计算机模拟的问题，基于判断计算机模拟与非计算机模拟方法孰优孰劣的问题。归纳以下运用计算机模拟的情况：（1）在一个实际系统还没有建立起来之前，要对系统的行为或结果进行分析研究时，计算机仿真是一种行之有效的方法。（2）在有些真实系统上做实验会影响系统的正常运行，这时进行计算机模拟就是为了避免给实际系统带来不必要的损失。如在生产中任意改变工艺参数可能会导致废品，在经济活动中随意将一个决策付诸行动可能会引起经济混乱。（3）当人是系统的一部分时，他的行为往往会影响实验的效果，这时运用系统进行仿真研究，就是为了排除人的主观因素的影响。

(完整版)波束形成-Matlab程序

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仿真

电子科技大学通信与信息工程学院标准实验报告实验名称：随机数的产生及统计特性分析

电子科技大学实验报告学生姓名：吴振国学号：2011019190006 指导教师：周宁实验室名称：通信系统实验室实验项目名称：随机数的产生及统计特性分析【实验内容】 1、编写MATLAB 程序，产生正态分布或均匀分布或二项分布或泊松分布或你感兴趣的分布的随机数，完成以下工作：（1）、测量该序列的均值，方差，并与理论值进行比较，测量其误差大小，改变序列长度观察结果变化；（2）、分析其直方图、概率密度函数及分布函数，并与理论分布进行比较；（3）、计算其相关函数，检验是否满足 Rx(0)=mu^2+sigma2，观察均值mu 为0和不为0时的图形变化；（4）、用变换法产生正态分布随机数，或用逆变换法产生其他分布随机数，（5）、重新完成以上内容，并与matlab 函数产生的随机数的结果进行比较。 2、已知随机信号：仿真M 个样本，估计其自相关函数和样本的功率谱(用自相关法和周期图法)，并利用样本估计序列X （n ）的功率谱。【实验原理】本实验采用matlab 实验方法进行实验，相关采样方法，作图方法等均在matlab 的学习中有过使用！下面不作具体介绍！【实验程序】 1.程序1： clear; sigma=1; mu=1; N=100; X=normrnd(sigma,1,1,N); average=sigma; variable=sigma^2; 1212()cos(80)4cos(200)(),,~[0,2],()~(0,1)X n t t N t U N t N πφπφφφπ=++++白噪声

阵列天线方向图的MATLAB实现

阵列天线方向图的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仿真

阵列天线方向图及其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

Poisson 泊松方程的差分方法matlab实现

Poisson 泊松方程的差分方法问题：设G 是如下图所示的十字形区域，由5 个相等的正方形构成。试用五点差分格式求解下面的Possion 问题：解法分析：原方程用五点差分格式写出来就变成了：

源代码： function F=fivepointdiff(l,n) h=l/n; N=2*(n-1)*n+(3*n-1)*(n-1); XY=zeros(2,N);%分割xy轴后每一个节点的坐标 for i=1:n for j=1:n-1 XY(:,(n-1)*(i-1)+j)=[l+j*h;i*h]; end end for i=1:n-1 for j=1:3*n-1 XY(:,n*(n-1)+(3*n-1)*(i-1)+j)=[j*h;l+i*h]; end end for i=1:n for j=1:n-1 XY(:,n*(n-1)+(3*n-1)*(n-1)+(n-1)*(i-1)+j)=[l+j*h;2*l+(i-1)*h]; end end A=zeros(N,N); for i=1:N for j=1:N

if(i==j) A(i,j)=4; else if(((XY(1,i)-XY(1,j))^2+(XY(2,i)-XY(2,j))^2)<2*h*h)%若是相邻点择系数为-1 A(i,j)=-1; end end end end f=zeros(N,1);%就是等号右边F for i=1:N f(i,1)=h*h; end U=bicg(A,f,0.1,100);%求解Au=F F=[XY;U'];%输出命令框中输入： fivepointdiff(1,25); x=ans(1,:); y=ans(2,:); z=ans(3,:); plot3(x,y,z) 得到的结果：

MATLAB仿真天线阵代码17页

天线阵代码 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仿真天线阵的方向图及半功率宽度概要

山东科技大学通信工程系创新性实验设计报告实验项目名称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;

MATLAB中关于常见的概率分布密度函数的语句及格式

MATLAB中关于常见的概率分布密度函数的语句及格式 MATLAB中关于常见的概率分布密度函数的语句及格式 normpdf(x,mu,sigma) 正态分布密度函数。 uifpdf(x,a,b) 均匀分布（连续）密度函数 exppdf(x,a) 指数分布密度函数 geopdf(x,p) 几何分布密度函数 binopdf(x,n,p) 二项分布密度函数 poisspdf9x,n) 泊松分布密度函数 unidpdf(x,n) 均匀分布（离散）密度函数 chi2pdf(x,3) X^2分布密度函数 fpdf(x,m,n) F分布密度函数 tpdf(x,n) t分布密度函数一、常见连续分布的密度函数MATLAB实现 1正态分布 x=-8:0.1:8; >> y=normpdf(x,0,1); >> figure(1);plot(x,y); >> grid on; >> y1=normpdf(x,1,2); >> figure(2);plot(x,y,x,y1,':') >> grid on;

2均匀分布 >> clear all >> x=-10:0.1:10; >> r=1; >> y=unifpdf(x,0,2*pi*r); >> plot(x,y,'r*'); >> grid on;

x=0:0.1:30; >> y=exppdf(x,4); >> plot(x,y,'m-.') >> grid on 二、常见离散分布的密度函数1几何分布 x=0:30; >> y=geopdf(x,0.5); >> plot(x,y,'bo') >> grid on