数值分析韩旭里答案.doc

数值分析韩旭里答案 【篇一：数值分析上机题目】 1631110xxxx 材料科学与工程学院 一．第 2 章插值法 l2.7 给定数据表 2-15. 用 newton 插值公式计算 3 次插值多项式 n3(x). 表 2-15 f(x) 1.25 1.5 2.50 1.00 5.50 a. matlab 代码如下， two.m, %第二章， p45, 练习题 2 第七题 clear(); x=[1,1.5,0,2]; y(:,1)=[1.25,2.50,1.00,5.50];% 已知点集合 x 和 y syms t w; w(1)=1; %计算基函数序列 w 和差商表 y，以及函数序列的权数 算的牛顿三次多项式表述为 t 的函数 for j=2:length(x) for  diag(y),  计 i=j:length(x) y(i,j)=(y(i,j-1)-y(i-1,j-1))/(x(i)-x(i-j+1)); i=i+1; end w(j)=prod(t-x(1:j-1)); j=j+1; end disp( 三次牛顿插值多项式为 ); disp(collect(w*diag(y))); plot(x,y(:,1),*); hold on; fplot(collect(w*diag(y)),[-0.5,2.5]); legend({ 已知点集 ,三次牛顿插值多项式函 数},location,northwest,fontsize,14); xlabel(x,fontsize,16); ylabel(y,fontsize,16); hold off; 计算结果如下： 二．第 3 章函数逼近与数据拟合 a. matlab 代码， three.m ， %第三章函数逼近与数据拟合， p68 练习题，第 2 题 clear(); syms x; %所使用的非线性基函数序列，用符号表示 y=abs(x);% 被逼近函数 f=[1,x^2,x^4]; %求解法方程的系数矩阵 a*gn=b ，其中 a 和 b 均为行向量 gn=ones(length(f),length(f)); for i=1:length(f) for j=1:length(f) gn(i,j)=int(f(i)*f(j),-1,1);j=j+1; end b(i)=int(f(i)*y,-1,1); i=i+1; end a=b/gn;% 最佳平方逼近的系数行向量 disp( 逼近函数表达式 ); disp(vpa(f*a)); disp( 最佳函数逼近得平方误差 ); disp(vpa(int(y^2,-1,1)-a*b)); fplot(y,[-1,1]); hold on; fplot(a*f,[-1,1]); legend({ 被逼近函数 ,逼近函数 },location,north,orientation,horizontal,fontsize,16,fontweight,b old); xlabel(x,fontsize,20,fontweight,bold); ylabel(y,fontsize,20,fontweight,bold); hold off; 运行结果如下： .第 4 章数值积分与数值微分 例 4.9 用 romberg 求积法计算定积分 0 1sin?(??) ?? a. matlab 代码， four.m %romberg 求积公式，外推原理 clear(); clear(); format long; a=0; b=1; t(1,1)=(b-a)/2*(f(a)+f(b)); t(2,1)=1/2*t(1,1)+(b-a)/2*f((a+b)/2); t(1,2)=(4*t(2,1)-t(1,1))/(4-1); col=2; while abs(t(1,col)-t(1,col-1))0.5*10^-6 %t(1,col) 对应的计算的是多少步的值 ,col →coln 关系 col=col+1;% 此时求得是第 n+1 次均分后的结果，使用的是第 的结果，注意在矩阵 %计算的第 n 斜列是第 n-1 次均分的结果  n 次 for j=1:colif j==1 h=(b-a)/2^(col-2);%  使用  n+1  之前的第  n 次结果 【篇二：数值分析 a 教学】 一、课程基本信息 二、课程目的和任务 “数值分析 ”是理工科院校计算数学、力学、物理、计算机软件等专业的学生必须掌握的一门重要的基础课程。它研究用计算机求解各种数学问题的数值计算方法及其理论与软件实现，既有纯数学高度抽象性与严密科学性的特点，又有应用的广泛性与实际实验的高度技术性的特点，是一门与计算机使用密切结合的实用性很强的