- 博客(4)
- 收藏
- 关注
原创 相似上三角矩阵变换
function [ B, P ] = my_same_tr( A ) %求与A相似的上三角矩阵和变换矩阵P %B = P2^(-1)*P1^(-1)*A*P1*P2 A = sym(A); n = size(A,2); for i = 1:n-1 [V1,~] = eig(A(i:n,i:n)); %V1 = V(:,1); V2 = null(V1'); V = [V1,V2]; B(
2020-07-02 19:42:58 3386
原创 用Givens旋转进行矩阵QR分解
不多废话,直接贴代码function [A_, T_, T] = my_qr_givens( A )%利用givens旋转进行qr分解%输出%A_ 每次变换后的A矩阵%T_ 对应于A_的变换矩阵 A_ = sym([]); T_ = sym([]); A = sym(A); n = size(A,1); T = sym(eye(n)); sum = 1; for j = 1:n-1 %从第一列到第n-1列,全部变换到第j列第j个元素上
2020-07-02 19:41:03 4720
原创 利用Householder变换进行矩阵QR分解
利用Householder变换进行非奇异矩阵的QR分解非奇异矩阵QR分解Householder变换贴代码非奇异矩阵QR分解实(复)非奇异矩阵A\boldsymbol{A}A,化成正交(酉)矩阵的Q\boldsymbol{Q}Q与实(复)非奇异上三角矩阵的R\boldsymbol{R}R的乘积的形式,即A=QR\boldsymbol{A=QR}A=QRHouseholder变换又称作对称变换,将需要变换的向量以某一镜像向量为参照进行对称变换,两向量关于镜像向量对称。假设Householder变换
2020-07-02 19:36:07 9066 2
原创 利用MATLAB求解一阶线性常系数非齐次微分方程组
用矩阵函数求解一阶线性常系数齐次微分方程组主要步骤1.问题形式2.求矩阵函数3.代入矩阵A的指数函数得最终解主要步骤本来想用在矩阵论期中开卷考试验证计算结果的,结果一个解方程组的题也没考…在一些学习网站白嫖惯了,也该分享一下知识,希望这里的code能对搜索到的人有帮助。主要应用了matlab的符号矩阵,讲解不侧重推导,只侧重编程实现,详情可以看矩阵论相关的参考书。1.问题形式一阶线性常系...
2020-04-16 23:25:07 9370 2
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人