三维曲线拟合函数的MATLAB程序

matlab 矩形线圈磁场仿真，可实现图形化分析

基于Matlab编写压缩感知重建算法集，包括OMP，CoSaMP，IHT，IRLS，GBP，SP和ROMP.-Matlab codes for CS recovery algorithms, including OMP, CoSaMP, IHT, IRLS, GBP, SP and ROMP.

数字电子技术基础课后习题答案，第五版，张锁良编写

【实例简介】利用Clarke的SoS模型建立瑞利衰落信号，并且对信号的包络概率密度函数做了分析，并求了其自相关函数和多普勒功率谱

VINS 的功能模块可包括五个部分:数据预处理、初始化、后端非线性优化、闭环检测及闭环优化。代码中主要开启了四个线程,分别是:前端图像跟踪、后端非线性优化(其中初始化和 IMU 预积分在这个线程中)、闭环检测、闭环优化。、总体框架Measurement PreprocessingInitializationCamera(30hz)Feature Detectionnd rackerVisual-lnertialInitializedis- onlySfMAlignmentIMU (100hMU Pre-integrationLocal Visual-Inertial: OldestSliting WindowNewestNonli+、 Keyframe?OptimizationBundle Adjustment II Loop detectionwith RelocalizationStates from Loop ClosureFealure retrievel oop Deleted二二1---11------22===Global Pose Graph4-DoF Pose Graph OptimizationKeyframe DatabaseOptimization图1VINS框架ⅵINS的玏能模块可包括五个部分:数据预处理、初始化、后端非线性优化、闭环检测及闭环优化。代码中主要开启了四个线稈,分别是:前端图像跟踪、后端非线性伉化(其中初始化和IMU预积分在这个线程中)、闭环检测、闭环优化各个功能模块的作用上要有:1.I图像和MU预处理●图像:提取图像 Harris角点,利用金字塔光流跟踪相邻帧,通过 RANSAC去除异常点,最后将跟踪到的特征点push到图像队列中,并通知后端进行处理●IU:将IMU数据进行积分,得到当前时刻的位置、速度和旋转(PVQ),同时计算在后端优化中将用到的相邻帧的预积分增量,及预积分误差的 Jacobian矩阵和协方差项。1.2初始化首先,利用SFM进行纯视觉佔计滑窗內所有帧的位姿及3D点逆深度,最后与IMU预积分进行对齐求解初始化参数1.3后端滑窗优化将视觉约束、IMU约束和闭环约束放在·个大的目标函数中进行非线性优化,求解滑窗内所有帧的PVQ、bias等。L M States in the sliding windowIMU:k States from loop clos1Camera:冷 MU measurements>visual measurements★ Catur图2滑窗优化示意图14闭环检测和优化利用D)BoW进行闭环检测,当检测成功后进行重定位,最后对整个相机轨迹进行闭环优化。U预积分VisionIMUVision图3MU预积分示意图21当前时刻pVQ的连续形式将第k唢和第kl帧之间的所有IMU进行积分,可得第kHI帧的位置、速度和旋转(PVQ),作为视觉估计的初始值,这里的旋转采用的四元数。v△t+k+1∈[k,k+1]rW(at-ba ) -owletbk JtE[k, k+1]n(,-bdt∈[k,k+1]其中,a2和O为ⅠMU测量的加速度和角速度,是在Body自身坐标系, world坐标系是IMU所在的惯导系,上式的旋转公式推导可参考附录10.1。22当前时刻PVQ的中值法离散形式公式(1)给出的是连续吋刻的相机当前PVR的达代公式,为了跟代码致,下面给出基于中值法的公式,这与 Estimator:; processIMg(O函数中的Ps]、Rs]和Vs是一致的,IMU积分出来的第j时刻的物理量可以作为第j帧图像的初始值。tr t+a26t(2)ka,St其中q(a1-ba)-g"+q:+1(a+1-ba)(a;+o;+1)2.3两帧之间PVQ增量的连续形式通过观察公式(1)可知,IvU的预积分需要依赖与第k帧的ν和R,当我们在后端进行非线性优化时,需要迭代更新第κ唢的ν和R,这将导致我们需要根据每次迭代后值重新进行积分,这将非常耗吋。因此,我们考虑将优化变量从第k帧到第κ+1帧的IU预积分项中分离开来,通过对公式(1)左右两侧各乘Rb,可化简为:R(+p2k-=2△)+ak+1b其中DtElk, k+1t∈[k,k+1R k(at-bar)ldt)Wendtt∈[kk+1这样我们就得到了连续时刻的MU预积分公式,可以发现,上式得到的MU预积分的值只与不同时刻的a2和o相关。这里我们需要重新讨论下公式(5)预积分公式,以ab,为例,我们发现它是与MU的bias相关的,而bias也是我们需要优化的变量,这将导致的问题是,当每次迭代时,我们得到一个新的bias,又得根据公式(巧5)重新对第k帧和第k+1帧之间的IMU预积分,非常耗时。这里假设预积分的变化量与bias是线性关系,可以写成:ab,+/6n6ba+/16b+8 8ba +p(6)k+1sb24两帧之间PVQ增量的欧拉法离散形式面给出离散时刻的IMU预积分公式,首先按照论文中采用的欧拉法,给出第i个MU时刻与第i1个IMU时刻的变量关系为b+1k+的6t+元R(P)(1+R(P")(a2-bbn)δt25两帧之间PⅤQ增量的中值法离散形式卜面给出代码中采用的基」中值法的IMU预积分公式,这与 Estimator: processIMUO函数中的 Integration Base: push backo上是一致的。注意这里跟公式(2)是不一样的,这里积分出来的是前后两顿之间的IU增量信息,而公式(2)给出的当前帧时刻的物理量信息+1+B k St +=a, &tbb+1Bi + au1其中a,=slqilai-bai)+qiDi t aitl2.6连续形式下PVQ增量的误差、协方差及 JacobianIMU在每个吋刻积分出来的值是有误差的,下面我们对误差进行分析。首先我们直接给出在t时刻误差项的导数为:sa00016a000000-82(066hkR;0006|=00-(a-bh)0-1192k|+|000|mLL000016ba00101n000018b000F+ozk+ Gt其中:F25×15,G215×2,62x1,n12×,上式推导可参考附录102。下面我们讨论它的作用,将其可以简写为:6之k=F62z+Gtnt根据导数定义可知:62b=1m24-6262+8=62+628t=(+F6t)6z+(Gt6t)nt(11)这里我们对公式(1)的IMU误差运动方程再说明,将上式和EKF对比可知,上式恰好给出了如EKF一般对非线性系统线性化的过程,这里的意义是表示下一个时刻的IMU测量误差与上一个时刻的成线性关系,这样我们根据当前时刻的值,可以预测出下一个时刻的均值和协方差,而公式(1)给出的是均值预测,协方差预测公式如下Pb+6=(1+Ft)P(+Fl6t)7+(G,t)Q(G18t)ot(12)上式给出了协方差的选代公式,初始值Pk=0。其中,Q为表示噪声项的对角协方差矩阵000003000另外根据(11)式可获得诀差项的 Jacobian的迭代公式:(I+F26t)(14)其中 Jacobian的初始值为bk=12.7离散形式的PVQ增量误差分析我们首先直接给出PVQ增量误差在离散形式下的矩阵形式,为了与代码一致,我们修改下变量顺序,这和代码中 midPointIntegration(函数是一致的。(但不知为何计算的V中与前四个噪声项相关的差个负号?)1t fo660f106t‖loeBk+1=0f211f20016bδb0[6b102001rnot000kRkotk+1(15006t0n0000δt其中,推导可参考附录10.3:stE(ak-ba)02-4B+1(kk+121k+1b.)6t|6t2(Rr+ rk+18t2Stn=71=Rk+1(a+1-b)6tWr+ wf1=Ik+11+Gb。)δt-Rn+1(ak+121=-2配+1Stl st21(RK+Ruts)4rula1RrotstStR+1(a1R,+114/+11t28离散形式的PVQ增量误差的 Jacobian和协方差将公式(15)简写为:k+1F15×158215×1+V15×13Q则 Jacobian的迭代公式为k+15×15=F/k(16)其中, Jacobian的初始值为/k=l。这里计算出来的k+1只是为了给后面提供对bias的acoblar。协方差的迭代公式为P+15×15=FPFr+vQv(17)其中,初始值P=0。Q为表示噪声项的对角协方差矩阵:00000000aa000Q18×180a00(18)000000三、后端非线性优化31状态向量状态向量共包括滑动窗口内的n+l1个所有相机的状态(包括位置、朝向、速度、加速度计bias和陀螺仪bias)、 Camera到IMU的外参、m+1个3D点的逆深度X=[xr=pw,vb bpc,q3.2目标函数吗+(喻,2)+2(19)其中三个残差项即误差项分别为边缘化的先验信息、IMU测量残差、视觉的重投影残差。三种残差都是用马氏距离表示。根据《十四讲》中高斯牛顿法,若要计算目标函数的最小值,可以理解为,当优化变量有一个增量后,目标函数值最小,以IU残差为例,可写成如下所示:nin lre2bk, X+8Xrk x)+HSⅩDk+1oXk+1k+1其中HB,为B关于 XIK Jacobian,将上式展开并令关于6X的导数为0,可得增量δx的计算公式:H k 8X=k+1TB那么,公式(28)可写成+∑+∑Tk∑1rc上式中,B为MU预积分噪声项的协方差,P为vual观测的噪声协方差。当MU的噪声协方差P越大时,其信息矩阵Pk,将越小,意味着该MU观测越不可信,换句话说,因MU噪声较大,越不可信IMU预积分数据,而更加相信 visual观测。注意,这里的IMU和vsua协方差的绝对值没有意义,因为考虑得是两者的相对性可将上式继续简化为:(Ap+AB +Acox=bp +bB +bc其中,Ap,AB和Ac为 Hessian矩阵,上述方程称之为增量方程。33MU约束1)残差:两帧之间的PVQ和bias的变化量的差△tx+k+1bk qbk+1bR+1 xyz+g"△t)-Bk(20)sbbbb其中各增量关于bias的 Jacobian可从公式(16)的大 Jacobian中的相应位置获得。上面与代码中 Integration base: evaluateD对应,2)优化变量pb, 0W, Svb ,8ba:,bor Opb,, 80W ,Swb,, bakr, Sba3)Jacobian:计算 Jacobian时,残差对应求偏导对象分别为p6e,6vB,6h,ba],6b,6b

基于XY2-100协议的振镜控制转换板的设计与实现,介绍了XY2-100协议，以及如何使用。《自动化与仪器仪表》2014年12期(总第182期)结束语2」李志洪计算机在办公自动化中的应用才智,2008(02)经济社会的发展促进办公自动化的发展,大数据时代使得[3]杨庆万计算机网络技术与办公自动化[广东科技,20071办公宰自动化对于计算机技术的使用提出越来越高的要求,企4宁长明,刘洪刚应用计算机网络系统实现办公自动化自动化技事业单位和办公人员都有责任共同努力,促进办公室自动化术与应用,2001(04)智能化,提高办公效率。办公室工作人员才是系统的主体部5余小惠计算机网络与办公自动化引进与咨询000分,要保持对新计算机技术的敏感性和可接受性。[6]唐世梅计算机在办公自动化中的应用[科技信息科学教研),2007(22)参考文献[7王民川.计算机在办公系统中综合运用的探索J中小企业管理与科技(下旬刊,2010(07)[l]龙玲.高校办公自动化应用教学的推广价值探究[小科技信息,2010[8]贺铿.大力推进统计系统办公自动化建设U中国统计,2001(08)(03)(上接第4页)图10表明输入控制量与光斑位置变化成线性关系,说明了输出数字量的准确性,振镜控制转换板可以有效地对振镜进行控制。5结束语本文在分析XY2-100协议的数据结构及接口定义的基础激光器振镜控制R32ARM上,结合DSP控制技术设计了振镜控制转换板,实现了对振镜振镜2204转换被的控制。从XY2-100协议入手,设计了DsPF2812的硬件电图9测试实验系统路,制作了振镜控制转换板。软件上通过对XY2-100的数据分扫描平面上的坐标值与两振镜转角的关系式为:析设计了软件流程,解决了协议的转换和与ARM11嵌入式系统通信的问题。通过实验表明设计的振镜转换板卡完成了协议y=d tane(1)转换的功能,与外界通信良好,实现了对振镜的X、Y轴方向的控制,取得了良好的控制效果。综上所述,这种高性能的振(Id 2ty+e)*tan日(2)镜控制转换板的设计具有重婆的应用价值。其中d为Y轴振镜到标记面的距离,X、Y轴振镜转角分别为,e为两个转镜之间的距离为l8mm。在距离坐标平面的距离参考文献d为235mm时,当Y轴输入控制量为0时,对应的坐标为]陈苗海中国激光加工产业现状和发展前景门激光与红外200434(2600mm,100mm),当输入的控制量为2时,对应的坐标为(1):73~77.(260.0mm,260.0mm),将坐标代人公式(1)计算得Y轴扫描[2]潘涌骆公序折射式扫描系统设计及应用应用激光,20123范围约为(-28°,+28°)。同理X轴扫描坐标范围:(260.0mm,[3] Newson Engineering. XY2-100 technical datasheet[EB/OL]. (2007-0331).http://www.new.be/rhotho1350mm),(l350mm,135.0mm),395.0mm,135.0mm)代人公式[4]柳宁基于MCU和DSP的运动控制研究硬件平台设计叮微计算(2,计算X轴扫描范围约为(-28°,+28°)。振镜控制转换板输机信息2006、22)11-2出的是20位的数据信号,其中起到控制转角作用的是16位的5]TMs320F28 x Product Data Sheet. Texas Instruments. SPRS74.2003数字量,范围是0-2的数字量,扫描范围约为28°,输入每隔6]TPs767D318 Product Data Sheet , Texas Instruments sLvS209H2004096个控制量振镜对应转角约为3.5°,在(-28,0)的范围[7]TMs32060003 oard Design for JTAG.Texas Instruments. SPRA584C,测试每转过35°与Y轴对应靶标上的位置。2002Y/mm[8]尹勇欧光军关荣锋DSP集成开发环境CCS使用指南[M北京:北京航空航天大学出版社,2003年聊入数字量50010000150002000250003000035000图10输出控制量与光斑在Y轴的位置万方数据

关于GPS信号如何产生、如何捕获、如何追踪的所有matlab程序大集合！

资源中包含神经网络、遗传算法、支持向量机、退火算法、粒子群算法等算法matlab源码、工具包及数据集。

1 数字调制系统2ASK与2FSK ………………………………………………… 1.1二进制幅度键控（2-ASK） ……………………………………………...　　1.2二进制幅度键控（2-FSK） ………………………………………………… 2 数字调制系统的仿真设计　　　　………………………………………... 2. 1　数字调制系统各个环节分析　……………………………………………. 2.2 仿真系统　　　　　　　　　…………………………………………. 2.2.1 　2ASK 　　　　　　　　………………………………………

Broadcom SDK 5.6 平台指南，关于SDK编译方法、步骤的编程手册，是学习了解Broadcom SDK的很有用的参考手册