第3节差分微分matlab.pptVIP

微分方程与差分方程上机 合肥师范学院 MATLAB的符号运算 符号表达式的建立 创建符号变量和表达式 符号矩阵 符号表达式的代数运算 符号表达式的代数运算 【例】求矩阵的行列式值、非共轭转置和特征值。 【例】符号表达式f=2x2+3x+4与g=5x+6的代数运算。 符号表达式的操作和转换 1. 自由变量的确定原则 小写字母i和j不能作为自由变量。 符号表达式中如果有多个字符变量，则按照以下顺序选择自由变量：首先选择x作为自由变量；如果没有x，则选择在字母顺序中最接近x的字符变量；如果与x相同距离，则在x后面的优先。 大写字母比所有的小写字母都靠后。 2. findsym函数 如果不确定符号表达式中的自由符号变量，可以用findsym函数来自动确定。 语法： findsym(f,n) %确定自由符号变量 说明：f可以是符号表达式或符号矩阵；n为按顺序得出符号变量的个数，当n省略时，则不按顺序得出f中所有的符号变量。 3、符号表达式的替换 MATLAB 中，可以通过符号替换使表达式的形式简化。符号工具箱中提供了两个函数用于表达式的替换：1．subexpr 该函数自动将表达式中重复出现的比较长的子表达式或字符串用变量替换，该函数的调用格式为： subexpr(s,s1)，指定用符号变量 s1 来代替符号表达式s（可以是矩阵）中重复出现的字符串。替换后的结果由 ans 返回，被替换的字符串由 s1返回； [Y,s1] = subexpr(X,‘s1’)，该命令与上面的命令不同之处在于第二个参数为字符串，该命令用来替换表达式中重复出现的字符串。 2. subs 函数 subs 可以用指定符号替换表达式中的某一特定符号。 subs(s) subs(s, new) subs(s,old,new) 5、 符号表达式与多项式的转换 符号极限、微积分和级数求和 符号微分 符号积分 符号级数 【例】求级数  1+x+x2+…+xk+…的和。 syms x k s1=symsum(1/k^2,1,10) %计算级数的前10项和 s1 = 1968329/1270080 s2=symsum(1/k^2,1,inf) %计算级数和 s2 = 1/6*pi^2 s3=symsum(x^k,k,0,inf) %计算对k为自变量的级数和 s3 = -1/(x-1) 【例】求ex的泰勒展开式 syms x s1=taylor(exp(x),8) %展开前8项 s1 = 1+x+1/2*x^2+1/6*x^3+1/24*x^4+1/120*x^5+1/720*x^6+1/5040*x^7 s2=taylor(exp(x)) %默认展开前5项 s2 = 1+x+1/2*x^2+1/6*x^3+1/24*x^4+1/120*x^5 符号方程的求解  符号常微分方程 不同求解器Solver的特点 不同求解器Solver的特点 参数设置 参数设置 参数设置 参数设置 参数设置 参数设置 运行命令文件runf3.m [t,x]=ode45(function3,[0,20],[30;20]); plot(t,x); xlabel(time t0=0,tt=20); ylabel(x values x1(0)=30,x2(0)=20); 创建函数function4，存在function4.m中 function f = function4(t,x) global U f = [x(2);U*(1-x(1)^2)*x(2)-x(1)]; 使用于精度较低的情形 一步算法；2，3阶Runge-Kutta方程；累计截断误差达(△x)3 非 刚 性 ode23 大部分场合的首选算法 一步算法；4，5阶Runge-Kutta方程；累计截断误差达(△x)3 非 刚 性 ode45 说明 特点 ODE类型 求解器 Solver 若ode45失效时，可尝试使用 多步法；Gear’s反向数值微分；精度中等 刚性 ode15s 适度刚性情形 采用梯形算法 适度刚性 ode23t 计算时间比ode45短 多步法；Adams算法；高低精度均可到10-3～10-6 非刚性 ode113 说明 特点 ODE类型 求解器Solver 当精度较低时，计算时间比ode15s短 梯形算法；低精度 刚性 ode23tb 当精度较低时，计算时间比ode15s短 一步法；2阶Rosebrock算法；低精度 刚性 ode23s 说明 特点 ODE类型 求解器Solver 相对误差对应于解向量中的所有元素。在每步(第k步)