- 积分
- 23
- 注册时间
- 2010-11-7
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2010-12-19 09:22:57
|
显示全部楼层
来自 山东烟台
本帖最后由 liuyalong008 于 2010-12-19 09:39 编辑
以下是我自己写的关于离散序列求导的函数
x不一定是等间距的,k是求导次数
- function [f,x_out]=cal_diff(x,y,k)
- len=length(x);
- dy=diff(y,k);
- for i=1:k
- tempx=zeros(len-k,1);
- dis=tempx;
- for j=1:(len-i)
- tempx(j)=(x(j)+x(j+1))/2;
- dis(j)=x(j+1)-x(j);
- end
- x=tempx;
- end
- L=size(dy);
- if L(1)==1
- dy=dy';
- end
- f=dy./dis;
- x_out=x;
复制代码-
- x=0:0.1:2*pi;
- y=sin(x);
- [a,b]= cal_diff(x,y,1);
- [c,d]= cal_diff(x,y,2);
- plot(x,y,b,a,d,c)
- legend('origin','1st derivative','2nd derivative')
复制代码 |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
评分
-
1
查看全部评分
-
|