MATLAB的符号运算实习V2.0.docVIP

MATLAB语言基础实习指导 MATLAB的符号运算 目的：通过实习MATLAB符号运算，掌握MATLAB矩阵的符号运算。 要求：熟悉符号表达式，符号矩阵的创建，符号线性代数，因式分解、展开和简化，符号代数方程求解，符号微积分，符号微分方程等。 1、sym和syms建立符号对象，运用findsym()函数查找符号表达式中的符号变量 syms a b x X Y k=sym(3) z=sym(c*sqrt(delta)+y*sin(theta)) f=a*z*X+(b*x^2+k)*Y 试写出下面命令的结果，上机验证，并思考原因： findsym(f) findsym(f,1) findsym(f,2) findsym(f,3) findsym(f,7) findsym(f,8) 2、符号表达式的两种建立方法 方法一： syms x f1=sin(x)+cos(x) 方法二： f2=sym(’sin(x)+cos(x)’) 3、subs函数进行符号替换 syms x a k=sym(‘3’) f=k*x+a subs(f,x,a) syms x y f=2*x+y; x=3,y=4; subs(f) syms x y subs(f,x,’a’) syms x y a b f=2*x+y subs(f,[x,y],[3,4]) subs(f,{x,y},{3,4}) subs(f,x,[1:3]) subs(f,{x,y},{[1:3],[5:7]}) subs(f,{x,y},{a+b,a-b}) subs(f,{x,y},{x+y,x-y}) 4、符号矩阵的建立 (1)使用sym函数直接生成 A=sym(’[1+x, sin(x); 5, exp(x)]’) D= [1+x, sin(x); 5, exp(x)] (2)将数值矩阵转化成符号矩阵 B=[2/3, sqrt(2); 5.2, log(3)] C=sym(B) (3)符号矩阵中元素的引用和修改 A=sym(’[1+x, sin(x); 5, exp(x)]’) A(1,2) A(2,2)=sym(’cos(x)’) 5、设A是一个符号矩阵（定义如下），试指出findsym(A,1) 的输出结果，并由此得出结论： findsym 确定自由变量时，是对整个矩阵进行的，还是对各个矩阵元素分别进行？ syms a b t u v x y A=[a+b*x,sin(t)+u; x*exp(-t),log(y)+v] findsym(A,1) 6、因式分解函数factor syms x f=x^6 +1 s=factor(f) 7、设a、b定义如下，试上机输出factor(a)和factor(b)的结果，并指出哪个结果才是12345678901234567890的因式分解，为什么？ a=sym(‘12345678901234567890’) b=sym(12345678901234567890) factor(a) factor(b) 8、展开函数expand syms x y expand( (x+1)^3 ) expand( sin(x+y) ) 9、合并同类项函数collect syms x y collect (x^2*y+y*x-x^2-2*x ) collect (x^2*y+y*x-x^2-2*x,y ) collect (x*exp(-2*x)+3/16*exp(-2*x), exp(-2*x) ) 10、简化函数simple和simplify syms x f= x^3+3*x^2+3*x+1 [R,HOW]=simple(f) simple(f) simplify(sin(x)^2+cos(x)^2) 11、分式通分函数numden syms x y [n d]=numden(x/y+y/x) [n d]=numden(sym(112/1024)) [n d]=numden( sym(0.24) ) 12、嵌套形式的多项式函数horner syms x f=x^6+x^5+x^4+x^3+x^2+x+1 horner(f) 13、计算极限函数limit syms h n x N=limit((log(x+h)-log(x))/h,h,0) M=limit((1-x/n)^n,n,inf) L=limit(abs(x)/x,x,0,’left) R=limit(abs(x)/x,x,0,’right) 14、计算导数函数diff syms a x y=sin(a*x) A=diff(y,x) B=diff(y,a)