Matlab与化工数值计算-第一讲简介与基本数学运算分解.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Matlab关系和逻辑运算 MATLAB关系运算符 运 算 符 号 运 算 符 号 大于 小于 等于 == 不等于 ~= 大于等于 = 小于等于 = 关系运算符： 逻辑、关系运算优先级 运算优先级： 括号 转置，幂 一元的加减，逻辑否 乘除 加减 冒号操作符 关系运算 逻辑运算 先决与 先决否 高 低 脚本文件和函数文件 脚本编辑窗口 函数文件和脚本文件 脚本文件 脚本仅仅是一连串可执行的MATLAB命令，它具有全局性 脚本文件中不能定义函数 函数文件 函数定义的一般格式： function [y1,y2,…,yn] = FuncName(x1,x2,…,xn) %函数声明语句 y1 = … % （表达式1） y2 = … % （表达式2） … yn = … % （表达式n） 其中，输入参数为x1,x2,…,xn,输出参数为y1,y2,…,yn。各参数可以是标量、向量或矩阵。 函数文件的编写与运行 函数声明语句：function Y=function_name(x) 编写表达式：对x进行加减乘除等运算，然后赋值给Y。 保存在当前路径下，文件名与函数名要相同。 在命令窗口中输入：文件名（参数值）。按回车后开始运行。 Function out=test（in） out=in^2-3*in-4 编写一个函数文件，其输出参数是其输入参数的10倍。（写上姓名学号后上交） 课堂作业 函数的调用 在调用函数的主函数中，直接采用函数名调用 调用函数得到函数的返回参数值。 function avg=average(x) %求矢量中各数值的平均值 n=length(x) %该函数可以得到矢量含多少个元素。 avg=mean(x,n) function a=mean(u,d) %此为子函数 a=sum(u)/d %sum命令可以求得矢量各元素的和。 2.通过函数句柄调用 调用函数本身。 Function a=result(x) a=feval(@zwb,x) %feval命令的意思是获得zwb这个函数在x处的函数值。 function y=zwb(x) %此为子函数 y=x^3+x^2-3*x-3; 内联函数（inline function） 内联函数是Matlab提供的一个对象，它的表现和函数文件一样，但内联函数的创建比较容易 内联函数的创建 inline(‘CE’) y= inline(x^3-2) 例如：y(3)=25 inline(‘CE’,arg1,arg2,...) y=inline(sin(2*pi*f + theta), f, theta) 匿名函数（anonymous function） 匿名函数用于在命令行、函数文件或脚本文件中创建简单形式的函数，避免另外定义新的函数。 匿名函数的定义形式 f=@(arglist)expression f=@(x) x.^2 a=f(5) 结果：a＝25 f=@(x) x.^2; g=@(x) 3*x; h=@(x) g(f(x)); h(3) 结果：ans＝27 alpha=0.9; f=@(x) sin(alpha*x)； f(pi) 结果：ans＝0.3090 Matlab流程控制语句 for循环结构 for i=1:10 x(i)=i; end x for 循环变量 = 表达式1（初值）：表达式2（步长）：表达式3（终值） statements （语句组） end字符串：包含在一对单引号中的字符集合 为了得到高效代码，应尽量提高代码的向量化程度，避免使用循环结构 为了得到高效代码，在循环指令之前应尽量对数组进行预定义 while循环结构 while condition（表达式） statements（执行语句组） end Fibonacci数组的元素满足Fibonacci规则：ak+2=ak+ak+1，（k=1,2,...);且a1=a2=1。求该数组中第一个大于10000的元素。 a(1)=1;a(2)=1;i=2; while a(i)=10000 a(i+1)=a(i-1)+a(i); i=i+1; end i, a(i) if-else-end分支结构 if语句的一般格式： if condition1 statements %如果condition1的值为True，则执行该语句组 elseif condition2 statements