核辐射测量数据处理成都理工大学程序设计报告
成都理工大学数据库实验报告
本科生实验报告实验课程《数据库》学院名称专业名称学生姓名学生学号指导教师实验地点实验成绩二〇一五年三月二〇一五年四月学生成绩管理系统设计第一章系统规划1.1设计背景及目的传统的学生成绩管理工作主要是成绩登记、成绩查询、学生信息查询。
该工作主要由任课教师来完成。
由于随着教师授课班级学生人数规模的持续壮大，学生的各种信息也不断增加。
面对庞大的信息量，传统的人工式的学生管理导致工作繁重，人力物力过多的浪费，况且工作质量也难以得到可靠的保证，往往导致一系列错误的发生，从而影响了整个学生信息的管理。
时代的进步，计算机越来越广泛的运用，原先主要由手工操作的学生信息管理，正在慢慢的由计算机代替，并且日益完善。
学生信息管理系统就是为了解决和减轻繁琐的手工管理，使教师对学生信息可以上升到一个完善的电子化，信息化管理。
在全球信息化的今天，是同计算机进行生产和管理已经成为一种趋势。
计算机可以为我们选择最优的方案来达到我们的要求。
在面对越来越繁重的教学任务环境下，使用计算机管理可以最大程度减少操作难度，减轻教师的工作负担。
实现“以人为本，科学管理“的目标。
1.2可行性分析学生信息管理系统是利用计算机技术使学生个人信息管理、成绩查询、成绩登记更改等有机的结合在一起方便教师进行统一高效率的管理。
对于系统的结构特性设计，应遵循提高系统实用性和操作简便灵活性的设计的原则。
1.2.1经济性可行性研究目前计算机广泛应用于各个领域，尤其是在数据处理方面表现出了巨大的优势，利用计算机将各种复杂的数据，都制作成数据库，交由电脑来管理。
用电脑管理数据，可以最大程度减少操作难度，减轻工作人员的劳动负担，并且运算速度快，可靠性高。
这使得创建一个好的信息管理系统成为一个必要的工作。
但是，建立学生成绩管理系统在经济上是否可行呢？一方面，随着计算机的普及，硬件设备价格不断下降。
建立学生成绩管理系统的费用并不是十分的昂贵。
通常一个小型的学生成绩管理系统只需一台普通计算机。
成都理工辐射防护试验四
实验四：β射线和中子的辐射屏蔽模拟一、实验目的1、掌握辐射防护模拟软件的使用；2、通过不同屏蔽材料、不同屏蔽厚度的模拟，了解β射线防护和屏蔽的效果；3、通过不同能量中子、不同屏蔽材料、不同屏蔽厚度的模拟，了解中子的屏蔽特点和效果。
二、实验原理辐射防护实验模拟软件，通过调用蒙特卡罗模拟软件的方式实现程序的链接，不受辐射源、材料等实验条件的限定。
通过参数的调节，可以研究射线的防护。
实验一为主要模拟β射线的屏蔽效果，以本次设计所用的错误！未找到引用源。
的β源为例，其β射线能量在0-2.28MeV之间。
错误！未找到引用源。
的半衰期为28.6年，它发射的β粒子最大能量为0.546MeV，错误！未找到引用源。
衰变后成为错误！未找到引用源。
，的半衰期为64.1小时，它发射的β粒子最大能量是2.27MeV，衰变后成为错误！未找到引用源。
，因而在0至2.27MeV的范围内形成连续的能谱。
该模拟选择了比较简单的模型，减少其他参数带来的影响。
如图1所示，由于β射线很容易被屏蔽，所以整个设定是在真空中的空间中屏蔽层为可调厚度、可变材料的圆盘。
空间4为检测栅元，用于检测粒子注量，以检测屏蔽效果。
材料中建立了两种不同的材料分别是铝和锡，射线类型选择了束状射线，通过计数卡记录相关问题。
图1 β射线屏蔽模拟实验的模型实验二是对中子屏蔽的模拟实验。
可以选用不同的材料、不同的中子能量、不同的屏蔽厚度进行实验。
如用铜、钨、铅、水、硼砂这几种材料作为屏蔽材料对14MeV能量的快中子进行慢化模拟，每种材料的模拟厚度可以自行调整，得到屏蔽后的中子通量，比较这几种材料的慢化效果，并选取最佳慢化材料以及最佳厚度；对于1MeV能量的中子，选取聚乙烯、聚氯乙烯、水、石蜡、硼砂进行屏蔽模拟，每种材料的模拟厚度也可以自行调整。
基于研究材料屏蔽性能的需要，建立了一个半径为100cm球体（球心为坐标原点）空间内建立屏蔽层模型，该模型由4个空间组成（如图2）。
成都理工大学工程技术学院
成都理工大学工程技术学院核辐射探测实验课程教案教师姓名所在系部能源与环境工程系授课专业核工程与核技术课程代码总学时学分教材名称核辐射探测实验设备清单：1000ml直量筒、玻璃杯（加水用）4．实验步骤（1）固定放射源、探测器和量筒（物料容器）的几何位置，按照设备连接框图连接图3.1 高斯正态分布曲线计数器工作状态的调整计数器是气体探测器的一种。
在实验2中已经使用过。
在其阴极和阳极之间所加高电压的高低直接决定其工作状态。
其电压特性表示如图段的后部﹚，应立即降低高压，不必测量完整的坪曲线，避免损坏计数器，段的前部。
图3.2 G－M计数器电压特性曲线实验步骤1）熟悉所使用的设备和电路插件。
按设备连接框图连接好各个插件。
经教师检查后再接通电源。
2）将高压电源调节旋扭调到最低。
单道分析器处于积分状态。
）放好放射源﹙距探测器约5 cm﹚。
打开高压电源，逐渐升高高压电源电压；在不同高压点，用定标器记录单道分析器输出信号的计数率。
测定曲线”。
并选定其工作电压。
将高压调到工作电压位置。
）移开放射源，测量实验室本底，这时计数率应很低。
若仍然较高，说明有大量图3.3 计数散点图阴影部分面积占总面积60.6%图3.4 计数直方图）数据统计计算平均值N。
计算标准误差σ。
统计在N±σ范围之内所占面积与总面积之比。
二、教学方法和手段、师生互动通过基本理论讲授，使学生学习基本的理论知识和分析解决问题的方法，深学生对所学知识点的掌握。
鼓励学生提问题，并随时解答学生的问题。
让学生成立小组做实验，通过互帮使大部分同学达到本实验要求。
通过思考题使学生巩固所学知识。
三、时间安排：3学时。
成都理工大学电离辐射计量学 辐射剂量学实验3方案
实验一内照射剂量估算一、实验目的1、学习indo2012内照射估算软件的基本使用方法；2、通过软件估算加深内照辐射对器官形成损害的认识；3、了解影响内照辐射损伤程度的若干因素；4、通过对碘、钾、锶等放射性同位素的内照估算实例，学习不同放射性核素的器官分布特点；二、实验原理内照射是指放射性物质经呼吸道、消化道、皮肤、粘膜和伤口以及其他各种途径进入机体后，放射性核素发出的核射线由体内对机体进行的照射。
放射性物质进入人体内主要有三种途径：呼吸吸入、口腔食入、创伤皮入，造成体内污染。
上述几种途径中，以呼吸道入体的毒理学意义最大。
经消化道进入：放射性核素吸收率最高的是碱族元素（钠、钾、铯）和某些非金属元素（碘、碲），可达90％以上；其次是碱土族元素（锶、钡）为10％～40％；镧系和锕系元素的吸收率最低，约为0.01％～0.1％。
工作人员有时用口接触被放射性污染的器具、物品，或环境受到放射性污染，通过食物、饮用水等导致居民和工作人员长时间摄入放射性物质。
胃肠道对核素的吸收率主要取决于化学性质，碱族和卤族元素极易溶雨水，吸收率较高。
经呼吸道进入：放射性核素可以气态、气溶胶或微小粉尘的形式存在于空气中，气态放射性核素（氡、氙、氚）易经呼吸道粘膜或透过肺泡被吸收入血。
粉尘或气溶胶态的放射性核素在呼吸道内的吸收决定于粒径大小及化合物性质。
一般粒径愈大，附着在上呼吸道粘膜上愈多，进入肺泡内愈少，吸收率低。
难溶性化合物在肺内溶解度很低，多被吞噬；而可溶性化合物则易被肺泡吸收入血。
粒径大于1μm者，大部分被阻滞在鼻咽部、气管和支气管内；粒径在0.01～1μm的落下灰危害最大，大部分沉积在肺部（包括细支气管、肺胞管、肺泡、肺泡囊）。
部分吸收入血，部分被吞噬细胞吞噬后滞留在肺内成为放射灶。
沉积在鼻咽部，气管和支气管的放射性灰尘大部分通过咳痰排出体外或吞入胃内，仅少部分吸收入血。
经伤口和皮肤粘膜进入：当皮肤创伤时，放射性物质通过伤口进入，吸收率较高。
放射性核素的剂量计算设计
成都理工大学工程技术学院毕业论文放射性核素的剂量计算设计作者姓名：严俊专业名称：核工程与核技术指导教师：张艳丽讲师放射性核素的计量计算设计摘要辐射剂量的计算在辐射防护、环境监测、辐射调查等多方面都有非常重要的应用。
基于此，本论文根据探测器测得的粒子注量，粒子能量、以及周围介质的吸收系数等，设计出了一种计算监测点处辐射剂量大小的一种思路。
在论文中，我们将实测值与计算值作了比较，比较结果发现，两者的变化趋势吻合较好，但是两者在数值大小上略有差别，本论文就此做了详细的分析。
关键字：粒子注量能量注量吸收剂量剂量当量闪烁体探测器AbstractRadiation dose calculation has very important aspects of the application in radiation protection, environment monitoring, radioation investigation,etc. Based on this, the paper designs a train of thought about how to calculate the size of the radiation dose on the monitoring place according to the particle fluence measured by the probe, particle energy , and the absorption coefficient of the surrounding media. In the paper, we have compared the measured values with the calculated values. The results shows that the variation trend of the measured values with the calculated values are better matched, but they have a slightly difference in numerical size, this paper made a detailed analysis of this.Keywords:particle fluence, energy note amount, absorb dose, dose equivalent, scintillation detectors目录摘要 (I)Abstract (II)目录 (III)前言 (1)1探测器的分类及特点 (2)1.1闪烁体探测器的简介 (2)1.1.1闪烁体的基本组成和工作原理 (2)1.1.2闪烁体探测的分类 (2)1.2气体探测器的简介 (3)1.3半导体探测器的简介 (3)1.4各种探测器优缺点 (3)1.4.1气体探测器的优缺点 (3)1.4.2闪烁体探测器的优缺点 (4)1.4.3半导体探测器的优缺点 (4)2辐射防护与辐射剂量学中常用的物理量 (5)2.1辐射防护中常用的物理量 (5)2.1.1放射性活度 (5)2.1.2粒子注量与粒子注量率 (5)2.1.3能量注量与能量注量率 (6)2.1.4相互作用系数 (7)2.2辐射剂量学中的物理量 (8)2.2.1吸收剂量 (8)2.2.2授予能 (8)2.2.3比释动能与比释动能率 (9)2.3其它常用的物理量 (10)2.3.1照射量 (10)2.3.2剂量当量概念及单位 (10)2.4物理量之间的关系 (11)2.4.1粒子注量与能量注量的关系 (11)2.4.2吸收剂量与比释动能的关系 (11)2.4.3照射量与吸收剂量率之间的关系 (12)3辐射剂量计算设计思路 (14)3.1设计思路简述 (14)3.2实验测量某点的粒子注量 (14)3.2.1实验目的 (14)3.2.2实验仪器设备及连接框图 (15)3.2.3实验步骤 (15)3.2.4实验数据 (16)3.3辐射剂量计算 (16)3.3.1实测值与计算值的比较 (17)3.3.2误差分析 (18)总结 (19)致谢 (20)参考文献 (21)前言辐射与物质（包括生物体）发生作用，引起物理、化学和生物学等方面变化。
核辐射测量方法实验报告
实验二 γ射线的吸收一、实验目的：1、了解γ射线在物质中的吸收规律；2、测量γ射线在不同介质中的吸收系数。
二、实验器材：1、KZG03C 辐射检测仪一台；2、Cs137点放射源一个；3、铅准直器一个；4、40×40×dcm3的水泥、铝、铁、铜、铅吸收屏若干块（附屏支架）；5、手套、长钳夹子、尺子、绳子各一套。
三、实验原理：天然γ射线与物质相互作用的三种主要形式：光电效应、康普顿散射和形成电子对效应。
由于三种效应的结果，γ射线通过物质时发生衰减（吸收），其总衰减系数应为三者之和：实验证明，γ射线在介质中的衰减服从指数规律：de I I μ-=0，mm d e I I μ-=0μ=(- Ln(I/I O ))/d ， μm =(- Ln(I/I O ))/d m式中：I 为射线经过某一介质厚度的仪器净读数（减去本底）；I 0为起始射线未经过介质的仪器净读数（减去本底）； d 为介质厚度，单位为cm; d m 为介质面密度，单位为g/cm 2 ；μ 为γ射线经过介质的线吸收系数，单位为cm -1；κστμ++=μm 为γ射线经过介质的质量吸收系数，单位为g/cm 2 ； 半吸收厚度：为使射线强度减少一半时物质的厚度，即021I I =时，μ2ln 21=d 或 212ln d =μ四、实验内容：1. 选择良好的测量条件（窄束），测量 Cs 137源的γ射线在同一组吸收屏（水泥、铝、铁、铜、铅）中的吸收曲线，并由半厚度定出吸收系数；2. 用最小二乘拟合的方法计算出吸收系数与1中的结果进行比较；3. 测量不同散射介质时（同一角度，同一厚度）γ射线的强度。
五、实验步骤： 1. 吸收实验1) 调整装置，使放射源、准直孔、探测器的中心在一条直线上； 2) 测量本底I 0’；3) 将源放入准直器中，测量无吸收屏时γ射线强度I 0”；4) 逐渐增加吸收屏，并按相对误差在N ±δ的要求测出对应厚度计数I d ’，每个点测三次取平均植；5) 更换一种吸收屏，重复步骤4，测量时注意测量条件不变。
计数器设计-核辐射监测方法与仪器
成都理工大学《核辐射监测方法与仪器》实验（报告）脉冲计数器设计Pulse counter design2011 级核技术与自动化工程学院专业辐射防护与环境工程学号 201106080112学生姓名袁子程指导教师王广西完成日期 2015 年 01月 07日前言脉冲计数是核辐射探测的一项重要中间环节，本着严谨的态度，实验过程中，我多次借鉴教科书中内容。
书中将脉冲计数统一到脉冲幅度分析系统中，将脉冲幅度的甄别与计数作为统一的系统进行设计，也就是说脉冲计数至少就是一个单道脉冲仪。
设置独立的总道是在理论电路中可以体现出来的，我想主要讨论一下计数部分的实现方法与模拟手段。
设计思路：我设计的脉冲计数器的目的是记录显示盖革计数器的计数率（单位时间内的计数），并没有设计存储功能，是一种事实监测的计数器。
模拟波形部分采用NE555产生所需波形，经过整形、计数、译码、显示，此方案性比价高，工作所需环境低，可靠性好，硬件电路较为合理。
NE555定时器产生脉冲信号→施密特触发器整形→计数器计数→数字信号译码→数码管显示NE555定时器电路，可以通过调节电阻电容，用示波器观察，形成我们想要得核脉冲信号，并计算幅度与前置放大器的关系，好进一步匹配。
CC40106由六个斯密特触发器电路组成。
每个电路均为在两输入端具有斯密特触发器功能的反相器。
触发器在信号的上升和下降沿的不同点开、关。
上升电压(V T+)和下降电压(V T-)之差定义为滞后电压。
CD40106引脚2 4 6 8 10 12 数据输出端1 3 5 9 11 13 数据输入端14 电源正7 接地切换时间波形：74LS161是常用的四位二进制可预置的同步加法计数器<74ls161引脚图>管脚图介绍：时钟CP和四个数据输入端P0~P3清零/MR使能CEP，CET置数PE数据输出端Q0~Q3以及进位输出TC. (TC=Q0·Q1·Q2·Q3·CET)<74LS161功能表>从74LS161功能表功能表中可以知道，当清零端CR=“0”，计数器输出Q3、Q2、Q1、Q0立即为全“0”，这个时候为异步复位功能。
成都理工大学电离辐射计量学 辐射剂量学实验一方案
实验一环境中伽玛射线剂量测量与放射性测量统计规律认识一、实验目的1、 了解DETECTIVE‐EX核素甄别器工作原理，掌握基本使用方法；2、使用DETECTIVE‐EX测量环境中伽玛辐射剂量并验证放射性计数的统计规律；3、了解统计误差的意义，掌握计算统计误差的方法；4、学习探测器的应用实例——测量就地X-γ剂量率的方法；5、学习个人剂量计的标定。
二、实验原理当放射源的半衰期足够长（即在实验测量时间内可以认为其强度基本上不变的情况下）时，在作重复的放射性测量中，即使保持完全相同的实验条件和足够的测量精度，每次的测量结果也不会完全相同，而是围绕其平均值上下涨落，有时甚至有很大的差别，这种现象叫做放射性计数的统计性。
放射性计数的统计性是放射性原子核衰变本身固有特性，与使用的测量仪器及技术无关。
对放射源所测数据进行处理时，可以作频率直方图。
频率直方图可以形象地表明数据的分布情况。
作频率直方图的方法如下：首先将测得的数据分组。
分组时可以将数据由小到大排列起来，用最小值（或比最小值还小一点的数）作为分组的起始点；也可以在算出平均值后，将平均值置于组中央来分组。
组数一般不做硬性规定。
但如果组数太少，会掩盖各组内数据随分布的变动情况；组数太多，由于增加了随机因素的影响，会使各组高度参差不齐而得不出简明的规律。
组距的大小通常以将数据分为十几组左右为宜。
我们这里建议以σ/2为组距，并将平均值置于组中央来分组。
显然，用这样的方法分组时，各组的分界点分别为而各组的中间值则为分好组后，统计好数据落在各组中的频数（k i），算出相应的组频率（k i/A,A总的测量次数），并以k i/A为纵坐标，画出一系列小矩形，即成为表示数据统计分布的频率直方图。
三、主要仪器及器材1、DETECTIVE‐EX核素甄别器一台2、个人剂量计5台3、137Cs、 60Co放射源四、实验方法与步骤1、根据仪器使用说明书和操作规程，学会DETECTIVE‐EX核素甄别器的基本使用方法。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
程序设计课程设计(论文)
设计（论文）题目
谱的显示、谱光滑、定量分析（峰面积）程序学院名称核技术与自动化工程学院
专业名称辐射防护与环境工程
学生姓名袁子程
学生学号2
任课教师马永红
设计（论文）成绩
教务处制
2015年 1 月10 日
一、实习目的
学习使用各种编程软件，利用计算机工具开展专业信息处理工作。
使用面向对象编程思想实现谱数据处理功能以及谱显示功能。
二、人员组成及分工
本人题目为：谱的显示、谱光滑、定量分析（峰面积）程序
具体要求：
多项式最小二乘拟合实现谱光滑
对称零面积法实现寻峰
Wasson（瓦森、沃森）峰面积法实现峰面积计算
寻峰应在5、7、9、11之间可选，讨论最佳参数
三、实习计划
1.首先进行MFC绘图区域的学习和了解，实现基本绘图功能
2.建立单文档MFC程序，设计一个合理美观的工作区
3.实现文件打开读取功能，并能打开不同文件
4.进行谱数据处理的算法编写
5.编写成谱图形功能并给出处理结论
四、主要工作介绍（详细介绍实习中各功能的实现方法、理论公式和计
算原理，程序流程图和程序运行界面截图等成果信息）
1.建立单文档类MFC可执行程序
2.建立绘图区函数，初始化绘图设备
3.添加文件打开虚函数
4.编写读取文件及初步寻峰代码
5.添加绘制谱图代码
6. 尝试读取mca 文件
7. 添加谱数据处理代码
①多项式最小二乘拟合法：
∑-=+=
m
m
j j
i j
b
i data
A K data 1
，b=2m+1，为平滑宽度。
谱光滑
5
7
9
11
实际代码：
//多项式最小二乘法谱光滑（参数为5） for( i=2;i<1023;i++)
{
data[i]=1/35*((-3)*data[i-2]+12*data[i-1]+17*data[i]+12*data[i+1]+(-3)*d
ata[i+2]); }
②零面积对称法：
基本思想：面积为零的“窗”函数与实验谱数据进行褶积变换，且要求“窗”函数为对称函数。
数学表达式：j j m
m
j j
m
m
j j
i j
i C C C
data
C y --=-=+===∑∑0
~
其中，i y ~
为变换后的谱数据，data i+j 为原始实验谱数据，C j 为对称零面积变换函数，W=2m +1为窗宽（变换宽度）。
对称零面积变换函数：∑-==
-=m
m
j j
j j G
W
d d
G C 1
G j ：[]
⎪⎪
⎪⎩⎪⎪⎪⎨⎧==+=-=)634.2(sec ))2(cos )4))(2ln 4exp )22
2
2
2H j h G d H j G c j H H G b H j G a j j j j 双曲正割函数：
余弦平方函数：柯西函数：高斯函数：π 实际代码：
//零面积对称法寻峰(m=16,W=11,H=4) for( i=16;i<1007;i++) { for(j=0;j<16;j++) G[i]+=G[j]; d=1/11*G[i]; C[i]=G[i]-d;
}
for( i=16;i<1007;i++) { if((yi/dy)>f) feng[k++]=data[i];
}
③Wasson （瓦森、沃森）峰面积法： 计算总面积：∑-=+=
n
n
i i
i data
S 0
计算本底面积：
L L
R n data n L i L R data data b +----=
-)(0
L L
R n data n L i L
R data data b ++---=)(0
)21)((2
1
2)(++=++=--n b b n b b B n n n n
计算净峰面积：B S A -=
实际代码：
// Wasson （瓦森、沃森）峰面积法：
for (i=L;i<=R;i++)
{
RL=right-left;
n=(right-left)/2;
im=(right+left)/2;
S+=data[i];
BL=(data[L]-data[R])/RL*(im-L-n)+data[L];
BR=(data[L]-data[R])/RL*(im-L+n)+data[L];
B=(BL+BR)*(n+1/2);
A=S-B;
}
计算峰位能量
//根据能量刻度系数求峰位能量
/* for( i=0;i<1023;i++)
{
E[i]=0.017584*i-0.138133;
}
for( j=0;j<1024;j++)
{
if(E[dao[j]]<(6.403-6.403*0.1)||E[dao[j]]>(6.403+6.403*0.1)) ch[dao[j]]="Fe";
else if(E[dao[j]]<(7.477-7.477*0.1)||E[dao[j]]>(7.477+7.477*0.1))
ch[dao[j]]="Ni";
else if(E[dao[j]]<(8.047-8.047*0.1)||E[dao[j]]>(8.047+8.047*0.1))
ch[dao[j]]="Cu";
else if(E[dao[j]]<(8.638-8.638*0.1)||E[dao[j]]>(8.638+8.638*0.1))
ch[dao[j]]="Zn";
else if(E[dao[j]]<(10.5-10.5*0.1)||E[dao[j]]>(10.5+10.5*0.1))
ch[dao[j]]="Ap";
else if(E[dao[j]]<(11.725-11.725*0.1)||E[dao[j]]>(11.725+11.725*0.1)) ch[dao[j]]="As";
else if(E[dao[j]]<(12.611-12.611*0.1)||E[dao[j]]>(12.611+12.611*0.1)) ch[dao[j]]="Ni";
} */
8.调试运行
五、体会
在本次程序设计实习中，我再次学习了MFC编程方法。
过程中，我遇到了许多曾经都没有遇到的问题，比如文件的打开，文件按钮
虚函数的添加，绘图无效等。
在一步一步解决这些问题的同时，渐渐进一步的熟悉掌握了谱数据处理的方法和编程思想。
但是依旧有一个问题尚未解决，那就是数据类型的变化。
这个使得我无法在数据处理中正确地获取应有的结果。
令我费解的是成
谱图像却是与原始数据非常的吻合。
对于这个问题我依然百思不得
其解。
六、建议
我对本门课程的掌握程度也非常有限，提出以下建议供老师和学弟学妹们参考。
1.MFC编程时一个系统过程，需要掌握的细节可以随时网上学
习，但是总体思路要牢记，不容易临时掌握。
2.对于本专业同学们而言，这次实习是扩展我专业同学技能的大
好机会，建议不要轻易放弃，认真编写，即使最后没有达到自己预
期的计划也是收获颇丰的。
3.在编写的过程中要有良好的格式及注释习惯，以免一段时间甚
至几分钟后就看不懂自己的思路。
4.在编写算法过程中，一定要回归理论，弄明白实际数学过程，
对算法的编写很有帮助。
5.希望老师能够对非算法部分给予同学更多的帮助，教授同学建立软件的过程，指导同学完成非专业的部分，让同学着眼于数学导算法的翻译。
最后，衷心感谢在本次程序设计实习过程中，马永红老师细心的指导，您不厌其烦地解答我的各种疑问，并且帮助我寻找原始数据的数据格式和数据类型。
我还要感谢芦瑞琪、蒋卓辰同学对我在思路上的重要启发，谢谢你们给予我的大力支持！
附：程序包附在附件中，项目名称为37gai。