数值计算方法丁丽娟课后习题答案

数值计算方法丁丽娟课后习题答案

【篇一：北京理工大学数值计算方法大作业数值实验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 上述方程组可用矩阵表示为：