Question

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

已知某次实验测试中测得一个支点在几个特定点的速度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分.。。。

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

Answer 1

绝对能用的，正好我也正学
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 4.0000

>> y=[0.1000 0.3000 0.6000 0.7000 0.9000 1.2000 1.4000]

y =

0.1000 0.3000 0.6000 0.7000 0.9000 1.2000 1.4000

>> newton_interpolation(x,y,3)

d =

0.4000 0.2000 -0.4000 0.4000 -0.2400 0.0889
0.6000 -0.4000 0.4000 -0.2000 0.0267 0
0.2000 0.2000 -0.0000 -0.1333 0 0
0.4000 0.2000 -0.2667 0 0 0
0.6000 -0.2000 0 0 0 0
0.4000 0 0 0 0 0

ans =

0.9000

Answer 2

0.4000 0.2000 -0.4000 0.4000 -0.2400 0.0889
0.6000 -0.4000 0.4000 -0.2000 0.0267 0
0.2000 0.2000 -0.0000 -0.1333 0 0
0.4000 0.2000 -0.2667 0 0 0
0.6000 -0.2000 0 0 0 0
0.4000 0 0 0 0 0

Answer 3

呵呵 就是要你分析问题嘛！

Answer 4

问什么都不说？？！