08-函数文件及应用.pptVIP

MATLAB程序设计 MATLAB函数文件及应用 8.1 M文件的建立 1）进入程序编辑器（MATLAB Editor/Debug） 2）输入程序 3）保存程序 8.2命令M文件及其运行 命令M文件：一些命令行的组合。 例如：% 文件名ex1.m x=4;y=6;z=2; items=x+y+z cost=x*25+y*22+z*99 average_cost=cost/items 1）在程序编辑窗口输入, 以名为ex1.m的M文 件存磁盘 2）在MATLAB命令编辑窗口键入ex1，便会顺序执行该文件中的命令。 8.3函数M文件及其调用 函数M文件的形式为 function [输出变量列表]=函数名（输入变量列表） 函数体语句； 例如： function f=f1(x) f=100*(x(2)-x(1)^2)^2-(1-x(1))^2; 调用格式： x=[2 3]; f=f1(x) 结果为：f = 99 8.3函数M文件及其调用 又如： function [F,G]=f2(x) F=2*x(1)^2+2*x(2)^2-2*x(1)*x(2)-4*x(1)-6*x(2); G=[x(1)+5*x(2)-5 2*x(1)^2-x(2); -x(1) -x(2)]; 调用格式： x1=[4 5]; [F1,G1]=f2(x1) 结果为：F1 = -4 G1 = 24 27 -4 -5 无输出参数时，如： function printresults(x) 或 function [ ]=printresults(x) 8.4匿名函数 7.0以后引入的新函数类型 替代内联函数(inline) 8.4匿名函数 在程序编写中，经常需要计算某些数学表达式的值，如何方便构造函数？写成M文件调用当然可以，但不够简洁 有时通过符号计算得到了一个函数表达式，接下来想通过数值计算对这个表达式进行进一步的操作，如求值，数值积分等，如何进行 很多时候要对含参变量的函数进行操作，譬如f(x)=xa,如何在得到a的具体值后再得到f(x)的具体函数表达式？ 8.4匿名函数 8.4.1匿名函数的具体定义 fhandle=@(arglist) expr 如 f=@(x) x.^2 %单变量匿名函数 f = @(x)x.^2 fx=f(1:10) fx = 1 4 9 16 25 36 49 64 81 100 8.4匿名函数 8.4.1匿名函数的具体定义 fhandle=@(arglist) expr 如 g=@(x,y) x.^2+y.^2 %多变量匿名函数 g = @(x,y)x.^2+y.^2 gxy=g(1:10,2:11) gxy = 5 13 25 41 61 85 113 145 181 221 8.4匿名函数 8.4.2多重匿名函数 fhandle=@(arglist1) @(arglist2) expr 如 f=@(a,b) @(x) a*x+b f = @(a,b)@(x)a*x+b a,b是外层变量，x是内层变量 gab=f(a,b)是一个以x为变量的单变量匿名函数. 8.4匿名函数 8.4.3匿名函数在求解方程中的应用 1 求下列方程的根 。 f=@(x) exp(x)+x^2+x^(sqrt(x))-100 f = @(x)exp(x)+x^2+x^(sqrt(x))-100 x0=fzero(f,3) x0 = 4.1635 f(x0) ans = 2.8422e-014 求下列方程的根 。 8.4匿名函数 2 对于α=[0,0.01,0.02,…,2]，求下列方程相应的x值，并画出a和相应x的图像。 f=@(a) @(x) exp(x)+x^a+x^(sqrt(x))-100 f = @(a)@(x)exp(x)+x^a+x^(sqrt(x))-100 aa=0:0.01:2; plot(aa,arrayfun(@(a) fzero(f(a),4),aa),*-) 求下列方程的根 。 8.4匿名函数 8.4匿名函数 8.4.4匿名函数在显示表示隐函数方面的应用 1 显示表示