matlab对离散数据求导

loveblue

试验数据，一组时间，一组位移，可以直接绘时间-位移曲线，如果想绘时间速度曲线下面的程序问题在哪儿呢？
time=a(:,1);
disp=a(:,2);
rate=diff(disp,time,1);
plot(rate);

qibbxxt

请lz仔细看看帮助文件，不要按照自己的想法去写程序，应该遵守matlab的语法规则

liuyalong008

以下是我自己写的关于离散序列求导的函数
x不一定是等间距的，k是求导次数

1. function [f,x_out]=cal_diff(x,y,k)
   
2. len=length(x);
   
3. dy=diff(y,k);
   
4. for i=1:k
   
5. tempx=zeros(len-k,1);
   
6. dis=tempx;
   
7. for j=1:(len-i)
   
8. tempx(j)=(x(j)+x(j+1))/2;
   
9. dis(j)=x(j+1)-x(j);
   
10. end
    
11. x=tempx;
    
12. end
    
13. L=size(dy);
    
14. if L(1)==1
    
15. dy=dy';
    
16. end
    
17. f=dy./dis;
    
18. x_out=x;

复制代码

1. 
   
2. x=0:0.1:2*pi;
   
3. y=sin(x);
   
4. [a,b]= cal_diff(x,y,1);
   
5. [c,d]= cal_diff(x,y,2);
   
6. plot(x,y,b,a,d,c)
   
7. legend('origin','1st derivative','2nd derivative')
   

复制代码

gfl

也可以先写成样条函数spline(t, u)，再用样条函数求导的命令求导，这样做应该比直接用中心差分/差分求导效果好。

loveblue

谢谢两位的详细指导，两种方法都试一下。谢谢！

liushuangq05

同意gfl 的方法,对于离散数据先进行曲面拟合,在对拟合曲面求导作为原离散数据点的导数,这是通用的方法,也是许多专家推荐的方法.