Matlab 第四讲：符号运算-数学软件 Matlab.ppt

数学软件 Matlab 主要内容 符号运算 Matlab 符号运算 符号运算举例 内容提要 符号对象 符号对象的建立 符号对象的建立 符号表达式 基本符号运算 符号矩阵 内容提要 symvar subs subs 举例 vpa 内容提要 因式分解 函数展开 合并同类项 函数简化 函数简化 函数简化 horner 多项式 求反函数 计算极限 计算导数 计算积分 符号级数求和 代数方程求解 微分方程求解 dsolve 求解析解 dsolve 的使用 dsolve 的使用 dsolve 的使用 上机作业 上机要求 求微分方程数值解*（自学，选学） Matlab的ODE求解器*（自学，选学） 参数说明*（自学，选学） 数值求解举例*（自学，选学） 数值求解举例*（自学，选学） solve(f,v) % 求方程关于指定自变量的解 详细用法见Matlab08：多项式运算与代数方程求解器 代数方程求解 这里 f 可以用字符串表示或符号表达式 solve 也可解方程组（通常是非线性的） 得不到解析解时，给出数值解 例： solve(2*x-3) % 或 solve(2*x-3=0) syms x; solve(2*x-3) % 不能写成 solve(2*x-3=0) syms x; solve(2*x-sin(x)+1) 用 Maltab自带函数 解初值问题 求微分方程解析解：dsolve 求微分方程数值解*（自学，选学）：  ode45、ode23、 ode113、ode23t、ode15s、 ode23s、ode23tb dsolve 的使用 y=dsolve(eq1,eq2, ... ,cond1,cond2, ... ,v) 其中 y 为输出， eq1、eq2、...为微分方程，cond1、cond2、...为初值条件，v 为自变量。 例 1：求微分方程 的通解，并验证。 y=dsolve(Dy+2*x*y=x*exp(-x^2),x) syms x; diff(y)+2*x*y - x*exp(-x^2) 几点说明 如果省略初值条件，则表示求通解； 如果省略自变量，则默认自变量为 t dsolve(Dy=2*x,x); ％ dy/dx = 2x dsolve(Dy=2*x); ％ dy/dt = 2x 若找不到解析解，则返回其积分形式。 微分方程中用 D 表示对 自变量 的导数，如： Dy y； D2y y； D3y y 例 2：求微分方程 在初值条件 下的特解，并画出解函数的图形。 y=dsolve(x*Dy+y-exp(x)=0,y(1)=2*exp(1),x) ezplot(y); dsolve 的使用 例3：求微分方程组 在初值条件 下的特解，并画出解函数的图形。 [x,y]=dsolve(Dx+5*x+y=exp(t),Dy-x-3*y=0, ... x(0)=1, y(0)=0, t) ezplot(x,y,[0,1.3]); 注：解微分方程组时，如果所给的输出个数与方程个数相同，则方程组的解按词典顺序输出；如果只给一个输出，则输出的是一个包含解的结构（structure）类型的数据。 例： [x,y]=dsolve(Dx+5*x=0,Dy-3*y=0, ... x(0)=1, y(0)=1,t) r = dsolve(Dx+5*x=0,Dy-3*y=0, ... x(0)=1, y(0)=1,t) 这里返回的 r 是一个 结构类型 的数据 r.x %查看解函数 x(t) r.y %查看解函数 y(t) 只有很少一部分微分方程（组）能求出解析解。大部分微分方程（组）只能利用数值方法求数值解。 dsolve的输出个数只能为一个 或 与方程个数相等 教材P33的习题6~10，写入文件m04_1.m； 2. 教材P32的习题4，写入文件m04_2.m； 3. 求解下列微分方程（组），写入文件m04_3.m，并且对于初值问题还要求画出解函数的图形；(5)(6)选做。 将完成每题所用的命令写入一个规定文件名的文件中 然后将这些文件作为附件，通过 foxmail 以邮件形式发给 admin@system.mail 邮件主题为：机号-学号-姓名 其中机号为 两位数 三个字段之间用英文状态下的减号链接 上机