求助,数值分析问题!

xuxiaoshuai

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

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

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

邮箱xuhongyan198461@163.com

37wing

....

我也想要

我们也要用.....

ljack

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

 

[em12]

niumeng

[em06]

灵幽二娃

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

西雅图3173

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

wgbdarling

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

梦天涯M

手算也太累了吧！这是matlab程序：
function [f,f0]=Newton(x,y,x0)
%求已知数据点的均差形式的牛顿插值多项式
$ S4 j% C1 X% e$ t+ K! L%已知数据点的x坐标向量：x
%已知数据点的y坐标向量：y
%插值点的x坐标：x0
4 `+ n& i3 M+ K  Z%求得的均差形式的牛顿插值多项式:f
4 e! v* i) q: S$ [' n1 X%x0处的插值:f0
/ i/ f& R1 t& p6 {, o, Gsyms t;
9 U$ ?- g. Q% p4 t0 Y) Nif(length(x)==length(y))
" V* g9 B, t% R3 \; P$ ]    n=length(x);
" [! I4 V; ~4 X0 M3 z( G    c(1:n)=0.0;
else
    disp('x,y的维数不相得！');
7 K  u% a8 f4 H2 @- g3 N' W5 B    return;
  Y$ p  h' F1 X! d2 m( {6 q- Yend
% ]8 m1 d* _$ ?' Z+ Qf=y(1);
9 {( N1 C, J! A! v2 Xy1=0;
3 w/ H; G$ ~0 a1 G8 z+ x5 Z2 K3 Kl=1;
for i=1:n-1
4 M; c: |3 d( e    for j=i+1:n
& s' W: e% V2 C) R% ?        y1(j)=(y(j)-y(i))/(x(j)-x(i));
    end
    c(i)=y1(i+1);
    l=l*(t-x(i));
    f=f+c(i)*l;
& n$ Q1 U# y7 x2 e9 ]    simplify(f);
    y=y1;
end
; T+ t  J3 l3 d' C8 B. D& g) nf0=subs(f,'t',x0);
. |3 w7 o# Y" ~- E1 ~( H