MATLAB求参数方程2阶偏导数.ppt

制作人：倪亮军MATLAB求参数方程2阶偏导数时的问题已知MATLAB参数方程偏导数求解方法设参数方程为x=g(t)和y=h(t)所确定的函数y=f(x),则y=f(x)的导数为计算参数方程k阶导数的函数调用格式为：这样k阶导数求解是否正确?构建偏导数求解函数functionvarargout=text_diff2(x,y,t)T=pi/4:pi/50:3*pi/4;n=length(T);X=zeros(1,n);Y1=zeros(1,n);Y2=zeros(1,n);%求1阶导数dydx1=diff(y,t)/diff(x,t);disp('参数方程的1阶导数的简化符号表达式：');s_dydx1=simplify(dydx1)dydx2=diff(y,t,2)/diff(x,t,2);%算法可取？disp('参数方程的2阶导数的简化符号表达式：');s_dydx2=simplify(dydx2)%按照高等数学中的算法重新求解该参数方程的3阶导数dydx_2=diff(s_dydx1,t)/diff(x,t);s_dydx_2=simplify(dydx_2)reply=input('是否需要图形显示？Y/N[Y]:','s')ifisempty(reply)||reply=='Y'fork=1:nX(k)=vpa(subs(x,t,T(k)),6);%subs函数将函数dydx1中的t替换%vpa函数将表达式运算结果近似到第6位有效数字Y1(k)=vpa(subs(s_dydx2,t,T(k)),6);Y2(k)=vpa(subs(s_dydx_2,t,T(k)),6);endplot(X,Y1,X,Y2);xlabel('\pi/4\leqx\leq3\pi/4');legend('初始2阶导','高数2阶导');title('两种2阶导数求解的直观图');elseifreply=='N'disp('不需要图形显示');elsedisp('输入有误！')endend参数方程输入例1参数方程例2参数方程参数方程MATLAB输入例1例2symst>>x=cos(t);>>y=sin(t);symst>>x=log(cos(t));>>y=cos(t)-sin(t);text_diff2(x,y,t)参数方程的1阶导数的简化符号表达式：s_dydx1=-cos(t)/sin(t)参数方程的2阶导数的简化符号表达式：s_dydx2=sin(t)/cos(t)s_dydx_2=-1/sin(t)^3参数方程的1阶导数的简化符号表达式：s_dydx1=(cos(t)*(cos(t)+sin(t)))/sin(t)参数方程的2阶导数的简化符号表达式：s_dydx2=cos(t)^2*(cos(t)-sin(t))s_dydx_2=(cos(t)*(2*cos(t)*sin(t)-2*cos(t)^2+cos(t)^4-cos(t)^3*sin(t)+1))/(cos(t)^2-1)^2绘制例1二阶导数图像例2二阶导数图形结论使用求解k>1的参数方程偏导数并不可取。在进行求解时，需注意实际操作与理论知识的衔接。谢谢！