求助,数值分析问题!

xuxiaoshuai

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

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

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

邮箱xuhongyan198461@163.com

37wing

....

我也想要

我们也要用.....

ljack

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

 

[em12]

niumeng

[em06]

灵幽二娃

我手算过。。。可是现在书没在身边。。。。
1 i  t/ H. i) b看清华的数学实验吧

西雅图3173

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

wgbdarling

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

梦天涯M

手算也太累了吧！这是matlab程序：
function [f,f0]=Newton(x,y,x0)
! N* D' [7 D7 O- I# v/ l%求已知数据点的均差形式的牛顿插值多项式
9 M8 K. c- \6 P% r4 K8 D" K%已知数据点的x坐标向量：x
%已知数据点的y坐标向量：y
! w5 Q" q+ H$ {6 W: b8 [" v%插值点的x坐标：x0
%求得的均差形式的牛顿插值多项式:f
%x0处的插值:f0
5 k# H+ a" Z; Wsyms t;
  [: ?9 l$ ^3 ]1 M- z- E% Zif(length(x)==length(y))
2 N  z. @8 }/ D; l    n=length(x);
- c# |; W. K& a6 ^9 c. }( {; _( z    c(1:n)=0.0;
else
    disp('x,y的维数不相得！');
( f. u0 S6 Q& S  @$ {% O    return;
0 R* S2 _8 x9 [end
* Z4 V3 [/ w3 F3 a6 If=y(1);
y1=0;
) m4 E3 d" k% A: l. [" {. Kl=1;
8 L' `  e0 x: ffor i=1:n-1
    for j=i+1:n
: q* X* |6 W) |% c1 m4 ^        y1(j)=(y(j)-y(i))/(x(j)-x(i));
8 J" q! V: [: _! [5 m6 k5 {    end
$ F8 w$ z1 \( {: I- ]    c(i)=y1(i+1);
    l=l*(t-x(i));
    f=f+c(i)*l;
    simplify(f);
    y=y1;
$ U* S1 [. g2 xend
f0=subs(f,'t',x0);
& f" A( A, r9 Q( l