时隙 Aloha 及伪贝叶斯算法性能仿真
设一个时隙 Aloha 系统的时隙长度为 1,所有节点的数据包均等长且等于时隙长度。网络中的节点数为 m,各节点数据包以泊松过程到达。 1 假定每个节点的数据包到达强度均为 λ /m,在不同的 λ 下,仿真时隙Aloha 数据包传送的成功概率,绘制呼入强度和成功概率的曲线,和理论结果进行对照。 仿真思路: 1) 生成一个二项分布列来模拟数据包的到达过程 2) 因为数据包以泊松过程到达,所以二项分布的 P 定为(1- m eλ− ) 3) 对生成的数列求和,只有当其和恰等于 1 即有且仅有一个数据包到达时,才可以成功发送,这时成功个数计数+1 4) 2.选取合理的引,,qa,m,采用延时的下界,仿真时隙Aoha系统数据传输过程,统计在不冋同η下,到达率及离开率,绘制它们随n的分布情况,和理论值进行对照qn:等待重传的节点在每一时隙内重传数据包的概率qa:每个发送节点有新数据包到达的概率m:系统内总的节点数n:每个时隙开始时等待重传的节点数仿真思路:1)用二项分布模拟数据包的到达及发送过程2)生成两个数列:一个表示等待重传的节点以q,重传的情况;一个表示新到达的数据包情况因为题日说明采用延时的下界,即不缓冲,每个节点最多容纳一个数据包,有包则扔。所以第一个数列前n项令为1,后一个前n项令为0,之后两个数列可以进行简单加和3)发送成功率:对两个数列相加之后求和,如果sum等于1,说明此时隙内到达和发送的总数为1,只有在这种情况下发送才有可能成功,计数加1到达率:在每N次实验中,对“表示到达的数列”求和,统计4)对n做循环以表示到达率和离开率随n的变化情况;每个n下进行N次实验,数理统计3/8仿真结果0.40.35*0.30.250.20.150.10.0550607080901001/曲线为理论曲线:Ps=G exp(-G)and G=(m-n)a+n gr2/仿真值基木与理论曲线吻合在仿真的过程中,合理选取个参数值对能否得到埋想的曲线起了重要的作用下图分别为qr=0.02,0.05,0.08s时的曲线。可以看到,随着qr的增加,曲线向左移,导致第二个交叉点也左移,这个时候重传的延时将会减小。反之,曲线右移。当q,增加到一定程度的时候,系统只有一个稳定点了。4/840.350.30.20.1501020304050607080901c03仿真时隙Aoha系统下的伪贝叶斯算法,通过仿真结果眼正在n的估计误差较大的情况下的收敛特性及到达率小于1/e下的稳定性。仿真思路:1、伪贝叶斯算法的主要思路是对新数据包和积压节点等同对待:当有新数据包到达的时候,暂不发送,下一时刻与以前的积压节点一起以4r发送。所以修改2中的仿真模型:1)依旧是一列表到达,一列表上一时隙的积压节点2)对两列加和,统计其中为1的个数,设为d3)以qr为概率,d为长度,生成又一个二项分布数列 depart,表示发送的情况4)对depa求和,如果 depart的和为一,说明恰发送成功,n(k+1)=d-1,否则n(k+1)=n(k)5)循环,进行数理统计2、仿真收敛特性和稳定特性哩论值:根据给岀的伪贝叶斯算法的具体步骤,由给出的n(k),不断模拟生成n(k+1)5/8仿真值:由仿真模型及给出的n(k),生成n(k+1)观察两种方式得到曲线的走向3、给出不同的值,观察n(k+1)随时间变化的情况判断标准如果要保持系统的稳定,至少n(k+1)应该保持在一个恒定的状态,或者逐渐趋于零。如果n(k+1)不断增加,则系统最终将趋于饱和,无法再接纳新的数据包,此时系统不稳定。仿真结果:1、验证在n的佔计误差较大的情况下的收敛特性:1)n=170;估计nt=20;m=100:20.2:N=100016030040050060070080g001CC08002)n=50;估计nt=180;m=1000=-1-02:N=80结果说明可以看出,当估计值与系统本身的积压数据包数有很大差别的时候,无论是大还是小,最终都可以趋于实际值,从而收敛特性得到验证。1)同时可以看到,改变的值:当λ增大的时候,收敛地更快;2)当n不变的时候,改变m的值,如果n/m变大,那么发生碰撞的几率就变大,也会导致估计的n值更快地趋向理论n值这些都是于课堂分析的理论情况相吻合的6/82、验证系统的稳定性下图分别为A=02:=10.1:=1:4=1+02:=0.3时候的情况。可以看到,当λ
- 2020-12-09下载
- 积分:1
船舶动力定位参数辨识
对船舶数学模型的各个推进器参数进行系统辨识,具有实际价值李文华,等:船舶动力定位系统数学模型参数辨识方法研究针对动力定位技术的发展,我国研究人员也进表1离散时间摘要扩展卡尔曼滤波行了积极有益的探索。文献[]用固定增益的卡尔f(k+1)=F((k),(k)+vw(k)曼滤波估计低频运动,而高频运动则用一个参数模系统模型量测|/(k)~N(O,Q(k》)x(k)=H((k)();型来模拟,并用递推增广最小二乘法来估计参数,从u(k)-N(0,Q(k)而估计出船舶的髙频运动。通过控制计算和模拟试初始情况1(0)+5np(0)=验取得了良好的效果。文献[12]提出了水面舰船动力定位控制系统模型参数的离线最速下降寻优的状态估计传递|(k+1)=F((k),()误差协方差传递辨识方法,提高了动力定位系统研制过程的工作效P(h+1)=(h)P()(k)+r(kQ(k )r(k)率。文献[13]在建立船舶三维几何模型基础上,对K()=P(k)HT(k)H()(kH()+R()]-满载船舶从浅水40m到深水500m的水动力系状态估计更新)=()1)((数进行数值计算。利用三维线性势流理论在频域误差协方差更新P)=kk)H(FLK)H()y里研究船舶在浅水中的辐射问题,应用三维源汇分K()R(k)K(k)布法对不同水深下船舶运动的水动力系数,包括附定义φ(k)=0()JH(-)加质量和阻尼系数进行数值计算与分析,得出了有(k)H(k)=0(k)限深水域的附加质量和阻尼系数的渐进特性。文献14]考虑具有修正PM波谱的长峰不规则浪,基于尾部隧道式侧推m,艏部隧道式侧推,艄部方海浪幅值响应算子(RAO研究了船舶在海浪中的六角式推进器。质量阵M可利用文献9]里介绍自由度运动预报模型。为了有效地量化海洋环境对的 Strip Theory计算得到:动力定位船舶的作用,文献[5提出了海洋环境负从/1127400018902-00744载(包括风、海浪和海流)的建模方法,并运用00.07440.1278MATLAB的M文件和SIMUⅠNK分别编制了风干为了得到需要辨识的量需重复进行3项(每项扰力和力矩计算及随机海浪的仿真程序。在三级海2次,共6次)海上试验,以此提高参数估计器的收况下,实现了对海洋环境的仿真,得到了合理的仿真敛性和表现。具体如下结果。文献[1]考虑到船舶的动态特性存在固有的第1项:解耦了的纵荡运动。船舶仅依靠主螺强非线性以及非线性控制改善系统性能和鲁棒性的旋桨山和实现恒速前进,艏向通过艏侧推控制。能力,将非线性控制理论应用到船舶动力定位控制第2项:结耦了的横荡与艏摇运动。通过三个隧道系统的设计中,对某供应船的计算机模型进行仿真,式推进器砌、4、实现两次结耦了的横荡与艏摇运验证了非线性控制系统是有效的。文献[17]提出并动。第3项:在结耦的横荡与艏摇运动中得到方位验证了基于线性核函数在线支持向量回归的模型预角式推进器u的推力系数K6测控制方案。在线支持向量回归算法的引入可以通第1项是为了计算主螺旋桨的推力系数K1和过在线调整,确保预测模型的精确性。Xa,需要的输入量是X本文中X的计算方法是利22船舶数学模型参数辨识用文献[19]里介绍的切片法。第2项是为了计算结文献[18]讨论了使用两个并行测量序列来估计耦了的横荡与艏摇运动的参数数值,可以辨识出的动力定位船舶模型参数的离线并行扩展卡尔曼滤波向量为[ YYNNK3K5]第3项是为了计算全方器算法(O- line parallel extended Kalman filter位推进器的推力系数K6( EKF) Algorithm),见表1。最后采用一项以供给船使用动量方程来代替标准动力学方程,不仅可为对象的全尺度的海上试验来验证提出的参数估计以显著提高状态和参数估计器的性能还具有以下器的收敛性和鲁棒性。优点:实验对象以挪威ABB公司的“ Far Scandia”号供(1)增加数据冗余度;给船为原型。该船总长762m,船宽18:8m,型深(2)降低量测噪声;825m,吃水625m,净吨位4200t,主发动机功率(3)降低环境干扰;3533kW。推进器配置左右舷两个主推进器u1、l2,(4)增加数据记录长度第23卷第3期(总第135期)船羔vd.23N.3012年6月shiP boatJune, 2012(5)以对参数分批进行辨识等手段提高参数辨风。将风速分量定义为识的精确度。L=v,cos(ψ图2显示了实验辨识得到的A和。其中A包(5)W V sin(B-0)含的待求未知量[XyYM而R包含的待求末式中和v分别为风速在X轴和Y轴的分量;v知量是[kk2k3k4k5k6]。和月分别表示风速和风向。如图1所示。假设风速远大于船速,风在纵荡、横荡和艏摇方向的负荷向量可表述为As elements.pAcM(o)V, IV0.5p.A_C-(r )V,V6)0.SpA,Lo C(rm)VV.式中,风的相对角为y=ψp为空气密度,单位e号为kgm3;Lm为船舶总长,单位为m;V为相对风速,103K elements单位为kn;A-和A为正投影面积和侧投影面积,2015单位均为m2;C(y)Cn(y)和C(4)分别为纵荡横荡和艄摇方向的无因次风系数,是通过 Isherwood半经验公式得到的。00003,波浪扰动数学模型波浪干扰力一般分为两种:一种是一阶波浪干图2实验辨识得到的参数曲线扰力,也称高频波浪干扰力。这是在假设波浪为微幅波,未引起船舶大幅摇荡的情况下,认为船舶受到经实验辨识出的动量方程中的量:与波高成线性关系并且与波浪同频率的波浪力。另0.03180种是二阶波浪力,也称波浪漂移力该波浪力与波A000602006l8高平方成比例。0.0075_0.2454这种具有高频率小振幅振荡特性的波浪所产生K=103ding([93,93,20,2.0,28,26]的一阶波浪干扰力最主要是引发船舶的纵摇和垂荡经过计算公式D-M得到运动,对横摇的影响稍次之,而对横荡及艄摇运动的002820影响相对来说就小一些。至于具有慢时变特性的二00.0130475900.081419676阶波浪干扰力,本身同时又是非线性的,它仍然和波写成动力定位模式下的状态空间表达式为:浪的频率有关。波浪的二阶漂移力不但会改变船舶元=AU+Bx(4)疔的航向和航迹,尤其对于在锚泊状态下船舶位置的移动及钻井平台的动力定位系统的工作等均有式中A=MA4M,并且B=MTK。其数值表达式为:重要影响00318000.062800030下面介绍一种估算二阶波浪漂移力方法。19740.0046-0.2428年, Newman提出一种应用频域波浪漂移力系数的0008200082000估算方法。通过把波谱(通常选用PM谱)分为N0∞505-069000108等份,每份有相对应的波浪频率m和波幅A。这样波浪漂移力对横荡、纵荡、艏摇运动的作用力计算公3环境扰动数学模型式为131风扰动数学模型A,(T(W,B=-y)1(W+)风的作用可分为平缓变化的风和快速变化的李文华,等:船舶动力定位系统数学模型参数辨识方法研究[5] Fossen T L. Handbook of Marine Craft Hydrodynamics and式中,T()x0是频域波浪偏移力公式fB是平均Motion Control[M]. Wiley Sons Ltd, 2011: 81-83.波浪方向:是随机的相角。[6] Balchen J G, Jenssen N A, Saelid S Dynamic Positioning可以通过对本估算式进行改变,以避免在数值Using Kalman Filtering and Optimal Control Theory[C]/上产生无物理意义的高频分量。还可对本式进行扩Proceedings of IFAC/IFIP Symposium on Automation in展,用来包括波浪蔓延( wave spreading)。Offshore Oil Field Operation Norway 1976: 183-18633海流扰动数学模型[7]Balchen J G, Jenssen N A Mathisen E, et al. Dynamic作用在海上动力定位船舶上的海流具有方向和Positioning System Based on Kalman Filtering and OptimalControl[J]Modeling, Identification and ControL 1980, 1(3)速度的特征,研究中一般不考虑在大地坐标系下铅135-163垂方向运动。海流分为恒定流和潮汐流。恒定流一般[8] Strand JP, Fossen t inonlinear Passive Observer Design为固定方向和速度的海流,如洋流。潮汐流指海洋for Ships with Adaptive Wave Filtering, In: New Directions因为潮汐运动而引起的海水流动,其典型的表现为in Nonlinear Observer Design(Nijmeijer H, Fossen T L)海流方向的缓慢变化。但对于动力定位来说,海流[M].London: Springer-Verlag London Ld, 1999: 113-134的大小与方向可以认为是确定的,所以海流的模型[9] Guttorm t, Jerome J, Fosset I. Nonlinear Dynamic可以统一按照大小和方向恒定来确立。流的速度分Positioning of Ships with Gain-Scheduled Wave Filtering量表示为5:[C]//The Proceedings of 43rd IEEE Conference orL=V2cos(ψ)Decision and Control, Atlantis, Paradise Island, BahamasDecemher2004:5340-5347ve=y sin(8-n)式中:和v分别为流速在X轴和y轴的分量;V10 i Do K d. Global Robust and Adaptive Output FeedbackDynamic Positioning of Surface Ships[C]/The Proceedings和月分别代表流速和流向。如图1所示。of 2007 IEEE Internati在此没有考虑第摇方向的流速,而海流对水面Automation. Roma, April 2007: 10-14船舶的作用可以通过将各海流速度分量引人到船的1]王晓声船舶动力定位系统设计及试验研究门J国造运动方程中由相对速度向量v=[u-,-a,r丁体现。船,1991(3):12-21[12]边信黔,严渐平,施小成船舶动力定位系统参数辨识4结论方法的研究[J]船舶工程,19994):36-38[13]姜哲,石珦,王磊动力定位船舶水动力参数数值试验本文讨论了船舶及推进器动力学数学模型与船研究[门]实验室研究与搡索,2005(12):14-17.舶外界环境干扰因素数学模型的建模策略。通过对14]李文魁张博田蔚风等.一种波浪中的船舶动力定位已有研究方法的分析研究与总结,有助于建立适用运动建模方法研究[]仪器仪表学报,2007(6):1051于各种海况和操作模式的船舶动力定位系统非线性数学模型。[15]施小成王元慧船舶动力定位海洋环境的建模与仿真J,计算机仿真,2006(11):237-239[16]刘芙蓉陈辉基于非线性控制理论的船舶动力定位控[参考文献制系统的数学模型[〕船海工程,209(5):92-95[1]杜佳璐,张显库汪思源,等船舶动力定位系统的自适[17]邓志良,胡寿松,张军峰船舶动力定位系统的在线模应非线性控制器设计[ C]/proceedings of the2 g chinese型预测控制[门中国造船,2009(6):879Control Conference. Beijing, 2010: 585-589.[2]周利,王磊,陈恒动力定位控制系统研究[船海[18] Fossen T I.Identification of Dynamically Positioned Shipe[].Control Engineering Practice, Volume 4, Issue 3, March程,008,37(2)86-911996:369-376[3]马超庄亚锋陈俊英船舶动力定位系统技术[J中国[19] FaltinsenO M Sea Loads on Ships and Oishore Structures造船,2009,50(增刊):52-57[4]贾欣乐,杨盐生船舶运动数学模型机理建模与数学建[M].Cambridge University Press, 1990:41-45模[M]大连大连海事大学出版社,199:294-356船舶动力定位系统数学模型参数辨识方法研究旧WANFANG DATA文献链接作者:李文华,杜佳璐,张银东,宋健,孙玉清,陈海泉, LI Wen-hua, DU Jia-luZHANG Yin-dong, SONG Jian, SUN Yu-ging, CHEN Hai-quan作者单位李文华,张银东,宋健,孙玉清,陈海泉, LI Wen-hua, ZHANG Yin-dong, SONG Jian, suN Yu-qing, chen Hai-quan(大连海事大学轮机工程学院大连116026),杜佳璐, DU Jia-lu(大连海事大学信息科学技术学院大连116026)刊名:船舶英文刊名:Ship boat年,卷(期):2012,23(3)参考文献(19条1. Balchen J G; Jenssen N A; Mathisen E Dynamic Positioning System Based on Kalmon Filtering andOptimal Control 1980(03)2. Balchen J G; Jenssen N A; Saelid S Dynamic Positioning Using Kalman Filtering and Optimal ControlTheory 19763. Fossen T I Handbook of Marine Craft Hydrodynamics and Motion Control 20114贾欣乐;杨盐生船舶运动数学模型机理建模与数学建模19995.马超;庄亚锋;陈俊英船舶动力定位系统技术2009(增刊)6.周利;王磊;陈恒动力定位控制系统研究[期刊论文]船海工程2008(02)7. Faltinsen 0 M Sea Loads on Ships and Offshore Structures 19908. Fossen t I Identification of Dynamically Positioned Ships 19969.邓志良;胡寿松;张军峰船舶动力定位系统的在线模型预测控制2009(06)10.刘芙蓉;陈辉基于非线性控制理论的船舶动力定位控制系统的数学模型[期刊论文]船海工程2009(05)11.施小成;王元慧船舶动力定位海洋环境的建模与仿真[期刊论文]计算机仿真2006(11)12.李文魁;张博;田蔚风一种波浪中的船舶动力定位运动建模方法硏究[期刊论文]仪器仪表学报2007(06)13.姜哲;石珣;王磊动力定位船舶水动力参数数值试验硏究[期刊论文]实验室硏究与探索2005(12)14.边信黔;严浙平;施小成船舶动力定位系统参数辨识方法的硏究[期刊论文]船舶工程1999(01)15.王晓声船舶动力定位系统设计及试验研究1991(03)Do K d Global robust and Adaptive Output Feedback Dynamic Positioning of Surface Ships 200717. Guttorm T; Jer(o)me J; Fossen T I Nonlinear Dynamic Positioning of Ships with Gain-Scheduled WaveFiltering 200418. Strand J P; Fossen T I Nonlinear Passive Observer Design for Ships with Adaptive Wave Filtering19.杜佳璐;张显库;汪思源船舶动力定位系统的自适应非线性控制器设计2010本文链接http://d.g.wanfangdata.comcn/periodiCalcb201203011.aspx
- 2020-12-03下载
- 积分:1
利用Matlab提取图片中的数据
从事科研或者工程的人员在文档撰写过程中,常需要将文献中的曲线与自己的结果对比,为获取原始数据,一种常用的办法是手动描点,即将原始曲线放大然后打印出来,选取一定数量的点,读出其横纵坐标,然后重绘。对于较为平坦的曲线,这种方法当然可行,但当曲线数量增加,曲线变化复杂,这种方法工作量可想而知。前段时间由于原始数据丢失,仅剩几十幅图片,本人尝试过手动描点,经历几个小时奋战,实在无法继续,索性转向matlab,借助其强大的数据处理能力,编写了两个GUI的小软件image2data、data_poly提取数据,如今大功告成,遂于大家分享。坐标轴标定按下面板上的按钮进行图像放大,按钮恢复初始显示大小,钮采用鼠标拖动图像按钮退出放大或者拖动的鼠标操作模式,空格键表示取点操作,键表示删除上一次取点操作,状态栏的显示当前鼠标取点总数目(注意,初始点数为,然后存处个坐标轴标定坐标,剩余用来存储曲线坐标)。值得指山的是,每次放大或者拖动桨作完毕后,必须按下按钮,才能用空格键进行取点操作。匹回图00.400.351 MHz足0.300.2530 MHz0.20frequency>100 MHz型30681012141618REVERSE VOLTAGE (VFigure 1. RF Capacitance vs ReverseBias. HSMP-3810 Series图坐标轴定标曲线描点按照上述操作反复进行图像放大、拖动、取点,状态栏的和用于显示当前坐标(注紊,这个坐标图像坐标,轴方向向下,后面坐标变换需要考虑),下图给出了描点完毕后的曲线,可以看出取点基本代表了曲线的全部信息。0.45E0400351 MHz0.3002530 MHZ.frequency> 100 MHzpr/。no0.1568101214161820REVERSE VOLTAGE (V)Figure 1. RF Capacitance vs Reversem[1Bias. hsMP-3810 Series图曲线描点数据处理及存储。按下按钮,描点后的曲线会显示在图像当中,按下按钮,程序自动进行坐标转换,得到所有描点的真实坐标,按卜按钮,使会生成一个文木,数据记录其和按钮下都有文本输入框,本别代表输入文本和存储文本的名称,不带后缀)。045"z55sd0.351M3002530 MHafrequency>100 MHz0.1502468101214161820REVERSE VOLTAGE (VFigure l. RF Capacitance vs Reversese4Bias. HsMP-3810 Series图拟合曲线效果记事本巴回囟文件〕编辑巴)格式迫!查看y1帮助)6,6ag15-92h432562gP-6811.18153B60-91488g2c-8912.15288c-E914.81Q68c-691.93177一Bg1571882c-6E-8Ube-出14.168y!:e-Uk1面=363211g-9g1py9py6-6817.99日969e-B914。394g55e-919.6696599e-B913.941218e-0911.119913-Bg3.E786699g-6511.36361243.779192-6811,5邵3777e十3.71813e-6117699899e+053.6164171e-6911_97了5775巴十B035小B725-6612-3111uB!g318286-B12.了7阳7P??5"24"9236P-R13.12096E5e+gs2502的9e-6613.526859e-E日31756187e-651391427562B9.1158646e-6814_287193.806E87ge-651h。47785F7eB98,38E17e-61↓816F877e+2.9877229e-6615.342g3B6e-Eg2.92049e-6615,6428g1eB9g2.B86ge6168697eBgg28616e-816652032cE9g78786ggc-691图数据记录文本数据后处理由于以上数据是于动选取的,故分布不够均匀,下面我们通过数据拟合,然后重新采样得到等间距的数据(可能大家会问,为什么两个功能不做在起呢?数据拟合是个比较味烦的事情,本组曲线采用多项式拟合即可,可对于更多的曲线采用指数函数、正弦函数等才能得到比较好的结果中的工具箱就包含了很多的拟合函数,为避免重复工作,仅绵写了这个小软件用来数据拟合再采样,其他的拟合就靠了)。运行代表多项式拟合的阶数(一般就够用了),代表重新釆样的数据个数,其他几个就不用解释了,默认输入文本为输出文本为,数据拟合结果如下图所示Xmax 23 ymax045fing图数据拟合再采样结果至此,数据提取过程完毕,可能操作上有些不便(毕竞只是小T具而已),但比起于动描点的速度和精度,可谓小巫见大巫。软件编写要点这两个小软件从构思到完成大概用了两大,功能的完成绝大部分归功」丰富的函数库和方便的,其中用到的主要函数上要有列衣如下,更详的介绍请参阅文档。表主要函数列表除了上述函数的掌握之外,还需要对的数据结构和函数响应有一定的理解,在此就不多讲了(多看相关例程就明白了)。编写程序之前,首先心中要有一个框架,做些什么,怎么做,顺序如何等。本软件的结构如下图所示:导入图像创建数据结构坐枟轴定标曲线描点坐标变換数据后处彐:拟合、重新采样、存储图程序沇程示意导入图像程序段:创建数据结构程序段坐标轴定标及坐标变换:数据后处理代码段附录3.1描点数据列表数据拟合结果3.3、重新采样数据列表
- 2020-12-04下载
- 积分:1