求助,数值分析问题!
<p>数值分析老师留了一个大作业,让结合一些工程数据作一个数学模型,用数值分析的牛顿或者插值或者线形方程组解出来,之后分析误差!</p><p></p><p>虽然我觉得自己学的不错,可是我知道这个作业不是要求做的特别复杂,就是因为这样才不知道怎么下手!老师说不需要用软件,只用手算就可以了!</p><p>哪位有做过这个作业的帮助一下了!谢谢了!或者有这个方面的例子发给我看一下!只要知道怎么下手就行!呵呵</p><p>邮箱<a href="mailto:xuhongyan198461@163.com">xuhongyan198461@163.com</a></p> <p>....</p><p>我也想要</p><p>我们也要用.....</p> <p>不会啊!应该在那个网站可以看到这方面的东西呢?</p><p> </p> 我手算过。。。可是现在书没在身边。。。。看清华的数学实验吧 我们也用啊,是不是一个学校的啊 数值分析,无非就是迭代,拟合,最小二乘 手算也太累了吧!这是matlab程序:
function =Newton(x,y,x0)
%求已知数据点的均差形式的牛顿插值多项式
%已知数据点的x坐标向量:x
%已知数据点的y坐标向量:y
%插值点的x坐标:x0
%求得的均差形式的牛顿插值多项式:f
%x0处的插值:f0
syms t;
if(length(x)==length(y))
n=length(x);
c(1:n)=0.0;
else
disp('x,y的维数不相得!');
return;
end
f=y(1);
y1=0;
l=1;
for i=1:n-1
for j=i+1:n
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;
simplify(f);
y=y1;
end
f0=subs(f,'t',x0);
页:
[1]