当前位置：文档之家 › 西安交通大学计算方法B上机报告

西安交通大学计算方法B上机报告

  • 格式：doc
  • 大小：747.50 KB
  • 文档页数：33

计 算 方 法 上 机 报 告

姓 名： 学 号： 班 级： 能动 上课班级： 计算方法（B）上机报告 1 题目及求解： 一、对以下和式计算：0681581482184161nnnnSn，要求： ① 若只需保留11个有效数字，该如何进行计算； ② 若要保留30个有效数字，则又将如何进行计算； 1 算法思想 （1）根据精度要求估计所加的项数，可以使用后验误差估计，通项为：

142111416818485861681nnnannnnn

；

（2）为了保证计算结果的准确性，写程序时，从后向前计算； （3）使用Matlab时，可以使用以下函数控制位数： digits(位数)或vpa(变量，精度为数)

2 算法结构

;0s



681581482184161nnnnt

n；

for 0,1,2,,ni if 10mt

end; for ,1,2,,0niii ;sst

3 Matlab源程序 clear; %清除工作空间变量 clc; %清除命令窗口命令 m=input('请输入有效数字的位数m='); %输入有效数字的位数 s=0; for n=0:50 计算方法（B）上机报告 2 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; end end; fprintf('需要将n值加到n=%d\n',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; %求和运算 end s=vpa(s,m) %控制s的精度 4 结果与分析 若保留11位有效数字，则n=7，此时求解得： s =3.1415926536； 若保留30位有效数字时，则n=22, 此时求解得： s =3.14159265358979323846264338328。 通过上面的实验结果可以看出，通过从后往前计算，这种算法很好的保证了计算结果要求保留的准确数字位数的要求。

二、某通信公司在一次施工中，需要在水面宽度为20米的河沟底部沿直线走向铺设一条沟底光缆。在铺设光缆之前需要对沟底的地形进行初步探测，从而估计所需光缆的长度，为工程预算提供依据。已探测到一组等分点位置的深度数据(单位:米)如下表所示： 计算方法（B）上机报告

3 分点 0 1 2 3 4 5 6 深度 9.01 8.96 7.96 7.97 8.02 9.05 10.13 分点 7 8 9 10 11 12 13 深度 11.18 12.26 13.28 13.32 12.61 11.29 10.22 分点 14 15 16 17 18 19 20 深度 9.15 7.90 7.95 8.86 9.81 10.80 10.93 ① 请用合适的曲线拟合所测数据点； ② 预测所需光缆长度的近似值，作出铺设河底光缆的曲线图； 1 算法思想 如果使用多项式差值，则由于龙格现象，误差较大，因此，用相对较少的插值数据点作插值，可以避免大的误差，但是如果又希望将所得数据点都用上，且所用数据点越多越好，可以采用分段插值方式，即用分段多项式代替单个多项式作插值。分段多项式是由一些在相互连接的区间上的不同多项式连接而成的一条连续曲线，其中三次样条插值方法是一种具有较好“光滑性”的分段插值方法。 在本题中，假设所铺设的光缆足够柔软，在铺设过程中光缆触地走势光滑，紧贴地面，并且忽略水流对光缆的冲击。海底光缆线的长度预测模型如图2-1所示，光缆从A

点铺至B点，在某点处的深度为ih。

图2-1 海底光缆线的长度预测模型 计算光缆长度时，用如下公式： 200()Lfxds

20'2

0()1()fxfxdx

191'20()1()kkkfxfxdx

 计算方法（B）上机报告 4 22()()xy

2 算法结构 1） For ni,,2,1,0 1.1 iiMy 2） For 2,1k 2.1 For knni,,1, 2.1.1 ikiiiiMxxMM)/()(1 3）101hxx 4）For 1-,,2,1ni 4.1 11iiihxx 4.2 bacchhhiiiiii2;1;)/(11 4.3 iidM16 5）0000;;cMdMdnn nnnbab2;;20

6）1111,db 7）获取M的矩阵元素个数，存入m 8）For mk,,3,2

8.1 kkkla1/ 8.2 kkkkclb1- 8.3 kkkkld1- 9）mmmM/ 10）For 1,,2,1mmk 10.1 kkkkkMMc/)(1 计算方法（B）上机报告 5 11）获取x的元素个数存入s 12）k1

13） For 1,,2,1si

13.1 if ixx~ then ki；break else ki1 14）xxxxxxhxxkkkkˆ~;~;11

yhxhMyxhMyxMxMkkkkkk~/]ˆ)6()6(6ˆ6[2211331

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=y; for k=2:3; %计算二阶差商并存放在M中 for i=n:-1:k; M(i)=(M(i)-M(i-1))/(x(i)-x(i-k+1)); end end 计算方法（B）上机报告 6 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); 计算方法（B）上机报告 7 end M(n)=y1(n)/u(n); %追赶法求解样条参数M(i) for k=n-1:-1:1; M(k)=(y1(k)-c(k)*M(k+1))/u(k); end s=zeros(1,N); for m=1:N; k=1; for i=2:n-1 if X(m)<=x(i); k=i-1; break; else k=i; end end 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 计算方法（B）上机报告 8 %% 计算所需光缆长度 L=0; %计算所需光缆长度 for i=2:N L=L+sqrt((X(i)-X(i-1))^2+(s(i)-s(i-1))^2); end disp('所需光缆长度为 L='); disp(L); figure plot(x,y,'*',X,s,'-') %绘制铺设河底光缆的曲线图 xlabel('位置','fontsize',16); %标注坐标轴含义 ylabel('深度/m','fontsize',16); title('铺设河底光缆的曲线图','fontsize',16); grid; 4 结果与分析 铺设海底光缆的曲线图如图2-2所示：

图2-2 铺设河底光缆曲线图

合集下载

成本核算的会计上机报告

成本核算的会计上机报告

成本核算的会计上机报告在当今竞争激烈的商业环境中，准确的成本核算对于企业的决策制定、资源分配和盈利能力评估至关重要。

作为一名会计专业的学生，通过本次成本核算的会计上机实践，我对成本核算的理论知识有了更深入的理解，并掌握了相关软件的操作技能。

以下是我对本次实践的详细报告。

一、实践目的本次成本核算会计上机实践的主要目的是让我们熟悉成本核算的流程和方法，熟练掌握会计软件在成本核算中的应用，提高我们的实际操作能力和问题解决能力。

同时，通过实践加深对成本会计理论知识的理解，为今后从事相关工作打下坚实的基础。

二、实践内容（一）基础设置在开始成本核算之前，我们需要进行一系列的基础设置。

这包括建立账套、设置会计科目、定义成本核算对象、设置成本项目和费用项目等。

在建立账套时，我们需要根据企业的实际情况选择合适的会计制度和核算方法。

会计科目的设置要符合企业的业务特点和成本核算要求，确保能够准确地记录和反映成本信息。

（二）原始数据录入完成基础设置后，我们开始录入原始数据。

这包括材料采购数据、生产领料数据、人工费用数据、制造费用数据等。

在录入数据时，要确保数据的准确性和完整性，避免出现错误和遗漏。

同时，要注意数据的格式和单位，按照软件的要求进行录入。

（三）成本计算与分配录入原始数据后，我们使用会计软件进行成本计算和分配。

成本计算的方法包括品种法、分批法、分步法等，根据企业的生产特点和管理要求选择合适的方法。

在成本分配过程中，要合理确定分配标准，如材料消耗定额、工时定额等，确保成本在各产品之间、各部门之间进行准确分配。

（四）成本报表生成成本计算和分配完成后，我们利用会计软件生成各种成本报表，如成本计算单、成本汇总表、成本分析表等。

通过对成本报表的分析，我们可以了解企业的成本构成和成本变动情况，为成本控制和决策提供依据。

三、实践过程中遇到的问题及解决方法（一）数据录入错误在录入原始数据时，由于粗心大意，出现了数据录入错误的情况。

石大 计算方法实验报告1

石大 计算方法实验报告1

《计算方法》实验报告二级学院：计算机学院专业：指导教师：班级学号：姓名：实验一1、 实验目的：1) 熟悉用牛顿法解非线性方程的过程；熟悉用弦截法求解非线性方程的过程 2) 编程实现牛顿法、弦截法求非线性方程的根。

2、 实验要求:1) 认真分析题目的条件和要求，复习相关的理论知识，选择适当的解决方案和算法；2) 编写上机实验程序，作好上机前的准备工作；3) 上机调试程序，并试算各种方案，记录计算的结果（包括必要的中间结果）；4) 分析和解释计算结果； 5) 按照要求书写实验报告； 3、 实验内容：1)用牛顿法求解01553=-x 的根，取初始值为10。

2) 用弦截法求解数学方程。

010*15.110*4.181.9*002.0)(255.15=--=--x x x f 4、题目、原理、程序及结论 题目：非线性方程的迭代数值解1.牛顿法。

2.单点弦截法，双点弦截法。

原理： 牛顿法：设r 是方程f(x) = 0的根，选取x0作为r 初始近似值，过点（x0,f(x0））做曲线y = f(x ）的切线L ，L 的方程为y = f(x0)+f'(x0)(x-x0），求出L 与x 轴交点的横坐标 x1 = x0-f(x0)/f'(x0），称x1为r 的一次近似值。

过点（x1,f(x1））做曲线y = f(x ）的切线，并求该切线与x 轴交点的横坐标 x2 = x1-f(x1）/f'(x1），称x2为r 的二次近似值。

重复以上过程，得r 的近似值序列，其中x(n+1）=x(n ）－f(x(n))/f'(x(n)），称为r 的n+1次近似值，这就是牛顿迭代法的原理。

弦截法：设方程f(x) = 0的根在两个数x1，x2之间，连接这两点所对应的函数值，连线与x 轴的交点为新的x ，若f(x)与f(x1)同号，则把x 当作新的x1，将新的x1与x2连接，如此循环，如果f(x)与f(x1)异号，则把把x 当作新的x2，将x1与新的x2连接，如此循环，直到满足误差要求求出解，这就是弦截法的原理。

电力系统潮流上机计算报告

电力系统潮流上机计算报告

电力系统潮流上机计算报告
系别：电力工程系
班级：
姓名：
学号：
课号课序号：0110361-2
选课顺序号：22
程序说明
包括：程序设计思想、程序流程图、程序使用说明。

给定题目的手算过程（迭代两次）
包括：原题目、节点导纳矩阵、雅克比矩阵、第一次和第二次迭代结果。

给定题目的程序计算结果
包括：原题目、节点导纳矩阵、雅克比矩阵、程序输入和输出文件（误差0.0001）。

编程特色与创新
包括：程序能够完成的基本功能；程序能够完成的高级功能（如：是否包括平行支路、接地支路、非标准变比变压器支路，是否采用了稀疏矩阵技术，是否增加了人机对话界面，程序的通用性和实用性如何）。

五、总结
包括：手算结果与程序计算结果的分析比较；本次上机体会，如：独立编程体会、跟踪调试技能的掌握情况，C语言中结构体、指针、文件输入输出的掌握情况等。

报告要求：
报告必须手写（最好使用黑色水笔）。

报告统一采用A4打印纸书写（留出页边距: 1.5~2厘米）。

不使用实验报告纸。

封面按上述格式书写。

装订统一在左侧1厘米，三个钉。

上述五部分内容必须齐全，各部分内容可以扩充。

报告书写要求字迹清楚，不得潦草。

报告必须与本人提交程序吻合，否则取消成绩。

报告不得有雷同，否则全部取消成绩。

数学的上机实验报告

数学的上机实验报告

实验题目：线性代数求解方程组一、实验目的1. 理解线性代数中方程组的求解方法。

2. 掌握利用计算机求解线性方程组的算法。

3. 熟悉数学软件（如MATLAB、Python等）在数学问题中的应用。

二、实验内容本次实验主要利用数学软件求解线性方程组。

线性方程组是线性代数中的一个基本问题，其求解方法有很多种，如高斯消元法、矩阵求逆法等。

本实验以高斯消元法为例，利用MATLAB软件求解线性方程组。

三、实验步骤1. 编写高斯消元法算法程序。

2. 输入方程组的系数矩阵和常数项。

3. 调用程序求解方程组。

4. 输出解向量。

四、实验代码及分析1. 高斯消元法算法程序```matlabfunction x = gaussElimination(A, b)[n, m] = size(A);assert(n == m, 'The matrix A must be square.');assert(n == length(b), 'The length of b must be equal to the number of rows in A.');% 初始化解向量x = zeros(n, 1);% 高斯消元for i = 1:n-1% 寻找最大元素[~, maxIdx] = max(abs(A(i:n, i)));maxIdx = maxIdx + i - 1;% 交换行A([i, maxIdx], :) = A([maxIdx, i], :);b([i, maxIdx]) = b([maxIdx, i]);% 消元for j = i+1:nfactor = A(j, i) / A(i, i);A(j, i:n) = A(j, i:n) - factor A(i, i:n); b(j) = b(j) - factor b(i);endend% 回代求解for i = n:-1:1x(i) = (b(i) - A(i, i+1:n) x(i+1:n)) / A(i, i); endend```2. 输入方程组的系数矩阵和常数项```matlabA = [2, 1, -1; 1, 2, 1; -1, 1, 2];b = [8; 5; 2];```3. 调用程序求解方程组```matlabx = gaussElimination(A, b);```4. 输出解向量```matlabdisp('解向量为：');disp(x);```五、实验结果与分析实验结果：```解向量为：2-13```实验分析：通过高斯消元法，我们成功求解了给定的线性方程组。