第三章差分微分matlab.pptVIP

微分方程与差分方程上机; MATLAB的符号运算;符号表达式的建立;创建符号变量和表达式;2.使用syms命令创建符号变量和符号表达式 ——syms用于创建多个符号变量 语法： syms(‘arg1’, ‘arg2’, …,参数)  %把字符变量定义为符号变量 syms arg1 arg2 …,参数  %把字符变量定义为符号变量的简洁形式 ;符号矩阵;符号表达式的代数运算;符号表达式的代数运算;2. 函数运算 (1) 三角函数和双曲函数 三角函数包括sin、cos、tan；双曲函数包括sinh、cosh、tanh；三角反函数除了atan2函数仅能用于数值计算外，其余的asin、acos、atan函数在符号运算中与数值计算的使用方法相同。 (2) 指数和对数函数 指数函数sqrt、exp的使用方法与数值计算的完全相同；对数函数在符号计算中只有自然对数log(表示ln)，而没有数值计算中的log2和log10。 (3) 复数函数 复数的共轭conj、求实部real、求虚部imag和求模abs函数与数值计算中的使用方法相同。但注意，在符号计算中，MATLAB没有提供求相角的命令。 (4) 矩阵代数命令 MATLAB提供的常用矩阵代数命令有diag，triu，tril，inv，det，rank， poly，eig 、expm等，它们的用法几乎与数值计算中的情况完全一样。;【例】求矩阵的行列式值、非共轭转置和特征值。;【例】符号表达式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中所有的符号变量。;2、符号表达式的化简;3、符号表达式的替换;;5、 符号表达式与多项式的转换;符号极限、微积分和级数求和;符号微分;符号积分;符号级数;【例】求级数  1+x+x2+…+xk+…的和。;【例】求ex的泰勒展开式;符号方程的求解; 符号常微分方程;一、利用Matlab求微分方程的解析解; 解 输入命令: y=dsolve(D2y+4*Dy+29*y=0,y(0)=0,Dy(0)=15,x);解 输入命令 ： [x,y,z]=dsolve(Dx=2*x-3*y+3*z,Dy=4*x-5*y+3*z,Dz=4*x-4*y+2*z, t)； x=simple(x) % 将x化简 y=simple(y) z=simple(z);建立数值解法的一些途径;2、使用数值积分;3、使用泰勒公式;;不同求解器Solver的特点;;不同求解器Solver的特点;参数设置;参数设置;参数设置;参数设置;参数设置;参数设置;创建函数function2, 保存在function2.m中function f=function2(t,x) f=-x.^2;;;创建函数function3,保存在function3.m中： function f=function3(t,x) f=[x(1)-0.1*x(1)*x(2)+0.01*t;... -x(2)+0.02*x(1)*x(2)+0.04*t];;;;创建函数function4，存在function4.m中 function f = function4(t,x) global U f = [x(2);U*(1-x(1)^2)*x(2)-x(1)];;