数模培训-matlab3MATLAB符号计算功能课件.pptVIP

例13 求极限 syms x n; limit((1 + x/n)^n,n,inf) limit(1/x,x,0,left) ans = exp(x) ans = -inf 例14 级数求和 s1=symsum(x^2,x,0,10) s2 = symsum(x^k,k,0,inf) s1 = 385 s2 = -1/(x-1) 四.符号代数方程求解 MATLAB符号运算能够解一般的线性方程、非线性方程及一般的代数方程、代数方程组。当方程组不存在符号解时，又无其他自由参数，则给出数值解。 命令格式： solve(f) —— 求一个方程的解 Solve(f1,f2, …fn) —— 求n个方程的解 使用solve时，要求未知数的个数和方程的个数一致。 例15. f = ax2+bx+c 求解 A= solve(a*x^2+b*x+c) —— 对缺省变量x求解 A = [1/2/a*(-b+(b^2-4*a*c)^(1/2))] [1/2/a*(-b-(b^2-4*a*c)^(1/2))] 计算机 格式 一般格式 例16. 分别对符号方程cos(x)=sin(x) tan(2*x)=sin(x)求解 f1=solve(cos(x)=sin(x)), f1 = 1/4*pi solve( f, b ) —— 对指定变量b求解 ans = -(a*x^2+c)/x f2=solve(tan(2*x)=sin(x)) f2 = [ 0] [acos(1/2+1/2*3^(1/2))] [acos(1/2 -1/2*3^(1/2))] numeric(f2) ans = 0 0 + 0.8314i 1.9455 例17 用solve解方程 s=solve(2*x^2=( x-1)), n=double(s), s=vpa(n,10) s = [ 1/4+1/4*i*7^(1/2)] [ 1/4-1/4*i*7^(1/2)] n = 0.2500 + 0.6614i 0.2500 - 0.6614i s = [ .25000+.6614378278*i] [ .25000-.6614378278*i] 例18 解方程组 x+y+z=1 x-y+z=2 2x-y-z=1 f=solve(x+y+z=1,x-y+z=2,2*x-y-z=1) f.x,f.y,f.z f =x: [1x1 sym] y: [1x1 sym] z: [1x1 sym] ans =2/3 ans =-1/2 ans =5/6 [x,y,z]=solve(‘x+y+z=1’,… x-y+z=2,2*x-y-z=1) x = 2/3 y =-1/2 z =5/6 五 符号微分方程求解 —— 用一个函数可以方便地得到微 分方程的符号解 符号微分方程求解指令：dsolve 命令格式：dsolve(f,g) f —— 微分方程，可多至12个微分方程的求解；g为初始条件 默认自变量为 ‘t,可任意指定自变量‘x, u等 微分方程的各阶导数项以大写字母D表示 或 或 或 y的一阶导数—— Dy y的二阶导数—— D2y y的 n 阶导数—— Dny [y1,y2…]=dsolve(x1,x2,…xn) —— 返回微分方程的解。 一阶微分方程 例19 dsolve(Dx=y,Dy=x,x(0)=0,y(0)=1) ans = x(t) = sin(t), y(t) = cos(t) 二阶微分方程 例20 dsolve(D2y=-a^2*y,y(0)=1,Dy(pi/a)=0‘,’x’) ans = cos(a*x) 例21 y=dsolve(D2y+2*Dy+2*y=0,y(0)=1,Dy(0)=0) ans = exp(-x)*cos(x)+exp(-x)*sin(x) 求该方程的解 例22 微分方程组 [f,g]=dsolve(Df=3*f+4*g,... Dg=-4*f+3*g,... f(0)=0,g(0)=1) f