实验三 常微分方程的差分方法实验
一. 实验目的
（1）深入理解常微分方程的差分方法的原理，学会用差分方法解决某些实际的常微分方程问题，比较这些方法解题的不同之处。
（2）熟悉Matlab 编程环境，利用Matlab 实现具体的常微分方程。
二. 实验要求
用Matlab 软件实现欧拉方法、改进的欧拉方法、龙格-库塔方法和亚当姆斯方法，并用实例在计算机上计算。
三. 实验内容
1. 实验题目
（1）取h =0.1，用欧拉方法、改进的欧拉方法、四阶龙格-库塔方法求解初值问题：
?????=-+='0
)0(21122y y x y 40≤≤x 并与精确解2
1x x y +=比较计算结果 A.欧拉方法：
a ．编写文件Euler.m,内容如下所示：
b.B
b.编写文件
f3.m:
c.编写文件
solvef3(x):
d.运行如下所示：
可见精度是比较粗糙的，具有两位有效数字；
B．改进的欧拉方法
a．编写文件MendEuler.m,内容如下所示：
c．运行如下：
可见部分值具有3为有效数字了
C．四阶龙格-库塔方法
a．编写文件Rungkuta4.m,内容如下所示：
b．运行如下所示：
（2）分别用二阶亚当姆斯预估校正系统、改进的四阶亚当姆斯预估校正系统求解初
值问题0)0(1=-='y y y ，，取181.0)2.0(2.0==y h ，，计算)0.1(y 。
a.编写文件Adams2PC.m,即为二阶亚当姆斯预估校正系统的程序：
b.编写改进的四阶亚当姆斯预估校正系统的程序，如下所示：
C ．编写文件f5.m,即为微分方程的程序；编写文件solvef5.m ，为解函数的程序，用以比较结果：
d．运行
从上面的输出可以看到CA4(:,2)和CA4的输出时不同的，前面一个值输入矩阵的第二列；
2. 设计思想
要求针对上述题目，详细分析每种算法的设计思想。
首先差分方法是一类离散化方法，这类方法将寻求解y（x）的分析问题转化为计算离散值y的代数问题，使问题得到实质性的简化；差分方法在解初值问题时采取步进式，进一步将计算模型化归为仅含一个变元的代数方程——差分格式
A．欧拉方法：确定步长，施行离散化手续，用差商代替倒数；
B．改进的欧拉方法：在欧拉方法的基础上，改变倒数项，通过首尾两点的倒数的一半来作为平均倒数，从而改善精度；
C．龙格-库塔方法：在每一步上多预报几个点的斜率，然后将它们加权平均作为平均斜率K*，即该平均斜率将比上面两种方法更加精确，则有可能构造出更高精
度的计算格式；
D．二阶亚当姆斯预估校正系统：充分利用先前计算出来的老信息来为新信息服务，期望来达到减少计算量的目的；预报过程采用松弛手续，来校正产生的值，使
值的精度更高；
E．改进的四阶亚当姆斯预估校正系统：在二阶的基础上，更近一步地改善精度，更加充分地利用来信息，又不失二阶时的预报校正的优点。
四.实验体会
对实验过程进行分析总结，对比求解常微分方程的不同方法，指出每种算法的设计要点及应注意的事项，以及自己通过实验所获得的对常微分方程的差分方法的理解。
在实验的过程也是一次对书本进行复习的过程，通过实验可以更加地理解和学会利用已知的结论来求解初值问题。算法的计算是建立在使用价值上的，如上面的设计思想中可以看出，每一种算法的建立都是在满足一定的实用价值上的，既要有一定的精度，也要控制在一定的计算范围之内。利用松弛手续来提高精度，这种设计是优秀的。
编写程序时，以前并不是很注意结尾要不要分号，但这个细节有时必须要注意，正如在运行窗口下时，要想得出计算结果是不能带分号的，假如你带了分号就会自动换行，
提示继续输入，如，而当去了分号后就可以运行了；当然这个道理前面几个实验已经知道；这次实验的情况是微分方程是另外存在一个函数文件中的，那么这个文件的结尾到底要不要分号呢？从习惯上还是提倡写上分号，假如在这次中忘了写的话，在其它“.m”文件中引用该函数时，虽然只是中间结果，但每用一次就好输出一次，使输出内容与预知的多了些；
还有一点是，像上面的内容假如直接复制，粘贴再运行会出错，因为把“>>”也复制上了，粘贴的“>>”和换行时自动出现的符号在意义上是不同的，故会出错，不过在最新版本的matlab上就没有这个问题，在粘贴时自动把多出的符号删去。
实验中二阶亚当姆斯预估校正系统题目中给了两个初值，从而使得程序更加地简单，不必利用其它的方法来求得初值再启动运行。
《计算方法》实验报告 姓名： 班级： 学号： 实验日期: 2011年10月26日
一、实验题目： 数值积分 二、实验目的： 1．熟悉matlab 编写及运行数值计算程序的方法。 2．进一步理解数值积分的基础理论。 3．进一步掌握应用不同的数值积分方法求解给定的积分并给出数据结果及误差分析。 三、实验内容： 1．分别用复合梯形求积公式及复合辛普森求积公式计算积分xdx x ln 10 ? , 要求计算精度达到410-,给出计算结果并比较两种方法的计算节点数. 2．用龙贝格求积方法计算积分dx x x ?+3 021,使误差不超过510-. 3．用3=n 的高斯-勒让德公式计算积分?3 1 sin x e x ,给出计算结果. 4．用辛普森公式(取2==M N ) 计算二重积分.5 .00 5 .00 dydx e x y ? ? - 四、实验结果： 1.（1）复合梯形法: 将区间[a,b]划分为n 等份，分点n k n a b h kh a x k ,2,1,0,,=-=+=在每个区间[1,+k k x x ](k=0,1,2,···n-1)上采用梯形公式，则得 )()]()([2)()(1 11 1 f R x f x f h dx x f dx x f I n n k k k b a n k x x k k ++===∑?∑? -=+-=+ 故)]()(2)([21 1 b f x f a f h T n k k n ++=∑-=称为复合梯形公式 计算步长和划分的区间 Eps=1E-4 h1=sqrt(Eps/abs(-(1-0)/12*1/(2+1))) h1 =0.0600 N1=ceil(1/h1) N1 =17 用复合梯形需要计算17个结点。 复合梯形： function T=trap(f,a,b,n) h=(b-a)/n;
《计算方法》上机实验报告 班级：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 四、运行结果： 实验二：
DES 差分分析实验报告 四大队四队五班 刘杰 一、实验目的 差分密码分析是一种选择明文攻击，是现代分组密码分析的重要方法之一，也是理论分析密码算法和算法抗攻击测试的重要依据之一。本实验通过3轮DES 简化算法的差分分析来达到加深学员对差分分析方法原理的理解和利用该原理分析实际问题的操作能力。 二、实验内容 （1）3轮DES 简化算法的差分分析； （2）通过三组明密文对（每组两个相关明文和相应密文），利用差分原理提取密钥。 明 文 密 文 748502CD38451097 03C70306D8A09F10 3874756438451097 78560A0960E6D4CB 486911026ACDFF31 45FA285BE5ADC730 375BD31F6ACDFF31 134F7915AC253457 357418DA013FEC86 D8A31B2F28BBC5CF 12549847013FEC86 0F317AC2B23CB944 三、实验原理 设DES 两个明密文对:=00m L R ***=00m L R =33c L R *** =33c L R 计算过程: (,)(,)(,)(,)=⊕=⊕=⊕⊕322312300123R L f R k R f R k L f R k f R k
(,)(,)****=⊕⊕300123R L f R k f R k 令：*'=⊕000L L L (,)(,)(,)(* **''=⊕=⊕⊕⊕⊕333001012323R R R L f R k f R k f R k f R k 观察得：在本次实验原始数据中，明文对*=00R R ，即* '=⊕=00000000000R R R 则(,)(,)** ''=⊕=⊕⊕33302323R R R L f R k f R k 同时有：=00m L R ***=00m L R =23R L ** =23R L 则可计算出：*'=⊕000L L L *'=⊕333R R R (,)(,)* ''⊕=⊕232330f R k f R k R L 则可得出： S 盒输入差：(())(())()()* *⊕⊕⊕=⊕232333E R k E R k E L E L S 盒输出差：()*-''⊕=⊕13 0D D P R L 分析过程： 令：()()*⊕=3312345678E L E L B B B B B B B B ()-''⊕=13 012345678P R L C C C C C C C C ()=312345678E L A A A A A A A A =312345678 k J J J J J J J J ()⊕=3312345678E L k X X X X X X X X *()⊕=3312345678E L k Y Y Y Y Y Y Y Y 基本思路：（分别计算12345678J J J J J J J J ） {|,()()∈=⊕⊕=⊕=i i i i i i i J T e s t x A x y B S x S y C ,,,,,,,=12345678i 对于本次实验的3个具有明文差（*,0）的明密文对，则可构造上面的3个 Test 集合，显然 ()()( )∈12 i i i i J Test Test Test t ,,,,,,,=12345678i 一种确定Ji 的直接方法： 1.建立26=64长度的数组J[64]={0}; 2.对Testi(r)，r = 1,2,…,t ，若a ∈Testi(r)，则 J[a] = J[a] + 1。 3.若J[b] =3，则6比特串b 就是可能的密钥比特 Ji 。 四、实验环境 Microsoft visual c++ 五、实验步骤 （1）计算简化算法第3轮S 盒输入差
本科实验报告 课程名称：计算机数值方法 实验项目：方程求根、线性方程组的直接解 法、线性方程组的迭代解法、代数插值和最 小二乘拟合多项式 实验地点：行勉楼 专业班级： ******** 学号： ********* 学生姓名： ******** 指导教师：李誌，崔冬华 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. 二力平衡的条件：作用在同一个物体上的两个力，如果大小相等，方向相反，并且在同一直线上，这两个力就平衡。 2. 在平衡力的作用下，静止的物体保持静止状态，运动的物体保持匀速直线运动状态。 3. 两个相互接触的物体，当它们做相对运动时或有相对运动的趋势时，在接触面上会产生一种阻碍相对运动的力，这种力就叫摩擦力。 4. 弹簧测力计拉着木块在水平面上做匀速直线运动时，拉力的大小就等于摩擦力的大小，拉力的数值可从弹簧测力计上读出，这样就测出了木块与水平面之间的摩擦力。
探究导引 探究指导： 关闭发动机的列车会停下来，自由摆动的秋千会停下来，踢出去的足球会停下来，运动的物体之所以会停下来，是因为受到了摩擦力。 运动物体产生摩擦力必须具备以下三个条件：1.物体间要相互接触，且挤压;2.接触面要粗糙;3.两物体间要发生相对运动或有相对运动的趋势。三个条件缺一不可。 摩擦力的作用点在接触面上，方向与物体相对运动的方向相反。由力的三要素可知：摩擦力除了有作用点、方向外，还有大小。 提出问题：摩擦力大小与什么因素有关? 猜想1：摩擦力的大小可能与接触面所受的压力有关。 猜想2：摩擦力的大小可能与接触面的粗糙程度有关。 猜想3：摩擦力的大小可能与产生摩擦力的两种物体间接触面积的大小有关。 探究方案： 用弹簧测力计匀速拉动木块，使它沿长木板滑动，从而测出木块与长木板之间的摩擦力;改变放在木块上的砝码，从而改变木块与长木板之间的压力;把棉布铺在长木板上，从而改变接触面的粗糙程度;改变木块与长木板的接触面，从而改变接触面积。 物理实验报告 .化学实验报告 .生物实验报告 .实验报告格式 .实验报告模板 探究过程： 1. 用弹簧测力计匀速拉动木块，测出此时木块与长木板之间的摩擦力：0.7N 2. 在木块上加50g的砝码，测出此时木块与长木板之间的摩擦力：0.8N 3. 在木块上加200g的砝码，测出此时木块与长木板之间的摩擦力：1.2N 4. 在木板上铺上棉布，测出此时木块与长木板之间的摩擦力：1.1N 5. 加快匀速拉动木块的速度，测出此时木块与长木板之间的摩擦力：0.7N 6. 将木块翻转，使另一个面积更小的面与长木板接触，测出此时木块与长木板之间的摩擦力：0.7N 探究结论：
工程电磁场 实验报告 ——有限差分法
用超松弛迭代法求解 接地金属槽内电位的分布 一、实验要求 按对称场差分格式求解电位的分布 已知： 给定边值：如图1-7示 图1-7接地金属槽内半场域的网格 给定初值)()(.1j 40 100 1j p 1 2j i -= --= ??? 误范围差: 510-=ε 计算：迭代次数N ，j i ,?，将计算结果保存到文件中 二、实验思想 有限差分法 有限差分法（Finite Differential Method ）是基于差分原理的一种数值计算法。其基本思想：将场域离散为许多小网格，应用差分原理，将求解连续函数?的泊松方程的问题转换为求解网格节点上? =?= V 100 ? 0 =?0 =?
的差分方程组的问题。 泊松方程的五点差分格式 )(4 1 4243210204321Fh Fh -+++=?=-+++?????????? 当场域中,0=ρ得到拉普拉斯方程的五点差分格式 )(4 1 044321004321??????????+++=?=-+++ 差分方程组的求解方法(1) 高斯——赛德尔迭代法 ][)(,)(,)(,)(,)(,2 k 1j i k j 1i 1k 1j i 1k j 1i 1k j i Fh 4 1 -+++=+++-+-+????? (1-14) 式中：??????=??????=,2,1,0,2,1,k j i ， ? 迭代顺序可按先行后列，或先列后行进行。 ? 迭代过程遇到边界节点时，代入边界值或边界差分 格式，直到所有节点电位满足ε??<-+)(,)(,k j i l k j i 为止。 （2）超松弛迭代法 ][) (,)(,)(,)(,)(,)(,)(,k j i 2k 1j i k j 1i 1k 1j i 1k j 1i k j i 1k j i 4Fh 4 ?????α??--++++=+++-+-+ (1-15) 式中：α——加速收敛因子)21(<<α 可见：迭代收敛的速度与α有明显关系 三、程序源代码 #include
上海交通大学致远学院 《常微分方程和偏微分方程的数值解法》教学大纲 一、课程基本信息 课程名称（中文）：常微分方程和偏微分方程的数值解法 课程名称（英文）：Numerical Methods for Ordinary and Partial Differential Equations 课程代码：MA300 学分 / 学时：4学分 / 68学时 适用专业：致远学院与数学系相关专业 先修课程：偏微分方程，数值分析 后续课程：相关课程 开课单位：理学院数学系计算与运筹教研室 Office hours: 每周二19：00—21：00，地点：数学楼1204 二、课程性质和任务 本课程是致远学院和数学系应用数学和计算数学方向的一门重要专业基础课程，其主要任务是通过数学建模、算法设计、理论分析和上机实算“四位一体”的教学方法，使学生掌握常微分方程与偏微分方程数值解的基本方法、基本原理和基本理论，进一步提升同学们利用计算机解决实际问题的能力。在常微分方程部分，将着重介绍常微分方程初值问题的单步法，含各类Euler方法和Runge-Kutta方法，以及线性多步法。将简介常微分方程组和高阶常微分方程的数值方法。在偏微分方程部分，将系统介绍求解椭圆、双曲、抛物型方程的差分方法的构造方法和理论分析技巧，对于椭圆型方程的边值问题将介绍相应变分原理与有限元方法。将在课堂上实时演示讲授的核心算法的计算效果，以强调其直观效果与应用性。本课程重视实践环节建设，学生要做一定数量的大作业。 三、教学内容和基本要求 第一部分：常微分方程数值解法 1 引论 1.1回顾：一阶常微分方程初值问题及解的存在唯一性定理
简单计算器 姓名: 周吉祥 实验目的：模仿日常生活中所用的计算器，自行设计一个简单的计算器程序，实现简单的计算功能。 实验内容： (1)体系设计： 程序是一个简单的计算器，能正确输入数据，能实现加、减、乘、除等算术运算，运算结果能正确显示，可以清楚数据等。 (2)设计思路： 1)先在Visual C++ 6.0中建立一个MFC工程文件，名为 calculator. 2)在对话框中添加适当的编辑框、按钮、静态文件、复选框和单 选框 3)设计按钮，并修改其相应的ID与Caption. 4)选择和设置各控件的单击鼠标事件。 5)为编辑框添加double类型的关联变量m_edit1. 6)在calculatorDlg.h中添加math.h头文件，然后添加public成 员。 7)打开calculatorDlg.cpp文件，在构造函数中，进行成员初始 化和完善各控件的响应函数代码。 (3)程序清单：
●添加的public成员： double tempvalue; //存储中间变量 double result; //存储显示结果的值 int sort; //判断后面是何种运算：1.加法2.减法3. 乘法 4.除法 int append; //判断后面是否添加数字 ●成员初始化： CCalculatorDlg::CCalculatorDlg(CWnd* pParent /*=NULL*/) : CDialog(CCalculatorDlg::IDD, pParent) { //{{AFX_DATA_INIT(CCalculatorDlg) m_edit1 = 0.0; //}}AFX_DATA_INIT // Note that LoadIcon does not require a subsequent DestroyIcon in Win32 m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME); tempvalue=0; result=0; sort=0; append=0; }
重庆工商大学 《统计学》实验报告 实验课程：统计学 _ 指导教师：陈正伟 _ 专业班级： 08 经济学 学生姓名：程剑波 学生学号： 2008011133 __
实验项目 实验日期实验地点80608 实验目的掌握统计学的基本计算方法和分析方法。 实验内容一、统计图绘制；二、动差、偏度系数、峰度系数的计算；三、趋势性的绘制； 四、相关分析与回归分析；五、时间数列的动态指标分析；六、循环变动的测 算分析。 通过统计学（2009.9.10-2009.12.15）实验报告如下： 一、统计图绘制； （一）过程： （二）结果： （三）分析： 二、动差、偏度系数、峰度系数的计算； （一）过程： （二）结果： （三）分析： 三、趋势性的绘制； （一）过程： （二）结果： （三）分析： 四、相关分析与回归分析； （一）过程： （二）结果： （三）分析：
五、时间数列的动态指标分析 （一）过程： （二）结果： （三）分析： 六、循环变动的测算分析。 （一）过程： （二）结果： （三）分析： 体会： 参考实验报告： 重庆工商大学数学与统计学院 综合评价方法及应用 实验报告
实验课程：非参数统计 _ 指导教师：陈正伟 _ 专业班级： 06市调2班 学生姓名：何春 学生学号： 2006004151 _
实验报告一 实验项目变异系数法相关系数法熵值发坎蒂雷法 实验日期2009-4-30 实验地点80608 实验目的 通过本实验本要求掌握综合评价指标体系中各个指标重要性权数的重要意义；掌握权数确定的定性和定量技术和技能；解决实际综合评价中重要性权数确定的处理技能。 实验内容 根据资料使用变异系数法、相关系数法、熵值法和坎蒂雷方法分别确定各个指标的权数。并进行权数比较分析。 检验方法的选择及实验步骤及结果： 1用变异系数求各个指标的权数： 基本步骤：（1）先求各个指标的均值Xi 和标准差 Si (2)接着求各个指标的变异系数Vi=Si/Xi (3)对Vi作作归一化处理，及得各个指标的权数 结果如下： 从这个表中可以看到最后一列的权数最大，即人均创造总收入这个指标在这项评价上的分辨信息丰富，这个指标的数值能明确区分开各个评价被评价对象差异。同理，第四列的权数最小，也就是说各个被评价对象在某项指标上的数值差异较小，那么这项指标区分开各评价对象的能力较弱。 2 用相关系数法求各个指标的权数： 基本步骤：（1）计算各个指标之间的相关系数矩阵 （2）构造分块矩阵 R1(去掉相关系数矩阵的第一行和第一列)R2 R3 R4 R5 R6 同理可得
实验十三差分编译码实验 一、实验目的 掌握差分编码/译码原理 二、实验内容 1、学习差分编译码原理 2、用示波器观察差分编码结果和译码结果 三、基本原理 差分码是一种把符号‘0’和‘1’反映在相邻码元的相对变化上的波形。比如，若以相邻码元的电位改变表示符号‘1’，而以电位不改变表示符号‘0’，如图13－1所示。当然，上述规定也可以反过来。由图可见，这种码波形在形式上与单极性或双极性码波形相同，但它代表的信息符号与码元本身电位或极性无关，而仅与相邻码元的电位变化有关。差分波形也称相对码波形，而相应地称单极性或双极性波形为绝对码波形。差分码波形常在相位调制系统的码变换器中使用。 图13-1差分码波形 组成模块如下图所示： cclk d_out 端口说明： CCLK：编码时钟输入端 DIN：编码数据输入端 Diff－OUT：差分编码结果输出端 DCLK：译码时钟输入端
Diff－IN：差分译码数据输入端 DOUT：译码结果输出端 四、实验步骤 1、实验所用模块：数字编解码模块、数字时钟信号源模块。 实验连线： CCLK：从数字时钟信号源模块引入一高频时钟，如512K。 DIN：从数字时钟信号源模块引入一低频时钟，如16K。 DIFF－OUT与DIFF－IN短接。 DCLK与CCLK短接。 2、用示波器两探头同时观测DIN与DIFF－OUT端，分析差分编码规则。 3、用示波器两探头同时观测DIN与DOUT端，分析差分译码结果。 五、实验报告要求 设信息代码为1001101，码速率为128K，差分码的编码时钟为码速率的四倍，根据实验观察得到的规律，画出差分码波形。
计算方法实验报告格式 小组名称: 组长姓名(班号)： 小组成员姓名(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一个完整的实验，应包括数据准备、理论基础、实验内容及方法，最终对实验结果进行分析，以达到对理论知识的感性认识，进一步加深对相关算法的理解，数值实验以实验报告形式完成，实验报告格式如下： 一、实验名称 实验者可根据报告形式需要适当写出． 二、实验目的及要求 首先要求做实验者明确，为什么要做某个实验，实验目的是什么，做完该实验应达到什么结果，在实验过程中的注意事项，实验方法对结果的影响也可以以实验目的的形式列出． 三、算法描述（实验原理与基础理论） 数值实验本身就是为了加深对基础理论及方法的理解而设置的，所以要求将实验涉及到的理论基础，算法原理详尽列出． 四、实验内容 实验内容主要包括实验的实施方案、步骤、实验数据准备、实验的算法以及可能用到的仪器设备． 五、程序流程图 画出程序实现过程的流程图，以便更好的对程序执行的过程有清楚的认识，在程序调试过程中更容易发现问题． 六、实验结果 实验结果应包括实验的原始数据、中间结果及实验的最终结果，复杂的结果可以用表格
形式列出，较为简单的结果可以与实验结果分析合并出现． 七、实验结果分析 实验结果分析包括对对算法的理解与分析、改进与建议． 数值实验报告范例 为了更好地做好数值实验并写出规范的数值实验报告，下面给出一简单范例供读者参考． 数值实验报告 小组名称: 小组成员(班号): 按贡献排序情况: 指导教师评语: 小组所得分数: 一、实验名称 误差传播与算法稳定性． 二、实验目的 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，则
计算方法实验报告（四） 方程和方程组的迭代解法 一、实验问题 利用简单迭代法，两种加速技术，牛顿法，改进牛顿法，弦割法求解习题5-1，5-2，5-3中的一题，并尽可能准确。 选取5-3：求在x=1.5附近的根。 二、问题的分析（描述算法的步骤等） （1）简单迭代法算法： 给定初始近似值,求的解。 Step 1 令i=0； Step 2 令(计算)； Step 3 如果,则迭代终止，否则重复Step 2。 （2）Aitken加速法算法 Step 1 令k=0，利用简单迭代算法得到迭代序列； Step 2 令-(计算得到一个新的序列，其中k=0,1,2…)；Step 3 如果,则迭代终止，否则重复Step 2。 (3)插值加速法算法 Step 1 令k=0，利用简单迭代算法得到迭代序列； Step 2 令+(计算得到一个新的序列，其中k=1,2,3…)； Step 3 如果,则迭代终止，否则重复Step 2。 (4)牛顿法算法
Step 1给定初始近似值； Step 2令,其中k计算得到的序列； Step 3如果,则迭代终止，否则重复Step 2。 （5）改进牛顿法的算法 Step 1给定初始近似值； Step 2令,其中k迭代计算得到的序列； Step 3如果,则迭代终止，否则重复Step 2。 （6）弦割法算法（双点弦割法） Step 1给定初始近似值，； Step 2令其中k计算得到的序列； Step 3如果,则迭代终止，否则重复Step 2。 三、程序设计 (1)简单迭代法 利用迭代公式进行迭代运算。 #include
实验报告 课程名称：计算方法 院系：数学科学系 专业班级：数应1001 学号：1031110139 学生姓名：姚海保 指导教师：沈林 开课时间：2012至2013学年第一学期
一、学生撰写要求 按照实验课程培养方案的要求，每门实验课程中的每一个实验项目完成后，每位参加实验的学生均须在实验教师规定的时间内独立完成一份实验报告，不得抄袭，不得缺交。 学生撰写实验报告时应严格按照本实验报告规定的内容和要求填写。字迹工整，文字简练，数据齐全，图表规范，计算正确，分析充分、具体、定量。 二、教师评阅与装订要求 1.实验报告批改要深入细致，批改过程中要发现和纠正学生实验报告中的问题，给出评语和实验报告成绩，签名并注明批改日期。实验报告批改完成后，应采用适当的形式将学生实验报告中存在的问题及时反馈给学生。 2.实验报告成绩用百分制评定，并给出成绩评定的依据或评分标准（附于实验报告成绩登记表后）。对迟交实验报告的学生要酌情扣分，对缺交和抄袭实验报告的学生应及时批评教育，并对该次实验报告的分数以零分处理。对单独设课的实验课程，如学生抄袭或缺交实验报告达该课程全学期实验报告总次数三分之一以上，不得同意其参加本课程的考核。 3.各实验项目的实验报告成绩登记在实验报告成绩登记表中。本学期实验项目全部完成后，给定实验报告综合成绩。 4.实验报告综合成绩应按课程教学大纲规定比例（一般为10-15%）计入实验课总评成绩；实验总评成绩原则上应包括考勤、实验报告、考核（操作、理论）等多方面成绩； 5.实验教师每学期负责对拟存档的学生实验报告按课程、学生收齐并装订，按如下顺序装订成册：实验报告封面、实验报告成绩登记表、实验报告成绩评定依据、实验报告（按教学进度表规定的实验项目顺序排序）。装订时统一靠左侧按“两钉三等分”原则装订。
南京信息工程大学实验（实习）报告 一、实验目的: 用最小二乘法将给定的十个点拟合成三次多项式。 二、实验步骤： 用matlab编制以函数为基的多项式最小二乘拟合程序，并用于对下列数据作三次多项式最小二乘拟合(取权函数wi=1) x -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 y -2.30 -1 -0.14 -0.25 0.61 1.03 1.75 2.75 4.42 6.94 给定直线方程为：y=1/4*x3+1/2*x2+x+1 三、实验结论： 最小二乘法：通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。最小二乘法还可用于曲线拟合。 一般地。当测量数据的散布图无明显的规律时，习惯上取n次代数多项式。 程序运行结果为： a = 0.9731 1.1023 0.4862 0.2238 即拟合的三次方程为：y=0.9731+1.1023x+0.4862*x2+0.2238*x3
-2.5 -2-1.5-1-0.5 00.51 1.52 2.5 -4-20246 81012 x 轴 y 轴 拟合图 离散点 y=a(1)+a(2)*x+a(3)*x.2+a(4)*x.3 结论： 一般情况下，拟合函数使得所有的残差为零是不可能的。由图形可以看出最小二乘解决了残差的正负相互抵消的问题，使得拟合函数更加密合实验数据。 优点：曲线拟合是使拟合函数和一系列的离散点与观测值的偏差平方和达到最小。 缺点：由于计算方法简单，若要保证数据的精确度，需要大量的数据代入计算。
实验报告的书写格式模版 有关实验报告的书写格式 一、完整实验报告的书写 完整的一份实验报告一般包括以下项目：实验名称： 实验目的： 实验器材： 实验原理： 实验步骤： 实验数据记录（表格）及处理： 实验结论（结果推导）： 实验讨论或分析等。 二、实验报告书写方法 1、实验名称：就是这个实验是做什么的。 2、实验目的：一般都写掌握什么方法啊；了解什么啊；知道什么啊；会什么啊；…… 等。 3、实验器材：就是做这个实验需要的所有器材（仪器）。 4、实验原理：就是这个实验是根据什么来做的，一般书上会写，抄一下也就可以啦。 5、实验步骤：就是你做实验的过程，开始操作时，（1）做什么； （2）做什么；（3）做什么；……
6、实验数据记录（表格）及处理：根据实验中涉及以及实验得到的数据，设计表格，将有关数据填在表格相应的位置；数据处理，就是该计算的，按要求计算后填入表格对应位置。 7、实验结论（结果推导）：就是做这个实验要得到的结果。 8、分析于讨论：写你的实验结果是否适合真实值？如果有误差要分析产生误差的原因，还有实验的一些比较关键的步骤的注意事项等。 对于初中生或小学生来说，书写的实验报告也可简单一点，有时也可不要分析于讨论，也可不写实验原理等。 三、探究实验书写一般有七个环节 1．提出问题：就是在生活中发现、提出问题。 2．猜想与假设：发现问题，就要弄清楚问题，在没有搞清楚之前总有基本的猜测和设想，这就是猜想与假设。 3．制定计划与设计实验：有了猜想，就有了实验的目的，再根据实验的目的设计实验方案，制定实验计划，包括取得证据的途径和方法，确定收集证据的范围。包括实验的理论依据（实验原理）、实验器材、实验步骤等。 4．进行实验与收集证据：上一步是动脑、思维活动，这一步是手脑并用的实验过程。 5．分析与论证：通过上面的实验，收集到一些数据，观察到一些现象，对其分析，得出事实与假设的关系，通过归纳、概括等方法，得到结论。
学生实验报告实验课程名称偏微分方程数值解 开课实验室数统学院 学院数统年级2013 专业班信计02班 学生姓名______________ 学号 开课时间2015 至2016 学年第 2 学期
数学与统计学院制 开课学院、实验室：数统学院实验时间：2016年6月20日
1、三层显格式建立 由于题中h 0.1, 0.1h，x 0,1 ,t 0,2，取N 10, M 200，故令网比r 0.1，h X j j h, j 0,1,2,L 10，t k k ,k O,1L ,200 ,在内网个点处，利用二阶中心差商得到如下格式: k 1 k U J 2U J 2- k 1 U j k k U j 1 2U j h2 k U j 1 o h2 略去误差项得到: k 1 U j 其中j 1,2丄9,k 对于初始条件 2 k r U J1 1,2,L ,199，局部截断误差为 U x,0 sin U J k U j k r U j 2 o k 1 U J h2。 (3) 对于初始条件-u x,0 t x，建立差分格式为： sin x j sin Jh , J 利用中心差商，建立差分格式为: 0,1,2,L 10 (4) 对于边界条件将差分格式延拓使综上（3 ）、 （4 ）、 k 1 u j 其中r山o.1 1 U J 2 1 U j 0,即U1二U j1, J 0,1,2,L 10 (5) 0,t 0,2 ，建立差分格式为： U N 0,k 0,1,L ,200 k 0为内点，代入（3）得到的式子再与（5）联立消去 1 1 2 0 ’ 2 0 1 5 r U, 1 1 r U, r J 2 J J 2 （7 ）得到三层显格式如下： U 0,t U 1,t k U0 (6 ) 、 2 k r U j 1 2 1 r2k 2 k U J r U J 1 k 1? U j , J U j （6） 1后整理得到: U j 1 (7) （局部截断误差为 1,2,L 9,k 1,2,L ,199 h2) 1 U j U J sin 1 2 0 2r U J 1 k U o X j k U N sin 2 0 r U j 0,k 0,1,2,L 10 Jh ,J 1 2r2u01, J 1,2,L 9 0,1L ,200 (8) 四?实验环境（所用软件、硬件等）及实验数据文件Matlab
第8章 常微分方程边值问题的数值解法 引 言 第7章介绍了求解常微分方程初值问题的常用的数值方法；本章将介绍常微分方程的边值问题的数值方法。 只含边界条件(boundary-value condition)作为定解条件的常微分方程求解问题称为常微分方程的边值问题(boundary-value problem). 为简明起见，我们以二阶边值问题为 则边值问题(8.1.1)有唯一解。 推论 若线性边值问题 ()()()()()(),, (),()y x p x y x q x y x f x a x b y a y b αβ'''=++≤≤?? ==? (8.1.2) 满足 (1) (),()p x q x 和()f x 在[,]a b 上连续； (2) 在[,]a b 上, ()0q x >, 则边值问题(8.1.1)有唯一解。 求边值问题的近似解，有三类基本方法： (1) 差分法(difference method)，也就是用差商代替微分方程及边界条件中的导数，最终化为代数方程求解； (2) 有限元法(finite element method)；
(3) 把边值问题转化为初值问题，然后用求初值问题的方法求解。 差分法 8.2.1 一类特殊类型二阶线性常微分方程的边值问题的差分法 设二阶线性常微分方程的边值问题为 (8.2.1)(8.2.2) ()()()(),,(),(), y x q x y x f x a x b y a y b αβ''-=<
实验报告 一、求方程f(x)=x^3-sinx-12x+1的全部根, ε=1e -6 1、 用一般迭代法; 2、 用牛顿迭代法; 并比较两种迭代的收敛速度。 一、首先，由题可求得：12cos 3)(2 ' --=x x x f . 其次，分析得到其根所在的区间。 ① 令()0=x f ，可得到x x x sin 1123 =+-. ② 用一阶导数分析得到1123 +-x x 和x sin 两个函数的增减区间；再用二阶导数分析得到 两个函数的拐点以及凹凸区间. ③ 在直角坐标轴上描摹出01123 =+-x x 和0sin =x 的图，在图上可以看到他们的交点，然后估计交点所在的区间，即是所要求的根的区间。经过估计，得到根所在的区间为 []3,4--，[]1,0和[]4,3. 1、 一般迭代法 （1）算法步骤： 设ε为给定的允许精度，迭代法的计算步骤为： ① 选定初值0x .由()0=x f 确定函数()x g ，得等价形式()x g x =. ② 计算()0x g .由迭代公式得()01x g x =. ③ 如果ε≤-01x x ，则迭代结束，取1x 为解的近似值；否则，用1x 代替0x ，重复步骤②和步骤③. （2）程序代码： ① 在区间[]3,4--内， 代码： clc
x0=-3.5; %初值0x iter_max=100; %迭代的最大次数 ep=1e-6; %允许精度 ε k=0; while k<=iter_max %k 从0开始到iter_max 循环 x1=(sin(x0)+12*x0-1).^(1/3); %代入0x ，算出1x 的值 if abs(x1-x0) 设计课题 设计一个具有恒流偏置的单端输入-单端输出差分放大器。 学校：延安大学 一： 已知条件 正负电源电压V V V V EE cc 12,12-=-+=+；负载Ω=k R L 20；输入差 模信号mV V id 20=。 二：性能指标要求 差模输入电阻Ω>k R id 10；差模电压增益15≥vd A ；共模抑制 比dB K CMR 50>。 三：方案设计及论证 方案一： 方案二 方案论证： 在放大电路中，任何元件参数的变化，都将产生输出电压的漂移，由温度变化所引起的半导体参数的变化是产生零点漂移的主要原因。采用特性相同的管子使它们产生的温漂相互抵消，故构成差分放大电路。差分放大电路的基本性能是放大差模信号，抑制共模信号好，采用恒流源代替稳流电阻，从而尽可能的提高共模抑制比。 论证方案一：用电阻R6来抑制温漂 ?优点：R6 越大抑制温漂的能力越强； ?缺点：<1>在集成电路中难以制作大电阻； <2> R6的增大也会导致Vee的增大（实际中Vee不 可能随意变化） 论证方案二 优点：（1）引入恒流源来代替R6,理想的恒流源内阻趋于无穷，直流压降不会太高，符合实际情况； （2）电路中恒流源部分增加了两个电位器，其中47R的用来调整电路对称性，10K的用来控制Ic的大小，从而调节静态工作点。 通过分析最终选择方案二。 四：实验工作原理及元器件参数确定 ?静态分析：当输入信号为0时， ?I EQ≈(Vee-U BEQ)/2Re ?I BQ= I EQ /(1+β) ?U CEQ=U CQ-U EQ≈Vcc-I CQ Rc+U BEQ 动态分析 ?已知：R1=R4,R2=R3(完整word版)差分放大器设计的实验报告
相关主题
文本预览