matlab7第8章.pptVIP

MATLAB 7.0从入门到精通 主要讲述内容 第1章 MATLAB简介 第2章 数值运算 第3章 单元数组和结构 第4章 字符串 第5章 符号运算 第6章 MATLAB绘图基础 第7章 程序设计 第8章 计算方法的MATLAB实现 第9章 优化设计 第10章 Simulink仿真初探 第8章 计算方法的MATLAB实现 8.1 方程求根 roots见多项式求根；roots(多项式系数矩阵) fzero可求解非线性方程 fzero(‘function’,x0)其中function为求解的方程，x0为估计的根，x0可为标量或长度为2的向量，为向量时函数的两端的值应该符号相反，此时求区间上的解。只能求解x0附近的一个解。即使在某个区间内有多个解，但是区间端点符号相同的话仍然出错。 fzero(x^3-3*x-1,2) ans = 1.8794 fzero(x^3-3*x-1,[2,4]) ??? Error using == fzero The function values at the interval endpoints must differ in sign. fzero(x^3-3*x-1,[1,4]) ans = 1.8794 fzero(x^2-3*x+2,0) ans = 1.0000 fzero(x^2-3*x+2,3) ans = 2.0000 fzero(x^2-3*x+2,[0,4]) ??? Error using == fzero The function values at the interval endpoints must differ in sign. 8.2 线性方程组数值解法 1、直接解法 A*X=B X=A\B或x=inv(A)*B 例： 12x(1)-3x(2)+3x(3)=15 -18x(1)+3x(2)-x(3)=-15 X(1)+x(2)+x(3)=6 a=[12 -3 3;-18 3 -1;1 1 1]; b=[15;-15;6]; x=a\b x = 1.0000 2.0000 3.0000 inv(a)*b ans = 1 2 3 2、线性方程组求解中的变换 上三角变换 U=triu(x)返回矩阵x的上三角部分； U=triu(x,k)返回第k条对角线以上部分的元素。 a=[12 -3 3;-18 3 -1;1 1 1]; triu(a) ans = 12 -3 3 0 3 -1 0 0 1 triu(a,1) ans = 0 -3 3 0 0 -1 0 0 0 triu(a,-1) ans = 12 -3 3 -18 3 -1 0 1 1 对角变换 U=diag(x)返回矩阵x主对角线上的元素，返回结果是一列向量形式； U=diag(x,k)返回第k条对角线上的元素值。 当x为向量时生成矩阵。 a=[12 -3 3;-18 3 -1;1 1 1]; diag(a) ans = 12 3 1 a=[12 -3 3;-18 3 -1;1 1 1]; diag(a,1) ans = -3 -1 a=[12 -3 3;-18 3 -1;1 1 1]; diag(a,-1) ans = -18 1 下三角变换 U=tril(x)返回矩阵x的下三角部分； U=tril(x,k)返回第k条对角线以上下部分的元素。 a=[12 -3 3;-18 3 -1;1 1 1]; tril(a) ans = 12 0 0 -18 3 0 1 1 1 tril(a,1) ans = 12 -3 0 -18 3 -1 1 1 1 tril(a,-1) ans = 0 0 0 -18 0 0 1 1 0 3、迭代解法 迭代解法非常适用于求解大型稀疏系数矩阵的方程组，在线性方程组常用的迭代解法主要有Jacobi迭代法、Gauss-Seidel迭代法。 Jacobi迭代法 编写jacobi函数，并使之计算 a=[10 -2 -1;-2 10 -1;-1 -2