- 博客(4)
- 收藏
- 关注
RSS订阅原创 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
1604
原创 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
1140
原创 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
1365
原创 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
329
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人