求助,数值分析问题!

xuxiaoshuai

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

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

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

邮箱xuhongyan198461@163.com

37wing

....

我也想要

我们也要用.....

ljack

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

 

[em12]

niumeng

[em06]

灵幽二娃

我手算过。。。可是现在书没在身边。。。。
看清华的数学实验吧

西雅图3173

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

wgbdarling

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

梦天涯M

手算也太累了吧！这是matlab程序：
6 P5 h0 D( u3 C( g- c/ xfunction [f,f0]=Newton(x,y,x0)
%求已知数据点的均差形式的牛顿插值多项式
%已知数据点的x坐标向量：x
+ m: _1 f% \1 y; p5 V%已知数据点的y坐标向量：y
( l. y  C5 R6 r  W$ b. f%插值点的x坐标：x0
9 I$ L: o0 x" z) n; N- \5 Q%求得的均差形式的牛顿插值多项式:f
# f* q; a! ?) |6 e$ U' u%x0处的插值:f0
syms t;
+ X+ Y7 \( K) ?1 M  K1 {6 B$ {if(length(x)==length(y))
" R/ x' X2 W9 `8 b    n=length(x);
    c(1:n)=0.0;
; X% c. w6 v8 s* D( a) Felse
, ?4 r5 p  [' i. R    disp('x,y的维数不相得！');
$ `7 R$ y0 J" P4 X    return;
end
f=y(1);
7 }8 }: e2 F( g& C% Q& f1 V1 ly1=0;
l=1;
3 y* [( I5 m1 A0 ]0 q, wfor i=1:n-1
    for j=i+1:n
5 D6 k1 \" n, l- b$ M, R8 Y        y1(j)=(y(j)-y(i))/(x(j)-x(i));
    end
# L6 I1 @- k1 o, A' Z' V. g' V, l    c(i)=y1(i+1);
( x* [- I* ]' v+ e* Q4 U1 [, Z    l=l*(t-x(i));
2 B( o' E' e8 M  f$ O- k" t) |. W( O2 i    f=f+c(i)*l;
    simplify(f);
    y=y1;
end
* \% j+ C) ?- ff0=subs(f,'t',x0);