Matlab
文章平均质量分 74
itheta
这个作者很懒,什么都没留下…
展开
-
关于ode和不同写法
% clc, clear, close all%%%%%% 离散写法clc, clear, close allx0=[-10;0.5];t_final=20;dt = 0.01;tmax = t_final/dt;KK = -1*[395 396 -18 -18 -18];x1(1)=-10;x2(1)=0.5; for k=1:(tmax) %xx(k)=[x1(k) x2(k) x1(k)*x1(k) x1(k)*x2(k) x2(k)*x2(k)...原创 2021-04-02 15:36:30 · 302 阅读 · 0 评论 -
状态反馈及Matlab实现
状态反馈及Matlab实现可控性验证使用Matlab计算可控性cont = ctrb(A,B)ANS = rank(cont);实例设系统的状态方程为x˙=[1101−10013]x+[100]u\dot{x}=\left[\begin{array}{ccc}1 & 1 & 0 \\1 & -1 & 0 \\0 & 1 & 3\end{array}\right] x+\left[\begin{array}{l}1 \\0 \\0原创 2021-01-24 10:32:17 · 8960 阅读 · 0 评论 -
YALMIP安装、教程与示例
Yalmip使用1. 创建决策变量: >> x = sdpvar(m, n [, option]):创建m*n的连续型决策变量矩阵,option是对矩阵的一些参数指定。 相应的,如果要创建整型或二值型决策变量,matlab语句分别为: >> x = intvar(m, n, [option]) >> x = binvar(m, n, [option])2. 添加约束: >> F = set(constraint [,原创 2020-11-29 23:49:50 · 40212 阅读 · 5 评论 -
Matlab画图
dout画图t=dout.time;y1=dout.data(:,1);y2=dout.data(:,2);%figure('name','Figd','Position',[300 200 380 280]) figure('name','Figd','Position',[800 200 380 280])subplot(2,1,1,'align');plot(t,y2,'black','linewidth',1);ylabel('(Nm)','Fontsize',15,'Fontname原创 2020-06-17 11:22:21 · 281 阅读 · 0 评论 -
matlab疑难问题笔记
mat文件的存储和读取先介绍下Matlab存储和读取数据的一种方法Matlab保存数据● save:将工作区中的所有变量保存在当前工作区中的文件中,文件名为 matlab.mat,MAT文件可以通过load函数再次导入工作区,MAT函数可以被不同的机器导入,甚至可以通过其他的程序调用。● save('filename'):将工作区中的所有变量保存为文件,文件名由filename指定。如果filename中包含路径,则将文件保存在相应目录下,否则默认路径为当前路径。...原创 2020-06-11 20:58:09 · 1754 阅读 · 0 评论 -
Duplication matrix and elimination matrix
Wiki上的解释Duplication matrix DnDnDn 是一个线性变换,用于将半线性化矩阵(half-vectorization)转换为线性化矩阵(vectorization).对于一个对称阵AAA而言,我们将其线性化矩阵表示为vech(A)vech(A)vech(A),线性化矩阵表示为vec(A)vec(A)vec(A).。例如,一个2×22 \times 22×2的对称矩阵 ...原创 2020-04-28 12:12:53 · 889 阅读 · 0 评论 -
Matlab求解微分方程
引言很多时候是在simulink搭建模块,这里写写在m文件中求微分方程的做法。采用Matlab画相轨迹y¨+μ(y2−1)y˙+y=0\ddot{y}+\mu\left(y^{2}-1\right) \dot{y}+y=0y¨+μ(y2−1)y˙+y=0状态方程为x˙1=x2,x˙2=−μ(x12−1)x2−x1\dot{x}_{1}=x_{2}, \dot{x}_{2}=-\mu...原创 2020-01-05 16:55:38 · 1411 阅读 · 0 评论