数学建模社区-数学中国
标题: 求助,数值分析问题! [打印本页]
作者: xuxiaoshuai 时间: 2006-12-3 10:39
标题: 求助,数值分析问题!
数值分析老师留了一个大作业,让结合一些工程数据作一个数学模型,用数值分析的牛顿或者插值或者线形方程组解出来,之后分析误差!
虽然我觉得自己学的不错,可是我知道这个作业不是要求做的特别复杂,就是因为这样才不知道怎么下手!老师说不需要用软件,只用手算就可以了!
哪位有做过这个作业的帮助一下了!谢谢了!或者有这个方面的例子发给我看一下!只要知道怎么下手就行!呵呵
邮箱xuhongyan198461@163.com
作者: ljack 时间: 2006-12-3 13:48
不会啊!应该在那个网站可以看到这方面的东西呢?
[em12]
作者: 37wing 时间: 2006-12-15 18:56
....
我也想要
我们也要用.....
作者: niumeng 时间: 2006-12-16 00:19
[em06]
作者: 灵幽二娃 时间: 2010-6-16 21:58
我手算过。。。可是现在书没在身边。。。。
! P% {( A. y5 P7 G看清华的数学实验吧
作者: 西雅图3173 时间: 2010-12-12 17:29
我们也用啊,是不是一个学校的啊
作者: matlabansys 时间: 2010-12-19 09:31
提示: 作者被禁止或删除 内容自动屏蔽
作者: wgbdarling 时间: 2012-5-11 21:35
数值分析,无非就是迭代,拟合,最小二乘
作者: 梦天涯M 时间: 2012-7-24 17:16
手算也太累了吧!这是matlab程序:. R6 _& U; C2 b0 G
function [f,f0]=Newton(x,y,x0)
u# v' h9 ]* L+ V( O; z%求已知数据点的均差形式的牛顿插值多项式
5 U' W0 `8 v3 M, g%已知数据点的x坐标向量:x
7 ~+ z; p g3 u%已知数据点的y坐标向量:y
% k7 O+ y: ^1 D% {+ T5 s+ p" k%插值点的x坐标:x0
4 R' q5 S- O9 f1 q% q% |%求得的均差形式的牛顿插值多项式:f
2 i" H& M! Z& l! m; O%x0处的插值:f0
) Q! L/ q% \1 ksyms t;
+ b' W! t0 Y8 p1 [6 p1 eif(length(x)==length(y)) M% V3 ?) i, Q: e/ R* M
n=length(x);
- ?8 H4 s- b. W c(1:n)=0.0;0 p/ u1 ]& d- T$ D4 m
else
4 u. m( x1 W2 N( A+ M) b disp('x,y的维数不相得!');
( Z2 L( o- k) `. Y' W return;
' q" n" G7 j! Y8 r" X5 D8 X1 _end
" L( u6 e/ ^: c6 z; ~8 W8 ^f=y(1);) B* `2 B; x+ }# b
y1=0;
, F7 j$ E) @" c9 A S o$ Ql=1;% k: s0 A4 p7 V9 d
for i=1:n-1
& ^- U2 F: s4 A' ^ for j=i+1:n* b0 k/ n! D: m" S
y1(j)=(y(j)-y(i))/(x(j)-x(i));9 G4 @; }! s8 N
end
& _9 C9 o5 x/ X- B c(i)=y1(i+1);
5 M6 j# w7 A" H; ?% N l=l*(t-x(i));
+ ?# F7 ~ P6 H% y# m f=f+c(i)*l;! @- d. j' @; x& A, @, V
simplify(f);
2 \4 ~& o+ F5 [. H& {' u y=y1; N7 X4 d6 Q4 |8 R& |
end
- t& R; a v" x! o f+ [f0=subs(f,'t',x0);
D" A( s( ~$ J+ b/ b
| 欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) |
Powered by Discuz! X2.5 |