![](https://img-blog.csdnimg.cn/20201014180756918.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
Matlab—矩阵论
刘志祥-兰州铁道设计院有限公司
刘志祥,男,1990.10-,2014-2017兰州交通大学硕士,2017-兰州理工大学博士在读,2017-兰州铁道设计院有限公司站场工程师
展开
-
《矩阵论引论》田振际——P30-1.4
clear clc %《矩阵论引论》田振际 王永铎... %P30 第4题 format rat syms a a1=[1 0 0 0]'; a2=[0 1 0 0]'; a3=[0 0 1 0]'; a4=[0 0 0 1]'; disp('解:基I:') I=[a1 a2 a3 a4] b1=[2 1 -1 1]'; b2=[0 3 1 0]'; b3=[5 3 2 1]'; b4=[6 6 ...原创 2018-04-17 11:50:55 · 1382 阅读 · 0 评论 -
《矩阵论引论》田振际——P40-例2.4+标准正交化
clear clc %《矩阵论引论》田振际 王永铎... %P40 例2.4 %% 施密特正交化方法 format rat a1=[1 1 0 0]; a2=[1 0 1 0]; a3=[-1 0 0 1]; a4=[1 -1 -1 1]; a1 = sym( a1,'f');%定义浮点类型 b1=a1 e1=simplify(b1/norm(b1))%单位化 b2=a2-dot(a2,e1)*e...原创 2018-04-17 15:56:35 · 948 阅读 · 0 评论 -
《矩阵论引论》田振际——状态方程的约当规范形
%% 求约当规范形 clear clc %% 给定一个线性时不变系统的状态方程为 disp('变换前系数矩阵:') A=[2 -1 -1;0 -1 0;0 2 1] B=[7 2 3]' %% (1)定出系统的特征值和特征向量。 [V D]=eig(A); for i=1:length(V) V(:,i)=V(:,i)/max(V(:,i)); end disp('解:特征向量:') V...原创 2018-03-29 16:51:56 · 1370 阅读 · 0 评论 -
《矩阵论引论》田振际——P78-例3.6-矩阵的奇异值分解2
%% 《奇异值分解MATLAB计算程序》 clear clc disp('矩阵') % format rat %为了便于观察,可以采用分数形式,‘0’的分数形式是‘*’ % 例1 % A=[1 0 0;2 0 0]; % 例2 % A=[1 0;0 1;1 1] %例3 % A=[1 0 1;0 1 1;1 1 2] %例4 A=[1 0 1;0 1 1;0 0 0] A=sym(A,'f');...原创 2018-04-19 11:37:27 · 760 阅读 · 0 评论 -
《矩阵论引论》田振际——P78-例3.6-矩阵的奇异值分解1
1 问题描述奇异值分解涉及大量重复性计算,为了提高效率,将计算部分用Matlab实现,具体见代码:2 程序及算例clear clc disp('矩阵') % format rat %为了便于观察,可以采用分数形式,‘0’的分数形式是‘*’ % 例1 % A=[1 0 0;2 0 0]; % 例2 % A=[1 0;0 1;1 1] %例3 % A=[1 0 1;0 1 1;1 1 2] %例4 ...原创 2017-04-11 00:41:15 · 1332 阅读 · 2 评论