03微分及积分方程等分解.ppt

六、M文件调试 　应用程序的错误有两类：一类是语法错误；另一类是逻辑错误，指程序运行过程中得到的结果不是用户需要的结果。 与调试(Debug)有关的菜单项及按钮： Set/Clear Breakpoint：设置/清除光标处的断点； Clear All Breakpoint：清除程序中的所有断点； Stop if Error/ Warning…：运行至…结束(自动断点)； Run：执行程序至结束遇到断点； step：单步执行程序； Step In：深入下层局部工作区； Continue：恢复程序运行至结束或另一个断点。 1、调试器的使用 指令 说明 dbclear 清除已经设置好的断点 dbcont 在断点后继续执行，等同于工具栏中的 按钮 dbdown/dbup 修改当前工作空间的上、下文关系 dbquit 退出调试状态 dbstack 显示谁调用谁 dbstatus 显示所有的已经设置的断点 dbstep 执行应用程序的一行或者多行代码 dbstop 设置断点 dbtype 显示M文件代码和相应的行号 2、调试命令 dbstop：设置断点；调用格式有：①dbstop at 行号 in 文件名：在指定文件的指定行处设置断点；② dbstop in 文件名：在指定文件的第一个可执行处设置断点； ③ dbstop if error：在任何运行出错时产生中断（不包括语法错误）。出错中断后用户不能恢复程序； ④ dbstop if naninf/infnan：当出现Inf或nan值时产生中断。 dbcont：用来恢复对于执行dbstop指令而导致中断的程序。一旦遇到其它已设置的断点或出现错误，就返回基本工作空间； dbstep：用于执行一行或多行代码。在调试模式下，dbstep允许用户实现逐行跟踪； dbstack：列出调用关系。在调试模式下，此指令显示当前的文件名及中断点的行号，其调用者的、调用者的调用者的…文件名和断点行号，直到最外层的文件。 dbstatus：用来列出全部断点。dbstatus （filename）可旬出指定文件的全部断点。 dbtype：用来显示带行号的文件内容，以协助用户设置断点，它有两种格式：①dbtype (文件名)：列出指定文件的全部内容，每行前带行号；② dbtype 文件名 始行号: 末行号：列出指定文件自起始行号到末行行号部分的内容。 dbquit：用来退出调试模式，将控制转向基本工作空间。此时，文件执行还没有完成，也不产生返回结果。 常用警示信息： error：显示预先设置的错误信息，程序停止运行； warning：显示预先设置的警告信息，程序继续运行； lasterr：显示MATLAB自动判断的最新出错原因并终止程序运行； lastwarn：显示MATLAB自动给出的最新警告提示，程序继续运行。 1.5.4 符号方程求解 一、符号代数方程求解 在Matlab中，利用solve函数进行对方程组求解。其调用格式 有以下几种： solve (eq)：计算单一方程，eq为输入的方程。 solve (eq, var)：计算单一方程，eq为输入的方程，var用来指明待求变量。 solve (eq1, eq2, …, eqn)：计算方程组，eq1,…eqn为输入的方程组。变量为默认变量。 solve (eq1, eq2, …, eqn, var1, …, varn)：计算方程组，eq1,…eqn为输入的方程组；var1,…varn用来指明每个方程待求的变量。 syms x a; f=3*x^2-3*a^2=0; x=solve(f,x) %程序运行结果为： x = -a a 例：求方程 的解。 [x,y]=solve(3*x+y=a,x-y=a,x,y) %程序运行结果为： x = 1/2*a y = -1/2*a 例：求方程组 的解。 f=sin(x)+y^2+log(z)-7 g=3*x+2^y-z^3+1 h=x+y+z-5 [x,y,z]=solve(f,g,h) %程序运行结果为： x = 0.5991 y = 2.3959 z = 2.0050 例：求方程组 的解。 二、微分方程求解 dsolve (‘eq1’, ‘eq2’, ‘…’ , ‘cond1’,‘cond2’,‘…’, ‘ v1’, ‘v2’ , ‘... ’ ) 求解常微分方程组eq1, eq2,…在初始条件cond1, cond2, …下的特解。参数v描述方程中的自变量，默认变量为t。若没有给出初