数学实验MATLAB第三章.ppt

*/16 符号计算与微积分 符号表达式创建 微积分符号计算 级数求和与泰勒展式 常微分方程符号解 ? ? ? ? 符号表达式——符号变量、运算符、函数、数字组成 表达式语句: 变量名=符号表达式 符号表达式以符号变量为基础;符号变量需声明后才可以投入使用。 MATLAB的符号计算主要功能是处理符号表达式，即针对常用数学表达式进行代数符号计算。 声明符号变量的语句格式: syms 开始, 空格后跟变量列表, 变量名以英文字母为首 syms x1 x2 y1 y2 z1 z2 syms x ; f = exp(-0.2*x)*sin(0.5*x); ezplot(f,[0,8*pi]) 例3.2 用符号表达式定义 f(x)= e-0.2xsin0.5x 并绘图. Name Size Bytes Class f 1x1 168 sym object x 1x1 126 sym object Whos ? 例3.4 用符号计算验证三角恒等式 syms x1 x2; y1= sin(x1)*cos(x2)-cos(x1)*sin(x2); y2=simple(y1) expand(y2) y2 = sin(x1-x2) ans = sin(x1)*cos(x2)-cos(x1)*sin(x2) 练习:分别用simple和simplify对符号表达式 cos(x)^2-sin(x)^2 进行化简，观察结果 化简符号表达式方法——simple() 或simplify() 微积分符号计算 diff(f) — 对缺省变量求导数 diff(f,v) — 对指定变量 v 求导数 diff(f,v,n) —对指定变量 v 求n阶导数 int(f) — 对f表达式的缺省变量求积分 int(f,v) — 对f表达式的v变量求积分 int(f,v,a,b) — 对f表达式的v变量在[a, b] 区间求定积分 例 计算 f = arctan(x) 的导数并绘图 syms x y=atan(x);y1=diff(y,x), y2=diff(y1,x), figure,ezplot(y) figure, ezplot(y1) figure, ezplot(y2) y1=1/(1+x^2) y2=-2/(1+x^2)^2*x 例3.9 计算不定积分 syms a b x F=int(exp(a*x)*sin(b*x));F=simplify(F) F = exp(a*x)*(-b*cos(b*x)+a*sin(b*x))/(a^2+b^2) f=diff(F);f=simplify(f) f = exp(a*x)*sin(b*x) 并求原函数导数，验证不定积分正确性 例 加油站油库深度12米呈啤酒桶形状,顶部和底部半径为2米,中心半径为3米。纵截面右侧是抛物线 t=linspace(-6,6,30);r=3-t.^2/36; cylinder(r,30),axis off syms y D f=3-y*y/36;V=pi*int(f*f,y,-6,D) V12=subs(V,D,6) V =pi*(1/6480*D^5+216/5-1/18*D^3+9*D) V12 = 271.4336 利用公式 计算储油量 符号表达式中变量替换方法——subs() 例3.11 计算曲线段 f(x)=exp(a x)sin(b x), 绕X轴旋转的旋转曲面体积 数据转化为数值数据: double(A) syms a b x f=exp(a*x)*sin(b*x); f1=subs(f,a,-0.2); f2=subs(f1,b,0.5); V=pi*int(f2*f2,x,0,2*pi) double(V) V =pi*(-125/116*exp(-4/5*pi)+125/116) ans = 3.1111 theta=linspace(2*pi,0,20); r=exp(-0.2*theta).*sin(0.5*theta); figure(1),cylinder(r) axis off 旋转曲面的绘图注记： [X,Y,Z]=cylinder(r); figure(2),mesh(-Z,X,Y) axis off colormap([0 0 0]) 例 卫星信号对地球的覆盖面积 通信卫星的轨道位于地球赤道平面内,卫星运行的角速度与地球自转的角速度相同.卫星距地面