Matlab函数线性方程组解法和矩阵特征值特征向量.ppt

1. Matlab function-函数 2.线性方程组的解法3.矩阵的特征值和特征向量 张恩维 2007-5-24 函 数-Function 概述 函数文件是另一类m文件。MatLab中函数文件是能够创建局部独立空间的程序文件。函数文件内定立的所有变量名对函数而言是局部的，既不影响其它程序文件和函数文件中同样的变量名又不被其所影响。 所有的MatLab内置函数都是这样。 函 数-Function 使用函数的优点 避免代码重复 限制因改变特定代码段而带来的影响 减少整体代码的复杂性，增加可读性和可管理性 隔离复杂操作 易于调试和隔离错误 使每段程序被最优化，提高整体性能 函 数-Function 编写函数的原则 即一个函数具有特定的功能。也就是说，它做一件事情而不是其它任何事情。 函数的基本格式 function 返回变量=函数名（输入变量） 函数的命名原则 直观： 与程序功能有关系 长度： 5~20个字符之间 函 数-Function 函数文件的编写 格式： function OutputVariable=FunctionName(InputVariables) % comments 注释段 expression(s) 函数m文件第一行必须以单词function作为引导词 m文件的文件名必须是函数名 .m. （FunctionName.m） 程序中的变量均为局部变量，不保存在工作空间中。其变量只在函数运行期间有效。 MatLab用紧跟函数接口语句之后的注释创建此函数的帮助信息—当在MatLab命令行中输入help FunctionName时，所有相邻的注释行都将出现在MatLab命令窗口中。 例如：计算如下等式 x=cos(at)+b, y=|x|+c. function zanswer=ComputeXY(t,a,b,c) % Computation of x=cos(at)+b, y=|x|+c. % Scalars: a, b, c . Vector: t % Returned Matrix: zanswer x=cos(a*t)+b; zanswer=[x; abs(x)+c]; 此时当输入 help ComputeXY Computation of x=cos(at)+b, y=|x|+c. Scalars: a, b, c . Vector: t Returned Matrix: zanswer 调用 vick=ComputeXY(0:pi/4:pi,1.4,2.0,0.75); 返回 vick = 3.0000 2.4540 1.4122 1.0123 1.6910 3.7500 3.2040 2.1622 1.7623 2.4410 两种特殊情况 只有输入没有输出 function FunctionName(InputVariables) 2. 没有输入只有输出 function OutputVariable=FunctionName 函 数-INLINE INLINE 创建局部函数 inline函数的优点：不必存储为单独文件。 格式： FunctionName=inline(‘expression’,’p1’,’p2’,…) Expression – 函数中的表达式 p1, p2… 参数，即表达式中变量的名称。 inline的规则： 1.它不能再调用另一个inline函数； 2. 只能由一个MatLab表达式组成 并只能返回一个变量。 例如创建FofX函数求解下式 f(x)=x2cos(ax)-b FofX=inline(x.^2.*cos(a*x)-b,x,a,b) FofX = Inline function: FofX(x,a,b) = x.^2.*cos(a*x)-b G=FofX([pi/3 pi/3.5],4,1) G = -1.5483 -1.7259 例 求满足下式的a 值 function z=suma(a) z=sum(1./([1:1000].^2-a)); 初始值设为 pi/2 options=optimset(display,off); fofa=fzero(‘