求助,数值分析问题!

xuxiaoshuai

数值分析老师留了一个大作业，让结合一些工程数据作一个数学模型，用数值分析的牛顿或者插值或者线形方程组解出来，之后分析误差！

虽然我觉得自己学的不错，可是我知道这个作业不是要求做的特别复杂，就是因为这样才不知道怎么下手！老师说不需要用软件，只用手算就可以了！

哪位有做过这个作业的帮助一下了！谢谢了！或者有这个方面的例子发给我看一下！只要知道怎么下手就行！呵呵

邮箱xuhongyan198461@163.com

37wing

....

我也想要

我们也要用.....

ljack

不会啊!应该在那个网站可以看到这方面的东西呢？

 

[em12]

niumeng

[em06]

灵幽二娃

我手算过。。。可是现在书没在身边。。。。
' d4 N+ y, h2 N) ~* h' ~) B看清华的数学实验吧

西雅图3173

我们也用啊，是不是一个学校的啊

wgbdarling

数值分析，无非就是迭代，拟合，最小二乘

梦天涯M

手算也太累了吧！这是matlab程序：
2 g" N) T! f# E, ~+ q5 t! Mfunction [f,f0]=Newton(x,y,x0)
& v) W2 z* V4 Z- `( C( i5 ~%求已知数据点的均差形式的牛顿插值多项式
$ e3 ~1 @5 k! T: k& T4 d%已知数据点的x坐标向量：x
%已知数据点的y坐标向量：y
& ~" q+ T4 {$ p8 y3 v+ m( ?%插值点的x坐标：x0
. X5 u: T1 F- h' T8 i. z0 Y! `  d%求得的均差形式的牛顿插值多项式:f
2 s/ X- p' D. W%x0处的插值:f0
: R* ?, v, _* ^5 D8 vsyms t;
if(length(x)==length(y))
# {; Q2 i' b9 O6 q4 k7 D% [2 y    n=length(x);
; p& g. Y7 S5 l9 w1 t- d    c(1:n)=0.0;
$ v4 a6 @7 i8 p: x( belse
% T% t$ t" I; c8 B1 ]; B+ M" Y. h    disp('x,y的维数不相得！');
    return;
end
f=y(1);
y1=0;
  J8 @6 ]. _) |4 Al=1;
for i=1:n-1
+ c4 h/ M  ?- i& s% h    for j=i+1:n
7 h1 F& v$ [3 A8 Q  z* `, t  ~        y1(j)=(y(j)-y(i))/(x(j)-x(i));
* B) s" p. b3 k2 j! Y% I/ D# _& I6 |    end
; _) A4 r+ a5 w7 \% @! y3 e    c(i)=y1(i+1);
$ u: N, u/ B9 F7 P4 }    l=l*(t-x(i));
    f=f+c(i)*l;
% d! b3 f* F; J7 w( X, Y    simplify(f);
& y( \+ M/ ?, I# N: S( x    y=y1;
; f( v1 M& \5 d5 _. Q; Eend
f0=subs(f,'t',x0);