现代控制理论(机器人方向)考核要求与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)简明教程的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!


原文地址:https://blog.csdn.net/ZhangRelay/article/details/88654172
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.chinasem.cn/article/715454

相关文章

Nexus安装和启动的实现教程

《Nexus安装和启动的实现教程》:本文主要介绍Nexus安装和启动的实现教程,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、Nexus下载二、Nexus安装和启动三、关闭Nexus总结一、Nexus下载官方下载链接:DownloadWindows系统根

CnPlugin是PL/SQL Developer工具插件使用教程

《CnPlugin是PL/SQLDeveloper工具插件使用教程》:本文主要介绍CnPlugin是PL/SQLDeveloper工具插件使用教程,具有很好的参考价值,希望对大家有所帮助,如有错... 目录PL/SQL Developer工具插件使用安装拷贝文件配置总结PL/SQL Developer工具插

Python远程控制MySQL的完整指南

《Python远程控制MySQL的完整指南》MySQL是最流行的关系型数据库之一,Python通过多种方式可以与MySQL进行交互,下面小编就为大家详细介绍一下Python操作MySQL的常用方法和最... 目录1. 准备工作2. 连接mysql数据库使用mysql-connector使用PyMySQL3.

Java中的登录技术保姆级详细教程

《Java中的登录技术保姆级详细教程》:本文主要介绍Java中登录技术保姆级详细教程的相关资料,在Java中我们可以使用各种技术和框架来实现这些功能,文中通过代码介绍的非常详细,需要的朋友可以参考... 目录1.登录思路2.登录标记1.会话技术2.会话跟踪1.Cookie技术2.Session技术3.令牌技

如何搭建并配置HTTPD文件服务及访问权限控制

《如何搭建并配置HTTPD文件服务及访问权限控制》:本文主要介绍如何搭建并配置HTTPD文件服务及访问权限控制的问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、安装HTTPD服务二、HTTPD服务目录结构三、配置修改四、服务启动五、基于用户访问权限控制六、

Python使用Code2flow将代码转化为流程图的操作教程

《Python使用Code2flow将代码转化为流程图的操作教程》Code2flow是一款开源工具,能够将代码自动转换为流程图,该工具对于代码审查、调试和理解大型代码库非常有用,在这篇博客中,我们将深... 目录引言1nVflRA、为什么选择 Code2flow?2、安装 Code2flow3、基本功能演示

Java Spring 中的监听器Listener详解与实战教程

《JavaSpring中的监听器Listener详解与实战教程》Spring提供了多种监听器机制,可以用于监听应用生命周期、会话生命周期和请求处理过程中的事件,:本文主要介绍JavaSprin... 目录一、监听器的作用1.1 应用生命周期管理1.2 会话管理1.3 请求处理监控二、创建监听器2.1 Ser

MySQL 安装配置超完整教程

《MySQL安装配置超完整教程》MySQL是一款广泛使用的开源关系型数据库管理系统(RDBMS),由瑞典MySQLAB公司开发,目前属于Oracle公司旗下产品,:本文主要介绍MySQL安装配置... 目录一、mysql 简介二、下载 MySQL三、安装 MySQL四、配置环境变量五、配置 MySQL5.1

MQTT SpringBoot整合实战教程

《MQTTSpringBoot整合实战教程》:本文主要介绍MQTTSpringBoot整合实战教程,本文通过实例代码给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考... 目录MQTT-SpringBoot创建简单 SpringBoot 项目导入必须依赖增加MQTT相关配置编写

在Java中基于Geotools对PostGIS数据库的空间查询实践教程

《在Java中基于Geotools对PostGIS数据库的空间查询实践教程》本文将深入探讨这一实践,从连接配置到复杂空间查询操作,包括点查询、区域范围查询以及空间关系判断等,全方位展示如何在Java环... 目录前言一、相关技术背景介绍1、评价对象AOI2、数据处理流程二、对AOI空间范围查询实践1、空间查