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

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

【篇一：北京理工大学数值计算方法大作业数值实验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

voidmain()

{

floatm=log(6.0)-log(5.0),n;

inti;

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(

);}

getch();

}

（ii）c语言程序—稳定解≈??[??+??

+????+??按公式=??(??)#includestdio.h

#includeconio.h

#includemath.h

voidmain()

{

floatm=(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;

fori=-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第一版）.丁丽