工程应用数学 第2版 教学课件 ppt 作者 万金保 5-6.pptVIP

* * 工程应用数学 §5.6 MATLAB与导数运算 §5.6 MATLAB与导数运算 1、求导运算在MATLAB里由命令函数diff()来 完成，其具体形式为： diff(function,’variable’,n) 参数function为需要进行求导运算的函数， variable为求导运算的独立变量，n为求导的 阶次。命令函数diff()默认求导的阶次为1阶； 如果表达式里有多个符号变量，并且没有在 参数里说明，则按人们习惯的独立变量顺序 确定进行求导的变量。 §5.6 MATLAB与导数运算 2、求N个方程N个未知数的方程组解的命令格 式为： [var1,var2,….,varN] =solve(‘eqn1’,’eqn2’,….,’eqnN’,’var1’,’var2’,… .,’varN’) 其中’eqnN’表示第N个方程，’varN’为第N个变 量. §5.6 MATLAB与导数运算 3、MATLAB也提供了另一个功能强大的画图 函数ezplot，格式为: ezplot(‘fun’,[a,b])， 其中[a,b]可省略，缺省状态下为[ ]. §5.6 MATLAB与导数运算 例5-60 求函数y = x5 + 4sinx – cosx + 7的导数。 解：在命令窗口中输入 syms x y=x^5+4*sin(x)-cos(x)+7; diff(y,x) 所以 y’ = 5x4 + 4cosx + sinx 按回车键，显示结果为： ans = 5*x^4+4*cos(x)+sin(x) §5.6 MATLAB与导数运算 例5-61 求函数 的导数。 解：在命令窗口中输入 clear syms x y = ((x+1)^2/((x-1)*(x+2)))^(1/3); diff(y,x) §5.6 MATLAB与导数运算 按回车键，显示结果为： ans = 1/3/((x+1)^2/(x-1)/(x+2))^(2/3)*(2*(x+1)/(x-1)/(x+2)-(x+1)^2/(x-1)^2/(x+2)-(x+1)^2/(x-1)/(x+2)^2) §5.6 MATLAB与导数运算 上式显示结果较为复杂，我们可以用pretty函数，令结果显示得更为直观 在命令窗口中输入： pretty(ans); 所以把结果进行简化可以写成 §5.6 MATLAB与导数运算 例5-62 求s = e–tcost的二阶导数。 解：在命令窗口中输入 clear syms t ; s=exp(-t)*cos(t); diff(s,t,2) 按回车键，显示结果： ans = 2*exp(-t)*sin(t) 所以 s??=2e–tsint §5.6 MATLAB与导数运算 例5–63 以初速度v0，发射角?发射炮弹， 其运动方程为 求炮弹在任何时刻的运动速度的大小和方向。 解：在MATLAB中输入以下命令： syms a v0 t g ; x=v0*cos(a)*t; §5.6 MATLAB与导数运算 y=v0*sin(a)*t-1/2*g*t^2; vx=diff(x,t); vy=diff(y,t); v=sqrt(vx^2+vy^2); %求炮弹的运动速度 v=simplify(v) %simplify表示对函数v进行化简 tanb=vy/vx ; %求炮弹的方向 pretty(tanb); §5.6 MATLAB与导数运算 按回车键，显示结果为： v =(v0^2-2*v0*sin(a)*g*t+g^2*t^2)^(1/2) v0 sin(a) - g t --------------- v0 cos(a) 所以 * * *