[曾德学]

diff的原代码有人知道吗？
我编了一个，感觉有错误，大家帮忙该改，
意思是，已知y了，求一个x=(1/1*pi*1000)dy/dt+y,各位帮帮忙，谢谢了


t=-0.001:1/5000:0.001;
y=0.3*sin(2000*2*pi*t);
m=0.002*5000;
for i=1:m
n=i+1;
y1(i)=(y(n)-y(i))/(1/5000);
i=i+1;
end
y1(m+1)=0;
x=(1/(2*pi*1000))*y1+y;
plot(x)