急求“用MATLAB编程数值分析的一道题目（题目见补充）”

来源：百度知道 编辑：UC知道 时间：2024/06/24 09:29:26
已知某次实验测试中测得一个支点在几个特定点的速度y为
t（s） 1.0000 1.5000 2.0000 2.5000 3.0000 3.5000 4.0000
y（m/s） 0.1000 0.3000 0.6000 0.7000 0.9000 1.2000 1.4000

最好今天晚上十点之前做好，如果可以做好，追加100分。。。。。。

明天早上八点前做好追加50分.。。。

过时不候
忘记问问题了，求在这段时间内质点的位移！

绝对能用的，正好我也正学
function fp = newton_interpolation(x,y,p)
% Script for Newton's Interpolation.
% Muhammad Rafiullah Arain
% Mathematics & Basic Sciences Department
% NED University of Engineering & Technology - Karachi
% Pakistan.
% ---------
% x and y are two Row Matrices and p is point of interpolation
%

n = length(x);
a(1) = y(1);
for k = 1 : n - 1
d(k, 1) = (y(k+1) - y(k))/(x(k+1) - x(k));
end
for j = 2 : n - 1
for k = 1 : n - j
d(k, j) = (d(k+1, j - 1) - d(k, j - 1))/(x(k+j) - x(k));
end
end
d
for j = 2 : n
a(j) = d(1, j-1);
end
Df(1) = 1;
c(1) = a(1);
for j = 2 : n
Df(j)=(p - x(j-1)) .* Df(j-1);
c(j) = a(j) .* Df(j);
end
fp=sum(c);

例如：
>> x=[1.0000 1.5000 2.0000 2.5000 3.0000 3.5000 4.0000]

x =

1.0000 1.5000 2.0000 2.5000 3.0000 3.5000