matlab程序总结.docVIP

程序总结 1、简单计算A:x255=x·x···x B:x255=x·x2·x4·x8·x16·x32·x64·x128 2. 算法1 (输入a(i)(i=0,1，…,n),x;输出 算法2 （秦九韶算法） 3. s=0; for i=1:n s=s+abs(x(i)); End 4. s=0; s=0; for i=1:n for i=1:n s=s+x(i)^2; s=s+x(i)*x(i); end end s=sqrt(s) s=sqrt(s) 5. s=0; for i=1:n if abs(x(i))s,s=abs(x(i)); end end 6. LU分解的matlap程序 function A=lud(A) %功能：对方阵A作三角分解A=LU,其中， % L为单位下三角阵，U为上三角阵， %输入：方阵A。 %输出：紧凑存储A=[L\U]. %注意：当A的主元=0时退出Matlab for k=1:n-1 for i=k+1:n if A(k,k) ==0 quit; end A(i,k) =A(i,k)/ A(k,k); A(i,k+1:n)= A(i,k+1:n)- A(i,k) *A(k,k+1:n); end end 7. 列主元Gauss消元法 Lupd.m %功能：对方阵A作列主元三角分解PA=LU,其中， % L为单位下三角阵，U为上三角阵，排列阵P %用向量p表示。 %输入：方阵A。 %输出：紧凑存储LU=[L\U]，以及p。 %注意：当A奇异时退出Matlab. function [LU,p]=lupd(A) %初始化 n=length(A); p=1:n; LU=A; %分解过程 for k=1:n %搜索列主元ik [s,i]=max(abs(LU(k:n,k))); ik=i+k-1; %判断矩阵的奇异性 if s==0 quit; end %行交换 if ik~=k m=p(k); p(k)=p(ik); p(ik)=m; lk=LU(k,:); LU(k,:)=LU(ik,:); LU(ik,:)=lk; end %用消元法计算LU=[L\U] if k==n break; end LU(k+1:n,k)=LU(k+1:n,k)/LU(k,k); LU(k+1:n,k+1:n)=LU(k+1:n,k+1:n)-LU(k+1:n,k)*… LU(k,k+1:n); End 8. Householder矩阵变换 function [H,y]=holder1(x) n=length(x); M=max(abs(x)); if M==0, disp(‘M=0); return; else z=x/M; end; s=norm(z); if z(1)0 s=-s; end p=s*(s+z(1)); u=z; u(1)= s+z(1); H=eye(n,n)-p\u*u; y=zeros(n,1); y(1)=-M*s; 9、解上三角方程 function X=backsub(A,b) %Input—A is an n×n upper- triangular nonsingullar matrix % ---b is an n×1 matrix %Output—X is the solution to the system AX=b n=length(b); X=zeros(n,1); X(n)=b(n)/A(n,n); for i=n-1:-1:1 X(i)=(b(i)-A(i,i+1:n)* X(i+1:n))/A(i,i); End 10、matlap中高斯消元法 function X=gauss(A,b) %Input—A is an n×n nonsingullar matrix % ---b is an n×1 matrix %Output—X is the solution to the system AX=b [n n]=size(A);