上机实习2：超星数字图书馆

一、超星数字图书馆

我馆读者可以使用的超星资源库包括“超星中文电子图书”、“超星读秀学术搜索”、“超星发现系统”、“超星名师讲坛”。请通过上机实习了解四个资源库的资源收录情况和特色功能。———————————————————————————————————————1、请通过图书馆主页进入“超星数字图书馆”，超星数字图书馆是超星公司最早推出的产

品，收录上百万册的已获得版权许可的电子图书，读者均可以进行在线阅读，请浏览了解电子图书的收录情况。

2、请在超星数字图书馆中搜索覃彪喜著的《读大学，究竟读什么》一书，并尝试进行在线

阅读。

3、请通过图书馆主页进入“读秀学术搜索”数据库，它是超星公司旗下的产品，它的价值

在于为用户提供深入到图书章节和内容的全文检索。请您点击“知识”选项，并在输入框中输入“网络水军”，看看能找到哪些资料。注意在结果列别的右侧，按文献类型对检索结果进行了分类。您可以根据需要分类浏览。请在检索输入框中输入“案例”，然后点击“在结果中搜索”，看看系统能找到哪些关于网络水军的案例。

4、通过第3题，相信您已经体会到全文检索的好处，当查找概念、图表、注释、案例、公

式、数据、结论等隐含在文献中的信息片段时，采取全文检索是十分有效的方法。请您再次使用读秀的“知识”检索，查找一个您想了解的知识条目。

5、读秀的另一项服务是提供原文传递，针对图书馆没有购买全文的文献，读者可以申请“图

书文献传递服务”。请搜索《电子交易教程》，并尝试使用图书文献传递服务获取原文。

注意：该服务需要使用邮箱传递，请按照指引使用。

6、超星发现是超星公司的新产品，我们知道，如果我们要使用多个数据库的资源，要进行

多次检索，很费时。超星发现的价值在于整合了多个数据库系统，为读者提供统一检索平台。请使用超星发现系统，检索“网络水军”，然后点击任意一篇文献，系统将会自动向图书馆已经订购的数据库发出下载请求。

7、超星发现系统的一个重要特色是对检索结果提供可视化浏览。其中之一是发文量趋势

图，定量地反映研究主题的研究情况，此外还有其他视图，请查看“网络水军”检索结果的可视化视图。这些可视化视图都能从一个侧面反映出一个课题的研究现状，在做论文时，您可以根据需要在论文中使用，有助于您分析当前的研究现状。

8、超星名师讲坛是一个蕴藏了国内外众多知名学者多年的学术精华视频的资源库。请进入

超星名师讲坛，了解该资源库的收录情况。

9、尝试利用超星名师讲坛，查找你感兴趣的主题或名师的讲课视频。这些视频中，有单场

的专题讲座，有系统的课程讲座。如果想快餐式学习，请观看专题讲座，如果想系统深入地学习，请挑选课程讲座观看。例如，程焕文教授的“中国图书文化”就是一个有趣的书文化系列讲座。

大连理工程序设计 第二次上机作业

1 题目描述：4个数排序 从键盘任意输入4个整数，将其按从小到大顺序排序后输出。 题目分析：定义a,b,c,d,t.为a,b,c,d赋上四个值。用if语句，借用t，将a,b,c,d通过多次调换数值使得a,b,c,d对应的数值大小依次减小。再倒序输出。 流程图： 程序代码： #include main() { int a,b,c,d,t; printf("please enter four number"); scanf("%d%d%d%d",&a,&b,&c,&d); if(b>a) {t=a;a=b;b=t;} if(c>a) {t=a;a=c;c=t;} if(d>a) {t=a;a=d;d=t;} if(c>b) {t=b;b=c;c=t;} if(d>b) {t=b;b=d;d=t;} if(d>c) {t=c;c=d;d=t;} printf("After ascending sorting:%d %d %d %d\n",d,c,b,a); }

2. 题目描述：求1 + 1/2! +....+ 1/n! 题目分析：定义i，n，t，s。运用循环语句，让t的值为n！，同时s累计求和。因为最后要输出原式，所以要将公式的部分包含在循环语句中。 流程图：

#include main() { int i,n; float s,t; i=1; s=0; t=1; printf("please enter a number"); scanf("%d",&n); printf("1"); while(i<=n) {t=t*i; s+=1/t; if(i!=1) printf("+1/%d!",i); i++;} printf("=%.4f\n",s); } 程序截屏：

计算方法上机实验报告

. / 《计算方法》上机实验报告 班级：XXXXXX 小组成员：XXXXXXX XXXXXXX XXXXXXX XXXXXXX 任课教师：XXX 二〇一八年五月二十五日

前言 通过进行多次的上机实验，我们结合课本上的内容以及老师对我们的指导，能够较为熟练地掌握Newton 迭代法、Jacobi 迭代法、Gauss-Seidel 迭代法、Newton 插值法、Lagrange 插值法和Gauss 求积公式等六种算法的原理和使用方法，并参考课本例题进行了MATLAB 程序的编写。 以下为本次上机实验报告，按照实验内容共分为六部分。 实验一： 一、实验名称及题目： Newton 迭代法 例2.7(P38):应用Newton 迭代法求在附近的数 值解，并使其满足. 二、解题思路： 设'x 是0)(=x f 的根,选取0x 作为'x 初始近似值,过点())(,00x f x 做曲线)(x f y =的切线L ,L 的方程为))((')(000x x x f x f y -+=,求出L 与x 轴交

点的横坐标) (') (0001x f x f x x - =,称1x 为'x 的一次近似值,过点))(,(11x f x 做曲线)(x f y =的切线,求该切线与x 轴的横坐标) (') (1112x f x f x x - =称2x 为'x 的二次近似值,重复以上过程,得'x 的近似值序列{}n x ,把) (') (1n n n n x f x f x x - =+称为'x 的1+n 次近似值，这种求解方法就是牛顿迭代法。 三、Matlab 程序代码： function newton_iteration(x0,tol) syms z %定义自变量 format long %定义精度 f=z*z*z-z-1; f1=diff(f);%求导 y=subs(f,z,x0); y1=subs(f1,z,x0);%向函数中代值 x1=x0-y/y1; k=1; while abs(x1-x0)>=tol x0=x1; y=subs(f,z,x0); y1=subs(f1,z,x0); x1=x0-y/y1;k=k+1; end x=double(x1) K 四、运行结果：

计算方法上机实验报告

《计算方法》上机实验报告 班级：XXXXXX 小组成员：XXXXXXX XXXXXXX XXXXXXX XXXXXXX 任课教师：XXX 二〇一八年五月二十五日

前言 通过进行多次的上机实验，我们结合课本上的内容以及老师对我们的指导，能够较为熟练地掌握Newton 迭代法、Jacobi 迭代法、Gauss-Seidel 迭代法、Newton 插值法、Lagrange 插值法和Gauss 求积公式等六种算法的原理和使用方法，并参考课本例题进行了MATLAB 程序的编写。 以下为本次上机实验报告，按照实验内容共分为六部分。 实验一： 一、实验名称及题目： Newton 迭代法 例2.7(P38):应用Newton 迭代法求 在 附近的数值解 ，并使其满足 . 二、解题思路： 设'x 是0)(=x f 的根,选取0x 作为'x 初始近似值,过点())(,00x f x 做曲线)(x f y =的切线L ,L 的方程为))((')(000x x x f x f y -+=,求出L 与x 轴交点的横坐标) (') (0001x f x f x x - =,称1x 为'x 的一次近似值,过点))(,(11x f x 做曲线)(x f y =的切线,求该切线与x 轴的横坐标) (') (1112x f x f x x - =称2x 为'x

的二次近似值,重复以上过程,得'x 的近似值序列{}n x ,把 ) (') (1n n n n x f x f x x - =+称为'x 的1+n 次近似值，这种求解方法就是牛顿迭代法。 三、Matlab 程序代码： function newton_iteration(x0,tol) syms z %定义自变量 format long %定义精度 f=z*z*z-z-1; f1=diff(f);%求导 y=subs(f,z,x0); y1=subs(f1,z,x0);%向函数中代值 x1=x0-y/y1; k=1; while abs(x1-x0)>=tol x0=x1; y=subs(f,z,x0); y1=subs(f1,z,x0); x1=x0-y/y1;k=k+1; end x=double(x1) K 四、运行结果： 实验二：

linux上机实习报告

目录 Linux 第一次上机(运行环境RedHat) (1) 调试运行CD唱片应用程序 (1) 运行过程 (1) 心得体会 (3) 编写shell脚本，求1到100的和 (3) 运行过程 (3) 源代码 (4) 编写shell脚本,从键盘输入两个数，求这两个数的和 (4) 运行过程 (4) 源代码 (5) 等待特定用户登录，每30秒确认一次 (5) 运行过程 (5) 源代码 (6) 找出系统中当前消耗磁盘容量最大的前10个用户，并向他们发送警告邮件 (6) 运行过程 (6) 源代码 (8) 查找输入文件的路径 (8) 运行过程 (8) 源代码 (9) Linux 第二次上机(运行环境RedHat) (9) 定制自己的ls命令。提供至少三种带参数的执行方式 (9) 运行过程 (9) 源代码 (11) 调试编译串行口通信程序p6.5.c (16) 运行过程 (16) Linux 第三次上机(运行环境RedHat，CentOS) (16) 创建一个系统监听守护进程，一旦接收到其他进程发来的信号，马上给出报告 (16) 运行过程 (16) 源代码 (18) 分别利用本地socket套接字和INTENET套接字实现进程间文件传输 (20) 本地socket套接字运行过程 (20) 本地socket套接字服务器server1.c源代码 (22) 本地socket套接字客户端client1.c源代码 (23) INTENET套接字运行过程 (25) INTENET套接字服务器server2.c源代码 (28) INTENET套接字客户端client2.c源代码 (29) 感谢与收获 (30)

文献检索第二次上机实习作业 (2)

1．检索有关钒钛中锰含量的测定的国家标准，给出标准号、起草单位及检索流程。 标准号YB/T 505.6-2007 起草单位：冶金工业信息标准研究院攀研院 以钒钛中锰含量的测定和国家标准为收索主题进行百度收索 2.在cnki的《中国期刊全文数据库》中检索篇名为《青梅果酒发酵工艺优化》的论文。下载并浏览全文。在阅读该篇文章的过程中有些地方该文章是引用其他文章的内容，记下该篇论文的第三篇参考文献。到维普中文科技期刊数据库中检索你记下的该篇参考文献，浏览该篇文章的全文，给出该文章作者的单位。并说明这种查找文献的方法是什么方法。同样的内容在万方数字化资源库中在进行一次检索，比较二者有何不同。 第三个参考文献《果梅资源的开发和利用》作者：曾凡骏张月天 作者单位四川大学轻工与食品工程学院: 查找文献方法：找出该文献，然后再上知网搜索。万方的文献带有超级链接，搜索更加方便。 3.检索一篇华南理工大学的硕士学位论文，题名为《具有修护作用的洗发和护发产品的配方研究》，说明可以到哪些数据库中检索，其中哪些数据库可以下载全文。写出该论文作者是谁。 中国知网中可以下载全文，作者：王新全 4.到中国重要报纸全文数据库中检索发表在人民日报上的题名含有攀枝花的文章。写出检索策略及结果，下载一篇caj格式的全文，并说明文章共有多少页。其实网上也可以免费检索和下载该文章的全文，写出检索步骤。 进中国知网，点击‘报纸栏’再关键词栏输入‘攀枝花’。文章：<< 目光里,始终有攀枝花传奇>> 5.请在维普信息资源系统中检索2008年涉及“数字图书馆”的相关文章。（1）全文中含有“数字图书馆”的文章，共有__1145_____篇。 （2）标题中含有“数字图书馆”，有_893_篇。 6.请在万方数据库中的“学术会议”数据库中检索涉及“室内空气质量”的会议文献，请摘录其中一篇. 常州城市空气质量变化特征及其与气象条件的关系 利用常州环境监测中心站空气质量资料和常规气象观测资料，对2001-2006年常州空气质量的变化特点进行了统计分析，结果表明：(1)常州最严重污染年份为2001、2002年；空气污染呈显著下降趋势。(2)污染最高发季节为冬季，其次为春季，再次为秋季，最少发生的季节为夏季。(3)污染出现最多的月份是1、3、10、11月，而6、7、8、9月则是污染出现较少的月份。(4)雾会加重空气污染程度，浓雾或强浓雾加重程度更深。由于常州空气质量呈逐渐好转的趋势，雾对API指数大小的贡献也呈下降趋势。(5)强降水会显著改善空气质量，

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

计算方法上机报告

姓名： 学号： 班级：能动上课班级：

题目及求解： 一、对以下和式计算： ∑ ∞ ? ?? ??+-+-+-+=0681581482184161n n n n S n ，要求： ① 若只需保留11个有效数字，该如何进行计算； ② 若要保留30个有效数字，则又将如何进行计算； 1 算法思想 （1）根据精度要求估计所加的项数，可以使用后验误差估计，通项为： 1421114 16818485861681 n n n a n n n n n ε??= ---<< ?+++++??； （2）为了保证计算结果的准确性，写程序时，从后向前计算； （3）使用Matlab 时，可以使用以下函数控制位数： digits(位数)或vpa(变量，精度为数) 2 算法结构 ;0=s ?? ? ??+-+-+-+= 681581482184161n n n n t n ； for 0,1,2,,n i =??? if 10m t -≤ end; for ,1,2,,0n i i i =--??? ;s s t =+ 3 Matlab 源程序 clear; %清除工作空间变量 clc; %清除命令窗口命令 m=input('请输入有效数字的位数m='); %输入有效数字的位数 s=0;

for 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; 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.8。 通过上面的实验结果可以看出，通过从后往前计算，这种算法很好的保证了计算结果要求保留的准确数字位数的要求。 二、某通信公司在一次施工中，需要在水面宽度为20米的河沟底部沿直线走向铺设一条沟底光缆。在铺设光缆之前需要对沟底的地形进行初步探测，从而估计所需光缆的长度，为工程预算提供依据。已探测到一组等分点位置的深度数据(单位:米)如下表所示：

华中科技大学c第2次上机作业

第二章上机实验 一、实验目的 ●学会编写选择结构的C++程序。 ●掌握if 语句的使用 ●掌握switch 语句的使用 ●学会编写简单的循环结构C++程序 二、实验内容 1．在vc环境下编写下述程序代码，编译调试其语法和逻辑错误，直至得到完全正确的结果。 #include using namespace std; int main(){ int a,b,c; int max,min; cout<<"请输入三个不等整数:\n"; cin>>a>>b>>c; if(a>b) max=a; min=b; else ; max=b; min=a; if(c>max) cout<

发系统给出很多条error信息，因此，发现一个错误并修改后最好重新编译一次，以便提高工作效率。 警告错误处理： 一般是触发了C/C++的自动规则，如将一个单精度（浮点）型数据赋给整型变量，需要系统将单精度型数据自动转换为整型，此时小数部分会丢失，因而系统给出警告信息，警告信息不会影响程序执行，本例可以通过强制转换去掉警告信息。 连接错误（linking error） 产生连接错误的原因可能是缺少所需的库文件或目标文件，或程序中调用的外部函数没有定义等，只要补充相应文档再重新建立即可。 逻辑错误 运行程序时，可能发现程序没有编译错误，而且也能执行，但执行结果不对，此时存在的就是逻辑错误。 ①编译上述程序，按照上面介绍的方法找出并修改其语法错误。 ②语法错误修改完毕后用下列四组数组分别测试其结果，看是否正确，如果不正确，请 使用下述调试手段修改。 数据一： 33 60 50 数据二： 33 22 10 数据三： 88 45 67 数据三： 88 45 30 调试过程：+ 首先，在源程序中可能出现错误的行上设置断点，方法是将光标移至该行，然后按F9键，或单击工具栏上的手形按钮（再按一次F9键或单击手形按钮将取消断点），此时该行左侧出现一个红色圆点，断点设置成功，如图就是在A行设置了断点。 然后选择【组建】|【开始调试】|【GO】菜单命令（也可直接单击工具栏上的图标），程序开始执行，但执行到断点处停止，这时选择【查看】|调试窗口|子菜单的Watch 和Variables两个菜单项，打开监视窗口和变量窗口来观察变量的值，分析并查找出错原因。 在Watch窗口加入max和min两个变量，进行监视。Watch窗口的每一行可以显示一个变量，其中左栏显示变量名，双击它可进行编辑；右栏显示变量的值。接下来可按F10键从断点位置处单步执行。通过单步执行（一条一条语句执行）可以看出我们所观察的变量的变化以及程序执行流程是否正确，如果不正确是由哪条语句引起的，因而有针对性的检查错误原因。例如在本例中，尽管a>b,在执行了if后面的｛max=a; min=b; ｝后，通过单步执行发现仍然执行了else后面的｛max=b; min=a; ｝所以当程序执行到箭头所指位置时，max=45, min=48,如图所示，与预期结果不相符，说明程序的流程有问题。此时再仔细分析源程序，发现问题在于else后多余的分号。

计算方法第二章方程求根上机报告

实验报告名称 班级：学号：姓名：成绩： 1实验目的 1）通过对二分法与牛顿迭代法作编程练习与上级运算，进一步体会二分法与牛顿迭代法的不同特点。 2）编写割线迭代法的程序，求非线性迭代法的解，并与牛顿迭代法。 2 实验内容 用牛顿法和割线法求下列方程的根 x^2-e^x=0; x*e^x-1=0; lgx+x-2=0; 3实验步骤 1)根据二分法和牛顿迭代法，割线法的算法编写相应的求根函数； 2)将题中所给参数带入二分法函数，确定大致区间； 3)用牛顿迭代法和割线法分别对方程进行求解； 3 程序设计 牛顿迭代法x0=1.0; N=100; k=0; eps=5e-6; delta=1e-6; while(1) x1=x0-fc1(x0)/fc2(x0); k=k+1; if k>N disp('Newmethod failed')

break end if(abs(x1-x0)=delta) c=x1; x1=cutnext(x0,x1); x0=c; %x0 x1μYí?μ?μ?x1 x2 è?è?±￡′??úx0 x1 end k=k+1; if k>N disp('Cutline method failed') break; end if(abs(x1-x0)

计算方法第三次上机实习报告

实验报告 课程名称： 计算方法 指导老师： 太英 成绩： 实验名称： 第三次上机作业 实验类型： matlab 同组学生： 一、实验目的和要求（必填） 二、实验容和原理（必填） 三、主要仪器设备（必填） 四、操作方法和实验步骤 五、实验数据记录和处理 六、实验结果与分析（必填） 七、讨论、心得 一、实验目的 用龙贝格算法计算积分I =∫ssss s ss 01 ，要求误差不超过ε=12 ×105 二、实验原理 龙贝格算法是由递推算法得来的。由梯形公式得出辛普森公式得出柯特斯公式最后得到龙贝格公式。设将求积区间[a ，b]分为n 个等分，则一共有n+1个等分点,k x a kh =+,0,1,b a h k n -= =，n 。这里用n T 表示复化梯形法求得的积分 值，其下标n 表示等分数。 先考察下一个字段[1,k k x x +]，其中点()112 1 2 k k k x x x ++= +，在该子段上二分前后两个积分值 ()()112 k k h T f x f x += +???? ()()21124k k k h T f x f x f x ++? ??? =+ +?? ??????? 显然有下列关系 2112122k h T T f x +? ?=+ ??? 将关系式关于k 从0到n-1累加求和，即可得递推公式1 210 2122n n n k k h T T f x -+=?? =+ ??? ∑

需要强调指出的是，上式中的b a h n -= 代表二分前的步长，而12 12k x a k h +? ?=++ ??? 根据梯形法的误差公式，积分值n T 的截断误差大致与2h 成正比，因此步长减半后误差将减至四分之一，即有 211 14 n n T T -≈- 将上式移项整理，知 221 1()3 n n n T T T -≈- 按上式，积分值2n T 的误差大致等于21 ()3 n n T T -，如果用这个误差值作为2n T 的 一种补偿，可以期望，所得的()222141 333 n n n n n T T T T T T =+ -=-应当是更好的结果。 组合得到的近似值T 直接验证，用梯形二分前后的两个积分值n T 和2n T 按式组合， 结果得到辛普森法的积分值241 33 n n n S T T =- 再考察辛普森法。其截断误差与4h 成正比。因此，若将步长折半，则误差相应的减至十六分之一。既有 21 16n n I S I S -≈- 由此得 21611515 n n I S S ≈ - 不难验证，上式右端的值其实就等于n C ，就是说，用辛普森法二分前后的两个积分值n S 和2n S ，在按上式再做线性组合，结果得到柯特斯法的积分值n C ，既有 2161 1515n n n C S S ≈ - 重复同样的手续，依据斯科特法的误差公式可进一步导出龙贝格公式 2641 6363 n n n R C C ≈ - 在步长二分的过程中运用公式加工三次，就能将粗糙的积分值n T 逐步加工成精度较高的龙贝格n R ，或者说，将收敛缓慢的梯形值序列n T 加工成熟练迅速的龙贝 格 值序列n R ，这种加速方法称龙贝格算法。 三、实验过程 1.流程图

第二次上机作业参考答案

第二次上机作业参考答案： 1． 4个数排序 从键盘任意输入4个整数，将其按从小到大顺序排序后输出。 例：输入 5 12 -9 3 输出After ascending sorting: -9 3 5 12 #include void main() { int a,b,c,d,t; scanf("%d%d%d%d",&a,&b,&c,&d); if(a>b) { t=a; a=b; b=t; } if(a>c) { t=a; a=c; c=t; } if(a>d) { t=a; a=d; d=t; } if(b>c) { t=b; b=c; c=t; } if(b>d) { t=b; b=d; d=t; } if(c>d) { t=c; c=d; d=t; } printf("The sorted number is:%d %d %d %d\n",a,b,c,d); } 2．求1 + 1/2! +....+ 1/n! 输入正整数n，计算上式前n 项的和(保留 4 位小数)。 例：输入10 输出 1 + 1/2! +....+ 1/10!＝1.7183 #include int main( ) { int i,j,n; float s,t; scanf("%d",&n); for(s=0,i=1;i<=n;i++) { for(t=1,j=1;j<=i;j++) t*=j; s+=1/t; } printf("%0.4f\n",s); }

3. 编程计算表达式：data1 op data2 的值。（+ - * / ％） 例：输入23+12 输出23＋12＝35 #include int main( ) { int a,b,result,flag=1; char op; scanf("%d%c%d",&a,&op,&b); switch(op) { case '+': result=a+b; break; case '-': result=a-b; break; case '*': result=a*b; break; case '/': printf("%d%c%d=%.2f\n",a,op,b,(float)a/b); flag=0; break; case '%': result=a%b; break; } if(flag) printf("%d%c%d=%d\n",a,op,b,result); } 4. 用二分法求方程2x3－4x2+3x-6=0在（a，b）之间的根。 例：输入Please input the lower and upper boundaries: 3 10 输出No root in this boundary! Please input the lower and upper boundaries: 1 5 The root is 2.00 #include void main() { double x,fx,a,fa,b,fb; do { printf("Please input the lower and upper boundaries:"); scanf("%lf%lf",&a,&b); fa=2*a*a*a-4*a*a+3*a-6; fb=2*b*b*b-4*b*b+3*b-6; }while(fa*fb>0); while(1) { x=(a+b)/2; fx=2*x*x*x-4*x*x+3*x-6;

数值分析上机实验报告

数值分析上机实验报告

《数值分析》上机实验报告 1.用Newton 法求方程 X 7-X 4+14=0 在（0.1，1.9）中的近似根（初始近似值取为区间端点，迭代6次或误差小于0.00001）。 1.1 理论依据： 设函数在有限区间[a ，b]上二阶导数存在，且满足条件 {}α?上的惟一解在区间平方收敛于方程所生的迭代序列 迭代过程由则对任意初始近似值达到的一个中使是其中上不变号 在区间],[0)(3,2,1,0,) (') ()(],,[x |))(),((|,|,)(||)(|.4;0)(.3],[)(.20 )()(.110......b a x f x k x f x f x x x Newton b a b f a f mir b a c x f a b c f x f b a x f b f x f k k k k k k ==- ==∈≤-≠>+ 令 )9.1()9.1(0)8(4233642)(0)16(71127)(0)9.1(,0)1.0(,1428)(3 2 2 5 333647>?''<-=-=''<-=-='<>+-=f f x x x x x f x x x x x f f f x x x f 故以1.9为起点 ?? ?? ? ='- =+9.1)()(01x x f x f x x k k k k 如此一次一次的迭代，逼近x 的真实根。当前后两个的差<=ε时，就认为求出了近似的根。本程序用Newton 法求代数方程(最高次数不大于10)在（a,b ）区间的根。

1.2 C语言程序原代码： #include #include main() {double x2,f,f1; double x1=1.9; //取初值为1.9 do {x2=x1; f=pow(x2,7)-28*pow(x2,4)+14; f1=7*pow(x2,6)-4*28*pow(x2,3); x1=x2-f/f1;} while(fabs(x1-x2)>=0.00001||x1<0.1); //限制循环次数printf("计算结果：x=%f\n",x1);} 1.3 运行结果： 1.4 MATLAB上机程序 function y=Newton(f,df,x0,eps,M) d=0; for k=1:M if feval(df,x0)==0 d=2;break else x1=x0-feval(f,x0)/feval(df,x0); end e=abs(x1-x0); x0=x1; if e<=eps&&abs(feval(f,x1))<=eps d=1;break end end

文献检索第二次上机作业

第二次计算机检索实习题目 本次实习目的是使学生对以下三个外文数据库：美国EI COMPENDEX数据库、荷兰SDOS 数据库、德国Springer数据库的使用方法有所了解。通过题目，使学生分别掌握每个数据库所提供的各种获取文献的方法，包括浏览和检索的方法，并通过各种算符的组配，以达到获取正确检索结果的目的。 ”） 一、EI COMPENDEX 数据库检索题目 题目1检索2000以前的有关“计算机网络”方面的文献，检索要求： （1）文献类型为会议，文献处理类型为理论，语种为英文。 （2）记下检索到的数量及第一条记录的篇名、作者、刊名、年、卷、期，出版地、出版社等。 （3 检索到的数量：31931； 第一条记录的篇名：Fault-tolerant analysis of the spiral computer network topology； 作者：Mitchell, Tony L.、Nilsson, Arne A.； 刊名：Conference Record - International Conference on Communications； 年：1989；卷：1；期：1989 出版地：Boston, MA, USA； 出版社：Publ by IEEE, Piscataway, NJ, United States； 该篇文章三种显示格式： 第一种题录格式：Fault-tolerant analysis of the spiral computer network topology Mitchell, Tony L. (USAF Acad, Dep of Math Sci, Colorado, Springs, CO, USA); Nilsson, Arne A. Source: Conference Record - International Conference on Communications, v 1, p 363-367, 1989

太原理工大学数值计算方法实验报告

本科实验报告 课程名称：计算机数值方法 实验项目：方程求根、线性方程组的直接解 法、线性方程组的迭代解法、代数插值和最 小二乘拟合多项式 实验地点：行勉楼 专业班级： ******** 学号： ********* 学生姓名： ******** 指导教师：李誌，崔冬华 2016年 4 月 8 日

y = x*x*x + 4 * x*x - 10; return y; } float Calculate(float a,float b) { c = (a + b) / 2; n++; if (GetY(c) == 0 || ((b - a) / 2) < 0.000005) { cout << c <<"为方程的解"<< endl; return 0; } if (GetY(a)*GetY(c) < 0) { return Calculate(a,c); } if (GetY(c)*GetY(b)< 0) { return Calculate(c,b); } } }; int main() { cout << "方程组为：f(x)=x^3+4x^2-10=0" << endl; float a, b; Text text; text.Getab(); a = text.a; b = text.b; text.Calculate(a, b); return 0; } 2.割线法： // 方程求根（割线法）.cpp : 定义控制台应用程序的入口点。// #include "stdafx.h" #include"iostream"

心得体会 使用不同的方法，可以不同程度的求得方程的解，通过二分法计算的程序实现更加了解二分法的特点，二分法过程简单，程序容易实现，但该方法收敛比较慢一般用于求根的初始近似值，不同的方法速度不同。面对一个复杂的问题，要学会简化处理步骤，分步骤一点一点的循序处理，只有这样，才能高效的解决一个复杂问题。

计算方法上机实习题大作业(实验报告).

计算方法实验报告 班级： 学号： 姓名： 成绩： 1 舍入误差及稳定性 一、实验目的 （1）通过上机编程，复习巩固以前所学程序设计语言及上机操作指令； （2）通过上机计算，了解舍入误差所引起的数值不稳定性 二、实验内容 1、用两种不同的顺序计算10000 21n n -=∑，分析其误差的变化 2、已知连分数() 1 01223//(.../)n n a f b b a b a a b =+ +++，利用下面的算法计算f ： 1 1 ,i n n i i i a d b d b d ++==+ (1,2,...,0 i n n =-- 0f d = 写一程序，读入011,,,...,,,...,,n n n b b b a a 计算并打印f 3、给出一个有效的算法和一个无效的算法计算积分 1 041 n n x y dx x =+? (0,1,...,1 n = 4、设2 2 11N N j S j == -∑ ，已知其精确值为1311221N N ?? -- ?+?? （1）编制按从大到小的顺序计算N S 的程序 （2）编制按从小到大的顺序计算N S 的程序 （3）按两种顺序分别计算10001000030000,,,S S S 并指出有效位数 三、实验步骤、程序设计、实验结果及分析 1、用两种不同的顺序计算10000 2 1n n -=∑，分析其误差的变化 （1）实验步骤： 分别从1~10000和从10000~1两种顺序进行计算，应包含的头文件有stdio.h 和math.h （2）程序设计： a.顺序计算

#include #include void main() { double sum=0; int n=1; while(1) { sum=sum+(1/pow(n,2)); if(n%1000==0)printf("sun[%d]=%-30f",n,sum); if(n>=10000)break; n++; } printf("sum[%d]=%f\n",n,sum); } b.逆序计算 #include #include void main() { double sum=0; int n=10000; while(1) { sum=sum+(1/pow(n,2)); if(n%1000==0) printf("sum[%d]=%-30f",n,sum); if(n<=1)break; n--; } printf("sum[%d]=%f\n",n,sum); } （3）实验结果及分析： 程序运行结果： a.顺序计算

CAD上机实习实验报告.

CAD上机实习实验报告 班级：****1101 姓名：王超 学号：******** 2014年6月

一、作业题目 利用CAD做下图： 要求：全部作图过程截图10张，每张图片配有简要的操作步骤说明。提交上机报告。 二、步骤 1.第一步 绘制中心线、中心线上侧大体轮廓。 双击： 打开CAD2010。单击“文件”——选择“新建”— —新建acadiso.dwt文件。然后在工作面板中进行图 形的绘制。

A.在工具栏，单击“直线”工具，画出 齿轮结构正视图和左视图的中心线，单击工具栏 ，中心线的线型设定为ACAD_ISO08W100，。 B.单击“直线”工具，做齿轮中心线上部构件图。做 一高度61.5mm，宽度30mm的矩形框。在工具栏 ，选择“偏移”工具→输入偏移距离 1.5mm→空格键确定→选择作齿轮最顶部直线→单击偏移侧区域，按照该步骤（下面不在赘述），依次向下偏移1.5mm。 C.第一步完成。

第一步完成后，绘制所得图形如下： 图1 2.第二步 修正齿轮倒角，完成中心孔的绘制。 A.在快捷工具栏选择“倒角”，在底部命令栏中 输入“A”,设置倒角长度为3mm,倒角角度45o,选 中齿轮左外轮廓线和最顶部直线，做出45o倒角， 倒角被顶部轮廓线下方的直线穿过，被切成两个倒 角。选择“修剪”工具，修剪多余线头。 B.选择“直线”，作齿轮中心孔。中心孔半径14mm, 选 “偏移”工具，作直线和中心线平行，偏离 14mm, 选择“倒角”,在底部命令栏中输入“A”

设置倒角长度为1mm,倒角角度45o。沿水平线和倒 角线的交点作垂线。 C.作主视图键槽。键槽深2.4mm , 选“偏移”工具 ，作中心线的平行线，向上偏离16.4mm 。 D.选择“修剪”工具，清除多余线头。 第二步完成后，绘制所得图形如下： 图2 3.第三步 完成齿轮构件板上沟槽和通孔的绘制，修正倒角。 A.作齿轮通孔（对称线左侧部分，镜像得到完整通

C语言上机作业试题5套含答案

第一次上机作业（2014）： 要求： ?独立完成，上机调试通过后上交 ?作业提交邮箱： ?邮件主题命名：班级学号姓名（第1次作业）， 例如：电1301班2王鹏（第1次作业） ?以附件形式将每次作业的所有程序源程序压缩打包后提交，压缩文件命名格式同上。 ?程序必须采用缩进格式 1．大写字母转换成小写字母 从键盘输入一个大写英文字母，输出相应的小写字母。 例：输入G 输出g 2．求平方根 输入一个实数x，计算并输出其平方根（保留1 位小数）。 例：输入17 输出The square root of 17.0 is 4.1 3．温度转换 设计一个程序将华氏温度转换成摄氏温度c = 5/9(f-32) a)输入华氏温度（实型） b)输出的摄氏温度结果保留两位小数 例：Please input Fahrenheit temperature: 76.8 The corresponding Celsius temperature is 24.89 4. 计算旅途时间 输入二个整数time1 和time2，表示火车的出发时间和到达时间，计算并输出旅途时间。 (输入的有效的时间范围是0000 到2359，不需要考虑出发时间晚于到达时间的情况。) 例：输入712 1411 （出发时间是7：12，到达时间是14：11） 输出The train journey time is 6 hours 59 minutes. 5. 数字加密 输入一个四位数，将其加密后输出。方法是将该数每一位上的数字加9，然后除以10 取余，做为该位上的新数字，最后将第1 位和第3 位上的数字互换，第2 位和第4 位上的数字互换，组成加密后的新数。 例：输入1257 输出The encrypted number is 4601 思考题：你能否编程找出谁做的好事？ 有四位同学中的一位做了好事，不留名，表扬信来了之后，校长问这四位是谁做的好事。 ?A说：不是我。 ?B说：是C。 ?C说：是D。 ?D说：他胡说。 已知三个人说的是真话，一个人说的是假话。现在要根据这些信息，找出做了好事的人。 第二次上机作业（2014）： 要求： ?独立完成，程序必须采用缩进格式，上机调试通过后上交

计算方法实验报告格式

计算方法实验报告格式 小组名称: 组长姓名(班号)： 小组成员姓名(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一个完整的实验，应包括数据准备、理论基础、实验内容及方法，最终对实验结果进行分析，以达到对理论知识的感性认识，进一步加深对相关算法的理解，数值实验以实验报告形式完成，实验报告格式如下： 一、实验名称 实验者可根据报告形式需要适当写出． 二、实验目的及要求 首先要求做实验者明确，为什么要做某个实验，实验目的是什么，做完该实验应达到什么结果，在实验过程中的注意事项，实验方法对结果的影响也可以以实验目的的形式列出． 三、算法描述（实验原理与基础理论） 数值实验本身就是为了加深对基础理论及方法的理解而设置的，所以要求将实验涉及到的理论基础，算法原理详尽列出． 四、实验内容 实验内容主要包括实验的实施方案、步骤、实验数据准备、实验的算法以及可能用到的仪器设备． 五、程序流程图 画出程序实现过程的流程图，以便更好的对程序执行的过程有清楚的认识，在程序调试过程中更容易发现问题． 六、实验结果 实验结果应包括实验的原始数据、中间结果及实验的最终结果，复杂的结果可以用表格

形式列出，较为简单的结果可以与实验结果分析合并出现． 七、实验结果分析 实验结果分析包括对对算法的理解与分析、改进与建议． 数值实验报告范例 为了更好地做好数值实验并写出规范的数值实验报告，下面给出一简单范例供读者参考． 数值实验报告 小组名称: 小组成员(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一、实验名称 误差传播与算法稳定性． 二、实验目的 1．理解数值计算稳定性的概念． 2．了解数值计算方法的必要性． 3．体会数值计算的收敛性与收敛速度． 三、实验内容 计算dx x x I n n ? += 1 10 ，1,2,,10n = ． 四、算法描述 由 dx x x I n n ? += 1 10 ，知 dx x x I n n ?+=--101110，则

计算方法B上机报告

计算方法B 上机报告 第1题 某通信公司在一次施工中，需要在水面宽度为20米的河沟底部沿直线走向铺设一条沟底光缆。在铺设光缆之前需要对沟底的地形进行初步探测，从而估计所需光缆的长度，为工程预算提供依据。已探测到一组等分点位置的深度数据(单位:米)如下表所示： (1)请用合适的曲线拟合所测数据点； (2)估算所需光缆长度的近似值，并作出铺设河底光缆的曲线图； 问题分析和算法思想： 本题的主要目的是对21个测量数据进行拟合，同时对拟合曲线进行线积分即可得到河底光缆长度的近似值，可以用的插值方法很多：多项式插值、Lagrange 插值、Newton 插值、三次样条插值等。由于数值点较多时，采用高次多项式插值将产生很大的误差，用拉格朗日插值多项式会出现龙格现象。故为了将所有的数据点都用上，且题中光缆为柔性，可光滑铺设于水底，鉴于此特性，采用三次样条插值的方法较为合适。 计算光缆长度近似值，只需将每两点之间的距离算出，然后依次相加，所得的折线长度，即为光缆长度的近似值。 光缆长度计算公式： 19 1 k k k l +===∑? ? ? 算法结构： 三次样条算法结构见《计算方法教程》P110。 源程序： clear;clc; x=0: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]; d=y; plot(x,y,'k.','markersize',15) hold on %%%计算二阶差商 for k=1:2 for i=21:-1:(k+1) d(i)=(d(i)-d(i-1))/(x(i)-x(i-k)); end end %%%假定d的边界条件，采用自然三次样条 for i=2:20 d(i)=6*d(i+1); end d(1)=0; d(21)=0; %%%追赶法求解带状矩阵的m值 a=0.5*ones(1,21); b=2*ones(1,21); c=0.5*ones(1,21); a(1)=0;c(21)=0; u=ones(1,21); u(1)=b(1); r=c; yy(1)=d(1); %%%追的过程 for k=2:21 l(k)=a(k)/u(k-1); u(k)=b(k)-l(k)*r(k-1); yy(k)=d(k)-l(k)*yy(k-1); end %%%赶的过程 m(21)=yy(21)/u(21); for k=20:-1:1 m(k)=(yy(k)-r(k)*m(k+1))/u(k); end %%%利用插值点画出拟合曲线 k=1; nn=100; xx=linspace(0,20,nn); l=0; for j=1:nn for i=2:20 if xx(j)<=x(i) k=i;