![](https://img-blog.csdnimg.cn/20201014180756927.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
Matlab
雪shan飞hu
这个作者很懒,什么都没留下…
展开
-
矩阵范数
1、常用的矩阵范数主要是p=1,2,∞2、谱半径3、谱半径与矩阵范数的关系原创 2020-03-28 18:19:03 · 16561 阅读 · 0 评论 -
数值线性代数之cholesky平方根法(MATLAB)
1. 该方法只适用于对称正定矩阵所以我们可以解得L代入原方程组即可2. MATLAB代码function [L]=mycholesky (A)n=length(A);for k=1:n A(k,k)=sqrt(A(k,k)); A(k+1:n,k)=A(k+1:n,k)/A(k,k); for j=k+1:n ...原创 2020-03-27 21:32:06 · 3191 阅读 · 0 评论 -
数值线性代数之列主元Gauss消去法(MATLAB)
1. 这是啥?在使用Gauss变换的过程中,如果矩阵对角线的元素太小,会干扰计算结果,所以我们尽量把每列较大的元素通过行变换挪到对角线上。2. 第二个红框只是普通的Guass变换,第一个红框就是选择最大列主元的算法。3. 下图中,左边的L代表的是常规的Gauss变换,P代表行变换4. 于是我们要解的方程就从Ax=b变成了LUx=Pb,接下来使用上下三角矩阵的前代法和回代法解...原创 2020-03-27 21:11:52 · 5114 阅读 · 0 评论 -
数值线性代数之Gauss分解(MATLAB)
1. 什么是Gauss变换对于一个向量x,左乘一个矩阵Lk,使得x变成向量y,该y向量第k个元素起皆为0。 比如x=(2,2,6)', 将其转化为向量y=(2,0,0)'。其中的Lk为向量ek=(0,0,...,0,1,0,...,0),2. 你可能没怎么看懂或者懒得看,我们来举个栗子所以Gauss变换至少要保证求分母不为0,不是所有矩阵都可以用的。3. 上个例子是针...原创 2020-03-27 17:57:23 · 5069 阅读 · 1 评论 -
数值线性代数之前代法解下三角方程组和回代法解上三角(MATLAB)
下三角方程组形如解得即MATLAB代码:function b=myFM(L,b)n=length(b);for j=1:n-1 b(j)=b(j)/L(j,j); b(j+1:n)=b(j+1:n)-b(j).*L(j+1:n,j);endb(n)=b(n)/L(n,n);同理,对上三角矩阵MATLAB代码:...原创 2020-03-27 14:11:59 · 6384 阅读 · 1 评论 -
MATLAB实现矩阵given变换
一共三个function文件,分别如下function [H,Q,c,s,B]=givensqr(H)[n,m]=size(H);Q=eye(n);for k=1:ma=H(k,k);b=H(k+1,k);[c,s]=givens(a,b);H=gmatmat(H,c,s,k);M=eye(n);M(k,k)=c;M(k,k+1)=-s;M(k+1,k)=s;M(k+1,k...原创 2020-03-01 20:08:30 · 3444 阅读 · 1 评论 -
MATLAB实现矩阵雅可比(Jacobi)迭代
包含三个function,下列代码最后一段是主函数,其他都是function。三个function建议从下往上看。function [x,k,resvec,DD,ID,JD,D,Ab] = jacobis(AA,IA,JA,b,x,tol,kmax)%This function is an implementation of Jacobi’s%iterative method for so...原创 2020-03-01 20:13:08 · 8249 阅读 · 0 评论