MATLAB语言程序设计教程（第二版）张德喜赵磊生第3章MATLAB数值计算.pptVIP

第3章 MATLAB数值计算 学习目标 了解傅里叶变换、数据的数值运算； 理解曲线的拟合与插值； 掌握线性方程和线性方程组的数值求解、多项式的数值运算、数值的微分与积分及其微分方程的求解。 3.1 线性方程与线性方程组的数值求解 MATLAB可以很迅速地求出线性方程和线性方程组的数值解 求解方法大致分为两大类 直接法 叠代法 3.1.1 线性方程数值求解 线性方程的求根运算可以直接通过调用求根函数roots()来求解。 【例3.1】求解的根。 解 在MATLAB命令提示符下输入： clear p=[1 -1 -6]; roots(p) ans=3 -2 3.1.2 线性方程组数值求解 1．直接法 求解方程AX=B，只需输入A、B，然后做矩阵左除X=A\B，即可输出X为方程组的解。 【例3.2】用直接法线性方程组ax=b的解。 解 在MATLAB命令提示符下输入： clear a=[1 1 1 1;1 2 1 1;1 1 1 3;1 1 4 1] a=1 1 1 1 1 2 1 1 1 1 1 3 1 1 4 1 b=[1 1 1 1] b=1 1 1 1 x=a\b x=1 0 0 0 2．利用矩阵三角分解LU、正交分解QR和特征值分解求线性方程组的解 矩阵的三角分解LU分解、正交分解QR分解和特征值分解已经在第2章讲过，在此不再细讲。例如经过LU分解，方程组AX=B的系数矩阵A可以化为A=LU，则方程组的解为：X= U-1（L-1B） 【例3.3】用矩阵的三角分解求线性方程组ax=b的数值解。 解 在MATLAB命令提示符下输入： clear a=[4 5 8 7;8 2 1 3;9 1 0 2;6 5 7 3] a=4 5 8 7 8 2 1 3 9 1 0 2 6 5 7 3 [l,u]=lu(a) l=0.4444 1.0000 0 0 0.8889 0.2439 1.0000 0 1.0000 0 0 0 0.6667 0.9512 0.6410 1.0000 u=9.0000 1.0000 0 2.0000 0 4.5556 8.0000 6.1111 0 0 -0.9512 -0.2683 0 0 0 -3.9744 b=[ 7 6 2 1] b=7 6 2 1 x=inv(u)*(inv(l)*b) x=-0.7419 4.7032 -3.4323 1.9871 x=a\b x=-0.7419 4.7032 -3.4323 1.9871 由上例可以看出inv(u)*(inv(l)*b) 和a\b计算的结果是一样的。 3.2 非线性方程与非线性方程组的数值求解 与线性方程和线性方程组的求解相比，非线性方程与非线性方程组的求解情况要复杂很多。 通过调用函数solve()和fsolve()可以分别实现非线性方程与非线性方程组的数值求解。 3.2.1 非线性方程数值求解 通过调用函数solve()可以对非线性方程进行数值求解。 【例3.4】求解方程。 解 在MATLAB命令提示符下输入： clear solve(x-sin(x)-cos(x)=0) ans=1.2587281774926764586391391659652 3.2.2 非线性方程组数值求解 非线性方程组数值求解函数为fsolve()。调用格式为： X=fsolve(fun,x0) 即求函数fun()在x0处（或x0附近）的解。 【例3.5】求圆 和直线 的两个交点。 建立方程组函数文件fxyz2.m： function F=F(X) x=X(1);y=X(2);z=X(3); F(1)=x^2+y^2+z^2-12; F(2)=6*x+2*y+z; F(3)=2*x-5*y-z-3; 解 在MATLAB命令提示符下输入： clear X1=fsolve(fxyz2,[2,1,