PAGE
PAGE #
百度文库-
百度文库-让每个人平等地提升自我
PAGE
PAGE #
内江师范学院
数值分析\\
验证性实验报告
编制覃燕梅 审定牟廉明
专业： 数学与应用数学
班级：数学系2009级 班
学号：
姓名：
数学与信息科学学院
2011年8月
目录
TOC \o "1-5" \h \z \o "Current Document" 说明： 1
学生在做实验之前必须要准备实验，主要包括预习与本次实验相关的理论知识，熟练与本次实验相关 的软件操作，收集整理相关的实验参考资料， 要求学生在做实验时能带上充足的参考资料； 若准备不充分,
则学生不得参加本次实验，不得书写实验报告； 1
/ \
/ \
要求学生要认真做实验，主要是指不得迟到、早退和旷课，在做实验过程中要严格遵守实验室规章制
度，认真完成实验内容，极积主动地向实验教师提问等；若学生无故旷课，则本次实验计为 0分； 1
学生要认真工整地书写实验报告，实验报告的内容要紧扣实验的要求和目的，不得抄袭他人的实验报
告;
实验成绩=实验准备+实验态度+实验报告书写+实验报告内容，根据实验成绩来判定学生是否达到了实
验目的。
验目的。
2
验证性实验成绩汇总
\
实 验
实 验
实 验
实 验 四
实 验 五
实 验 六
实 验 七
验证性实验成绩-
（实验一 +实验二+实验三+实验 四+实验五+实验六+实验七）/7
分数
说明：
学生在做实验之前必须要准备实验，主要包括预习与本次实验相关的理论知识，熟练 与本次实验相关的软件操作，收集整理相关的实验参考资料，要求学生在做实验时能
带上充足的参考资料；若准备不充分，则学生不得参加本次实验，不得书写实验报告;
要求学生要认真做实验，主要是指不得迟到、早退和旷课，在做实验过程中要严格遵
PAGE
PAGE #
守实验室规章制度，认真完成实验内容，极积主动地向实验教师提问等；若学生无故 旷课，则本次实验计为0分；
学生要认真工整地书写实验报告，实验报告的内容要紧扣实验的要求和目的，不得抄 袭他人的实验报告；
实验成绩=实验准备+实验态度+实验报告书写+实验报告内容，根据实验成绩来判定学 生是否达到了实验目的。
PAGE
PAGE #
实验名称： Lagra nge插值（实验一） 指导教师： 覃燕梅
实验时数： 2 实验设备： 计算机
实验日期：2011年_月_日 实验地点： 第五教学楼北802
实验目的：
掌握Lagrange插值数值算法，能够根据给定的函数值表达求出插值多项式和函数在某一点的近似值。
实验准备：
在开始本实验之前，请回顾教科书的相关内容；
需要准备安装 Windows XP Professional操作系统和 Matlab、C++、VF等数学软件。
实验内容及要求
已知数据如下:
Xi
10
11
12
\ 13
yi
要求：
试用分别用线性插值、二次插值、三次插值分别计算在点 x 11.75时的函数近似值.
实验过程：
function [y0,N]=Lagrange_eval(X，丫,x0) m=le ngth(X);
N=zeros(m,1);
y0=0;\
for i=1:m
N(i)=1;
for j=1:m
if j~=i
N(i)=N(i)*(x0-X(j))/(X(i)-X(j))； end
end
old = digits;
digits (16);
y0=vpa(y0+Y(i)*N(i)); digits(old);
end
令 X=[11,12];Y=[,];x0=;
运行[y0,N]=Lagrange_eval(X,Y,x0)
PAGE
PAGE #
结果为 y0
结果为 y0 =
这是用线性插值所求得的结果
令 X=[10,11,12];Y=[,,];x0=;
运行[yO,N]=Lagrange_eval(X,Y,x0)
y0 = 这是用二次插值所求得的结果
令 X=[10,11,12,13];Y=[,,,];x0=;
运行[y0,N]=Lagrange_eval(X,Y,x0)
y0 = 这是用三次插值所求得的结果
实验总结（由学生填写）：
教师对本次实验的评价（下面的表格由教师填写）：
评价细目
\ 定性评价
量化分
实验准备（20分）
充分（15-20）□ 基本充分（0-15）□
实验态度（20分）
认真（15-20）□
比较认真（10-15） 口
/不认真（0-10）□
实验报告书写（20分）
规范（15-20）□
基本规范（10-15） 口
不规范（0-10）□
实验报告内容（40分）
符合实验要求
（30-40）□
基本符合实验要求/ r
（10-30） 口
不符合实验要求
（0-10）□
实验结论及实验成绩
达到实验目的
（80-100）□
基本达到实验目的
（60-80） 口
没有达到实验目的
（0-60）□
PAGE
PAGE #
实验名称： 曲线拟合的最小二乘方法(实验二) 指导教师： 覃燕梅
实验时数： 2 实验设备： 计算机
实验日期：2011年_月_日 实验地点： 第五教学楼北802
实验目的：
掌握最小二乘方法，并能根据给定数据求其最小二乘一次或二次多项式，然后进行曲线拟合。
实验准备：
在开始本实验之前，请回顾教科书的相关内容；
需要准备安装 Windows XP Professional操作系统和 Matlab、C++、VF等数学软件。
实验内容及要求
已知一组实验数据如下：.
i
1
2
3
4
Xi
2
4
6
8
yi
2
11
28
40
试用最小二乘法求一次拟合多项式，并将拟合曲线画出来.
实验过程：
MATLAB程序
x=2:2:8;
y=[2 11 28 40];
A=polyfit(x,y,1)
z=polyval(A,x);
plot(x,y, 'k*' ,x,z, 'r')
PAGE
PAGE #
40
实验总结（由学生填写）：
教师对本次实验的评价（下面的表格由教师填写）：
评价细目
\ 定性评价
量化分
实验准备（20分）
充分（15-20）□ 基本充分（0-15）□
实验态度（20分）
认真（15-20）□
比较认真（10-15） 口
/■不认真（0-10）□
实验报告书写（20分）
规范（15-20）□
基本规范（10-15） 口
不规范（0-10）□
实验报告内容（40分）
符合实验要求
（30-40）□
基本符合实验要求/ r
（10-30） 口
不符合实验要求
（0-10）□
实验结论及实验成绩
达到实验目的
（80-100）□
基本达到实验目的
（60-80） 口
没有达到实验目的
（0-60）□
PAGE
PAGE #
实验名称： Romberg积分法(实验三) 指导教师： 覃燕梅
实验时数： 2 实验设备： 计算机
实验日期：2011年_月_日 实验地点： 第五教学楼北802
实验目的：
掌握Romberg算法，并能根据给定的精度要求计算定积分。
实验准备：
在开始本实验之前，请回顾教科书的相关内容；
需要准备安装 Windows XP Professional操作系统和 Matlab、C++、VF等数学软件。
实验内容及要求
1 sin x
用Romberg法求函数积分| dx，精度为0.5e 6 .
0 x
实验过程：
fun cti on [quad,R]=Romberg(f,a,b,eps)
h=b-a;
R(1,1)=h*(feval(f,a)+feval(f,b))/2;
M=1;J=0;err=1;
while err>eps
J=J+1;
h=h/2;
S=0;
for p=1:M
x=a+h*(2*p-1);
S=S+feval(f,x);
end
R(J+1,1)=R(J,1)/2+h*S;
M=2*M;
for k=1:J
R(J+1,k+1)=R(J+1,k)+(R(J+1,k)-R(J,k))/(4Ak-1);
end
err=abs(R(J+1,J)-R(J+1,J+1));
end
quad=R(J+1,J+1);
format( 'long' , 'eng' )
fun cti on f=f1(x)
f=si n(x)/x;
运行
f=@f1;b=1;a=10A-100;eps=10A-4;
[quad,R]=Romberg(f,a,b,eps)
PAGE
PAGE #
quad =
R =
+000 +000
+000
f=@f1;b=1;a=10A-100;eps=10A-7;
[quad,R]=Romberg(f,a,b,eps)
quad =
R =
+000 +000 +000
+000 +000
+000
实验总结（由学生填写）：
教师对本次实验的评价（下面的表格由教师填写）：
评价细目
\ 定性评价 /
量化分
实验准备（20分）
充分（15-20）□ 基本充分（0-15）□
实验态度（20分）
认真（15-20）□
比较认真（10-15） 口
不认真（0-10）□
实验报告书写（20分）
规范（15-20）□
基本规范（10-15） 口
不规范（0-10）□
实验报告内容（40分）
符合实验要求 \
（30-40）□
基本符合实验要求
\（ 10-30） 口
不符合实验要求
（0-10）□
实验结论及实验成绩
达到实验目的
（80-100）□
基本达到实验目的
（60-80） 口
没有达到实验目的
（0-60）□
PAGE
PAGE #
实验名称： 改进的欧拉方法(实验四) 指导教师： 覃燕梅
实验时数： 2 实验设备： 计算机
实验日期：2011年 月 日实验地点： 第五教学楼北802
一 / 实验目的：
掌握欧拉方法与改进的欧拉方法，能够用改进的欧拉方法求解微分方程。
实验准备：
在开始本实验之前，请回顾教科书的相关内容；
需要准备安装 Windows XP Professional操作系统和 Matlab、C++、VF等数学软件。
实验内容及要求
用改进的欧拉方法求初值问题：
y' y 0 (Ox 1.5) \
y(o)=1
实验过程：fun cti on E=Me ndEuler(f,a,b,N,ya)
h=(b-a)/N;
y=zeros(1,N+1);
x=zeros(1,N+1);
y(1)=ya;
x=a:h:b;
for i=1:N
y1=y(i)+h*feval(f,x(i),y(i));
y2=y(i)+h*feval(f,x(i+1),y1);
y(i+1)=(y1+y2)/2;
end
E=[x',y']
fun cti on f=f1(x,y)
f=-y;
fun cti on y=solvef1(x)
y=exp(-x);
运行程序：
f=@f1;a=0;b=;N=15;ya=1;
E=MendEuler(f,a,b,N,ya);y=solvef1(a:(b-a)/N:b);m=[E,y'] \ /
PAGE
PAGE #
实验总结（由学生填写）：
教师对本次实验的评价（下面的表格由教师填写）：
评价细目
\ 定性评价 /
量化分
实验准备（20分）
充分（15-20）□ 基本充分（0-15）□
实验态度（20分）
认真（15-20）□
比较认真（10-15） 口
不认真（0-10）□
实验报告书写（20分）
规范（15-20）□
基本规范（10-15） 口
不规范（0-10）□
实验报告内容（40分）
符合实验要求 \
（30-40）□
基本符合实验要求
\（ 10-30） 口
不符合实验要求
（0-10）□
实验结论及实验成绩
达到实验目的
（80-100）□
基本达到实验目的
（60-80） 口
没有达到实验目的
（0-60）□
PAGE
PAGE #
实验名称： Newton方法(实验五) 指导教师： 覃燕梅
实验时数： 2 实验设备： 计算机
实验日期：2011年_月_日 实验地点： 第五教学楼北802
实验目的：
掌握Newt on迭代算法，能够根据所给方程求出在某一点附近的根。
实验准备：
在开始本实验之前，请回顾教科书的相关内容；
需要准备安装 Windows XP Professional操作系统和 Matlab、C++、VF等数学软件。
实验内容及要求
公元1225年，Lenardo宣布他求得方程 \
x3 2x2 10x 20 0
的一个根x 1.368808107 ?当时颇为轰动，但无人知道他是用什么方法得到的. 现在，请你试试用Newton
迭代法求解这个结果.
实验过程：
fun cti on [x,k]=Newt on _method(f,x0,emg)
[f1,d1]=feval(f,x0);
k=1;
x(1)=x0;
x(2)=x(1)_f1/d1;
while abs(f1)>emg
k=k+1;
[f1,d1]=feval(f,x(k)); x(k+1)=x(k)-f1/d1;
[f2,d2]=feval(f,x(k+1));
while abs(f2)>abs(f1)
x(k+1)=x(k)-f1/d1;
[f2,d2]=feval(f,x(k+1));
end
end
fun cti on [f,d]=fu nc3(x)
f=xA3+2*xA2+10*x-20;
syms y
d1=yA3+2*yA2+10*y-20;
d=subs(diff(d1),y,x);
PAGE
PAGE #
运行:
实验总结（由学生填写）：
教师对本次实验的评价（下面的表格由教师填写）：
评价细目
\ 定性评价 /
量化分
实验准备（20分）
充分（15-20）□ 基本充分（0-15）□
实验态度（20分）
认真（15-20）□
比较认真（10-15） 口
不认真（0-10）□
实验报告书写（20分）
规范（15-20）□
基本规范（10-15） 口
不规范（0-10）□
实验报告内容（40分）
符合实验要求 \
（30-40）□
基本符合实验要求/
（10-30） 口
不符合实验要求
（0-10）□
实验结论及实验成绩
达到实验目的
（80-100）□
基本达到实验目的
（60-80） 口
没有达到实验目的
（0-60）□
x=x1; PAGE
x=x1;
PAGE #
实验名称： 雅可比迭代方法(实验六) 指导教师： 覃燕梅
实验时数： 2 实验设备： 计算机
实验日期：2011年_月_日 实验地点： 第五教学楼北802
实验目的：
掌握列主元的高斯消去法思想，能够利用列主元的高斯消去法求解任意阶数的线性方程组。
实验准备：
在开始本实验之前，请回顾教科书的相关内容；
需要准备安装 Windows XP Professional操作系统和 Matlab、C++、VF等数学软件。
实验内容及要求
用的雅可比迭代方法求解下列方程组：
TOC \o "1-5" \h \z 10 -1 2 0 为 6 \
1 11 -1 3 x2 25 \
2 -1 10 -1 % 11 .
0 3 -1 8 x4 15
实验过程：
fun cti on [x,k]=Jacobimethod(A,b,xO,N,emg)
n=len gth(A);
x1=zeros (n ,1);
x2=zeros (n ,1);
x1=x0;k=0;
r=max(abs(b-A*x1));
while r>emg
for\ i=1: n
sum=0;
for j=1: n
if i~=j
sum=sum+A(i,j)*x1(j);
end
end
x2(i)=(b(i)-sum)/A(i,i);
end
r=max(abs(x2-x1));
x1=x2;
k=k+1;
if k>N
disp( ' fiu，ue§°u￡ ?;?□ '
return ;
end
end
PAGE
PAGE #
运行
A=[10,-1,2,0;-1,11,-1,3;2,-1,10,-1;0,3,-1,8];b=[6,25,-11,15]';x0=[0,0,0,0]';
[x,k]=Jacobimethod(A,b,x0,100,10A-5)
k =
16
实验总结（由学生填写）：
教师对本次实验的评价（下面的表格由教师填写）：
评价细目
\ 定性评价
量化分
实验准备（20分）
充分（15-20）□ 基本充分（0-15）□
实验态度（20分）
认真（15-20）□
比较认真（10-15） 口
/■不认真（0-10）□
实验报告书写（20分）
规范（15-20）□
基本规范（10-15） 口
不规范（0-10）□
实验报告内容（40分）
符合实验要求
（30-40）□
基本符合实验要求/ r
（10-30） 口
不符合实验要求
（0-10）□
实验结论及实验成绩
达到实验目的
（80-100）□
基本达到实验目的
（60-80） 口
没有达到实验目的
（0-60）□
PAGE
PAGE #
实验名称： 列主元的高斯消去法（实验七） 指导教师： 覃燕梅
实验时数： 2 实验设备： 计算机
实验日期：2011年_月_日 实验地点： 第五教学楼北802
实验目的：
掌握列主元的高斯消去法思想，能够利用列主元的高斯消去法求解任意阶数的线性方程组。
实验准备：
在开始本实验之前，请回顾教科书的相关内容；
需要准备安装 Windows XP Professional操作系统和 Matlab、C++、VF等数学软件。
实验内容及要求
解下列方程组
0.729
0.81
0.9
0.6867
1
1
1
X2
0.8338
1.331
1.21
1.1
X3
1.000
实验过程：
fun cti on x=Gauss_pivot(A,b)
n=len gth(b);
x=zeros( n,1);
c=zeros(1, n);
d1=0;
for i=1: n-1
max=abs(A(i,i));
m=i;
for j=i+1: n
if max<abs(A(j,i)) max=abs(A(j,i));
m=j;
end
end
if (m~=i)
for k=i: n
c(k)=A(i,k);
A(i,k)=A(m,k);
A(m,k)=c(k);
end
d1=b(i);
b(i)=b(m);
b(m)=d1;
end
for k=i+1: n
PAGE
PAGE #
for j=i+1: n
A(k,j)=A(k,j)-A(i,j)*A(k,i)/A(i,i); end b(k)=b(k)-b(i)*A(k,i)/A(i,i); A(k,i)=0;
end
end
%huidaiqiujie
x( n)=b( n)/A( n,n);
for i=n-1:-1:1
sum=0;
for j=i+1: n
sum=sum+A(i,j)*x(j);
end
x(i)=(b(i)-sum)/A(i,i);
end
运行
A=[,,;1,1,1;,,];b=[,,]';
>> x=Gauss_pivot(A,b)
得
x =
实验总结（由学生填写）：
教师对本次实验的评价（下面的表格由教师填写）：
评价细目
\ 定性评价 /
量化分
实验准备（20分）
充分（15-20）□ 基本充分（0-15）□
实验态度（20分）
认真（15-20）□
比较认真（10-15） 口
不认真（0-10）□
实验报告书写（20分）
规范（15-20）□
基本规范（10-15） 口
不规范（0-10）□
实验报告内容（40分）
符合实验要求 \
（30-40）□
基本符合实验要求/
（10-30） 口
不符合实验要求
（0-10）□
实验结论及实验成绩
达到实验目的
（80-100）□
基本达到实验目的
（60-80） 口
没有达到实验目的
（0-60）□
推荐访问:实验报告 数值 验证 实验 数值分析验证性实验报告