现代控制理论（机器人方向）考核要求与Matlab（Octave）简明教程

2024-02-16 19:20
文章标签 matlab 教程 方向 考核 机器人 理论 控制 现代 要求 简明 octave

本文主要是介绍现代控制理论（机器人方向）考核要求与Matlab（Octave）简明教程，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

课程全部资料请查阅课程分类：https://blog.csdn.net/zhangrelay/article/category/6161998

现代控制理论成绩构成为如下四个部分：

总成绩根据平时成绩（包括考勤、作业、课堂测试等占30%）、编程考核30%、创新实践报告10%、期末考试（占30%）综合评定。期末考试形式采用闭卷笔试。

创新实践报告模版：https://share.weiyun.com/5tHl9I6

Octave Online（Matlab）：https://octave-online.net/

编程示例：https://blog.csdn.net/ZhangRelay/article/details/51615389

成绩构成说明

Matlab简明教程（对应教材）：

可以使用电脑系统Windows/MacOS/Linux，也可以使用手机系统Android/ios等实现。

>> connector on
首次运行 MATLAB Connector 时，必须指定密码。
请输入在设备上设置 MATLAB Mobile 时所用的同一密码。
Password: *************
DNS 名称: 
IP 地址: 192.168.x.xxx
使用此链接可测试 MATLAB Connector:
http://192.168.x.xxx:31415/
如果测试成功，但 MATLAB Mobile 无法连接，
可能是因为您的计算机有多个 IP 地址。要确定
正确的地址，请参阅确定计算机的 DNS 名称或 IP 地址。

程序代码
执行结果
Python I
Python II

倒立摆案例：

A=[ 0 1 0 0; 0 0 0 0;0 0 0 1; 0 0 29.4 0]
B=[0 ; 1 ; 0 ;3 ]
C=[1 0 0 0; 0 0 1 0]
D=[0 ; 0]syms t
eAt=expm(A*t)ctrb(A,B)rank(ctrb(A,B))
rank(obsv(A,C))step(A, B ,C ,D)
flag=0;
[z,p,k]=ss2zp(A,B,C,D,1);
disp('系统零点，极点和增益为：');
z
p
k
n=length(A);
for ii=1:nif real(p(ii))>0flag=1;end
end
if flag==1disp('系统不稳定');
elsedisp('系统稳定');
end% Q=eye(4,4);
% P=lyap(A,Q);
% flag=0;
% n=length(A);
% for i=1:n
%     det(P(1:i,1:i))
%     if(det(P(1:i,1:i))<=0)
%         flag=1;
%     end
% end
% if flag==1
%     disp('系统不稳定');
% else
%     disp('系统稳定');
% endP=[-10 -10 -2-2*sqrt(3)*i -2+2*sqrt(3)*i]
K=acker(A,B,P)A =0    1.0000         0         00         0         0         00         0         0    1.00000         0   29.4000         0B =0103C =1     0     0     00     0     1     0D =00eAt =[ 1, t,                                                                                                          0,                                                                                                      0]
[ 0, 1,                                                                                                          0,                                                                                                      0]
[ 0, 0,                                           exp(-(7*3^(1/2)*5^(1/2)*t)/5)/2 + exp((7*3^(1/2)*5^(1/2)*t)/5)/2, (3^(1/2)*5^(1/2)*exp((7*3^(1/2)*5^(1/2)*t)/5))/42 - (3^(1/2)*5^(1/2)*exp(-(7*3^(1/2)*5^(1/2)*t)/5))/42]
[ 0, 0, (7*3^(1/2)*5^(1/2)*exp((7*3^(1/2)*5^(1/2)*t)/5))/10 - (7*3^(1/2)*5^(1/2)*exp(-(7*3^(1/2)*5^(1/2)*t)/5))/10,                                       exp(-(7*3^(1/2)*5^(1/2)*t)/5)/2 + exp((7*3^(1/2)*5^(1/2)*t)/5)/2]ans =0    1.0000         0         01.0000         0         0         00    3.0000         0   88.20003.0000         0   88.2000         0ans =4ans =4系统零点，极点和增益为：z =5.4222    0.0000-5.4222   -0.0000p =5.4222-5.422200k =1.00003.0000系统不稳定P =-10.0000 + 0.0000i -10.0000 + 0.0000i  -2.0000 - 3.4641i  -2.0000 + 3.4641iK =-54.4218  -24.4898   93.2739   16.1633

课后习题参考，编程示例不再重复列出。

第一章：状态空间表达式

num为传递函数分子参数，den为传递函数分母参数，tf为传递函数，ss为状态空间，

tf2ss传递函数转状态空间，ss2tf状态空间转传递函数。

掌握系统框图、模拟结构图、状态方程组、状态空间表达式（不唯一）、传递函数等。

% 1.6
num=[6];
den=[1 6 41 7];
[A B C D]=tf2ss(num,den)%1.7
num=[360 440];
den=[1 28 196 740];
[A B C D]=tf2ss(num,den)%1.8%1.9
A=[0 1 -1; -6 -11 6; -6 -11 5];
[P J]=eig(A)
inv(P)*[0;0;1]
[1 0 0]*P%1.10
[T J]=jordan(A)
inv(T)*[0;0;1]
[1 0 0]*T%1.11
A=[0 1 0; 0 0 1; 2 3 0];
[T J]=jordan(A)
%[P J]=eig(A)

第二章：表达式的解

step求解阶跃，plot画图。

%2.1 2.2 2.4 2.6
syms t
A=[0 1; -2 -3]
eAt=expm(A*t)%2.3 2.7
syms t
A=[0 1 0; 0 0 1; 2 -5 4]
eAt=expm(A*t)%2.8
syms t
A=[0 1; -2 -3]
B=[0;1]
x0=[0;0]
eAt=expm(A*t)
xt=eAt*x0+inv(A)*(eAt-1)*B*1

修订：

%2.8
syms t
A=[0 1; -2 -3]
B=[0;1]
x0=[0;0]
eAt=expm(A*t)
xt=eAt*x0+int(eAt*B*1,t)
%xt=eAt*x0+inv(A)*(eAt-1)*B*1%xt2.6 
syms t
A=[0 1; 0 0]
B=[0;1]
C=[1 0]
x0=[1;1]
eAt=expm(A*t)
xt=eAt*x0+int(eAt*B*1,t)
%xt=eAt*x0+inv(A)*(eAt-1)*B*1
yt=C*xt

第三章：能控性和能观性

%3.2
A=[-4 5;1 0]
B=[-5;1]
[T,J]=jordan(A)
inv(T)*B%3.5
M=[B, A*B]
rank(M)%3.8
A=[1 2 1; 0 1 0; 1 0 3]
B=[1 0; 0 1; 0 0]
M=[B A*B A*A*B]
rank(M)

第四章：稳定性和李雅普诺夫法

%4-1
A=[-1 0; 0 1];
B=[1;1];
C=[1 0];
D=[0];
flag=0;
[z,p,k]=ss2zp(A,B,C,D,1);
disp('系统零点，极点和增益为：');
z
p
k
n=length(A);
for i=1:nif real(p(i))>0flag=1;end
end
if flag==1disp('系统不稳定');
elsedisp('系统稳定');
end%4-5
A=[0 1; -1 -1];
Q=eye(2,2);
P=lyap(A,Q);
flag=0;
n=length(A);
for i=1:ndet(P(1:i,1:i))if(det(P(1:i,1:i))<=0)flag=1;end
end
if flag==1disp('系统不稳定');
elsedisp('系统稳定');
end%4-6
A=[0 1; -1 0];
Q=eye(2,2);
P=lyap(A,Q);
flag=0;
n=length(A);
for i=1:ndet(P(1:i,1:i))if(det(P(1:i,1:i))<=0)flag=1;end
end
if flag==1disp('系统不稳定');
elsedisp('系统稳定');
end%4-8
A=[1 1; -1 1];
Q=eye(2,2);
P=lyap(A,Q);
flag=0;
n=length(A);
for i=1:ndet(P(1:i,1:i))if(det(P(1:i,1:i))<=0)flag=1;end
end
if flag==1disp('系统不稳定');
elsedisp('系统稳定');
end%4-9
A=[0 1; -2 -3];
Q=eye(2,2);
P=lyap(A,Q);
flag=0;
n=length(A);
for i=1:ndet(P(1:i,1:i))if(det(P(1:i,1:i))<=0)flag=1;end
end
if flag==1disp('系统不稳定');
elsedisp('系统稳定');
end

第五章：线性定常系统综合

%5-2
A=[0 1 0;0 0 1;0 -2 -3];
B=[0;0;1];
P=[-2 -1+1i -1-1i];
M=[B,A*B,A*A*B];
n=length(A);
rank(M)
if rank(M)==ndisp('系统可控')disp('状态反馈')K=acker(A,B,P)
elsedisp('系统不可控')[Ac,Bc,Cc,T,K]=ctrbf(A,B,C)
end    
Ac=A-B*K
disp('配置后极点')
eig(Ac)

第六章：最优控制

-未完待续-

 

这篇关于现代控制理论（机器人方向）考核要求与Matlab（Octave）简明教程的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！


http://www.chinasem.cn/article/715454。

相关文章

Spring Security 基于表达式的权限控制

Spring Security 基于表达式的权限控制

前言 spring security 3.0已经可以使用spring el表达式来控制授权，允许在表达式中使用复杂的布尔逻辑来控制访问的权限。 常见的表达式 Spring Security可用表达式对象的基类是SecurityExpressionRoot。 表达式描述hasRole([role])用户拥有制定的角色时返回true （Spring security默认会带有ROLE_前缀）,去
阅读更多...
Spring Security 从入门到进阶系列教程

Spring Security 从入门到进阶系列教程

Spring Security 入门系列 《保护 Web 应用的安全》 《Spring-Security-入门（一）：登录与退出》 《Spring-Security-入门（二）：基于数据库验证》 《Spring-Security-入门（三）：密码加密》 《Spring-Security-入门（四）：自定义-Filter》 《Spring-Security-入门（五）：在 Sprin
阅读更多...
Makefile简明使用教程

Makefile简明使用教程

文章目录 规则makefile文件的基本语法：加在命令前的特殊符号：.PHONY伪目标: Makefilev1 直观写法v2 加上中间过程v3 伪目标v4 变量 make 选项-f-n-C Make 是一种流行的构建工具，常用于将源代码转换成可执行文件或者其他形式的输出文件(如库文件、文档等)。Make 可以自动化地执行编译、链接等一系列操作。 规则 makefile文件
阅读更多...
2024年流动式起重机司机证模拟考试题库及流动式起重机司机理论考试试题

2024年流动式起重机司机证模拟考试题库及流动式起重机司机理论考试试题

题库来源：安全生产模拟考试一点通公众号小程序 2024年流动式起重机司机证模拟考试题库及流动式起重机司机理论考试试题是由安全生产模拟考试一点通提供，流动式起重机司机证模拟考试题库是根据流动式起重机司机最新版教材，流动式起重机司机大纲整理而成（含2024年流动式起重机司机证模拟考试题库及流动式起重机司机理论考试试题参考答案和部分工种参考解析），掌握本资料和学校方法，考试容易。流动式起重机司机考试技
阅读更多...
系统架构师考试学习笔记第三篇——架构设计高级知识（20）通信系统架构设计理论与实践

系统架构师考试学习笔记第三篇——架构设计高级知识（20）通信系统架构设计理论与实践

本章知识考点：         第20课时主要学习通信系统架构设计的理论和工作中的实践。根据新版考试大纲,本课时知识点会涉及案例分析题(25分),而在历年考试中,案例题对该部分内容的考查并不多,虽在综合知识选择题目中经常考查,但分值也不高。本课时内容侧重于对知识点的记忆和理解,按照以往的出题规律,通信系统架构设计基础知识点多来源于教材内的基础网络设备、网络架构和教材外最新时事热点技术。本课时知识
阅读更多...
SWAP作物生长模型安装教程、数据制备、敏感性分析、气候变化影响、R模型敏感性分析与贝叶斯优化、Fortran源代码分析、气候数据降尺度与变化影响分析

SWAP作物生长模型安装教程、数据制备、敏感性分析、气候变化影响、R模型敏感性分析与贝叶斯优化、Fortran源代码分析、气候数据降尺度与变化影响分析

查看原文>>>全流程SWAP农业模型数据制备、敏感性分析及气候变化影响实践技术应用 SWAP模型是由荷兰瓦赫宁根大学开发的先进农作物模型，它综合考虑了土壤-水分-大气以及植被间的相互作用；是一种描述作物生长过程的一种机理性作物生长模型。它不但运用Richard方程，使其能够精确的模拟土壤中水分的运动，而且耦合了WOFOST作物模型使作物的生长描述更为科学。 本文让更多的科研人员和农业工作者
阅读更多...
深入理解RxJava：响应式编程的现代方式

深入理解RxJava：响应式编程的现代方式

在当今的软件开发世界中，异步编程和事件驱动的架构变得越来越重要。RxJava，作为响应式编程（Reactive Programming）的一个流行库，为Java和Android开发者提供了一种强大的方式来处理异步任务和事件流。本文将深入探讨RxJava的核心概念、优势以及如何在实际项目中应用它。 文章目录 💯 什么是RxJava？💯 响应式编程的优势💯 RxJava的核心概念
阅读更多...
沁恒CH32在MounRiver Studio上环境配置以及使用详细教程

沁恒CH32在MounRiver Studio上环境配置以及使用详细教程

目录 1.  RISC-V简介 2.  CPU架构现状 3.  MounRiver Studio软件下载 4.  MounRiver Studio软件安装 5.  MounRiver Studio软件介绍 6.  创建工程 7.  编译代码 1.  RISC-V简介         RISC就是精简指令集计算机（Reduced Instruction SetCom
阅读更多...
嵌入式方向的毕业生，找工作很迷茫

嵌入式方向的毕业生，找工作很迷茫

一个应届硕士生的问题： 虽然我明白想成为技术大牛需要日积月累的磨练，但我总感觉自己学习方法或者哪些方面有问题，时间一天天过去，自己也每天不停学习，但总感觉自己没有想象中那样进步，总感觉找不到一个很清晰的学习规划……眼看 9 月份就要参加秋招了，我想毕业了去大城市磨练几年，涨涨见识，拓开眼界多学点东西。但是感觉自己的实力还是很不够，内心慌得不行，总怕浪费了这人生唯一的校招机会，当然我也明白，毕业
阅读更多...
前端技术(七)——less 教程

前端技术(七)——less 教程

一、less简介 1. less是什么？ less是一种动态样式语言，属于css预处理器的范畴，它扩展了CSS语言，增加了变量、Mixin、函数等特性，使CSS 更易维护和扩展LESS 既可以在 客户端 上运行 ，也可以借助Node.js在服务端运行。 less的中文官网:https://lesscss.cn/ 2. less编译工具 koala 官网 http://koala-app.
阅读更多...

Copyright China编程 23002807@qq.com

沪ICP备2023033072号-2