计算方法上机作业

计算方法上机报告。

姓名：学号：

班级：上课班级：

1） 算法思想。

2）算法结构。

3）matlab源程序。

clear; %清除工作空间变量。

clc; %清除命令窗口命令。

m=input('请输入有效数字的位数m输入有效数字的位数。

s=0for n=0:50

t=(1/16^n)*(4/(8*n+1)-2/(8*n+4)-1/(8*n+5)-1/(8*n+6));

if t<=10^(-m判断通项与精度的关系。

break;

endend;

fprintf('需要将n值加到n=%d',n-1需要将n值加到的数值。

for i=n-1:-1:0

t=(1/16^i)*(4/(8*i+1)-2/(8*i+4)-1/(8*i+5)-1/(8*i+6));

s=s+t求和运算。

ends=vpa(s,m控制s的精度

4）结果与分析。

1）算法思想。

海底光缆线的长度**模型。

2）算法结构。

3）matlab源程序。

clear;

clc; x=0:1:20产生从0到20含21个等分点的数组。

x=0:0.2:20;

y=[9.01,8.96,7.

96,7.97,8.02,9.

05,10.13,11.18,12.

26,13.28,13.32,12.

61,11.29,10.22,9.

15,7.90,7.95,8.

86,9.81,10.80,10.

93等分点位置的深度数据。

n=length(x等分点的数目。

n=length(x);

% 求三次样条插值函数s(x)

m=yfor k=2:3计算二阶差商并存放在m中。

for i=n:-1:k;

m(i)=(m(i)-m(i-1))/x(i)-x(i-k+1));

endend

h(1)=x(2)-x(1计算三对角阵系数a,b,c及右端向量d

for i=2:n-1;

h(i)=x(i+1)-x(i);

c(i)=h(i)/(h(i)+h(i-1));

a(i)=1-c(i);

b(i)=2;

d(i)=6*m(i+1);

end m(1)=0选择自然边界条件。

m(n)=0;

b(1)=2;

b(n)=2;

c(1)=0;

a(n)=0;

d(1)=0;

d(n)=0;

u(1)=b(1对三对角阵进行lu分解。

y1(1)=d(1);

for k=2:n;

l(k)=a(k)/u(k-1);

u(k)=b(k)-l(k)*c(k-1);

y1(k)=d(k)-l(k)*y1(k-1);

endm(n)=y1(n)/u(n追赶法求解样条参数m(i)

for k=n-1:-1:1;

m(k)=(y1(k)-c(k)*m(k+1))/u(k);

ends=zeros(1,n);

for m=1:n;

k=1; for i=2:n-1

if x(m)<=x(i);

k=i-1;

break;

elsek=i;

endend

h=x(k+1)-x(k在各区间用三次样条插值函数计算x点处的值。

x1=x(k+1)-x(m);

x2=x(m)-x(k); s(m)=(m(k)*(x1^3)/6+m(k+1)*(x2^3)/6+(y(k)-(m(k)*(h^2)/6))*x1+(y(k+1)-(m(k+1)*(h^2)/6))*x2)/h;

end% 计算所需光缆长度。

l=0; %计算所需光缆长度。

计算方法上机作业

表2 新算法迭代结果1 初值 0，精度ep 初始选择因子p 修正因子 0.1，最大迭代次数n 0，局部修正次数m分别取1,2,3次，f 49881192646492934000000000000.000000 表3 新算法迭代法结果2 初值 0，精度ep 初始选择因子p 修正因子 0.05，最大迭代...

计算方法上机作业

计算方法第四次上机报告。2.用欧拉方法解初值。y 10x 1 y 0 x 1 y 0 0 取步长h 0.1,保留5位有效数字，并与准确解相比较。分析 该题目考察欧拉方法解初值问题。程序如下 function heun a,b,y0,n h b a n x a h b y y0 ones 1,n 1 ...

数值计算方法上机作业

地球科学与技术学院地球探测与信息技术研11 6班孙殿强s1 实验2.1迭代法 初始值与收敛性。实验目的 初步认识非线性问题的迭代法与线性问题的迭代法的差别，迭代法及初始值与迭代收敛性的关系。问题提出 迭代法是求解非线性方程 组 的基本方法，与线性方程的情况一样，其构造方法多种多样，但关键是怎样才能使...