笔记
Krarey
他很懒,懒到什么都没留下!?
展开
-
Matlab实现乘幂法
function [u,r,res]=myPM(A,u)res=1;while (res>1e-6) y=A*u; [a,b]=max(abs(y)); r=y(b,1); u=y/a; res=norm(A*u-r*u);end原创 2021-05-30 08:21:02 · 1467 阅读 · 0 评论 -
Matlab实现高斯消去解方程组
function [x,res]=MyLU(A,b)AA=A;[m,n]=size(A);x=zeros(n,1);y=zeros(n,1);for k=1:n-1 A(k+1:n,k)=A(k+1:n,k)/A(k,k); A(k+1:n,k+1:n)=A(k+1:n,k+1:n)-A(k+1:n,k)*A(k,k+1:n);endl=tril(A);for i=1:n l(i,i)=1;endu=triu(A);for i=1:n sum=0;原创 2021-05-30 08:19:12 · 1045 阅读 · 0 评论 -
Matlab实现Cholesky分解解方程组
function [x,res]=mych(A,b)AA=A;[m,n]=size(A);y=zeros(n,1);x=zeros(n,1);for k=1:1:n A(k,k)=sqrt(A(k,k)); A(1:k-1,k)=0; A(k+1:n,k)=A(k+1:n,k)/A(k,k); for j=k+1:n A(j:n,j)=A(j:n,j)-A(j:n,k)*A(j,k); endend U=A';L=A; y(1,原创 2021-05-30 08:13:50 · 1251 阅读 · 0 评论 -
Matlab实现列选主元高斯消去解方程组
function[x,res]=myPLU(A,b)AA=A;bb=b;[m,n]=size(A);A=[A,b];y=zeros(n,1);x=zeros(n,1);for k=1:n-1 [m,p]=max(abs(A(k:n,k))); p=p+k-1; c=A(p,:); A(p,:)=A(k,:); A(k,:)=c; if A(k,k)~=0 A(k+1:n,k)=A(k+1:n,k)/A(k,k); A(k+原创 2021-05-30 08:09:23 · 280 阅读 · 0 评论