在区间[-1,1]上对函数，选取不同的插值节点构造插值多项式，比较他们的误差。

二维稳态导热的数值计算主要采用了热平衡法。
用差分法建立节点的热平衡方程，将节点所在的单元体的四个方向传递的热流密度，内热源在单元体产生的热流密度，根据能量守恒的原则建立方程，可以得到每一个节点的离散化代数方程。
进行数值计算的方法是：先设定初值，在根据初值对每一个节点进行迭代可以求得节点的值。
再将初值与新值进行比较，判断迭代的敛散性。
比较常用的迭代方法有两种：Gauss-Seidel法和Jacobi法。
Gaus-Seidel法每次迭代计算，均是使用节点温度的最新值。
Jacobi迭代法每次迭代计算均用上一次迭代计算出的值。
对于一个代数方程组，若选用的迭代方式不合适有可能导致迭代过程发散，而对于常物性导热问题组成差分方程组，每一个方程都选用导出方程的中心节点温度作为迭代变量则迭代一定收敛。

利用MATLAB强大的图形处理功能、符号运算功能和数值计算功能，实现离散系统的Z域分析的仿真波形。
1、用MATLAB绘制离散系统极零图，根据极零图分布观察系统单位响应的时域特性并分析系统的稳定性，至少以六个例子说明。
2、用MATLAB实现离散系统的频域特性分析：以二个实例分别代表低通，高通滤波器，绘出极零图，幅频特性，相频特性；
3、用MATLAB实现巴特沃兹滤波器分析；
4、撰写MATLAB课程设计说明书。

1.1单项选择题1.数据结构是一门研究非数值计算的程序设计问题中,数据元素的①、数据信息在计算机中的②以及一组相关的运算等的课程。
①A．操作对象　　　Ｂ．计算方法　　Ｃ．逻辑结构　　Ｄ．数据映象②A．存储结构Ｂ．关系Ｃ．运算Ｄ．算法2.数据结构DS(DataStruct)可以被形式地定义为DS=（D，R），其中D是①的有限集合，R是D上的②有限集合。
①A．算法Ｂ．数据元素Ｃ．数据操作Ｄ．数据对象②A．操作Ｂ．映象Ｃ．存储Ｄ．关系3.在数据结构中，从逻辑上可以把数据结构分成。
A．动态结构和静态结构Ｂ．紧凑结构和非紧凑结构Ｃ．线性结构和非线性结构Ｄ．内部结构和外部结构4.算法分析的目的是①，算法分析的两个主要方面是②。
①A.找出数据结构的合理性B.研究算法中的输入和输出的关系C.分析算法的效率以求改进D.分析算法的易懂性和文档性②A.空间复杂性和时间复杂性B.正确性和简明性C.可读性和文档性D.数据复杂性和程序复杂性5.计算机算法指的是①，它必具备输入、输出和②等五个特性。
①A.计算方法B.排序方法C.解决问题的有限运算序列D.调度方法②A.可行性、可移植性和可扩充性B.可行性、确定性和有穷性C.确定性、有穷性和稳定性D.易读性、稳定性和安全性1.2填空题（将正确的答案填在相应的空中）1.数据逻辑结构包括、、和四种类型，树形结构和图形结构合称为。
2.在线性结构中，第一个结点前驱结点，其余每个结点有且只有个前驱结点；
最后一个结点后续结点，其余每个结点有且只有个后续结点。
3.在树形结构中，树根结点没有结点，其余每个结点有且只有个直接前驱结点，叶子结点没有结点，其余每个结点的直接后续结点可以。
4.在图形结构中，每个结点的前驱结点数和后续结点数可以。
5.线性结构中元素之间存在关系，树形结构中元素之间存在关系，图形结构中元素之间存在关系。
6.算法的五个重要特性是____,____,____,____,____。
7.分析下面算法（程序段），给出最大语句频度，该算法的时间复杂度是____。
for(i=0;i＜n;i++)for(j=0;j＜n;j++)A[i][j]=0;8.分析下面算法（程序段），给出最大语句频度，该算法的时间复杂度是____。
for(i=0;i＜n;i++)for(j=0;j＜i;j++)A[i][j]=0;9.分析下面算法（程序段），给出最大语句频度，该算法的时间复杂度是____。
s=0;for(i=0;i＜n;i++)for(j=0;j＜n;j++)for(k=0;k＜n;k++)s=s+B[i][j][k];sum=s;10.分析下面算法（程序段）给出最大语句频度，该算法的时间复杂度是____。
inti=0，s=0;while(s＜n){i++;s+=i;//s=s+i}11.分析下面算法（程序段）给出最大语句频度，该算法的时间复杂度是____。
i=1;while(i＜=n)i=i*2;

中国石油大学《数值计算方法》历年期末试卷（含答案）

基于龙贝格算法的积分，可计算一维和二维积分学数值计算的兄弟们可以看看精度可以自己设定具体文件里有说明

研究了面向海洋应用的光纤法布里-珀罗高压传感器，通过建立有限元数值模型对传感器满量程腔长变化量进行分析。
数值仿真显示，有限元模型的满量程腔长变化量处于固支模型和简支模型之间，且随着法布里-珀罗腔半径的减小和硅膜片厚度的增加而偏离固支模型。
引入固支边界条件偏离度β对偏离程度进行量化分析。
制作了三种不同规格的传感器进行压力实验研究。
实验结果显示，实际测量得到的传感器芯片满量程腔长变化量与有限元数值计算的结果基本吻合，使用该有限元模型设计传感器芯片可将满量程腔长变化量误差降低到13.4%以下。
传感器最大量程达到105MPa，满量程测量精度均优于0.100%。

通过软件，可以很方便的计算出采集AD，通过程序查表，就能得出温度。
超级方便。

《ANSYS_LS_DYNA模拟鸟撞飞机风挡的动态响应》鸟撞问题在飞机设计中至关重要，尤其是在飞机起飞和降落时，高速运动的飞机与鸟类相撞可能导致严重损伤，甚至造成机毁人亡的灾难。
特别是飞机的前风挡部分，由于迎风面积大，成为鸟撞概率较高的区域，而风挡玻璃的强度相对较低，因此对风挡受鸟撞冲击的模拟分析显得尤为必要，以提升飞行安全性。
早期的抗鸟撞设计主要依赖实验方法，但随着计算机技术和有限元数值计算理论的发展，现在越来越多地采用数值计算来分析鸟撞问题。
目前的有限元模型主要分为解耦解法和耦合解法。
解耦解法将鸟撞冲击力作为已知条件，单独求解风挡的动态响应，但鸟撞载荷模型的不确定性会影响求解精度。
耦合解法则考虑碰撞接触，通过协调鸟体与风挡接触部位的条件，联合求解，能更直观地模拟整个鸟撞过程。
本文采用ANSYS_LS_DYNA软件，建立鸟撞风挡的三维模型，研究鸟撞风挡的动态响应特征。
在建立有限元模型时，使用ANSYS软件，简化了计算过程，忽略了对风挡动态响应影响不大的结构因素，如机身、后弧框和铆钉等，将其替换为边界固定。
风挡结构为圆弧形，材料为特定型号的国产航空玻璃，鸟撞击点设在风挡中部，撞击角度为29°。
选用LS-DYNA材料库中的塑性动力学材料模型，破坏准则设定为最大塑性应变失效模式，当材料塑性应变达到5%时材料破坏。
鸟体的模拟是鸟撞分析的一大挑战，由于真实鸟体的本构特性难以准确描述，通常采取弹性体、弹塑性体或理想流体等简化模型。
本文中，鸟体被简化为质量1.8kg、直径14cm的圆柱体，材料选用弹性流体模型。
计算结果显示，当鸟撞速度达到540km/h（相对于风挡的绝对速度）时，风挡的后弧框处有效塑性应变达到5%，风挡破坏。
据此，计算得出风挡的安全临界速度为150m/s。
在这一速度下，风挡后弧框处首先发生破坏，成为结构弱点。
撞击时的最大应力主要集中在后弧框及其下方，而非撞击点。
此外，鸟撞还会导致风挡结构产生位移。
风挡下方通常布置有精密仪器，因此必须考虑鸟撞引起的位移情况。
鸟体撞击后在风挡上滑行，挤压风挡表面，产生较大位移。
计算表明，在150m/s的撞击速度下，最大位移可达38mm，位于撞击点和后弧框之间。
风挡表面位移随着时间呈现出先向下位移，然后因弯曲波反弹而振荡的行为。
总结来说，鸟撞风挡的最危险区域位于后弧框及其下方。
不同结构的风挡有不同的鸟撞安全临界速度、最大位移和撞击时间。
对于本文的风挡模型，临界速度为450km/h，最大位移为38mm，撞击时间约为7ms。
这些分析结果对于飞机设计改进和飞行安全性的提升具有重要指导意义。

用Python实现蒙特卡洛模拟(MonteCarlosimulation)，蒙特·卡罗方法（MonteCarlomethod），也称统计模拟方法，是二十世纪四十年代中期由于科学技术的发展和电子计算机的发明，而被提出的一种以概率统计理论为指导的一类非常重要的数值计算方法。

在日常工作中，钉钉打卡成了我生活中不可或缺的一部分。然而，有时候这个看似简单的任务却给我带来了不少烦恼。 每天早晚，我总是得牢记打开钉钉应用，点击"工作台"，再找到"考勤打卡"进行签到。有时候因为工作忙碌，会忘记打卡，导致考勤异常，影响当月的工作评价。而且，由于我使用的是苹果手机，有时候系统更新后，钉钉的某些功能会出现异常，使得打卡变得更加麻烦。 另外，我的家人使用的是安卓手机，他们也经常抱怨钉钉打卡的繁琐。尤其是对于那些不太熟悉手机操作的长辈来说，每次打卡都是一次挑战。他们总是担心自己会操作失误，导致打卡失败。 为了解决这些烦恼，我开始思考是否可以通过编写一个全自动化脚本来实现钉钉打卡。经过一段时间的摸索和学习，我终于成功编写出了一个适用于苹果和安卓系统的钉钉打卡脚本。