第8部分MATLAB代数方程求解学习课件.pptVIP

第8章MATLAB代数方程求解 * 8.1 MATLAB的线性方程求解 8.1.1 直接法 例 8-1 编制如下命令文件 clear;clc; A=[1, 1, 1, 1; 1, 2, 1, 1; 2, 2, 5, 1; 4, 2, 1, 6]; B=[4; 5; 10; 13]; X1=A\B X2=linsolve(A,B) 运行结果为 X1 =1.00001.00001.00001.0000 X2 =1.00001.00001.00001.0000 例 8-2 编制如计算程序 clear; clc; A=[1,1,2,4;1,2,2,2;1,1,5,1;1,1,1,6]; B=[4, 5, 10, 13]; X=B/A 运行结果为 X =1.0000 1.0000 1.0000 1.0000 8.1.2 迭代法 例 8-3 编制如下函数文件fun.m function [X,n] = fun(A,B,X0) D=diag(diag(A)); LU=A-D; X=D\B-D\LU*X0; n=1; while norm(X-X0)=1.0e-6X0=X;X=D\B-D\LU*X0;n=n+1; end 编制如下计算程序 clear;clc; A=[10,-2,0;-1,10,-2;0,-2,10]; B=[9;8;7]; X0=[