数学实验第三章.pptVIP

参 考 程 序 r=1:0.1:3; t=0:0.06:2*pi; x=r*cos(t); z=r*sin(t); y=sqrt(1+(1e-008)-(sqrt(x.^2+z.^2)-2).^2); mesh(x,y,z) hold on mesh(x,-y,z) colormap([0 0 0]) xlabel(x),ylabel(y),zlabel(z) axis equal 参 考 程 序 t=(0:60)*pi/30; theta=(0:40)*pi/20; r=2+cos(theta); z=r*sin(t); y=sin(theta)*ones(size(t)); x=r*cos(t); mesh(x,y,z) colormap([0 0 0]) xlabel(x),ylabel(y),zlabel(z) axis equal y=-1:0.1:1 r=sqrt(1-y.^2)+2 t=0:0.06:2*pi; x=r*cos(t); z=r*sin(t); mesh(x,y,z) colormap([0 0 0]) xlabel(x),ylabel(y),zlabel(z) axis equal */20 符号计算与微积分 符号表达式及其应用 微分和积分符号计算 台劳展开式符号计算 旋转曲面绘制方法 ? ? ? ? 在符号计算中, 符号表达式是主要操作对象. 符号表达式——符号变量、运算符、函数、数字组成 在定义符号表达式之前, 首先要创建符号变量. 符号计算又称为计算机代数，以符号形式处理数学表达式，关注准确的计算和公式推导。符号计算不仅用于数学研究，还可以用于工程计算 。1993年MathWorks公司购买Maple的使用权, 开发了符号计算工具箱——Symbolic Math Toolbox. syms 符号变量1 符号变量2 ······ 符号变量创建方法 syms x ; f = exp(-0.2*x)*sin(0.5*x); ezplot(f, [0,8*pi]) 例3.1 用符号表达式定义 f(x)= e-0.2xsin0.5x 并绘图. Name Size Bytes Class f 1x1 168 sym object x 1x1 126 sym object Whos ? S1=subs(S, ‘old’, ‘new’) 用‘new’置换符号表达式S中的’old’ 例3. 2 输入不同的参数a, b绘制函数图形 f(x)=exp(-a x)sin(b x) 符号表达式中变量替换方法 syms a b x f=exp(-a*x)*sin(b*x); f1=subs(f,a,0.8); f2=subs(f1,b,0.5); ezplot(f2,[0,2*pi]) 例3.3 用符号计算验证三角恒等式 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 进行化简，观察结果 化简符号表达式方法——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] 区间求定积分 例3.4 验证求导数公式 (k=1，2，3，4 ) syms x s=[diff(sin(x),1)==sin(x+pi/2); diff(sin(x),2)==sin(x+pi); diff(sin(x),3)==sin(x+3*pi/2); diff(sin(x),4)==sin(x+2*pi)] s = 1 1 1 1 Name Size Bytes Class s 1x4 32 double array x