Matlab符号运算LIU.pptVIP

* * syms a b x;f=[a*x,b*x^2;1/x,sin(x)]; disp(The integral of f is);pretty(int(f)) The integral of f is [ 2 3] [1/2 a x 1/3 b x ] [ ] [ log(x) -cos(x) ] * int（被积表达式，积分变量，积分上限， 积分下限）—— 定积分 ——缺省时为不定积分 s1=sym(6*x^3-4*x^2+b*x-5); s2=sym(sin(a)); int(s1) ans = 3/2*x^4-4/3*x^3+1/2*b*x^2-5*x int(s2) ans = -cos(a) * 例.计算二重不定积分 F=int(int(x*exp(-x*y),x),y) F= 1/y*exp(-x*y) syms x y F=int(int(x*exp(-x*y),x),y) F = 1/y*exp(-x*y) * 4.符号代数方程求解 matlab符号运算能够解一般的线性方程、非线性方程及一般的代数方程、代数方程组。当方程组不存在符号解时，又无其他自由参数，则给出数值解。 命令格式： solve(f) —— 求一个方程的解 Solve(f1,f2, …fn) —— 求n个方程的解 * 例1. f = ax2+bx+c 求解 f=sym(a*x^2+b*x+c‘); solve(f) —— 对缺省变量x求解 ans = [1/2/a*(-b+(b^2-4*a*c)^(1/2))] [1/2/a*(-b-(b^2-4*a*c)^(1/2))] 计算机 格式 一般格式 * 例. 符号方程cos(x)=sin(x) 求解 f1=solve(cos(x)=sin(x)), f1 = 1/4*pi solve(f , b ) —— 对指定变量b求解 ans = -(a*x^2+c)/x * 例. 解方程组 x+y+z=1 x-y+z=2 2x-y-z=1 g1=sym(x+y+z=1‘),g2=sym(x-y+z=2‘),g3=sym(2*x-y-z=1‘) [a,b,c]=solve(g1,g2,g3) f = a = 2/3 ,b = -1/2, c = 5/6, * 【例5.2.3.2-1】用简单算例演示subs的置换规则。 syms a x;f=a*sin(x)+5; f1=subs(f,sin(x),sym(y)) f1 = a*y+5 * f3=subs(f,{a,x},{2,pi/3}) f3 = 6.7321 * 求射程 syms v0 t theta g distancey=v0*t*sin(theta)-1/2*g*t^2; distancex=v0*t*cos(theta); impact_time=solve(distancey,t) impact_time = 0 2*v0*sin(theta)/g impact_distance=subs(distancex,t,impact_time(2)) impact_distance = 2*v0^2*sin(theta)/g*cos(theta) * 求0~pi/2的射程 impact_100=subs(impact_distance,{v0,g},{100,9.8}) impact_100 = 100000/49*sin(theta)*cos(theta) * ezplot(impact_100,[0,pi/2]) * 求射程最大时的发射角度和最大射程 max_angle=solve(diff(impact_100)) max_angle = -1/4*pi 1/4*pi max_distance=subs(impact_100,theta,max_angle(2)) max_distance = 50000/49 double(max_distance) ans = 1.0204e+003 * 【例5