求助,数值分析问题!

xuxiaoshuai

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

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

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

邮箱xuhongyan198461@163.com

37wing

....

我也想要

我们也要用.....

ljack

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

 

[em12]

niumeng

[em06]

灵幽二娃

我手算过。。。可是现在书没在身边。。。。
7 L' ^) `4 F4 Q# u( z0 ^看清华的数学实验吧

西雅图3173

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

wgbdarling

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

梦天涯M

手算也太累了吧！这是matlab程序：
: _2 ]3 I- j5 i+ ofunction [f,f0]=Newton(x,y,x0)
( b- ~- M2 f8 ?+ F% J3 G9 W%求已知数据点的均差形式的牛顿插值多项式
( l! k. Y* h2 s6 D%已知数据点的x坐标向量：x
. F, P  _- ^7 _$ [+ f. N9 b%已知数据点的y坐标向量：y
3 z* e$ }- u& p  S%插值点的x坐标：x0
%求得的均差形式的牛顿插值多项式:f
%x0处的插值:f0
$ S8 ]. ]) Y2 [, x' a9 xsyms t;
# h& L+ ^% g8 wif(length(x)==length(y))
    n=length(x);
( l7 j3 q5 d% _/ u# m. S    c(1:n)=0.0;
1 f$ v# g+ @1 E* L1 `: s2 v+ G' jelse
    disp('x,y的维数不相得！');
& r$ I% L. C9 g2 U    return;
end
: w% m) l3 Q! v/ c: C3 df=y(1);
& @$ ^0 s7 c2 m% \6 N* @( C0 G, ~y1=0;
! c& N5 Z* ]# P/ P# Al=1;
for i=1:n-1
- r  y( T, ~; {    for j=i+1:n
/ s! q! [, s& U$ ?' w        y1(j)=(y(j)-y(i))/(x(j)-x(i));
) e$ r+ B# z, {2 X1 [    end
    c(i)=y1(i+1);
6 r4 ~0 a! b( o) s  t# J2 C    l=l*(t-x(i));
    f=f+c(i)*l;
    simplify(f);
    y=y1;
end
" J6 A6 ~( w; F  W/ uf0=subs(f,'t',x0);
2 I/ e* r4 Y6 g4 ]2 |