数值计算课程设计
数值计算课程设计报告
系别信息与计算科学
专业数学与应用数学
学号
姓名
指导教师
成绩
教师评语：
指导教师签字：
2011年7月13日
1绪论
数值分析是研究分析用计算机求解数学计算问题的数值计算方法及其理论的学科，是数学的一个分支，它以数字计算机求解数学问题的理论和方法为研究对象。
MATLAB是主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案。
2用MATLAB进行插值计算
题目
下列数据点的插值
(1)用这9个点作8次多项式插值8()
L x.
(2)用三次样条（第一边界条件）程序求()
S x.
从得到的结果看在[0, 64]上，哪个插值更精确；在区间[0, 1]上，两种插值哪个更精确？
2.1 程序文件
建立新的m文件unit2.m：
代码如下：
clc; clear;
x = [0, 1, 4, 9, 16, 25, 36, 49, 64];
y = [0:8];
xi = [0:64];
p = polyfit(x, y, 8);
yi = polyval(p, xi);
subplot(2, 2, 1);
plot(x, y, 'x', xi, yi, 'k')
xlabel('x轴'); ylabel('y轴'); legend('插值点', '8次差值多项式') yi = interp1(x, y, xi, 'cubic'); subplot(2, 2, 2);
plot(x, y, 'x', xi, yi, 'k')
xlabel('x轴'); ylabel('y轴'); legend('插值点', '三次样条插值') xi = [0:0.01:1];
p = polyfit(x, y, 8);
yi = polyval(p, xi);
subplot(2, 2, 3);
plot(xi, yi, 'k')
xlabel('x轴'); ylabel('y轴'); legend('8次差值多项式')
yi = interp1(x, y, xi, 'cubic'); subplot(2, 2, 4);
plot(xi, yi, 'k')
xlabel('x轴'); ylabel('y轴'); legend('三次样条插值')
图1
2.2 图样
运行文件得到图1
2.3 分析
比较两个函数的图像可知，在区间[0, 64]上三次样条插值函数要更精确一些，在区间[0, 1]上拉格朗日插值函数仍然不如三次样条插值函数更精确。
3多项式曲线拟合
题目
由实验给出数据表
3.1程序文件
建立新的m文件unit3.m：
代码如下：
clc; clear;
x = [0.0, 0.1, 0.2, 0.3, 0.5, 0.8, 1.0];
y = [1.0, 0.41, 0.50, 0.61, 0.91, 2.02, 2.46];
a = polyfit(x, y, 3);
b = polyfit(x, y, 4);
xi = [0.0:0.01:1.0];
aa = polyval(a, xi);
bb = polyval(b, xi);
subplot(1, 2, 1);
plot(x, y, 'x', xi, aa, 'k')
xlabel('x轴'); ylabel('y轴');
legend('插值点', '3次曲线拟合')
subplot(1, 2, 2);
plot(x, y, 'x', xi, bb, 'k')
xlabel('x轴'); ylabel('y轴');
legend('插值点', '4次曲线拟合')
poly2str(a, 'x')
poly2str(b, 'x')
3.2 图形和结果
运行程序结果如下：
ans = - 6.6221 x^3 + 12.8147 x^2 - 4.6591 x + 0.92659
ans = 2.8853 x^4 - 12.3348 x^3 + 16.2747 x^2 - 5.2987 x + 0.94272 图像如下：
图2
3.2 5次多项式曲线拟合
下面进行的：
代码如下：
clc; clear;
x = [0.0, 0.1, 0.2, 0.3, 0.5, 0.8, 1.0];
y = [1.0, 0.41, 0.50, 0.61, 0.91, 2.02, 2.46];
a = polyfit(x, y, 5);
xi = [0.0:0.01:1.0];
aa = polyval(a, xi);
plot(x, y, 'x', xi, aa, 'k')
xlabel('x轴'); ylabel('y轴');
legend('插值点', '5次曲线拟合')
poly2str(a, 'x')运行程序结果如下：
ans = - 79.3261 x^5 + 195.4554 x^4 - 172.7104 x^3 + 69.0498 x^2 - 11.0044 x + 0.99547 图像如下：
图3
由图像可知，次数越高，拟合越准确。
4 数值积分
题目
用不同数值方法计算积分
94ln 1
-
=?
xdx x .
(1)取不同的步长h.分别用复合梯形及复合辛普森求积计算积分，给出误差中关于h 的函数，并与积分精确值比较两个公式的精度，是否存在一个最小的h ，使得精度不能呢个再被改善？
(2)用龙贝格求积计算完成问题(1). (3)用自适应辛普森积分，使其精度达到4
10- 4.1 复合梯形公式代码
建立新的m 文件unit41.m ： 代码如下： clc, clear;
x = 0.0000001:0.1:1;
y = x.^(1 / 2). * log(x);
Fx = trapz(x, y)
error = Fx + 4 / 9
运行文件结果为：
Fx = - 0.4123
error = 0.0321
取不同步长h:
h = 0.01，Fx = - 0.4431，error = - 0.0014；
h = 0.001，Fx = - 0.4444，error = - 5.4875e - 005；
h = 0.0001，Fx = - 0.4444，error = - 2.0338e - 006；
h = 0.00001，Fx = - 0.4444，error = - 6.4496e - 008.
可见，没有一个最小的h使精度不能再改变。
4.2 复合辛普森公式代码
建立新的m文件unit42.m：
代码如下：
Fx = quad(inline('x.^(1 / 2). * log(x)'), 0.000001, 1，0.0001)
error = Fx + 4 / 9
运行文件结果为：
Fx = - 0.4440
error = 4.3273e - 004
取不同的步长h：
h = 0.00001，Fx = - 0.4444，error = - 6.3537e - 005；
h = 0.000001，Fx = - 0.4444，error = - 2.9938e - 006；
h = 0.0000001，Fx = - 0.4444，error = - 3.0657e - 007；
h = 0.0000000000001，Fx = - 0.4444，error = - 9.6548e - 009；
h = 0.00000000000001，Fx = - 0.4444，error = - 9.6548e - 009。
可见，存在一个最小的h使精度不能在改善，且h在0.0000000000001与0.00000000000001之间。
4.3龙贝格算法
数值计算方法课程设计(C语言)
数值计算方法课程设计 姓名 学号 成绩
课程实际报告 实验一：秦九韶算法 题目 用选列主元高斯消去法解线性方程组 ???????=+- =-+-=-+-=--02 02 0 21 34343232121x x x x x x x x x x 算法语言： 利用c 语言的知识编写该算法程序 算法步骤叙述： 秦九昭算法的基思路是v[0]=a[0]*x+a[1] v[i]=v[i-1]*x+a[i+1];利用秦九昭算法计算多项式函数。 程序清单： #include
for(i=5;i>=1;i--) {sum=sum*x+a[i-1]; } printf("f(x)=%f/n",sum); } 输出结果计算：
实验总结： 通过运用C 语言，解决了秦九韶算法手写的复杂。为以后的雪地打下基础。 实验二：用选列主元高斯消去法解线性方程组 题目 用选列主元高斯消去法解线性方程组 ???????=+- =-+-=-+-=--02 0 2 0 21 34343232121x x x x x x x x x x 算法步骤叙述 第一步消元——在增广矩阵（A,b ）第一列中找到绝对值最大的元素，将其所在行与第一行交换，再对（A,b ）做初等行变换使原方程组的第一列元素除了第一行的全变为0； 第二步消元——在增广矩阵（A,b ）中第二列中（从第二行开始）找到绝对值最大的元素，将其所在行与第二行交换，再对（A,b ）做初等行变换使原方程组的第二列元素除了第一和第二行的全变为0； 第三步消元——在增广矩阵（A,b ）中第三列中（从第三行开始）找到绝对值最大的元素，将其所在行与第三行交换，再对（A,b ）做初等行变换使原方程组的第三列第四行元素为0； 第四，按x4-x3-x2-x1的顺序回代求解出方程组的解,x[n]=b[n]/a[n][n],x[i]=(b[i]-Σa[i][j]x[j])/a[i][i],i=n-1,…,2,1 程序清单： #include
计算器课程设计报告
课设报告 福建工程学院软件学院 题目：汇编计算器 班级： 1301 姓名 学号： 指导老师： 日期：
目录 1、设计目的 (3) 2、概要设计 (3) 2.1 系统总体分析 (3) 2.2 主模块框图及说明 (3) 3、详细设计 (4) 3.1 主模块及子模块概述 (4) 3.2各模块详运算 (4) 4、程序调试 (7) 4.1 运行界面分析 (7) 算法分析 (7) 4.2 调试过程与分析 (9) 5、心得体会 (11) 5.1 设计体会 (11) 5.2 系统改进 (11) 附录： (11)
1、设计目的 本课程设计是一次程序设计方法及技能的基本训练，通过实际程序的开发及调试，巩固课堂上学到的关于程序设计的基本知识和基本方法，进一步熟悉汇编语言的结构特点和使用，达到能独立阅读、设计编写和调试具有一定规模的汇编程序的水平。 2、概要设计 用8086汇编语言编写一个能实现四则混合运算、带括号功能的整数计算器程序。程序能实现键盘十进制运算表达式的输入和显示（例如输入：“1+2*(3-4)”），按“=”后输出十进制表示的运算结果。 2.1 系统总体分析 在8086的操作环境下，该计算器分成输入，数据存储，运算功能，输出几个大模块，实现了使用者使用该计算器时输入一个算式，能让系统进行计算。此计算器的实现功能是基本的数学的四则运算，结果范围在0~65535。 2.2 主模块框图及说明 此流程图简要的表现出了所要实现的功能以及一些功能的大概算法，同时也是我编写的一个总体的框架。 程序流程图说明：通过流程图，可以看出程序运行时，首先输出提示语气，当用户输入后，程序根据所输入内容进行判断，通过判断的结果来决定调用哪个功能模块，首要先要要判断的是否为0-9，“+”“-”“*”“/”这些字符，若不是就会报错，实则根据运算符号调用其功能模块完成运算。最后将运算的结果显示在主频幕上，返回主程序，使用户可以重新输入。
计算材料课程设计--计算BN的弹性常数
课程设计任务书 2011—2012 学年第1 学期 课程名称：计算材料学 一、设计题目：计算BN的弹性常数 二、完成期限：自2011 年12 月 4 日至2011 年12 月12 日共 2 周 内容及任务1.DFT基本理论，CASTEP使用方法 2.晶体模型的建立与几何优化，相关性质的计算。 3.计算BN的弹性常数 4.结果分析 5. 报告写作与修改 进度安排 起止日期工作内容 11-12-4-6 熟悉DTT理论，软件安装，认识界面，熟悉基本操作11-12-7 晶体模型建立，进行结构优化，计算物理性质 11-12-8 物理性质，力学性质的计算 11-12-9 计算BN的弹性常数 11-12-10-12 写出课程设计的总结实验报告.，修改成文 主要参考资料[1] Kohn W, Sham L J, Self-consistent equations including exchange and correlation effects [J]. Physical review, 1965, 140(4):A1133-A1338. [2] Hohenberg P, Kohn W. Inhomogeneous electron gas [J]. Physical review, 1964,136(3):B864- B871. [3] 谢希德, 陆栋.固体能带理论[M].上海：复旦大学出版，1998. [4] Perdew J P, Chevary J A, Vosko S H. Atoms, molecules, solids, and surfaces: Applications of the generalized gradient approximation for exchange and correlation[J]. Physical review B, 1992, 46(11): 6671-6687. 指导教师（签字）：年月日 系（教研室）主任（签字）：年月日 1
数值计算方法课程设计
重庆邮电大学 数学与应用数学 专业 《数值计算方法》课程设计 姓名： 李金徽 王莹 刘姝楠 班级： 1131001 1131002 1131002 学号： 2010213542 2010213570 2010213571 设计时间： 2012-6-4 指导教师： 朱伟
一、课程设计目的 在科学计算与工程设计中，我们常会遇到求解线性方程组的问题，对于系数矩阵为低阶稠密矩阵的线性方程组，可以用直接法进行消元，而对于系数矩阵为大型稀疏矩阵的情况，直接法就显得比较繁琐，而迭代法比较适用。比较常用的迭代法有Jacobi 迭代与Gauss － seidel 迭代。本文基于两种方法设计算法，并比较他们的优劣。 二、课程设计内容 给出Jacobi 迭代法和Gauss-Seidel 迭代法求解线性方程组的算法思想和MATLAB 程序实现，并对比分析这两种算法的优劣。 三、问题的分析（含涉及的理论知识、算法等） Jacobi 迭代法 方程组迭代法的基本思想和求根的迭代法思想类似，即对于线性 方程组Ax = b( 其中n n n R b R R A ∈?∈,），即方程组 )1(2211222221211 1212111?? ???? ?=+?++??=+?++=+?++n n nn n n n n n n b x a x a x a b x a x a x a b x a x a x a 将系数矩阵A 写为 )2(000000 21122 12122 11U L D a a a a a a a a a A n n n n nn --≡??? ?? ? ? ??---- ??????? ??----??????? ??= 若选取D M =，则U L A M N +=-=，方程组)1(转化为等价方程组 b x U L Dx ++=)(
数值分析课程设计
淮海工学院计算机工程学院课程设计报告书 课程名：《数值分析》 题目：数值分析课程设计 班级： 学号： 姓名：
数值分析课程设计 课程设计要求 1、研究第一导丝盘速度y与电流周波x的关系。 2、数据拟合问题运用样条差值方法求出温度变化的拟合曲线。 课程设计目的 1、通过编程加深对三次样条插值及曲线拟合的最小二乘法的理解； 2、学习用计算机解决工程问题，主要包括数据处理与分析。 课程设计环境 visual C++ 6.0 课程设计内容 课程设计题目1： 合成纤维抽丝工段中第一导丝盘的速度对丝的质量有很大的影响，第一丝盘的速度和电流周波有重要关系。下面是一组实例数据： 其中x代表电流周波，y代表第一导丝盘的速度 课程设计题目3： 在天气预报网站上获得你家乡所在城市当天24小时温度变化的数据，认真观察分析其变化趋势，在此基础上运用样条差值方法求出温度变化的拟合曲线。然后将该函数曲线打印出来并与原来的温度变化数据形成的曲线进行比较，给出结论。写出你研究的心得体会。 课程设计步骤 1、利用最小二乘法写出题1的公式和算法； 2、利用excel表格画出数据拟合后题1的图像； 3、在Visual C++ 6.0中编写出相应的代码； 4、搜索11月12日南通当地一天的温度变化数据； 5、在Visual C++ 6.0中编写出相应的代码； 6、利用excel表格画出数据拟合后题3的图像 课程设计结果 课程设计题目1 数值拟合
解：根据所给数据，在excel窗口运行： x=[49.2 50.0 49.3 49.0 49.0 49.5 49.8 49.9 50.2 50.2] y=[16.7 17.0 16.8 16.6 16.7 16.8 16.9 17.0 17.0 17.1] 课程设计题目3 数据为：X=[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23]; Y=[12,12,11,12,12,12,12,12,13,15,16,17,17,18,17,17,17,16,15,15,15,15,14,14]; 源代码为： 第一题： #include
单片机计算器的课程设计报告
目录 一、设计任务和性能指标 (1) 1.1设计任务 (2) 1.2性能指标 (2) 二、设计方案 (2) 3 3 4 5 5 6 6 7 7 20 20 20 20 21 参考文献 (21) 附录1、系统硬件电路图 (22) 附录2、硬件实物图 (23) 附录3、器件清单 (24)
一、设计任务和性能指标 1.1设计任务 自制一个单片机最小系统，包括复位电路，采用外部小键盘输入数据，能够实现加法、乘法及一个科学计算，计算结果显示在四位一体的数码管上。 要求用Protel 画出系统的电路原理图（要求以最少组件，实现系统设计所要 显 位 监测模块采用二极管和扬声器（实验室用二极管代替）组成电路。 键盘电路采用4*4矩阵键盘电路。 显示模块采用4枚共阳极数码管和74ls273锁存芯片构成等器件构成。 整个单片机的接口电路： P0用于显示输出； P1用于键扫描输入； P2用于数码管位选控制； P3用于键盘扩展（部分运算符输入）；
三.系统硬件设计 3.1单片机最小系统 单片机最小系统就是支持主芯片正常工作的最小电路部分，包括主控芯片、复位电路和晶振电路。 主控芯片选取STC89C52RC芯片，因其具有良好的性能及稳定性，价格便宜应用方便。 扩展键：“log”，“ln”，“x^2”“小数点”，“开方” 共计25个按键，采用4*4矩阵键盘，键盘的行和列之间都有公共端相连，四行和四列的8个公共端分别接P1.0~P1.7，这样扫描P1口就可以完成对矩阵键盘的扫描，通过对16个按键进行编码，从而得到键盘的口地址，对比P1口德扫描结果和各按键的地址，我们就可以得到是哪个键按下，从而完成键盘的功能。 以下为键盘接口电路的硬件电路图
吉林大学材料力学课程设计五种传动轴的静强度、变形及疲劳强度的计算
1．材料力学课程设计的目的 本课程设计的目的是在于系统学完材料力学之后，能结合工程中的实际问题，运用材料力学的基本理论和计算方法，独立地计算工程中的典型零部件，以达到综合运用材料力学的知识解决工程实际问题之目的。同时，可以使学生将材料力学的理论和现代计算方法及手段融为一体。既从整体上掌握了基本理论和现代的计算方法，又提高了分析问题，解决问题的能力；既把以前所学的知识（高等数学、工程图学、理论力学、算法语言、计算机和材料力学等）综合运用，又为后继课程（机械设计、专业课等）打下基础，并初步掌握工程中的设计思想和设计方法，对实际工作能力有所提高。具体的有以下六项： 1．使学生的材料力学知识系统化、完整化； 2．在系统全面复习的基础上，运用材料力学知识解决工程中的实际问题； 3．由于选题力求结合专业实际，因而课程设计可以把材料力学知识和专业需要结 合起来； 4．综合运用了以前所学的个门课程的知识（高数、制图、理力、算法语言、计算机等等）使相关学科的知识有机地联系起来； 5．初步了解和掌握工程实践中的设计思想和设计方法； 6．为后继课程的教学打下基础。 2．材料力学课程设计的任务和要求 要求参加设计者，要系统地复习材料力学的全部基本理论和方法，独立分析、判断、设计题目的已知条件和所求问题。画出受力分析计算简图和力图，列出理论依据和导出计算公式，独立编制计算程序，通过计算机给出计算结果，并完成设计计算说明书。
3．材料力学课程设计题目 传动轴的强度、变形及疲劳强度计算 3-1 设计题目 传动轴的材料为优质碳素结构钢（牌号45），许用应力[σ]=80MPa，经高频淬火处理，其=650MPa,=300MPa, =155MPa,磨削轴的表面，键槽均为端铣加工，阶梯轴过渡圆弧r均为2，疲劳安全系数n=2，要求： 1)绘出传动轴的受力简图; 2)作扭矩图及弯矩图； 3)根据强度条件设计等直轴的直径； 4)计算齿轮处轴的挠度；（按直径Φ1的等直杆计算） 5)对阶梯传动轴进行疲劳强度计算；（若不满足，采取改进措施使其满足疲 劳强度）； 6)对所取数据的理论根据作必要的说明。 说明： 1)坐标的选取均按下图6—1所示； 2)齿轮上的力F与节圆相切； 3)数据表中为直径D的皮带轮传递的功率，为直径为D1的皮带轮传递 的功率。
JAVA实现计算器课程设计
JAVA实现计算器课程设计 计算机科学系 计应0701班 指导老师:刘其昌 设计人员:陈秀桃 设计日期:2009年11月10日——2009年12月20日 计算器的设计 目录 第一章绪 论 ..................................................................... .. (2)
1.1 开发环 境 ..................................................................... . (2) 1.2 基本功能介 绍 ..................................................................... ......... 2 第二章系统设 计 ..................................................................... (3) 2.1 系统流程 图 ..................................................................... . (3) 2.2 系统功能框 图 ..................................................................... (3) 2.3 需求分 析 ..................................................................... ................ 4 第三章软件设 计 ..................................................................... (5) 3.1 界面设 计 ..................................................................... . (5) 3.2 代码设 计 .....................................................................
科学计算器课程设计报告C课程设计修订稿
科学计算器课程设计报告C课程设计 集团标准化工作小组 [Q8QX9QT-X8QQB8Q8-NQ8QJ8-M8QMN]
计算机科学与技术学部 C++课程设计 题目科学计算器 学部计算机科学与技术 班级计科1103 指导教师李军 姓名刘明 学号 2012年6月27日
摘要 计算器的产生和发展是建立在电子计算机基础之上的。硬件方面，自1946年第一台电子计算机诞生以来，计算机技术的发展可谓日新月异，从庞大的只能在实验室里供研究使用的计算机到如今能适应不同环境满足不同需求的各种各样的计算机；运算速度从每秒几千次到每秒几百亿次；处理器从焊有上百万个电子管的大的惊人的电子板到只有指甲大小的集成电路；现在计算机在硬件方面的发展已达到了每三个月更新换代一次的惊人速度。软件方面，也已从机器语言、汇编语言、高级语言发展到现如今的第四代语言——非结构化、面向对象、可视化的语言。 在这个计算器里面，我们实现了简单的四则运算以及更高功能的科学计算，它的外观简洁美观，使人们能快捷简单的操作。能准确的得到计算结果，大大减少了数字计算所需要的时间，为人们的生活带来便利。此系统在Windows 7环境下，使用VC++ 进行编写。 简单计算器包括双目运算和单目运算功能，双目运算符包含基本的四则运算及乘幂功能，单目运算符包含正余弦，对数，开方，阶乘，倒数，进制转换等运算。可对其输入任意操作数，包括小数和整数及正数和负数进行以上的所有运算并能连续运算。并且包含清除，退格功能等。我们所做的计算器其功能较Windows 7下的计算器还是很不够多，没有其菜单的实现功能项，没有其小巧的标准计算器。 关键词：计算器；运算；VC++等
材料力学课程设计
目录 一、 关于材料力学课程设计 (2) 二、 设计题目 (2) 三、 设计内容 (3) 3.1 柴油机曲轴的受力分析 (3) 3.2 设计曲轴颈直径d ，主轴颈直径D (6) 3.3 设计h 和b,校核曲柄臂强度 (6) 3. 4 校核主轴颈H —H 截面处的疲劳强度，取疲劳安全系数n=2。键 槽为端铣加工，主轴颈表面为车削加工 (6) 3.5 用能量法计算A —A 截面的转角y θ，x θ (7) 3.6对计算过程的几点必要说明 (9) 3.7 改进方案 (10) 四、 计算机程序设计 (10) 4.1程序框图 (10) 4.2计算机程序 (11) 4.3输出结果 (12) 五、 设计体会 (12) 六、 参考书目 (12) 一、 关于材料力学课程设计 1.材料力学课程设计的目的 本课程设计的目的是在于系统学完材料力学之后，能结合工程中的实际问题，运用材料力学的基本理论和计算方法，独立地计算工程中的典型零部件，以达到综合运用材料力学的知识解决工程实际问题之目的。同时，可以使学生将材料力学的理论和现代计算方法及手段融为一体，既从整体上掌握了基本理论和现代的计算方法，又提高了分析问题、解决问题的能力；既把以前所学的知识(高等数学、工程图学、理论力
学、算法语言、计算机和材料力学等)综合运用，又为后继课程(机械设计、专业课等)打下基础，并初步掌握工程中的设计思想和设计方法，对实际工作能力有所提高。具体的有以下六项: (1)使学生的材料力学知识系统化、完整化; (2)在系统全面复习的基础上.运用材料力学知识解决工程中的实际问题; (3)由于选题力求结合专业实际.因而课程设计可以把材料力学知识和专业需要结 合起来; (4)综合运用了以前所学的多门课程的知识(高数、制图、理力、算法语言、计算 机等等)使相关学科的知识有机地联系起来; (5)初步了解和掌握工程实践中的设计思想和设计方法; (6)为后继课程的教学打下基础 2.材料力学课程设计的任务和要求 参加设计者要系统地复习材料力学的全部基本理论和方法.独立分析、判断、设计题目的已知条件和所求问题.画出受力分析计算简图和内力图.列出理论依据和导出计算公式.独立编制计算程序.通过计算机给出计算结果.并完成设计计算说明书. 3.材料力学课程设计的一般过程 材料力学课程设计与工程中的一般设计过程相似.从分析设计方案开始到进行必要的计算并对结构的合理性进行分析.最后得出结论.材料力学设计过程可大致分为以下几个阶段: (1)设计准备阶段:认真阅读材料力学课程设计指导书.明确设计要求.结合设计题目复习材料力学课程设计的有关理论知识.制定设计步骤、方法以及时间分配方案等; (2)从外力变形分析入手,分析及算内力、应力及变形,绘制各种内力图及位移、转角曲线; (3)建立强度和刚度条件.并进行相应的设计计算及必要的公式推导; (4)编制计算机程序并调试; (5)上机计算,记录计算结果; (6)整理数据,按照要求制作出设计计算说明书; (7)分析讨论设计及计算的合理性和优缺点,以及相应的改进意见和措施; 二、设计题目 某柴油机曲轴可以简化为下图所示的结构，材料为球墨铸铁（QT450—5），弹性常数为E 、μ，许用应力[σ],G 处输入转矩为e M ，曲轴颈中点受切向力t F 、径向力r F 的作用，且2t r F F = 。曲柄臂简化为矩形截面，1.4≤h D ≤1.6，2.5≤h b ≤4，3l =1.2r ，有关数据如下表：
计算器课程设计报告
高级语言程序（JAVA）课程设计报告 系部名称：商学系专业班级：营销*** 学生姓名：墨璇 墨兰学号： ********** ********** 指导教师：王芬教师职称：讲师 2014年06月26日
目录 一、课程设计目的及意义 .................................... 错误!未定义书签。 二、课程设计任务 .......................................... 错误!未定义书签。 2.1程序设计要求....................................... 错误!未定义书签。 三、课程设计时间 .......................................... 错误!未定义书签。 四、课程设计地点 .......................................... 错误!未定义书签。 五、课程设计内容 .......................................... 错误!未定义书签。 5.1开发工具与平台..................................... 错误!未定义书签。 (1).开发工具 ...................................... 错误!未定义书签。 (2).开发平台 ...................................... 错误!未定义书签。 5.2设计思路........................................... 错误!未定义书签。 5.3 程序测试 .......................................... 错误!未定义书签。 5.4实验总结........................................... 错误!未定义书签。 六、课程设计感想 .......................................... 错误!未定义书签。 七、附录（程序代码） ...................................... 错误!未定义书签。
计算力学课程设计报告
计算力学课程设计说明书 作者: 学号： 学院: 专业(方向): 题目: ANSYS软件建模 指导者： (姓名) (专业技术职务) (姓名) (专业技术职务) 评阅者： (姓名) (专业技术职务) 2015 年 9 月
题目1： 下图所示为一个支座类零件的结构示意图，要求用ANSYS软件建立相应的实体模型和有限元离散模型，说明所用单元的种类、单元数和节点总数。（单位：mm） 解题步骤： １、定义工作文件名、标题，并更改工作目录 将文件名和标题均改为zuoye1，工作目录改为H:\ANSYS。 ２、选择计算类型 选择Main Manu>preferences,勾选Structural。 ３、定义单元类型 选择Main Manu>Preprocessor>Element Type>Add/Edit/Delete命令，在弹出的对话框内选择Solid Brick 8node 185单元，点击OK。如下图。 ４、自底向上创建几何模型 （1）绘制底座：选择Main Manu>Preprocessor>Modeling> Create>Volumes> Block>By Dimensions命令,输入相应的坐标参数，点击OK，生成底座。见下图。
（2）打圆孔：选择Main Manu>Preprocessor>Modeling>Create>Volumes> Cylinder>Solid Cylinder命令，输入相应参数，先后建立两个圆柱模型，然后选择Main Manu>Preprocessor>Modeling>Operate>Booleans>Subtract>Volumes 命令，减掉两个圆柱，从而打出两个圆孔。如下图。 （3）倒圆角：选择Main Manu> Preprocessor>Modeling>Create>Areas>Area Fillet命令，选择需要倒角的两个面，输入圆角半径，创建出两个圆弧面，然后选择Main Manu> Preprocessor>Modeling>Operate>Booleans>Divide>Volume By Area命令，切割出两个圆角，最后选择Main Manu>Preprocessor>Modeling> Delete>Volume and Below命令，删除圆角外的多余部分。如图所示。 （4）绘制底座凹槽：选择Main Manu> Preprocessor>Modeling>Create> Volumes>Block>By Dimensions命令，输入凹槽各点坐标，在凹槽位置创建出一个和底座重合的长方体块，然后选择Main Manu>Preprocessor>Modeling> Operate> Booleans>Subtract>Volumes命令，在底座中减去重合部分即可。如图所示。
计算方法课程设计
数理学院2014级信息与计算科学 课程设计 姓名：刘金玉 学号： 3141301240 班级： 1402 成绩：
实验要求 1．应用自己熟悉的算法语言编写程序，使之尽可能具有通用性。2．上机前充分准备，复习有关算法，写出计算步骤，反复检查，调试程序。（注：在练习本上写，不上交） 3．完成计算后写出实验报告，内容包括:算法步骤叙述，变量说明，程序清单，输出计算结果，结构分析和小结等。（注：具体题目 具体分析，并不是所有的题目的实验报告都包含上述内容！）4．独立完成，如有雷同，一律判为零分！ 5．上机期间不允许做其他任何与课程设计无关的事情，否则被发现一次扣10分，被发现三次判为不及格！非特殊情况，不能请 假。旷课3个半天及以上者，直接判为不及格。
目录 一、基本技能训练 (4) 1、误差分析 (4) 2、求解非线性方程 (6) 3、插值 (12) 4、数值积分 (12) 二、提高技能训练 (16) 1、 (16) 2、 (18) 三、本课程设计的心得体会（500字左右） (21)
一、基本技能训练 1、误差分析 实验1.3 求一元二次方程的根 实验目的： 研究误差传播的原因与解决对策。 问题提出：求解一元二次方程20ax bx c ++= 实验内容： 一元二次方程的求根公式为 1,22b x a -+= 用求根公式求解下面两个方程： 2210(1)320(2)1010 x x x x +-=-+= 实验要求： （1） 考察单精度计算结果（与真解对比）； （2） 若计算结果与真解相差很大，分析其原因，提出新的算法（如先求1x 再 根据根与系数关系求2x ）以改进计算结果。 实验步骤： 方程（1）： 根据求根公式，写出程序： format long a=1;b=3;c=-2; x1=((-1)*b+sqrt(b^2-4*a*c))/2*a x2=((-1)*b-sqrt(b^2-4*a*c))/2*a
《数值分析》课程设计报告
《数值分析》课程设计实验报告 龙格—库塔法分析Lorenz 方程 200820302033 胡涛 一、问题叙述 考虑著名的Lorenz 方程 () dx s y x dt dy rx y xz dt dz xy bz dt ?=-???=--???=-?? 其中s ，r ，b 为变化区域内有一定限制的实参数，该方程形式简单，表面上看并无惊人之处，但由该方程揭示出的许多现象，促使“混沌”成为数学研究的崭新领域，在实际应用中也产生了巨大的影响。 二、问题分析 Lorenz 方程实际上是一个四元一阶常微分方程，用解析法精确求解是不可能的，只能用数值计算，最主要的有欧拉法、亚当法和龙格- 库塔法等。为了得到较高精度的，我们采用经典四阶龙格—库塔方法求解该问题。 三、实验程序及注释 （1）算法程序 function [T]=Runge_Kutta(f,x0,y0,h,n) ％定义算法，其中f 为待解方程组， x0是初始自变量，y0是初始函数 值，h 是步长，n 为步数 if nargin<5 n=100; ％如果输入参数个数小于5，则步数 n=100 end r=size(y0);r=r(1); ％返回初始输出矩阵的行列数，并将 值赋给r(1) s=size(x0);s=s(1); ％返回初始输入矩阵的行列数，并 将值赋给s(1) r=r+s; T=zeros(r,n+1); T(:,1)=[y0;x0]; for t=2:n+1 ％以下是具体的求解过程 k1=feval(f,T(1:r-1,t-1)); k2=feval(f,[k1*(h/2)+T(1:r-1,t-1);x0+h/2]); k3=feval(f,[k2*(h/2)+T(1:r-1,t-1);x0+h/2]); k4=feval(f,[k3*h+T(1:r-1,t-1);x0+h]); x0=x0+h; T(:,t)=[T(1:r-1,t-1)+(k1+k2*2+k3*2+k4)*(h/6);x0]; end
计算机技术综合课程设计报告
计算机技术综合课程设计 设计题目锅炉液位控制系统学生姓名史婷艳 专业班级自动化1302班学号20134460203 指导老师洪镇南 2017年1 月3日
目录 前言 (2) 1 锅炉汽包水位控制对象与控制指标 (4) 1.1锅炉汽包水位的特征 (4) 1.2汽包水位动态特性 (4) 1.2.1汽包水位在给水流量W作用下的动态特性 (4) 1.2.2汽包水位在蒸汽流量D扰动下的动态特性 (5) 1.2.3燃料量B扰动下汽包水位的动态特性 (6) 2. 汽包水位控制方案 (7) 2.1单冲量控制方式 (7) 2.2 双冲量控制方式 (8) 2.3 三冲量控制方式 (9) 3. 三冲量串级PID控制 (11) 3.1 串级PID控制 (11) 3.2 智能整定PID控制 (12) 4 汽包水位模糊控制器设计及仿真 (12) 4.1 输入输出变量 (12) 4.2 隶属度函数 (15) 4.3基于MATLAB/Simulink 环境建立的系统仿真分析 (16) 4.3.1 基于MATLAB/Simulink 的系统模型 (16)
4.3.2 仿真结果分析 (18) 总结与体会 (18) 参考文献 (20) 前言 锅炉是典型的复杂热工系统，目前，中国各种类型的锅炉有几十万台，由于设备分散、管理不善或技术原因，使大多数锅炉难以处于良好工况，增加了锅炉的燃料消耗，降低了效率。同时，锅炉工作过程中各项指标的调节难以建立数学模型，具有非线性、不稳定性、时滞等特点，所以如何改善对锅炉的控制，保证其正常工作，提高效率一直是人们关注的焦点。而汽包液位是锅炉安全、稳定运行的重要指标，保证液位在给定范围内，对于高蒸汽品质、减少设备损耗和运行损耗、确保整个网络安全运行具有要意义。 现代锅炉的特点之一就是蒸发量显著提高，汽包容积相对变小，水位变化速度很快，稍不注意就容易造成汽包满水或者烧成干锅，这都对汽包液位控制系统提出了更高的要求。汽包液位过高，会影响汽包内汽液分离效果，使汽包出口的饱和蒸汽带水增多，蒸汽带水会使汽轮机产生水冲击，引起轴封破损、叶片断裂等事故。同时会使饱和蒸汽中含盐量增高，降低过热蒸汽品质，增加在过热器管壁和汽轮机叶片上的结垢。水位过低，则可能破坏自然循环锅炉汽水循环系统中某些薄弱环节，以致局部水冷管壁被烧坏，严重时会造成爆炸事故。 目前，对汽包液位位控制大多采用常规PID控制方式，从控制方式来看，它们要么系统结构简单成本低，不能有效的控制锅炉汽包“虚假水位”现象，要么能够在一定程度上控制“虚假现象”，系统却过于复杂，成本投入过大。常用的蒸汽锅炉液位调节系统有三种基本结构：单冲量调节系统结构、双冲量调节系统结构、串级三冲量调节系统结
水文分析计算课程设计
《水文分析与计算》课程设计指导书 ———设计年径流及设计洪水的计算 一、课程设计的目的 1．掌握PIII型频率曲线的制作方法 2. 掌握设计年径流及其年内分配的计算方法 3．掌握考虑历史特大洪水的设计洪水及其过程的计算方法 二、课程设计任务 1．根据所给资料推求设计年径流与设计年内分配过程 表1是某站1958～1976年各月径流量资料，根据所给资料推求P=10%的设计丰水年、P=50%的设计平水年、P=90%的设计枯水年的设计年径流量；并计算P=90%的设计枯水年径流年内分配过程。 要求：理论频率曲线采用PIII型分布，由矩法作参数无偏估计，并以估计值为初值，用目估适线法选配理想的理论频率曲线，注意比较验证均值X a、变差系数C V、偏态系数C S对频率曲线的影响效果。检查所选最终的理论频率曲线的合理性，并计算所求设计频率的相应设计年径流，年径流分配过程采用典型年同倍比放大法。 3
三、课程设计成果要求 要求提交设计成果：一份电子文档，一份打印文档。设计中的计算可采用采用excel 或编程计算，编程语言可采用FORTRAN 语言、C 语言、Basic 语言或同等功能的语言编程。要求程序正确、可靠、可运行，符合结构化程序设计思想，具有易读性、可修改性、可验证性、通用性，关键变量应作注释说明。计算结果要表格化，便于检查、保存和打印。设计设计报告，其重点是对计算成果的说明和合理性分析及其有关问题的讨论。要求文字流畅，简明扼要；图表整齐清楚，名称、编号齐全；封面统一，最后装订成册。 四、课程设计的考核 平日考勤、设计报告，加上抽查提问及上机操作，对成绩进行综合评定。 五、课程设计时间与地点 时间： 2013年5月9日星期四 地点： 学院 六、实验原理 1.经验频率计算 经验频率：P=m/(n+1)*100%，模比系数：Q Q Ki i = 2．线型选择 频率曲线一般应采用皮尔逊Ⅲ型。 3.频率曲线参数估计 平均值：n 1 ∑== n i i Q Q 变差系数：() 1 n 11 2 --= ∑=n i i v K C 4.偏态系数：Cs=2-3Cv 七、实验步骤 1、将测站所得数据年份及年平均流量数据复制与Excel 表格中，并列出序号，同时计算出年平均流量的均值。 2、另起一列，将年平均流量数据按从大到小排列。按数学期望公式计算出相应经验频率P=m/(n+1)*100%。在画图软件上绘制经验点距。再计算出各相应的模比系数Ki （Q Q Ki i =）和（Ki-1）2。 3、选定水文频率分布线型（选用皮尔逊Ⅲ型）。 表2 某站年径流量频率计算表
数值分析课程课程设计汇总
课 程 设 计 我再也回不到大二了， 大学是那么短暂 设计题目 数值分析 学生姓名 李飞吾 学 号 x x x x x x x x 专业班级 信息计x x x x x 班 指导教师 设 计 题 目 共15题如下 成绩
数值分析课程设计 1．1 水手、猴子和椰子问题：五个水手带了一只猴子来到南太平洋的一个荒岛上，发现那里有一大堆椰子。由于旅途的颠簸，大家都很疲惫，很快就入睡了。第一个水手醒来后，把椰子平分成五堆，将多余的一只给了猴子，他私藏了一堆后便又去睡了。第二、第三、第四、第五个水手也陆续起来，和第一个水手一样，把椰子分成五堆，恰多一只猴子，私藏一堆，再去入睡，天亮以后，大家把余下的椰子重新等分成五堆，每人分一堆，正好余一只再给猴子，试问原先共有几只椰子？（15621） 试分析椰子数目的变化规律，利用逆向递推的方法求解这一问题 解：算法分析：解该问题主要使用递推算法，关于椰子数目的变化规律可以设起初的椰子数为0p ，第一至五次猴子在夜里藏椰子后，椰子的数目分别为01234,,,,p p p p p 再设最后每个人分得x 个椰子,由题: 14 (1)5 k k p p +=- （k=0,1,2,3,4）51(1)5 x p =- 所以551p x =+,11k k p p +=+利用逆向递推方法求解 15 1,4 k k p p +=+ （k=0,1,2,3,4） MATLAB 代码： n=input('n= '); n= 15621 for x=1:n p=5*x+1; for k=1:5 p=5*p/4+1; end if p==fix(p), break end end disp([x,p]) 1．2 设，1 5n n x I dx x =+? （1）从0I 尽可能精确的近似值出发，利用递推公式： 11 5(1,2,20)n n I I n n -=-+= 计算机从1I 到20I 的近似值； （2）从30I 较粗糙的估计值出发，用递推公式：
java课程设计报告计算器
java课程设计报告 计算器 1 2020年4月19日
目录 一、课程设计目的 (5) 二、课程设计任务.................................................................................................... .. (2) 2.1、设计任务 (5) 2.2、课程设计要求： (6) 2.3、需求分析 (6) 三、开发工具与平台.................................................................................................... (3) 3.1、开发工具 (7) 3.2、开发平台 (7) 2 2020年4月19日
四、设计思路.................................................................................................... . (4) 4.1、界面设计.................................................................................................... . (4) 4.2.1、逻辑设计 (8) 4.2.2、程序流程图.................................................................................................... . (5) 4.2.3、主要代码展示及说明 (5) 4.3、程序测试 (18) 五、实验小结 (20) 3 2020年4月19日
计算材料课程设计
湖南工业大学 课程设计 资料袋 理学院（系、部）2011 ~ 2012 学年第一学期课程名称计算材料学指导教师雷军辉职称讲师 学生姓名张丽佳专业班级应用物理081班学号08411200125 题目计算BN的弹性常数 成绩起止日期2011年12月4日～2011年12 月12 日 目录清单 序号材料名称资料数量备注 1 课程设计任务书 2 课程设计说明书 3 课程设计图张 4 5 6 湖南工业大学 1
课程设计任务书 2011—2012 学年第1 学期 理学院学院（系、部）应用物理学专业081 班级 课程名称：计算材料学 一、设计题目：计算BN的弹性常数 二、完成期限：自2011 年12 月 4 日至2011 年12 月12 日共 2 周 内容及任务1.DFT基本理论，CASTEP使用方法 2.晶体模型的建立与几何优化，相关性质的计算。 3.计算BN的弹性常数 4.结果分析 5. 报告写作与修改 进度安排 起止日期工作内容 11-12-4-6 熟悉DTT理论，软件安装，认识界面，熟悉基本操作11-12-7 晶体模型建立，进行结构优化，计算物理性质 11-12-8 物理性质，力学性质的计算 11-12-9 计算BN的弹性常数 11-12-10-12 写出课程设计的总结实验报告.，修改成文 主要参考资料[1] Kohn W, Sham L J, Self-consistent equations including exchange and correlation effects [J]. Physical review, 1965, 140(4):A1133-A1338. [2] Hohenberg P, Kohn W. Inhomogeneous electron gas [J]. Physical review, 1964,136(3):B864- B871. [3] 谢希德, 陆栋.固体能带理论[M].上海：复旦大学出版，1998. [4] Perdew J P, Chevary J A, Vosko S H. Atoms, molecules, solids, and surfaces: Applications of the generalized gradient approximation for exchange and correlation[J]. Physical review B, 1992, 46(11): 6671-6687. 指导教师（签字）：年月日系（教研室）主任（签字）：年月日 2