主函数:

clear;
%三阶矩阵
A=[0,1,-1;-6,-11,6;-6,-11,5];
B=[0;0;1];
C=[1,0,0];
n=size(A,1);
Q1=zeros(n);
Q2=zeros(n);J=jordan(A)%约当阵输出%判断是否能控
for i=0:n-1Q1(:,i+1)=A^i*B;
endn1=rank(Q1);%判断是否能观测for j=0:n-1 Q2(j+1,:)=C*A^j;endn2=rank(Q2);%%转化成能控、能观标准型if n1==n[A1,B1,C1,T]=nengkong(A,B,C,n)endif n2==n[A2,B2,C2,K]=nengguan(A,B,C,n)end

能控标准型转换:

function [D,E,F,T] = nengkong(A,B,C,n)
%能控标准型转换
T=zeros(n,n);
T(:,1)=B;
%求其特征多项式
a=poly(A);for i=2:nT(:,i)=A*T(:,i-1)+a(1,i)*T(:,1);
end
T=fliplr(T);
%得到新的A矩阵
D=inv(T)*A*T;
%得到新的B矩阵
E=inv(T)*B;
%得到新的C矩阵
F=C*T;
end

能观标准型转换:

function [D,E,F,T] = nengkong(A,B,C,n)
%能控标准型转换
T=zeros(n,n);
T(:,1)=B;
%求其特征多项式
a=poly(A);for i=2:nT(:,i)=A*T(:,i-1)+a(1,i)*T(:,1);
end
T=fliplr(T);
%得到新的A矩阵
D=inv(T)*A*T;
%得到新的B矩阵
E=inv(T)*B;
%得到新的C矩阵
F=C*T;
end

MATLAB:三阶矩阵约当阵,能控,能观标准型转化相关推荐

  1. matlab将空间描述转变为约当阵,matlab求约当标准型

    -- 矩阵的转置 jordan -- 约当标准型 simple -- 符号矩阵简化 solve(f) -- 求一个方程的解 经济数学实验班 1101 班 数学软件应用 MATLAB 复习资料 .... ...

  2. matlab矩阵相乘没有负数,关于matlab中矩阵与逆矩阵为何相乘不等于单位阵的问题...

    当我们用matlab计算矩阵和逆矩阵时,结果往往不等于单位阵 如下: >> a=magic(4) a = 16     2     3    13 5    11    10     8 ...

  3. 关于matlab中矩阵与逆矩阵为何相乘不等于单位阵的问题

    当我们用matlab计算矩阵和逆矩阵时,结果往往不等于单位阵 如下: >> a=magic(4) a = 16     2     3    13      5    11    10   ...

  4. 【几种特殊的矩阵:对角矩阵、上下三角矩阵、正态分布随机矩阵、魔方矩阵、希尔伯特矩阵、托普利兹矩阵】

    (这两个最不特殊了,线代中学过不过有点忘) 1.对角矩阵:       对角矩阵(diagonal matrix)是一个主对角线之外的元素皆为0的矩阵,常写为diag(a1,a2,...,an) .对 ...

  5. matlab ctrb(),4.8Matlab问题能控能观解释.ppt

    Ch.4 线性系统的能控性和能观性 Matlab问题(1/1) 4.8 Matlab问题 本章涉及的计算问题主要有 状态能控性/能观性判定. 系统能控能观分解. 能控/能观规范形变换以及 能控/能观规 ...

  6. matlab将矩阵分解成lu,10行代码实现矩阵的LU分解(matlab)

    最近由于数值分析实验课要求,需要通过matlab实现矩阵的LU分解.但是看了很多网友写的程序,基本上都是通过循环嵌套循环来实现矩阵的LU分解.略感琐碎,因此最近两天便一直在思考能否利用矩阵的乘v法,来 ...

  7. matlab矩阵特征分解,用MATLAB实现矩阵分解

    MATLAB求解线性方程的过程基于三种分解法则: (1)Cholesky分解,针对对称正定矩阵: (2)高斯消元法,  针对一般矩阵: (3)正交化,      针对一般矩阵(行数≠列数) 这三种分解 ...

  8. 大数据矩阵计算基础(二):R、Excel、Matlab中矩阵操作

    一.R语言 构造矩阵 a=matrix(c(1,2,3,4),ncol=2,byrow=T) b=matrix(c(5,6,7,8),ncol=2,byrow=T) c=matrix(c(1,2,3, ...

  9. matlab中矩阵的表示与简单操作

    matlab中矩阵的表示与简单操作 原文地址为:matlab矩阵的表示和简单操作 一.矩阵的表示 在MATLAB中创建矩阵有以下规则: a.矩阵元素必须在"[ ]"内: b.矩阵的 ...

最新文章

  1. 欧阳自远:有个性的嫦娥12345,如何不重复美国探月路?
  2. Webform(Linq高级查、分页、组合查询)
  3. 缩进动画Scaleanimation的一个小示例
  4. 详解 nginx 配置文件及构建 web 虚拟主机
  5. 三星Galaxy A80首款保护壳曝光:配件厂商这样解难题
  6. css3转换图形展示,CSS3的常见transformation图形变化用法小结
  7. Centos7.4源码搭建zabbix3.4.11企业级监控
  8. 项目cookie优化之cookie数量限制和多页面共享污染数据问题
  9. APACHE 在windows下的配置
  10. 风控模型面试问题汇总
  11. lammps后处理:Python调用Ovito模块配置方法
  12. VUE输入 身份证 号,获取出生年月日
  13. Ubuntu 18.04 网口创建网络共享
  14. 【算法】递归(recursion)+经典例题个人分析
  15. 《缠中说禅108课》23:市场与人生
  16. Windows Server 2008(修改计算机名并加入工作组)
  17. 山西大同大学计算机科学与技术在哪个校区,山西大同大学有几个校区及校区地址 哪...
  18. 凹凸/法线/移位贴图的区别
  19. 干货|红外热成像摄像头拆解分析
  20. python 嵌套函数及其练习题

热门文章

  1. Canvas+Socket搞出一个多人游玩的“我画你猜”
  2. 电脑tcp协议设置成服务器,电脑tcp协议设置成服务器
  3. 神经网络预测指标是什么,神经网络怎么预测数据
  4. Uos统信系统 本地APT源配置
  5. CAD保存时,提示 警告 多重从属对象 XXX
  6. iOS 获取当前的UIViewController
  7. 【LittleVGL】Windows环境下利用CodeBlocks搭建LittleVGL的PC模拟器环境配置问题
  8. php 加载慢,解决php加载慢的一个办法_PHP教程
  9. 001 线性代数之行列式:定义、逆序数、余子式与代数余子式、n个易算行列式、范德蒙行列式
  10. excel表格末尾添加一行_教你使用Python批量读写excel文件