《数值分析计算实习题》由会员分享,可在线阅读,更多相关《数值分析计算实习题(17页珍藏版)》请在人人文库网上搜索。
1、数值分析计算实习题姓名: 学号: 班级: 第二章1、程序代码Clear;clc;x1=0.2 0.4 0.6 0.8 1.0;y1=0.98 0.92 0.81 0.64 0.38;n=length(y1);c=y1(:);for j=2:n %求差商for i=n:-1:jc(i)=(c(i)-c(i-1)/(x1(i)-x1(i-j+1);endendsyms x df d;df(1)=1;d(1)=y1(1);for i=2:n %求牛顿差值多项式df(i)=df(i-1)*(x-x1(i-1);d(i)=c(i-1)*df(i);endP4=vpa(sum(d),5) %P4即为4次牛。
2、顿插值多项式,并保留小数点后5位数pp=csape(x1,y1, variational);%调用三次样条函数q=pp.coefs;q1=q(1,:)*(x-.2)3;(x-.2)2;(x-.2);1;q1=vpa(collect(q1),5)q2=q(1,:)*(x-.4)3;(x-.4)2;(x-.4);1;q2=vpa(collect(q2),5)q3=q(1,:)*(x-.6)3;(x-.6)2;(x-.6);1;q3=vpa(collect(q3),5)q4=q(1,:)*(x-.8)3;(x-.8)2;(x-.8);1;q4=vpa(collect(q4),5)%求解并化简多项式2。
3、、运行结果P4 =0.98*x - 0.3*(x - 0.2)*(x - 0.4) - 0.625*(x - 0.2)*(x - 0.4)*(x - 0.6) - 0.20833*(x - 0.2)*(x - 0.4)*(x - 0.8)*(x - 0.6) +