技术规范,过几天发个压缩包,一次下载全部拥有 下载地址:https://download.csdn.net/download/appe1943/12006549压缩包包括以下内容:1.《建筑设计防火规范》GB50016;2.《小型火力发电厂设计规范》GB50049;3.《供配电系统设计规范》GB50052;4.《低压配电设计规范》GB50054;5.《工程建设标准强制性条文》(电力工程部分)6.《爆炸和火灾危险环境...
stk中天体坐标系的定义 stk中天体坐标系的定义中心天体坐标系统每个中心天体都有相应的坐标系统。坐标系原点均位于中心天体的质心,区别在于坐标轴的定义。每个中心天体均支持固定坐标系(Fixed coordinate system,即相对于大地无位移)、ICRF坐标系、J2000坐标系和惯性坐标系(其坐标轴相对于J2000坐标系为常值旋转)地球和月球具有最多的坐标系(月球固定坐标系的定义有不同的版本),太阳具有黄...
将UTCG转换为MJD(IAU) %--------------------------------------------------------------------------% Purpose:%% Modified Julian Date from calendar date and time%% Input(s):%% Year % Month% Day% Hour ...
IAU SOFA 软件包介绍 国际天文联合会(International Astronomical Union),简称IAU,其担负着讨论并制定国际上通用的天文标准的任务。常用的天文常数,有关地球自转轴运动的岁差,章动和极移等都是其制定的,并且被作为天文界中的标准。SOFA(Standards Of Fundamental Astronomy)软件包就是IAU编写的有关天文基础标准的子程序库。目前一共有161个子程序,主要...
GPS时与GPS周 GPS采用10bit定义域来对每个GPS时间信息中的周数进行编码,最多可以处理1024周(2^9+2^8+2^7+2^6+2^5+2^4+2^3+2^2+2^1+2^0=1023,0~2013),这个期间称为一个Epoch(历元)。在每个Epoch结束时,接收机将周数复位为第0周,然后再次开始计数。每个Epoch共1024周=7168天=19.62年。第一个GPS Epoch起始于GPS时198...
天文学中的时间单位的发展历程 Time ScalesA definition for the term "time scale":16th General Assembly of the International Astronomical UnionGrenoble, France (1976)Resolution No. 4 by Commissions 4 (Ephemerides) and 31 (Tim...
1620年-1972年历书时ET与世界时UT之间的关系 % (last updated 2018-12-14)% Delta-T 1620-1972% delta-T = ET - UT for the years 1620 - 1972% ET-UT, s%% Year +0.0 +1.0 +2.0 +3.0 +4.0% -----------------------------...
TAI、GPS与UTC时间的转换(以下数据经stk验证) % Table of time scales 1972-present, and some predictions% TT = TAI+32.184s ==> UT1-UTC = TAI-UTC - (TT-UT1) + 32.184s% Starting at TAI-UTC GPS-UTC TT-UT1 UT1-UTC%% 1972-01-01 ...
GPS时间与TAI时间的转换(经stk验证) %××××××××××××××××××××××××××××××××××××××××% GPS time = TAI - 19 seconds. GPS time matched UTC from 1980-01-01% to 1981-07-01. No leap seconds are inserted into GPS time, thus% GPS time...
天文学中的时间尺度单位 The different time scalesTAI = International Atomic Time (Temps Atomique International = TAI) is defined as the weighted average of the time kept by about 200 atomic clocks in over 50 n...
TAI与TDT的转换(经过stk检验) %××××××××××××××××××××××××××××××××××××××××××××××××% TDT = Terrestial Dynamical Time. Was used 1984-2000 as a time-scale of % ephemerides from the Earth's surface. TDT = TAI + 32.184. Repla...
儒略日的计算(经过stk验证) %儒略日的计算%输入日期:Gregorian TDT(TDTG)-----stk%输出日期:Julian Date(JDate)------stk%验证两组输入:%-----第一组输入:J2000.0-----%stk结果:2451544.99925713 JDate%程序计算结果:2451545 JDate%对应的UTCG时间为:1 Jan 2000 11:58:55.816 U...
关于matlab多目标优化函数fgoalattain的一点小提示 多个优化目标量值如果差的太大,会影响优化精度,例如要同时优化半长轴和偏心率,半长轴范围可能是几千千米,偏心率可能是千分乃至万分之几,量值差的就太大了,因此需要对量值范围进行处理,不要相差太大。...
matlab找不到编译器 装了两个版本的matlab:R2009a和R2014a,输入命令:mex -setup,结果都显示没有编译器。这个问题困扰了我一天,后来找到原因了:我的操作系统是64位的,安装matlab的时候自动装成了64位,自带的LCC编译器目前好像只支持32位的,装了Visual Studio 2008也找不到编译器。解决的办法是重新安装了32位的R2009a,再运行mex -setup命令,就出现...
matlab GUI界面自动更新? 如果运算时间短,matlab GUI界面还能保持之前的样子。像我在做优化计算时,需要长时间计算,此时如果切换到别的窗口再切换回来,界面就花了,必须等到给GUI中的对象重新赋值时界面才能恢复。为了克服这个问题,我现在的做法是经常给GUI控件赋值,但这样做效果不好。下一步准备引入定时器,隔段时间给界面中的某个控件赋值。控件可以是edit控件,将这个控件设置成透明的,赋值量是空格,这样应该能够实现界...
利用matlab fgoalattain函数实现多目标优化 之前利用matlab做多目标优化的时候查了很多书,没有找到matlab多目标优化函数,所以就对多目标进行组合加权变成了单目标,然后利用matlab单目标优化函数例如fmincon等进行优化。后来发现matlab里面集成了一个多目标优化函数,挺好用的,函数名字为fgoalattain,具体用法此处不具体介绍,有兴趣的可以help查询一下。...
在stk场景里面设置北京时 stk场景默认是UTCG时间,如果想改为北京时,可按照以下步骤设置:(1)首先启动stk,载入或者新建一个场景。为了方便后面介绍,此处场景名称为Scenario(2)双击Scenario,打开属性页面,点击左面的Units,将“DateFormat”设置为Gregorian LCL(LCLG)。此后场景中的对象、报表、图形等都以北京时方式呈现。...
伪谱法编程的两个例子,自己编写底层代码 以下两个例子底层代码均为自己编写,包括配点计算、拉格朗日插值等等,求解用的snopt工具箱。总体感觉,伪谱法编写非常麻烦。一、《基于Radau伪谱法的制导炸弹最优滑翔弹道研究_袁宴波》二、《最优控制问题的Legendre伪谱法求解及其应用_徐少兵》中的例1...
我回来了 好久没更新了,回来后看到很多朋友留言。有些代码我放在下载频道,当初都是免积分的,这次回来看到csdn改版后都需要积分了。有些代码是我帮一些朋友做的,就不能放到下载频道了。编算法不易,有些代码要收些钱,还请谅解。qq今后会用起来,2274579702,有需要的朋友就在上面联系吧...
伪谱法优化程序已经开发出来 可以用于弹道优化等场合。《最优控制问题的Legendre伪谱法求解及其应用_徐少兵》在编程时起到了很大的帮助作用。伪谱法程序编起来比我想象得要困难很多,编了整整两个月,克服了大大小小的困难。目前优化算法采用snopt优化工具箱,我觉得遗传算法等也可以,后续将进行试验。
分享snopt优化工具箱matlab版 下载地址:http://download.csdn.net/download/appe1943/9926928该工具箱为matlab版本,可以用于伪谱法优化。里面有几个例子,我已经调通了。经过测试,该版本在32位和64位windows下都可以运行,操作系统试验过xp和win 7.
航空航天领域的一些国内学术会议(非EI检索) 1、第四届高分辨率对地观测学术年会http://www.chreos.org/2、第二届中国空天安全会议http://saa.csu.edu.cn/info/1003/2619.htm3、第三十届全国空间探测学术研讨会征文通知http://www.cssr.org.cn/cms/readNewsViewlet.rt;jsessionid=DDA0B1FBCDB912
Isight unable to evaluate matlab commands的解决办法 这两天在利用Isight进行优化时遇到一个问题,困扰了我很久。用到了matlab模块,优化简单的matlab程序时可以直接运行,复杂的matlab程序就出现“Isight unable to evaluate matlab commands”的错误提示,但matlab软件已经打开了,那为什么不能执行matlab程序呢?matlab路径设置里面已经添加了当前目录,代码本身也没有问题,后来还是用这篇文
(源代码)四元数与欧拉角间的转换 clearclc vartheta_c=0.1;fea_c=0.1;gamma_c=0;q1=sin(fea_c/2)*sin(vartheta_c/2)*cos(gamma_c/2)+cos(fea_c/2)*cos(vartheta_c/2)*sin(gamma_c/2);q2=sin(fea_c/2)*cos(vartheta_c/2)*cos(gamma_c/2)+cos(f
(源代码)一类匹配不确定非线性系统的动态逆全程滑模变结构控制 参考文献:《一类匹配不确定非线性系统的动态逆全程滑模变结构控制》clear clcH=[5 2 5; 8 1 7; 1 1 2];invH=inv(H);a11=invH(1,1);a12=invH(1,2);a13=invH(1,3);a21=invH(2,1);a22=invH(2,2);a23=invH(2,3);a31=inv
(源代码)PWPF调制 1、仿真来源:《应用时标分离和动态逆方法设计飞行器的姿态控制系统》,武立军,贺有智,现代防御技术,2007clearclc[row,av,g]=OnAir(50000); %计算大气密度,当地声速,重力加速度m=44; %质量 Tay=60*g; %姿控发动机推力L=1.945; %抛盖后长度DD=0.37;%弹体直径Jx=0.3; %绕x轴转动惯量Jy=
(源代码)《应用时标分离和动态逆方法设计飞行器的姿态控制系统》 1、仿真来源:《应用时标分离和动态逆方法设计飞行器的姿态控制系统》,武立军,贺有智,现代防御技术,20072、在发动机推力的处理上,尝试了三种方法:连续推力、通过幅值限制得到的离散推力和通过PWPF调制得到的离散推力。连续推力很容易实现;后两者与飞行器的参数相关,发动机推力等参数需要仔细设置,由于原文没有给出飞行器的具体参数,发动机推力的大小与原文不同。一、离散推力
(源代码)最优控制与姿控喷流在导弹姿态控制中的应用 参考文献:《最优控制与姿控喷流在导弹姿态控制中的应用》%跟踪正弦结果可以看出,实际偏航角与指令值存在一定的跟踪误差%从仿真结果可以看出,时间最优bang-bang控制律是以牺牲跟踪精度来满足时间最短的要求clearclcA=2500;B=10.625;a=A/B;Mf=245;x1=0;%偏航角初始值x2=0;%偏航角速率初始值n=1;t=0;dt=0.001;
isight不能调用matlab的解决办法 最近安上isight后发现不能调用matlab,花了几天时间终于解决了,办法如下:双击matlab模块,在弹出来的菜单里选择options页面,调用方式选择java socket,然后点击browse按钮,导向matlab安装目录/bin/matlab.exe。关闭页面后重新运行,就可以调用matlab了。以前用isight 5.6.1的时候没遇到这个问题,isight自动识别了,不知
stk论文资料下载 从cnki下载的关于stk的论文,包括硕士论文、期刊论文,最近几年的基本包括全了,600M左右,如果从cnki上下载可是要花1K左右哦!http://download.csdn.net/detail/appe1943/9742899#comment
如何安装光驱 这两天从电脑城买回来了个DVD刻录光驱,安装的时候费了老鼻子劲,不知道怎么安装,把机箱电源卸了也装不上去,前面板那有铁片挡着。后来顿悟,把机箱前面板卸下来,就能够把铁片拆下来,然后把光驱塞进去就可以了。
iSIGHT中优化方法种类 iSIGHT中优化方法种类iSIGHT里面的优化方法大致可分为三类:1 数值优化方法数值优化方法通常假设设计空间是单峰值的,凸性的,连续的。iSIGHT中有以下几种:(1)外点罚函数法(EP): 外点罚函数法被广泛应用于约束优化问题。此方法非常很可靠, 通常能够在有最小值的情况下,相对容易地找到真正的目标值。外点罚函数法可以通过使罚函数的值达到无穷值,把设计变量从不可行
攻击时间控制的动态逆三维制导律源代码 %参考文献:《攻击时间控制的动态逆三维制导律》,哈尔滨工程大学学报,2010%编程难点:视线坐标系与参考坐标系间的转换关系clearclcV=300;% 导弹速度,m/sx=-5000;%进攻弹坐标,my=5000;z=6000;xt=0;%目标坐标,myt=0;zt=0;Rm=sqrt((x-xt)^2+(y-yt)^2+(z-zt)^2);DRm=0;theta_m=
带落角约束(变结构和比例)源代码 1、变结构制导律和比例导引均集成在同一个程序中,其表达式分别见daoyin.m第48行和53行;如果不需要某个导引律,请将其注释(前面加“%”)。2、比例导引和变结构导引结果图可以画在同一幅图中(每运行完一种导引律后,不要关闭图,继续运行下一种导引律)。clearclc%-----------导弹参数---------V_m=300;X_m=0;Y_m=3000;theta_
基于lmbp神经网络的效能评估(源代码) %参考论文:基于lmbp神经网络的效能评估clearclc%------------------训练数据------------------%Tin_train为输入数据,Tout_train为输出数据(由于原文没有定义输出数据,此处为自定义)Tin_train=[0.89 0.97 0.91 0.88 0.95 0.82 0.89 0.78 0.85 0.69 0.92 0.87 0.
四元数与欧拉角间的转换(源代码) clearclc vartheta_c=0.1;fea_c=0.1;gamma_c=0;q1=sin(fea_c/2)*sin(vartheta_c/2)*cos(gamma_c/2)+cos(fea_c/2)*cos(vartheta_c/2)*sin(gamma_c/2);q2=sin(fea_c/2)*cos(vartheta_c/2)*cos(gamma_c/2)+cos(f
三维纯比例制导源代码 %具有碰撞角约束的三维纯比例导引律研究%无人机空中加油自主会合控制器设计%说明:在仿真中,下标"_m"表示追踪器参数,下标"_t"表示目标参数clearclc%---------------------追踪器参数设置-----------------------V_m=300; % 单位:m/s,飞行速度theta_m=45*pi/180;%单位:rad,初始弹道倾角fea_m
神经网络制导律(不带角度约束)源代码 %对原始数据进行规范化处理[pn,meanp,stdp,tn,meant,stdt]=prestd(Tin',Tout);%[pn,meanp,stdp,tn,meant,stdt]=prestd(inputs',output);%乱序N=length(pn);b=randperm(N);B=zeros(2,677);B(1,:)=pn(1,b);B(2,:)=pn(2,b);
具有角度和时间约束的导弹最优全弹道设计(源代码) %具有角度和时间约束的导弹最优全弹道设计%根据论文,初制导功能在于控制导弹快速爬升%指令高度:Hc=1000+35*t,以35m/s快速爬升%指令偏航角为-45degclearclc%---------------------导弹参数设置-----------------------V=260; % 单位:m/s,导弹飞行速度theta=45*pi/180;%单位:rad,初始弹道
带落角约束的变结构末制导律(非线性)源代码 参考文献:周荻那本变结构制导律的设计%===================================================================================%--------------------------------情景说明------------------------------------------%目标正弦机动,At=Asin
带落角约束的变结构末制导律(线性) 源代码 参考文献:周荻那本变结构导引律的设计clear%===================================================================================%--------------------------------情景说明------------------------------------------%目标正弦机动,
比例导引+弹道成型导引源程序 参考文献:《几种增大空地导弹落角的制导方式比较》《Strategic and tactical missile guidance》clearclc%-----------导弹参数---------V_m=260;%飞行速度 X_m=0; Y_m=1000; %初始飞行高度 theta_m=0*pi/180; %弹道倾角 %----------目标参数-------
Zarchan弹道成型导引源程序 参考文献:《几种增大空地导弹落角的制导方式比较》《Strategic and tactical missile guidance》clearclc%-----------导弹参数---------V_m=260;%飞行速度 X_m=0; Y_m=1000; %初始飞行高度 theta_m=0*pi/180; %弹道倾角 %----------目标参数---------
变结构制导律源程序 参考文献:clearclc%-----------导弹参数---------V_m=300;X_m=0;Y_m=3000;theta_m=0*pi/180;%----------目标参数---------V_t=0; %静止目标X_t=3000;Y_t=0;theta_t=0;dtheta_t=0;n_t=0;A_t=0;%变结构制导律参数设置k=1;d
过重比例导引源程序 参考文献:《几种增大空地导弹落角的制导方式比较》clearclc%-----------导弹参数---------V_m=260;%飞行速度 X_m=0; Y_m=1000; %初始飞行高度 theta_m=0*pi/180; %弹道倾角 %----------目标参数---------V_t=0; %静止目标X_t=4800;Y_t=0;theta_t=0;d
利用Isight软件优化matlab程序的一点心得 一、保证matlab源代码是参数优化范围内是可以运行的;二、一定要设置matlab文件搜索目录,否则无法调用自己编写的子程序!!!!!!!此处可参考本人的另一篇博文:在Isight软件里如何调用matlab子程序?
关于落点偏差的疑惑 落点偏差经典方法是通过卯酉半径和角度偏差算纵向和横向偏差,如果落点和目标点近的话可不可以利用发射系三个坐标直接做差、平方、求和,再开方,即L=sqrt((x1-x2)^2+(y1-y2)^2+(z1-z2)^2)这样子算出来的偏差是不是纵向和横向偏差对应的斜距?我用上式算出来偏了50米,但利用卯酉法却差了几公里,很是疑惑。。。发射系是旋转的,上式误差会那么大吗?
在Isight软件里如何调用matlab子程序? 这两天在利用Isight调用matlab时老是报错,执行到while循环时就出现错误,刚开始以为是while循环的问题,后来发现是子程序调用出现问题,不能执行子程序。我琢磨,为什么Isight能识别plot这类matlab自带的指令,为什么就不能识别我自己编的子程序呢?这么一想,就知道解决办法了,得设置搜索路径,让Isight碰到不能识别的指令时到指定的路径下搜索,这样不就解决了吗?!这么一试,还
Isight软件调用matlab matlab m文件里面有while循环(跟优化无关),为什么isight软件在优化时不执行while循环?如果去掉循环可以得到优化结果。循环在这个程序里非常重要,如何处理?
【源代码】信息熵法算权重 function ww=xxs(RD,n,m)%1、无量纲规范矩阵RD%2、计算每个指标的熵值k=1/log(n);for i=1:m for j=1:n f(j,i)=RD(j,i)/sum(RD(:,i)); temp(j,i)=f(j,i)*log(f(j,i)); end H(i)=-k*sum(temp(:,i));
【源代码】《关于敏捷性评价指标权重的研究》 %此子程序来源参考:《关于敏捷性评价指标权重的研究》function [alpha,beta]=xishu(W)[W_ascend,Index]=sort(W,'ascend');nn=length(W_ascend);sum_w=0;for i=1:nn sum_w=sum_w+i*W_ascend(i);endTs=sum_w/nn-(nn-1)/nn;alpha=nn
【源代码】AHP法算权重 %此子程序来源参考:《层次分析法在确定企业家评价指标权重中的应用》function [W,CR]=AHP(P)%计算指标权重[m,n]=size(P);for i=1:m M_t(i)=prod(P(i,:));%计算每一行元素的积 W_t(i)=M_t(i)^(1/n);endW=W_t/sum(W_t);%计算判断矩阵的最大特征值PW=P*W';PW_t=P
【源代码】基于扩展贝叶斯方法融合的作战效能评估研究 %来源《基于扩展贝叶斯方法融合的作战效能评估研究》clearclc%定性指标T2 T7 T8%以下Kcap和Ccap数据对应表6,将表6拆分为两个矩阵%横坐标Kcap:知识度%纵坐标Ccap:满意度%Kcap和Ccap相同的位置组合起来对应信念图上的一个点%Kcap和Ccap:行对应专家(3个专家),列对应指标(3个指标)Kcap=[0.8 0.78 0.85;
【源代码】指标权重赋权 《武器电子系统质量评估》第33页和39页例子%《武器电子系统质量评估》第33页和39页例子clearclcP=[0.5 0.6 0.7; 0.4 0.5 0.6; 0.3 0.4 0.5;];tau1=0.8;tau2=0.5;tau3=0.1;fea1=jisuan(P,tau1)fea2=jisuan(P,tau2)fea3=jisuan(P,tau3)
【源代码】关于D_S证据理论和可信度的计算 参考文献:《基于D_S证据理论的组网雷达_四抗_能力评价》,请自行下载clearclc%第2层权重e1=[0.056 0.157 0.237 0.452 0.099]; % e11 e12 e13 e14 e15e2=[0.483 0.157 0.088 0.272]; % e21 e22 e23 e24 e3=[0.370 0.345 0.185 0.100]; % e31 e3
[源代码]基于D-S证据理论的雷达探测信息融合 (节选自《基于粗糙集与卡尔曼法的防空雷达信息融合技术》)参考文档:http://wenku.baidu.com/view/c52cce66bceb19e8b9f6bade源代码如下:clearclc%D-S证据理论%功能:融合三个雷达的观测数据,判别目标类型%T1、T2、T3为三个观测周期,三行对应三个雷达(m1,m2,m3),六列对应目标类型(A1侦察机、A2直升机、A3
matlab坐标轴axes对象句柄消失的解决办法 今天在GUI编程时,利用axes对象显示图形,在第一次作图完成后,axes句柄消失,其它对象句柄(例如popup对象)仍然存在,后来发现解决的办法:打开axes对象的属性页面,将“NextPlot”属性更改为“ReplaceChildren”即可(默认属性为Replace)。参考资料:http://www.mathworks.com/matlabcentral/answers/103055
Matlab界面中选项卡的切换 Matlab界面中选项卡的切换界面采用matlab/GUIDE设计,若直接采用actxcontrol函数创建TabStrip控件,则不仅要在GUI窗口显示到屏幕之后才能使用actxcontrol函数(可以在OutputFcn函数中使用actxcontrol函数,而不能在OpeningFcn函数中使用),而且,同样不能创建TabStrip控件的回调函数。因此,软件编制过程中用到了选项卡ActiveX控
坐标的3-2-1转动和2-3-1转动 以前我对矩阵的转动次序没有重视,这段时间弹道仿真时发现这两种转动顺序存在很多差别,说一下注意的地方:1、两者获得的欧拉角转换矩阵形式不一样,但代入数值后基本一致;2、如果同时用四元数和欧拉角描述姿态,必须要取同一种转动次序下的对应矩阵;3、3-2-1转动多用于弹道导弹,2-3-1多用于飞机、战术弹、巡航导弹。4、对于两种转动次序,获得的弹道倾角、弹道偏角、倾侧角、俯仰角、偏航角、滚
关于鱼雷高速空投弹道仿真 鱼雷高速空投数学模型见http://download.csdn.net/detail/appe1943/8737337。这篇文档是别人摘抄的,我做了下仿真,难点在于确定gamma_t, gamma_p和FB。文档中说得很轻松,对等式求二阶导就可以求出来,方法倒是没错,但基本不具有操作性。我试着求二阶导展开,非常麻烦。后来我想到一招,不是直接求,而是在解空间寻优,将gamma_t, gamma_
吐槽红米2 2月3日在小米官网上抢了红米2,过了两天到货,发现通话时对方听不到声音。当时还在淘宝上买了一个,那个正常。尝试了更新系统、恢复出厂设置、输入×#检测都不行,百度一查,发现小米很多型号手机都存在这个问题,最早的能搜到2011年的,过了三四年了还是存在这个问题,真是让人失望,连基本的通话功能都有问题。今天去了小米授权服务中心,一看还真热闹,排队的人就跟银行一样多。旁边三星、苹果的客服一个人也没有,华为
如何播放巧虎正版DVD光盘 巧虎DVD带有正版复制保护功能,普通家用dvd播放机一般可以直接播放,但用电脑dvd光驱播放时则会出现画面不清,或者干脆无法播放的现象。经过笔者多次尝试,终于找到了一种电脑播放巧虎dvd的方法,总体思路是先将其制作成镜像iso文件,然后再用dvd播放软件播放。具体如下:1、下载DVDFab 软件,然后利用其复制功能做成ISO文件(普通的镜像制作软件像alcohol 120好像无法制作正版保护的
高超声速飞行器轨迹跟踪控制仿真研究 %仿真来源:,系统仿真学报,2011,谭湘敏,易建强等%由于之前已经做过类似的工作,仿真本文时速度的三阶导数、高度的四阶导数没有采用本文的形式,%而是采用之前推导的模型。控制器采用本文的方法。飞行器初始参数见本文表1.
空间科学学报 网址:http://www.cjss.ac.cn/CN/volumn/home.shtml《空间科学学报》 Chinese Journal of Space Science 《空间科学学报》是中国科学院空间科学与应用研究中心和中国空间科学学会共同主办的综合性学术刊物。它立足于空间科学的前沿,积极反映该领域内的新发现、新创造、新见解、新理论和新方法,展示空间科学发展的进程
最优控制与姿控喷流在导弹姿态控制中的应用 一、仿真来源林枫等,航天控制,2004,《最优控制与姿控喷流在导弹姿态控制中的应用》二、开关线%1、只考虑偏航角fea初始值的不同,固定偏航角速率的初始值,令其一直为零,得到图2曲线结果%2、只考虑偏航角速率初始值的不同,固定偏航角的初始值,令其一直为零,得到曲线结果与图2相比,只是位置左右平移而已%这说明,偏航角、偏航角速率的初始值不影响曲线的形状,只影响曲线的位置
一种免费下载CNKI、维普等资源的方法 科技信息资源免费下载服务论坛http://qbs.hznet.com.cn/bbs/杭州本地用户可以免费下载cnki、维普等资源,我用起来感觉下载速度还可以,很不错,特向大家推荐
应用时标分离和动态逆方法设计飞行器的姿态控制系统(源代码) 1、仿真来源:《应用时标分离和动态逆方法设计飞行器的姿态控制系统》,武立军,贺有智,现代防御技术,20072、在发动机推力的处理上,尝试了三种方法:连续推力、通过幅值限制得到的离散推力和通过PWPF调制得到的离散推力。连续推力很容易实现;后两者与飞行器的参数相关,发动机推力等参数需要仔细设置,由于原文没有给出飞行器的具体参数,发动机推力的大小与原文不同。%滚动通道:推力阈值0.5N,大
高超声速滑翔飞行器摆动式机动突防弹道设计(源代码) %谢愈,刘鲁华等,《高超声速滑翔飞行器摆动式机动突防弹道设计》,航空学报,2011%算法有两个控制量:攻角和倾侧角,攻角只是起辅助作用,主要还是倾侧角的设计%由于倾侧的发生导致维持平衡的升力发生变化,由此攻角需要随之变化%可见,倾侧角是起主动控制的作用,而攻角则是被动响应%论文没有给出攻角的具体表达式(不影响仿真结果),此处仿真也没有体现%论文第四部分,约束条件的处理,在后面的
最优再入机动末制导方法(源代码) 源代码共3个独立的m文件,包括俯仰平面、偏航平面以及三维空间的最优导引。当初编程的时候直接编三维空间,程序老是得不到正确结果,于是化繁为简,先编两个平面程序然后综合。在综合得到三维空间程序的时候,又出现问题,主要是俯冲平面角度定义问题,不能直接使用俯仰平面角度定义%陈克俊,刘鲁华等,远程火箭飞行动力学与制导,P335,最优再入机动末制导方法%国防工业出版社,2014年
攻击时间控制的动态逆三维制导律(源代码) %参考文献:《攻击时间控制的动态逆三维制导律》,哈尔滨工程大学学报,2010%编程难点:视线坐标系与参考坐标系间的转换关系一、采用本文方法结果二、偏航通道采用增广比例导引律
大攻角导弹非线性动态自动驾驶仪研究(源代码) %参考文献:《大攻角导弹非线性动态自动驾驶仪研究》,现代电子技术,2011%原文式(21)中k6和k8的表达式写错了,修正后的表达式见本程序第76和78行%原文第(22)~(24)式错误,实际上(21)式才是指令舵偏角,根据式(6)计算实际舵偏角%原文式(6)中舵机固有角频率为150rad/s,此处改为1000rad/s%由于原文没有给出具体参数,仿真结果与原文有一些不同,但基本能
大气参数计算子程序 function [row,av,g]=OnAir(YH) H=6356.766*YH/(6356766+YH);if H H=0;elseif H Tb=288.15; Lb=-6.5; Pb=101325.3; Hb=0; elseif H Tb=216.65;