Matlab—矩阵论
刘志祥-兰州铁道设计院有限公司
刘志祥,男,1990.10-,2014-2017兰州交通大学硕士,2017-兰州理工大学博士在读,2017-兰州铁道设计院有限公司站场工程师
展开
-
《矩阵论引论》田振际——P30-1.4
clearclc%《矩阵论引论》田振际 王永铎...%P30 第4题format ratsyms aa1=[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 · 1450 阅读 · 0 评论 -
《矩阵论引论》田振际——P40-例2.4+标准正交化
clearclc%《矩阵论引论》田振际 王永铎...%P40 例2.4%% 施密特正交化方法format rata1=[1 1 0 0];a2=[1 0 1 0];a3=[-1 0 0 1];a4=[1 -1 -1 1];a1 = sym( a1,'f');%定义浮点类型b1=a1e1=simplify(b1/norm(b1))%单位化b2=a2-dot(a2,e1)*e...原创 2018-04-17 15:56:35 · 984 阅读 · 0 评论 -
《矩阵论引论》田振际——状态方程的约当规范形
%% 求约当规范形clearclc%% 给定一个线性时不变系统的状态方程为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));enddisp('解:特征向量:')V...原创 2018-03-29 16:51:56 · 1420 阅读 · 0 评论 -
《矩阵论引论》田振际——P78-例3.6-矩阵的奇异值分解2
%% 《奇异值分解MATLAB计算程序》clearclcdisp('矩阵')% 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]%例4A=[1 0 1;0 1 1;0 0 0]A=sym(A,'f');...原创 2018-04-19 11:37:27 · 794 阅读 · 0 评论 -
《矩阵论引论》田振际——P78-例3.6-矩阵的奇异值分解1
1 问题描述奇异值分解涉及大量重复性计算,为了提高效率,将计算部分用Matlab实现,具体见代码:2 程序及算例clearclcdisp('矩阵')% 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 · 1369 阅读 · 2 评论