2017年06月_马小_乐

原创自动控制小结——稳定裕度与校正网络

第四章：根轨迹法1.定义闭环系统特征方程，获得s=f(K),K值变化即得到s（根）的轨迹。2.零点：分子根极点：分母根3.根轨迹画图本质是研究闭环方程得到的法则，核心在P105-106自动控制原理书上。第五章频域分析法1.传递函数频率特性：让式子变成 a+bj.然后幅频特性即为 sqrt(a2+b2),相频特性为 act(b/a)。

2017-06-28 15:06:22 1651

原创 M插值之我见

问：1.Matlab读数据（C读数据） 2.超出边界如何处理自：按旁边斜率处理众：？见： 3.减少重复程序interp,单独拉出来 4.减少输入量，矩阵长度自己计算测试案例：clc;clear;%%%一维插值测试A1=[0;1;2;3;4;5];A2=[1.2;2.4;3.6;4.8;

2017-06-17 10:35:53 313

原创关于控制之理解我见

一：1.开环传递函数：主反馈量/输入量（拉氏变换） 2.闭环传递函数：1）输出/输入（拉氏变换，不考虑扰动） 2）冻结输入（=0），获得输出=f1（）*扰动 ,与输出=f2()*输入相加，即为扰动作用下的闭传 3.误差：E(S)=R-反馈（输入信号与反馈信号差值）稳态误差ess（t）:E(S)拉反后的稳态分量（另一部分是瞬态分量）。

2017-06-15 11:02:55 414

原创龙格库塔的计算写法

参考《导弹飞行力学》对部分参数的解释:dx/dt=f(t,x): 之前一直看不懂f（t,x）到底指的哪个式子，其实在开头就提出来了，f是dy K2=△t*f(tk+△t/2,xk+1/2*K1)：t在导弹飞行力学，是y(0)(M中是1)，所以可忽略。然后xk+1/2*K1即dx/dt的变量进行变换，在飞行力学没有直接加，需要通过y来改变dy而y的变换则是借助y+h*dy去改变

2017-06-12 17:12:57 1727

原创关于导引法

1）数值积分项为 R,q,xm,ym,x,y,ny,Eta(比例导引法涉及，并需要给出初值)2）关于matlab的ode45使用方法：[t,y,Tend,Yend,evennum] = ode45(@dif_missle_trace1,tspan,y0,op); 变步长，无法设置步长项输入：tspan:[t0,tf]初试时间及结束时间。y0:初始值dif_

2017-06-12 15:38:06 1173

原创关于多学科框架搭建

1.兰文博：系统级变量与目标函数，学科级变量与目标函数与约束。解释：道理我都懂，大系统是否需要加上代理模型从上图可以看出，大系统也加响应面模型。上述为学科框架。（过载传递是不太懂的）

2017-06-06 11:22:43 378

原创关于应用Isight做代理模型及EI采点+优化

1.代理模型应用MATLAB的DACE工具箱制作制作流程代码如下：load('shuju1.txt');X(:,1)=shuju1(:,1);X(:,2)=shuju1(:,2);YY(:,1)=shuju1(:,3);theta=[10 10];lob=[1e-1 1e-1];upb=[20 20];[dmodel,perf]=dacefit(X,YY,@re

2017-06-05 16:07:34 5502 9

转载 MATLAB中自带遗传算法函数GA的用法

ga用遗传算法寻找函数的最优解语法规则x = ga(fitnessfcn,nvars)x = ga(fitnessfcn,nvars,A,b)x = ga(fitnessfcn,nvars,A,b,Aeq,beq)x = ga(fitnessfcn,nvars,A,b,Aeq,beq,LB,UB)%其中fitnessfc为函数的句柄或者为匿名函数

2017-06-05 10:47:26 83842 2