MATLAB入门教程-updated20131008汇编.ppt

例 * 5.4.3 积分 int f 求表达式f对默认自变量的积分值； int f, t 求表达式f对自变量t的不定积分值； int f, a, b 求表达式f对默认自变量的定积分 值，积分区间为[a,b]； int f, t, a, b 求表达式f对自变量t的定积分值，积 分区间为[a,b] * 例 * 5.5 方程求解 5.5.1代数方程 代数方程的求解由函数solve实现： solve f 求解符号方程式f solve f1,…,fn 求解由f1,…,fn组成的代数方程组 5.5.2常微分方程 使用函数dsolve来求解常微分方程： dsolve eq1, eq2, ..., cond1, cond2, ..., v * 例 syms a b c x f sym a*x*x+b*x+c 0 solve f ans [ 1/2/a* -b+ b^2-4*c*a ^ 1/2 ] [ 1/2/a* -b- b^2-4*c*a ^ 1/2 ] solve 1+x sin x ans -1.座机电话号码0座机电话号码537689 dsolve Dy x ,x %求微分方程y x的通解，指定x为自变量。 ans 1/2*x^2+C1 dsolve D2y 1+Dy ,y 0 1,Dy 0 0 %求微分方程y 1+y的解，加初始条件 ans -t+exp t [x,y] dsolve Dx y+x,Dy 2*x %微分方程组的通解 x -1/2*C1*exp -t +C2*exp 2*t y C1*exp -t +C2*exp 2*t * Thank You * * * 应用举例 例 将表达式 x-4 x+5 x2-6x+9 展开为多项式形式，并求其对应的一元n次方程的根。 p conv [1 -4],conv [1 5],[1 -6 9] px poly2str p,’x’ x roots p * 四 脚本编程基础 对于简单问题，使用直接输入命令简单有效；对稍复杂和多次重复的应用，直接输入命令比较麻烦。MATLAB提供了逻辑解决方案，它允许用户把多个命令放在一个简单的文本文件中，如同在MATLAB中键入命令一般，这种文件称为脚本文件，由于脚本文件以m为扩展名，它常称为M文件。脚本文件为文本形式的，对跨平台处理十分有利。 使用脚本文件，可以把命令保存在磁盘上，便于以后的访问；同时对使用大的数组也带来的方便；增加注释可以为脚本中的命令作文挡以免以后忘记。 * 4.1 流程控制 为了便于应用，MATLAB提供了一些流程控制的命令。这些命令对脚本编写带来了一些方便，但是需要注意的是，尽量不要使用这些流程控制命令，尤其是循环控制命令。 * 4.1.1 for 循环 for循环允许一组命令以固定的次数重复，例如 for i 1:5; for j 1:4; A i,j 1/ i+j-1 ; end end * For循环应当注意的一些方面 for循环不能使用内部重新赋值循环变量而终止； for循环内部接受任何有效的MATLAB数组； for循环可以嵌套； 只要有矩阵形式可以解决的问题，不要使用for循环。使用for循环的算法执行很慢，一个好的MATLAB算法不应当出现循环语句。Tic/toc 循环可以使用break跳出，但只跳出所在的循环，不跳出整个嵌套结构。 * 4.1.2 while 循环 与for循环以固定的次数执行一组指令，while循环以不定的次数执行一组指令。While循环的一般形式为： while expression commands end 只要表达式expression里的所有元素为真，就执行命令串commands。通常表达式求值给一个标量值，单数组值也同样有效。 * 4.1.3 if-else-end 结构 很多情况下，命令的序列需根据关系的检验有条件的执行，它由if-else-end结构提供。它的结构如下： if expression1 commands1 elseif expression2 commands2 elseif … … else commands end * 4.1.3 if-else-end 结构（续） 在执行过程中，MATLAB依次检查各个表达式，只执行第一个表达式为真的命令串，接下来的关系表达式不检验，跳过其余的if-else-end结构，而且，最后的else命令可有可无。 * 4.1.4 try-catch 结构 可以进行错误的处理，提高计算的可靠性 try 可能出错的语句 catch 错误处理 end *