不需要用循环,直接使用逻辑数组运算更快。
我假设你的x和x(1)--x(n)是不同的变量,那么我用a来表示你的x,用x数组代表你的x(1)--x(n)
我随便给这些数赋值
a=10;
x=1:3:30;
k=5;
temp = true(size(x));
temp(k) = false;
y = prod(a - x) / prod(x(k) - x(temp));
a-x是生成数组[a-x(1), a-x(2), , a-x(n)],prod()是吧这个数组相乘
temp是一个逻辑数组,和x相同大小,并且只有第k个数是false。x(temp)就是生成[x(1), x(2),, x(k-1), x(k+1), , x(n)]
其他的应该容易理解
补充:拉格朗日插值多项式可也这么写:
%x0,y0是已知各点坐标,共n点
x0 = [1, 2, 3];
y0 = [1, 8, 27];
n = length(x0);
%x是要求的点
x = 4;
y = 0;
for k = 1:n
temp = true(size(x0));
temp(k) = false;
y = y + y0(k) prod(x - x0(temp)) / prod(x0(k) - x0(temp));
end
倒立摆是自动控制理论研究的典型实验设备,由于其具有强非线性与 强耦合性,一直是先进控制算法研究的热点文章采用了一种微分几何方法--李理论,对倒立摆系统进行精确反馈线性化,此种线性化方法使模型更多地保留了原 系统主要的非线性部分,更能逼近实际系统在此线性化模型的基础上,对倒立摆系统进行最优稳定控制设计仿真结果表明,所提出的精确反馈线性化方法对于倒 立摆系统的控制器设计是有效的
以上就是关于问一下大神matlab应该怎么编程全部的内容,包括:问一下大神matlab应该怎么编程、基于matlab/simulink的二级倒立摆系统最优控制有效性通过什么来验证、等相关内容解答,如果想了解更多相关内容,可以关注我们,你们的支持是我们更新的动力!
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)