MatLab在理工课程中的应用 6.pptVIP

* 求1！+2！+3！+ … +5！的值 clear sum=0; for i=1:5 pdr=1; for k=1:i pdr=pdr*k; end sum=sum+pdr; end sum 运行结果： sum = 153 §5 MATLAB的其他函数库 5.1.1 函数的零值 求解函数的零值的函数fzero fzero(function,x0,TOL) function为方程表达式所在的函数文件，x0是初始值，TOL为误差限。fzero采用的是叠代的方法来求解零值，使得初始估计值接近零值。 我们首先将函数表达式写入一个函数文件中，然后再调用函数 fzero，函数文件tzero.m为： function zexer=tzero(x) zexer=x.^3-3*x+3; 【例】求解函数 的零值。 在命令窗口输入 span=-5:0.1:5; plot(span,tzero(span)),grid result=fzero(tzero,0.3) 由于这种方法要运用试初始值的方法来求出零值点，必须先要对函数的零值点所在区间给出一个估计。故当估计不足或者给出的初始值范围不当时，有可能漏掉零值点。 下面介绍一个求零值点更方便的函数 roots 【例】求解一元三次方程 的根。 一元高次代数多项式 ，在MATLAB中可以用它的系数向量表示 注意：如果x的某次幂的系数为零，这个零必须列入系数向量中。 首先建立表达式的向量p=[1 0 -3 3]，然后直接输入函数命令： p=[1 0 -3 3]; results=roots(p) 5.1.2 最值问题 在MATLAB中求最大最小值的方法是进行迭代计算。 本节介绍两个函数命令fmin和fmins,前者是求单变量函数的最小值，后者是求多变量函数的最小值。开始时应该给出一个初始值，这个初始值为数组。 【例】求解函数 在[-2，3]的最小值。 fmin(‘函数名’，初值xo, 终值xf) 首先建立一个函数文件，然后再调用函数 fmin 这里应该注意最小值和极值的区别，函数给出的是区间中的极小值，需要求最值还应该比较驻点和非驻点以及端点的函数值大小。 MATLAB中没有求最大值的函数，但是可以利用加符号求最小值的方法求最大值。 function exer=ymini(x) exer=(4*x+3)./((x-1.4).^2+3)+(2*x+1)./(x.^3+2) +(4*x+5)./((x-0.9).^2+1); minnum=fmin(@ymini,-2,3) %fminbnd (@ymini,-2,3) 【例】求解函数 的最小值。 function temp=ymins(x) temp=x(1).^2+x(2).^2-2*x(1)+x(2)-2; fmins(@ymins,[1,0]) %fminsearch(@ymins,[1,0]) 5.1.3 绘制函数曲线 命令格式： fplot(‘函数名’，[初值xo, 终值xf]) 【例】画出humps函数在x=0~2之间的曲线。 fplot(humps,[0,2]),grid fplot函数对于快速了解一些复杂特殊函数的波形很有用处。例如画出第一类bessel函数的曲线图 fplot(besselj(alpha,x),[0,10]) 以上给出的是这些函数调用的典型格式，还有其他选项可作变元，例如 fplot(@tan,[-2*pi,2*pi,-2*pi,2*pi],*),grid 还有一个简便画函数图的命令ezplot，它连自变量范围都无需规定，其默认的自变量范围为[-2 π, 2 π]，因此只要键入 ezplot tan(x),grid ezplot tan(sin(x))-sin(tan(x)) 定义y轴的上下限 5.1.4 任意函数的数值积分 （a）定积分子程序quad 其格式为： quad(‘函数名’，初值xo,终值xf) 例如，要画出humps函数在 x=1与2之间的定积分，键入 s=quad(‘humps’,1,2) 我们可以用定积分函数来求不定积分的数值解。只要固定积分下限，用for循环，把积分上限逐步增加即可。 quad、quad8、quadl等调用格式相同，只是内部算法不同。