MATLAB科学计算4讲述.ppt

insert-calculation 插入新的运算提示符 了解按钮功能 了解书写格式 如果有误，给出蓝色错误信息 使用的字符不要与符号运算中的系统保留字相同 例如，大写D是系统中的微分符号 注意 limit(#f, #x=#a) #f用表达式替换 #x用自变量替换 #a用具体数字替换 一、极限运算 sum(#f, #n=#a..#b) #f用表达式替换， #n用自变量替换， #a和#b分别用求和的起止数字或符号表达式替换 二、级数求和 三、求导 diff(#f, #x) #f 用函数用表达式替换 #x 用自变量替换 一阶导数 diff(diff(#f, #x), #x) diff(#f, #x, #y, …) diff(#f, x$n) 高阶导数 四、积分 int(#f, #x=#a..#b) #f用表达式替换 #a和#b分别用上下限替换 五、解常微分方程 solve(ode(#, #(#))) 其中ode({eq, initials}, y(x)) eq: 微分方程， = 导数表达为 y’(x) ,y’’(x) 高阶diff(y(x), x$a) ode({eq, initials}, y(x)) initials-定解条件 缺省，求通解 ode({eq, initials}, y(x)) y(x)-待求函数 解常微分方程组 ode({eq1,eq2,...,initials},{y1(x),y2(x),... }) initials-定解条件 缺省求通解 例1 求 和 例4 例7 例10 符号运算 符号运算功能函数 Mupad Notebook MATLAB 之 Notebook 在Word环境中使用MATLAB 作业 1、分别用符号运算功能函数 和符号运算工具箱求解 2、求一阶、二阶差分 5.2 符号运算功能函数 函数名 功能 函数名 功能 limit 求极限 ode 解微分方程 diff 求导 fourier Fourier变换 int 求积分 laplace Laplace变换 taylor Taylor展开 jacobian Jacobian矩阵 symsum 级数求和 f—f(x)的符号表达式； x—自变量； a—缺省值为0，可以inf或-inf； right or left-表示右极限或左极限 limit(f, x, a, right or left ) 一、limit 例1 limit(f, x, a, right or left ) f —通项表达式； n —级数中的项数变量； h和l —项数取值范围,步长是1 symsum(f, n,h, k) 二、symsum f —待展开的函数 x —变量名； a —展开点，缺省值为0，为迈克劳林展开（Maclaurin series expansion） value —截断阶数 taylor(f, x, a, Order,value) 三、taylor 例2 用Taylor展开法, 将sinx在x0=0处展开为五次幂级数。 taylor(f, x, a, Order,value) 对f中的变量x求n阶导数 n=1可以省略不写 四、diff differential diff(f, x, n) 例3 求sinx的3阶导数 diff(f, x, n) 求 和 例4 diff(f, x, n) 求 例5 已知 diff(f, x, n) polyder 多项式求导 derivation f: 被积函数； x: 积分变量，必须是符号变量； a, b：积分限，可以是inf；缺省时，给出不定积分，输出省略了积分常数 int(f, x, a, b) 五、int int(f, x, a, b) 例6 例7 vpa(p,d) subs(f,old,new) 例7 int(f, x, a, b) 例8 eqn-待解微分方程； conds-定解条件，条件之间用逗号间隔。 y=dsolve(eqn, conds) 六、dsolve 解微分方程 eqns-待解微分方程组，方程之间逗号间隔； conds-定解条件，条件之间逗号间隔； y1,...,yN- 解得的若干函数的名称， 依据eqns中函数名的英文字母先后顺序输出 [y1,...,yN]=dsolve(eqns, conds) 解微分方程组 1、n 阶导数的表示 diff(y,m) 表示 2、常微分方程eqn的表达 y=dsolve(eqn, conds) 3、初始条件的表示 首先定义y的各阶导数 Dy=diff(y) D2y=diff(y,2) … D(m-1)y=diff(y,m-1) y=