求助,数值分析问题!

xuxiaoshuai

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

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

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

邮箱xuhongyan198461@163.com

37wing

....

我也想要

我们也要用.....

ljack

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

 

[em12]

niumeng

[em06]

灵幽二娃

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

西雅图3173

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

wgbdarling

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

梦天涯M

手算也太累了吧！这是matlab程序：
- f2 @% u; M3 g3 R, Afunction [f,f0]=Newton(x,y,x0)
: K6 T/ `6 W! l  {: g& k) S%求已知数据点的均差形式的牛顿插值多项式
%已知数据点的x坐标向量：x
2 P) N: w7 `* V( i/ y7 M%已知数据点的y坐标向量：y
# U' D7 K$ p4 z6 Q; ]$ {%插值点的x坐标：x0
; J- A- N. r) Q+ J- w2 z%求得的均差形式的牛顿插值多项式:f
2 h% s2 E% C( w1 j2 f# _! q/ M4 M%x0处的插值:f0
6 l' K- ^( |) V9 s- \syms t;
$ h! [: d( G/ E, I) z; sif(length(x)==length(y))
$ d: }9 ?) n5 l& x- x    n=length(x);
& L, W7 `+ c; u$ J: @! Q! d0 c    c(1:n)=0.0;
else
4 P3 p# S& p2 q( V    disp('x,y的维数不相得！');
    return;
end
f=y(1);
y1=0;
2 G1 d8 K) N5 f- \9 @# O2 tl=1;
! z- |) Q- ]# Q# O0 R9 v1 kfor i=1:n-1
    for j=i+1:n
' J  w  ?( t9 R        y1(j)=(y(j)-y(i))/(x(j)-x(i));
! ]  e* E& |4 \9 r; \$ }! q    end
1 A% K1 D: G- u( R' K* D+ b  e( l% k    c(i)=y1(i+1);
    l=l*(t-x(i));
    f=f+c(i)*l;
8 w' Z. ^$ r) Y3 R3 }% W# |    simplify(f);
; j9 G& d0 ]( Y7 H( m" k" W    y=y1;
& Z8 A4 U. e$ }0 Yend
4 k4 J3 q! r8 b0 \+ C, d. sf0=subs(f,'t',x0);