求助,数值分析问题!

xuxiaoshuai

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

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

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

邮箱xuhongyan198461@163.com

37wing

....

我也想要

我们也要用.....

ljack

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

 

[em12]

niumeng

[em06]

灵幽二娃

我手算过。。。可是现在书没在身边。。。。
) h( H% K& S, v7 ^* h看清华的数学实验吧

西雅图3173

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

wgbdarling

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

梦天涯M

手算也太累了吧！这是matlab程序：
" Q2 R5 W  D4 y! n" W: bfunction [f,f0]=Newton(x,y,x0)
/ ]7 O$ |! [, _3 X9 z/ s%求已知数据点的均差形式的牛顿插值多项式
- L+ H. G8 I- w3 N! K+ x* ^6 _9 ~: V%已知数据点的x坐标向量：x
* f2 b9 a9 i3 U( u1 L: }1 o%已知数据点的y坐标向量：y
%插值点的x坐标：x0
%求得的均差形式的牛顿插值多项式:f
7 t, W+ [" h1 P. ~* G3 \%x0处的插值:f0
, q2 l* T" s6 H" d. g9 vsyms t;
; k$ k& u8 }) l! R) `1 w, pif(length(x)==length(y))
    n=length(x);
    c(1:n)=0.0;
+ u5 J# u/ H0 _, B9 }: W- O' e  Belse
3 U6 [. R, ]; \$ b3 @6 t1 [    disp('x,y的维数不相得！');
, r) [3 Q% y: o5 U/ n! v# T    return;
& F! \) I# z' N" e% gend
; Y- e  ]- }6 Z7 S* |/ gf=y(1);
y1=0;
l=1;
. {8 n1 C) n, e/ G1 hfor i=1:n-1
    for j=i+1:n
        y1(j)=(y(j)-y(i))/(x(j)-x(i));
0 r, w  m1 @7 ~4 J/ K2 i    end
. l% N- R' y, Y    c(i)=y1(i+1);
' @: n0 o: i, N5 w! b: `+ I9 {    l=l*(t-x(i));
4 l, t' ?$ o! k' c    f=f+c(i)*l;
    simplify(f);
8 _$ c" |- ?+ I    y=y1;
end
f0=subs(f,'t',x0);
% d* k  Y5 O; C# \