Matlab求导Matlab求导.doc

怎么来对已经求出来的速度曲线进行求导？ 我现在通过状态空间法，求解出来一个结构的速度和位移曲线，现在想通过对速度求导的方式来得到加速度，请各位好手指点，具体的表达式是怎么样的，这是个动力系统，我希望能通过求导来画出加速度的曲线 中心差分法:0.5*(diff(y(1:end-1)+diff(2:end))/dt这样比原始数据的头尾少两个点,可以自己补一下 以试试根据数据用polyfit求出多项式拟合，然后再diff得到加速度表达式，再作图 本帖最近评分记录 lxq +18 2007-6-22 11:28 理由：启发引导性回答 我用matlab求导后，用plot做曲线做不出来，总提示错误，不知道为什么，谁知道怎么出图啊？具体步骤说一下 你直接对y数据求导得到的数据长度比原来少了一个 plot肯定报错，x和y尺度不符 对符号函数求导之后得到符号函数 要么ezplot，要么代入数据再画图matlab的求导命令与求导法 技术专区 ??2009-03-10 16:19 ??阅读543???评论1 ? 字号： 大大? 中中? 小小 １．matlab命令. 建立符号变量命令sym和syms调用格式: x=sym(x)， 建立符号变量x； syms x y z ， 建立多个符号变量x，y，z； matlab求导命令diff调用格式: diff(函数) ， 求的一阶导数； diff(函数， n) ， 求的n阶导数(n是具体整数)； diff(函数，变量名)， 求对的偏导数； diff(函数， 变量名，n) ，求对的n阶偏导数； matlab求雅可比矩阵命令jacobian，调用格式: jacobian([函数；函数； 函数]， [])给出矩阵: ２．导数概念. 导数是函数的变化率，几何意义是曲线在一点处的切线斜率. （１）点导数是一个极限值. 例1.设，用定义计算. 解:在某一点的导数定义为极限: 我们记，输入命令: syms h；limit((exp(0+h)-exp(0))/h，h，0) 得结果:ans=１．可知 （２）导数的几何意义是曲线的切线斜率. 例2.画出在处()的切线及若干条割线，观察割线的变化趋势. 解:在曲线上另取一点，则的方程是: .即 取，分别作出几条割线. h=[3，2，1，0.1，0.01]；a=(exp(h)-1)./h；x=-1:0.1:3； plot(x，exp(x)，’r’)；hold on for i=1:5； plot(h(i)，exp(h(i))，’r.’) plot(x，a(i)*x+1) end axis square 作出在处的切线 　　plot(x，x+1，’r’) 从图上看，随着与越来越接近，割线越来越接近曲线的割线． ３．求一元函数的导数. （１）的一阶导数. 例3.求的导数. 解:打开matlab指令窗，输入指令: dy_dx=diff(sin(x)/x). 得结果: dy_dx=cos(x)/x-sin(x)/x^2. matlab的函数名允许使用字母、空格、下划线及数字，不允许使用其他字符，在这里我们用dy_dx表示 例4.求的导数. 解: 输入命令: dy_dx=diff(log(sin(x))). 得结果: dy_dx=cos(x)/sin(x). 在matlab中，函数用log(x)表示，而log10(x)表示 例5.求的导数. 解: 输入命令:dy_dx=diff((x^2+2*x)^20). 得结果: dy_dx=20*(x^2+2*x)^19*(2*x+2). 注意输入时应为2*x. 例6.求的导数. 解: 输入命令: dy_dx=diff(x^x). 得结果: dy_dx =x^x*(log(x)+1). 利用matlab 命令diff一次可以求出若干个函数的导数. 例7.求下列函数的导数: 1. 2. 3. 4. 解: 输入命令: a=diff([sqrt(x^2- 2*x+5)，cos(x^2)+2*cos(2*x)，4^(sin(x))， log(log(x))]). 得结果: a= [1/2/(x^2-2*x+5)^(1/2)*(2*x-2)，-2*sin(x^2)*x-4*sin(2*x)， 4^sin(x)*cos(x)*log(4)， 1/x/log(x)]. dy1_dx=a(1) dy1_dx=1/2/(x^2-2*x+5)^(1/2)*(2*x-2). dy2_dx=a(2) dy2_dx=-2*sin(x^2)*x-4*sin(2*x). dy3_dx=a(3) dy3