数值计算方法丁丽娟课后习题答案
【篇一：北京理工大学数值计算方法大作业数值实验1】
）书p14/4
分别将区间[?10,10]分为100,200,400等份，利用mesh或surf命
令画出二元函数的三维图形。
z=???|??|+?????? ??+?? +??
????+????+??【matlab求解】
[x,y]=meshgrid(-10:0.1:10);
a=exp(-abs(x));
b=cos(x+y);
c=1./(x.^2+y.^2+1);
z=a+b+c;
mesh(x,y,z);
[x,y]=meshgrid(-10:0.05:10);
a=exp(-abs(x));
b=cos(x+y);
c=1./(x.^2+y.^2+1);
z=a+b+c;
mesh(x,y,z);
[x,y]=meshgrid(-10:0.025:10); a=exp(-abs(x));
b=cos(x+y);
c=1./(x.^2+y.^2+1);
z=a+b+c;
mesh(x,y,z);
（二）书p7/1.3.2数值计算的稳定性（i）
取????= ??c语言程序—不稳定解 ??????+??????=ln1.2，按公式????=?????????? （n=1,2，…） ????#includestdio.h
#includeconio.h
#includemath.h
void main()
{
float m=log(6.0)-log(5.0),n;
int i;
i=1;
printf(y[0]=%-20f,m); while(i20)
{
n=1/i-5*m;
printf(y[%d]=%-20f,i,n);m=n;
i++;
if (i%3==0) printf(\n); }
getch();
}
（ii） c语言程序—稳定解????≈??[???????????? ??+?? +?? ??+??按公式 ????=??(??????????)#includestdio.h
#includeconio.h
#includemath.h
void main()
{
float m=(1/105.0+1/126.0)/2,n; k=n,n-1,n-2，…）（
【篇二：北京理工大学数值计算方法大作业数值实验4】 p260/1
考纽螺线的形状像钟表的发条，也称回旋曲线，它在直角坐标系中
的参数方程为
?? ??
= ?????????????????? ?? ??????????= ?????????????? ??曲线关于原点对称，取a=1，参数s的变化范围[-5,5]，容许误差限
分别是???????，,和?????????。选取适当的节点个数，利用数
值积分方法计算曲线上点的坐标，并画出曲线的图形。
【matlab求解】 ??
x=zeros(100,1);
y=zeros(100,1);
func1=inline(cos(1/2*(t.^2)));
func2=inline(sin(1/2*(t.^2)));
n=1;
for i= -5:0.1:5
x(n,1)=quad(func1,0,i,1e-6);
y(n,1)=quad(func2,0,i,1e-10);
n=n+1;
end
plot(x,y,b);
xlabel(x);
ylabel(y);
（代码运行图）
2.书p261/3 二重积分的计算
?????????（1） i= ??????? ?????????????(????
????+????)dy ???????????????(??+??)（2）
i= ?????? ????dy ??+??+??
【matlab求解】
（1）首先，用matlab自带的求解二重积分的命令对两个二重积分
分别求解
【第一问解】
【篇三：北京理工大学徐特立学院数值分析大作业上机
实验】
ss=txt>机实验选做
教材：数值计算方法（2011第一版）.丁丽娟,程杞元.高等教育出版
社
^以下代码作者原创^
超链接：
1.2（，，，）
2.2（，，，）
3.1（，，，，）
5.1（，，，）
5.3（，，，）
第一章：数值计算中的误差
2、题目简介：利用pi/4=1-1/3+1/5-1/7。。。级数计算pi的近似值。输入：误差值
输出：求和项数，并输出pi值
工具：c语言
运行环境：vc-6.0
计算公式及原理：利用pi/4=1-1/3+1/5-1/7。。。级数计算pi的近
似值，由数学原理可知误差会小于首次舍弃的项,可以编写循环实现。程序运行结果截图：
程序代码：（c语言）
#includestdio.h
void main()
{
printf(第一章第2题求pi,欢迎使用，请按提示操作。\n); int
i=1,n=0,k=1; double e,pi,er; printf(请输入误差(例如1e-4):);
scanf(%lf,e);
} printf(请稍候。。。\n); er=e; pi=0; while(er=e) {} pi*=4;
printf(%d项求和后可以达到%.10lf精度，这时pi=%.10lf\n,n,e,pi); getchar(); getchar(); pi+=k*1.0/i;k=-k;er=1.0/i; i+=2;n++;
第二章：解线性方程组的直接方法
2、题目简介：用matlab软件编程实现追赶法求解三对角方程组的算法，并考虑梯形电阻电路问题，电路如下：
工具：m语言
运行环境：matlab r2012.b
计算公式与原理：vr?220?8.1481 27
上述方程组可用矩阵表示为：