第4单元 MATLAB符号计算.pptVIP

4.6 方程的解析解 求解代数方程 函数solve 其调用格式为： solve (eqn1,eqn2,…,eqnN,var1,var2,…,varN) 例4.21 [x,y] = solve(x^2 + x*y + y = 3,x^2 - 4*x + 3 = 0) 常微分方程的求解 函数dsolve dsolve(eqn1,condition,var) 该函数求解微分方程eqn1在初值条件condition下的特解。参数var描述方程中的自变量符号，省略时按缺省原则处理，若没有给出初值条件condition，则求方程的通解。 dsolve在求微分方程组时的调用格式为： dsolve(eqn1,eqn2,…,eqnN,condition1,…,conditionN,var1,…,varN) 函数求解微分方程组eqn1、…、eqnN在初值条件conditoion1、…、conditionN下的解，若不给出初值条件，则求方程组的通解 Dsolve 参数说明 (1)微分方程输入方法： 例： 命令：Y=dsolve(‘D2y-2*Dy-3*y=0’,’x’) (2)condition 初始条件 y(a)=b; Dy(a)=c; D2y(a)=d 等 例4.22 求微分方程的通解。 命令如下： y=dsolve(‘Dy-(x^2+y^2)/x^2/2’, x) %方程的右端为0时可以不写 y=dsolve(‘Dy*x^2+2*x*y-exp(x)’, x) y=dsolve(‘Dy-x/y/sqrt(1-x^2)’, ‘x) 例4.23 求微分方程的特解。 命令如下： y=dsolve(‘Dy=2*x*y^2’, ‘y(0)=1’, x) y=dsolve(‘Dy-x^2/(1+y^2)’, ‘y(2)=1’, x) 第4章 MATLAB符号运算 * 本章内容 4.1 符号对象及其表达方式 4.2 符号算术运算 4.3 独立变量与表达式化简 4.4 符号微积分运算 4.5 符号积分变换 4.6 方程的解析解 4.1 符号对象及其表达式化简 1. 建立符号变量和符号常量 (1)sym函数 S=sym(A) sym函数用来建立单个符号变量，例如，a=sym(‘a’)建立符号变量a，此后，用户可以在表达式中使用变量a进行各种运算。 定义符号常量 pi2=sym(pi) a=sym(‘log(2)’) 例4.1考察符号变量和数值变量的差别。 在 MATLAB命令窗口，输入命令： a=sym(a);b=sym(b);c=sym(c);d=sym(d); %定义4个符号变量 A=[a,b;c,d] %建立符号矩阵A det(A) %计算符号矩阵A的行列式 --------------------------------------------- w=10;x=5;y=-8;z=11; %定义4个数值变量 B=[w,x;y,z] %建立数值矩阵B det(B) %计算数值矩阵B的行列式 例4.2 比较符号常量与数值变量在代数运算时的差别。 在 MATLAB命令窗口，输入命令： pi1=sym(pi);k1=sym(8);k2=sym(2);k3=sym(3); % 定义符号常数 pi2=pi;r1=8;r2=2;r3=3; % 定义数值变量 sin(pi1/3) % 计算符号表达式值 sin(pi2/3) % 计算数值表达式值 sqrt(k1) % 计算符号表达式值 sqrt(r1) % 计算数值表达式值 sqrt(k3+sqrt(k2)) % 计算符号表达式值 sqrt(r3+sqrt(r2)) % 计算数值表达式值 (2)syms函数 1.用于定义多个符号变量。 syms函数的一般调用格式为： syms var1 var2 … varn 函数定义符号变量var1,var2,…,varn等。用这种格式定义符号变量时不要在变量名上加字符分界符(‘)，变量间用空格而不要用逗号分隔。 2. 建立符号表达式及符号矩阵 例4.3: 用2种方法建立符号表达式。 在MATLAB窗口，输入命令： U=sym(3*x^2+5*y+2*x*y+6) %定义符号表达式U syms x y; %建立符号变量x、y V=3*x^2+5*y+2*x*y+6 %定