第10章数学软件使用入门.pptVIP

* 中山大学计算机科学系 10.3.6 符号计算 三、符号代数运算 包括化简、展开、因式分解、合并同类项等运算 例10.37 试一试符号代数运算功能 syms r s x y a; simplify(((r+s)/s+s/(r-s))/(s/(r-s))) %化简 expand((x+y)^2) %展开 expand(sin(2*x)) %展开 factor(x^2 - 7*x + 12) %因式分解 factor(80) %因式分解 p=x^2 - a*y^2*x^2 + 3*y^2*x + x - x*y; collect(p,x) %按x合并同类项 collect(p,y) %按y合并同类项 * 中山大学计算机科学系 10.3.6 符号计算 四、符号微积分运算 包括求极限、求导、求积分、求级数和等运算。 例10.38 试一试符号微积分运算 syms n; limit((1+2/n)^n,n,Inf) %当n→∞时 (1+2/n)^n 极限 syms x; limit(sin(x)^tan(x),x,pi/2) %当x→π/2时函数极限 limit(abs(x)/x,x,0,‘right’) %当x→0+时函数右极限 limit(abs(x)/x,x,0,‘left’) %当x→0-时函数左极限 diff(sin(x),x) %以x自变量，求函数 sin(x) 的1阶导数 diff(cos(x)^2 * log(x),x) %以x为自变量，求 1阶导数 diff(x^5,x,2) %以x为自变量，求2阶导数 * 中山大学计算机科学系 10.3.6 符号计算 四、符号微积分运算 int(exp(-2*x)*sin(3*x),x) %求不定积分 int(x+sin(x),x,0,pi/2) %求函数在[0,π/2]区间的定积分， int(1/(exp(x+1)-2),x,0,Inf) %求 [0,∞) 区间的定积分 symsum(x/n^2,n,1,10) % 求级数x/n^2从第1项至第10项的级数和 symsum(x/n^2,n,1,Inf) %求无穷级数之和2 symsum(x^n/(3^n*sym(‘n!’)),n,0,Inf) %求无穷级数之和，函数 sym(n!) 生成表示求阶乘的数学表达式 n! * 中山大学计算机科学系 10.3.6 符号计算 五、解符号方程 使用函数 solve 求解符号代数方程。 例10.39 试一试符号代数方程求解方法 syms x a b c; %解方程 a*x^2+b*x+c=0 solve(a*x^2+b*x+c,x) %解方程 sin(6*x)=1 solve(sin(6*x)-1,x) %解方程组 x^2+y^2=m，x-y=n syms m n [X,Y]=solve(x^2+y^2-m,x-y-n,x,y) * 中山大学计算机科学系 10.3.6 符号计算 五、解符号方程 使用函数 dsolve 求解符号常微分方程，格式是： dsolve(‘eq1’, ‘ eq2’, …) ， eq1、ep2、… 是方程或条件。其中，大写字母 D 有着特殊含义，如 Dy 表示 y 关于自变量的一阶导数， D2y 表示 y 关于自变量的二阶导数，依此类推。 例10.40 试一试符号常微分方程求解方法 * 中山大学计算机科学系 10.3.6 符号计算 五、解符号方程 %求微分方程 y‘ = x 的通解 syms x y dsolve(‘Dy=x’,‘x’) %求微分方程 y‘ = a*y 的特解，其中 x 是自变量，y(0) = b 是条件。 syms a b dsolve(‘Dy = a*y’, ‘y(0) = b’,‘x’) % 求微分方程组 x‘ = y, y’=-x 的通解。 [x,y] = dsolve(‘Dx = y’,‘Dy = -x’) p=[5,0,-21,3]; %定义表示第1个多项式的向量 q=[2,7]; %定义表示第2个多项式的向量 k=conv(p,q) %函数 conv 用于计算两个多项式的乘积 * x=input(x=); %提示输入 x y=input(y=); %提示输入 y if xy maxValue=x; else maxValue=y; end disp(sprintf(最大值是:%d,maxValue)) %显示计算结果 %%%%%%%%%%%%%%%%%%%%%%%% sel = input(中国的首都在哪个城市?\n1.香港 2.广州 3.北京 4.上海\n); switch sel