《计算机数学基础》数值部分第五单元辅导 14 常微分方程的数值解法 一.重点内容 1. 欧拉公式： ）心知1）a 儿+1 =儿 + hfg ，儿） m 1、 伙=0丄2,…川一 1） I 无=x Q +kh 局部截断误差是0（*）。 2. 改进欧拉公式： 预报一校正公式： 预报值 _v*+1 =儿+ hf （x k ,儿） - h - 校正值 y M = y k +-［f （x kt y k ） + /（x A+1, y M ）］ 即 儿+1 =儿+ ￡ "（忑'儿）+心+「儿+ hfg ,儿））］ 或表成平均的形式： 儿=儿+ hfg ,儿） '儿=儿+"（无+】，儿） +K ） 改进欧拉法的局部截断误差是0（2） 3. 龙格一库塔法 二阶龙格一库塔法的局部截断误差是0（爪） 三阶龙格一库塔法的局部截断误差是0（护） 四阶龙格F 塔法公式：儿计=儿+ 2（匕+ 2心+ 2? + ?） 四阶龙格一库塔法的局部截断误差是0（爪）。 二实例 y' = — y — xv f2（0 < x < 0.6） 例1用欧拉法解初值问题｛ ' ? -取步长/匸02计算过程保留 b （o ）= 1 4位小数。 解/i=0.2. f （x ）= —y —xy 2<,首先建立欧拉迭代格式 y*+i =儿+ hf g,y k ） = y k -hy k -hx k y ； =0.2 儿（4 - x k y k ）（k = 0,1,2） K 2=f(x n +^h, yk+-hK\)t gg+舟人，>'n +y/?A3）；
当k=0, xi=0.2 时，已知x（）=0,y（）=l,有 y(0?2)今i=0?2X l(4-0X 1)=0.8000 当k=\. M=0?4时，已知“=0?2」尸0?8,有 y(0?4)今2=0.2 X 0.8X(4-0.2X0.8)=0.614 4 当k=2, xs=0.6 时，已知x2=0.4,y2=0.6144,有 y(0?6)今3=0.2 X0.6144X (4-0.4 X 0.4613)=0.8000 「J, ,2 ?_ ZX 例2用欧拉预报一校正公式求解初值问题\y + v +V sinx=,取步长/?=0.2,计算 .y ⑴=1 y(0.2),y(0.4)的近似值，计算过程保留5位小数。 解步长力=0.2,此时/(x,y)=—y—fsiiu 欧拉预报一校正公式为： 预报值兀I = y k + hfg y k) - I J_ 校正值)3=儿+尹(忑，儿)+ fg,儿+1)] 有迭代格式] 预报值儿+] = y k 4-h(-y k -y； sin x k) =y k (0?8-0?2儿sin x k) < h 、—— 2 校止值y如]=儿 +尸[(一片一力sinxJ + LN+i-yl sin.v I+1)] ——?> =儿(°?9一0?1儿sin心)一0?1(儿+| +y；j sin心利) 、"M=0.別=1」)=1 时，Xj=1.2> 有 儿=yo(°?8-O?2yo sinx0) = 1 x (0.8-02x lsin 1) = 0.63171 y(1.2) ?= lx(0.9-0.1xlxsinl)-0.1(0.63171+0.631712sinl.2) = 0.71549 当 T xi=1.2, yi=0.71549 时,x2=1.4,有 y2 =儿(0.8-0?2儿sinXj) = 0.71549x(0.8-02x0.71549sinl.2) =0.47697 y(14) z y2 = 0.71549x(0.9-0.1x0.71549xsin 1.2)-0.1(0.47697+ 0.476972 sin 1.4) =0.52608 V = 8 — 3y 例3写出用四阶龙格一库塔法求解初值问题^ ‘的计算公式，取步长/匸0.2计 b(0) = 2 算y(0.4)的近似值。讣算过程保留4位小数。 解此处.心,刃=8 —3”四阶龙格一库塔法公式为 艰=儿 + % + 2? + 2勺 + ?) 1 h, y n+ y/?A3)： 本例计算公式为: 0 2 呱严儿+三(32?+2?+心
三极管电路放大电路及其分析方法 一、教学要求 1. 重点掌握的内容 （1）放大、静态与动态、直流通路与交流通路、静态工作点、负载线、放大倍数、输入电阻与输出电阻的概念； （2）用近似计算法估算共射放大电路的静态工作点； （3）用微变等效电路法分析计算共射电路、分压式工作点稳定电路的电压放大倍数A u和A us,输入电阻R i和输出电阻R0。 2. 一般掌握的内容 （1）放大电路的频率响应的一般概念； （2）图解法确定共射放大电路的静态工作点，定性分析波形失真，观察电路参数对静态工作点的影响，估算最大不失真输出的动态范围； （3）三种不同组态（共射、共集、共基）放大电路的特点； （4）多级放大电路三种耦合方式的特点，放大倍数的计算规律。 3. 一般了解的内容 （1）共射放大电路f L、f H与电路参数间的定性关系，波特图的一般知识＜多级放大电路与共射放大电路频宽的定性分析； （2）用估算法估算场效应管放大电路静态工作点的方法。 二?内容提要 1. 共射接法的两个基本电路 共射放大电路和分压式工作点稳定电路是模拟电路中最基本的单元电路。学习这两种基本电路的分析方法是学习比较复杂的模拟电路的基础。 2. 两种基本分析方法——图解法和微变等效电路法 在“模拟电路”中，三极管是非线性元件，因此不能简单地采用“电路与磁路”课中线性电路地分析方法。图解法和微变等效电路法就是针对三极管非线性的特点而采用的分析方法。 3. 放大电路的三种组态——共射组态、共集组态和共基组态 由于放大电路输入、输出端取自三极管三个不同的电极，放大电路有三种组态——共射组态、共集组态和共基组态。由于组态的不同，其放大电路反映出的特性是不同的。在实际中，可根据要求选择相应组态的电路。 4. 两种放大元件组成的放大电路——双极型三极管放大电路和场效应管放大电路 一般来说，双极性三极管是一种电流控制元件，它通过基极电流i B的变化控制集电极电流I c的变化。而场效应管是一种电压控制元件，它通过改变栅源间的电压U GS来控制漏极电流i D的变化；其次，双极性三极管的输入电阻较小，而场效应管的输入电阻很高，静态时栅极几乎不取电流。由于它们性能和特点的不同，可根据要求选用不同元件组成的放大电路。 5. 多级放大电路的三种耪合方式一一阻容耦合、直接耦合和变压器耦合 将多级放大电辟连接起来的时候，就出现了级与级之间的耦合方式问题。通过电阻和电容将两级放大电路连接起来的方式称为阻容耦合。由于电容的作用，
上海交通大学致远学院 《常微分方程和偏微分方程的数值解法》教学大纲 一、课程基本信息 课程名称（中文）：常微分方程和偏微分方程的数值解法 课程名称（英文）：Numerical Methods for Ordinary and Partial Differential Equations 课程代码：MA300 学分 / 学时：4学分 / 68学时 适用专业：致远学院与数学系相关专业 先修课程：偏微分方程，数值分析 后续课程：相关课程 开课单位：理学院数学系计算与运筹教研室 Office hours: 每周二19：00—21：00，地点：数学楼1204 二、课程性质和任务 本课程是致远学院和数学系应用数学和计算数学方向的一门重要专业基础课程，其主要任务是通过数学建模、算法设计、理论分析和上机实算“四位一体”的教学方法，使学生掌握常微分方程与偏微分方程数值解的基本方法、基本原理和基本理论，进一步提升同学们利用计算机解决实际问题的能力。在常微分方程部分，将着重介绍常微分方程初值问题的单步法，含各类Euler方法和Runge-Kutta方法，以及线性多步法。将简介常微分方程组和高阶常微分方程的数值方法。在偏微分方程部分，将系统介绍求解椭圆、双曲、抛物型方程的差分方法的构造方法和理论分析技巧，对于椭圆型方程的边值问题将介绍相应变分原理与有限元方法。将在课堂上实时演示讲授的核心算法的计算效果，以强调其直观效果与应用性。本课程重视实践环节建设，学生要做一定数量的大作业。 三、教学内容和基本要求 第一部分：常微分方程数值解法 1 引论 1.1回顾：一阶常微分方程初值问题及解的存在唯一性定理
二、基本放大电路及其分析方法 一个放大器一般是由多个单级放大电路所组成，着重讨论双极型半导体三极管放大电路的三种组态，即共发射极，共集电极和共基极三种基本放大电路。从共发射极电路入手，推及其他二种电路，其中将图解分析法和微变等效电路分析法，作为分析基础来介绍。分析的步骤，首先是电路的静态工作点，然后分析其动态技术指标。对于放大器来说，主要的动态技术指标有电压放大倍数、输入阻抗和输出阻抗。 2.1.共射极基本放大电路的组成及放大作用 在实践中，放大器的用途是非常广泛的，它能够利用三极管的电流控制作用把微弱的电信号增强到所要求的数值，为了了解放大器的工作原理，先从最基本的放大电路学习： 图2.1称为共射极放大电路，要保证发射结正偏，集电极反偏Ib=（V BB-V BE）/Rb，对于硅管V BE约为0.7V左右,锗管约为0.2V左右,I B=(V BB-0.7)/Rb这个电路的偏流I B决定于V BB 和Rb的大小,V BB和Rb一经确定后,偏流I B就固定了,所以这种电路称为固定偏流电路,Rb又称为基极偏置电阻,电容Cb1和Cb2为隔直电容或耦合电容,在电路中的作用是“传送交流,隔离直流”,放大作用的实质是利用三极管的基极对集电极的控制作用来实现的. 上图是共射极放大电路的简化图,它在实际中用得比较多的一种电路组态,放大电路的主要性能指标，常用的有放大倍数、输入阻抗、输出阻抗、非线性失真、频率失真以及输出功率和效率等。对于不同的用途的电路，其指标各有侧重。 初步了解放大电路的组成及简单工作原理后，就可以对放大电路进行分析。主要方法有图解法和微变等效法。 2.2.图解分析法 2.2.1.静态工作情况分析 当放大电路没有输入信号时，电路中各处的电压，电流都是不变的直流，称为直流工作状态简称静态，在静态工作情况下，三极管各电极的直流电压和直流电流的数值，将在管子的特性曲线上确定一点，这点称为静态工作点，下面通过例题来说明怎样估算静态工作点。 解:Cb1与Cb2的隔直作用,对于静态下的直流通路,相当于开路,计算静态工作点时,只需考虑图中的Vcc、Rb、Rc及三极管所组成的直流通路就可以了，I B=（Vcc－0.7）/Rb （I C=βI B+I CEO ） I C=βI B，V CE=V CC－I C R C 如已知β，利用上式可近似估算放大电路的静态工作点。 2.2.2.用图解法确定静态工作点 在分析静态工作情况时,只需研究由V CC、R C、V BB、Rb及半导体三极管所组成的直
i.常微分方程初值问题数值解法 常微分方程初值问题的真解可以看成是从给定初始点出发的一条连续曲线。差分法是常微分方程初值问题的主要数值解法，其目的是得到若干个离散点来逼近这条解曲线。有两个基本途径。一个是用离散点上的差商近似替代微商。另一个是先对微分方程积分得到积分方程，再利用离散点作数值积分。 i.1 常微分方程差分法 考虑常微分方程初值问题：求函数()u t 满足 (,), 0du f t u t T dt =<≤ （i.1a ） 0(0)u u = (i.1b) 其中(,)f t u 是定义在区域G : 0t T ≤≤, u <∞上的连续函数，0u 和T 是给定的常数。我们假设(,)f t u 对u 满足Lipschitz 条件，即存在常数L 使得 121212(,)(,), [0,]; ,(,)f t u f t u L u u t T u u -≤-?∈∈-∞∞ （i.2） 这一条件保证了（i.1）的解是适定的，即存在，唯一，而且连续依赖于初值0u 。 通常情况下，（i.1）的精确解不可能用简单的解析表达式给出，只能求近似解。本章讨论常微分方程最常用的近似数值解法-差分方法。先来讨论最简单的Euler 法。为此，首先将求解区域[0,]T 离散化为若干个离散点： 0110N N t t t t T -=<< <<= （i.3） 其中n t hn =，0h >称为步长。 在微积分课程中我们熟知，微商（即导数）是差商的极限。反过来，差商就是微商的近似。在0t t =处，在（i.1a ）中用向前差商 10()()u t u t h -代替微商du dt ，便得 10000()()(,())u t u t hf t u t ε=++ 如果忽略误差项0ε，再换个记号，用i u 代替()i u t 便得到 1000(,)u u hf t u -= 一般地，我们有 1Euler (,), 0,1, ,1n n n n u u hf t u n N +=+=-方法： （i.4） 从(i.1b) 给出的初始值0u 出发，由上式可以依次算出1,,N t t 上的差分解1,,N u u 。
3.2 基本放大电路的分析方法 3.2.1 放大电路的静态分析 放大电路的静态分析有计算法和图解分析法两种。 （1）静态工作状态的计算分析法 根据直流通路可对放大电路的静态进行计算 (03.08) I C= I B (03.09) V CE=V CC－I C R c (03.10) I B、I C和V CE这些量代表的工作状态称为静态工作点，用Q表示。 在测试基本放大电路时，往往测量三个电极对地的电位V B、V E和V C即可确定三极管的工作状态。 （2）静态工作状态的图解分析法 放大电路静态工作状态的图解分析如图03.08所示。 图03.08 放大电路静态工作状态的图解分析 直流负载线的确定方法：
1. 由直流负载列出方程式V CE=V CC－I C R c 2. 在输出特性曲线X轴及Y轴上确定两个特殊点 V CC和V CC/R c,即可画出直流负载线。 3. 在输入回路列方程式V BE =V CC－I B R b 4. 在输入特性曲线上，作出输入负载线，两线的交点即是Q。 5. 得到Q点的参数I BQ、I CQ和V CEQ。 例3.1：测量三极管三个电极对地电位如图03.09所示，试判断三极管的工作状态。 图03.09 三极管工作状态判断 例3.2：用数字电压表测得V B=4.5V 、V E=3.8V 、V C=8V，试判断三极管的工作状态。 电路如图03.10所示 图03.10 例3.2电路图 3.2.2 放大电路的动态图解分析 (1) 交流负载线 交流负载线确定方法：
1.通过输出特性曲线上的Q点做一条直线，其斜率为1/R L'。 2.R L'= R L∥R c，是交流负载电阻。 3.交流负载线是有交流输入信号时，工作点Q的运动轨迹。 4.交流负载线与直流负载线相交，通过Q点。 图03.11 放大电路的动态工作状态的图解分析 (2) 交流工作状态的图解分析 动画 图03.12 放大电路的动态图解分析（动画3-1）通过图03.12所示动态图解分析，可得出如下结论： 1. v i→↑ v BE→↑ i B→↑ i C→↑ v CE→↓ |-v o|↑； 2. v o与v i相位相反； 3.可以测量出放大电路的电压放大倍数； 4.可以确定最大不失真输出幅度。 (3) 最大不失真输出幅度 ①波形的失真
第一章 绪 论 1．设0x >,x 的相对误差为δ，求ln x 的误差。 解：近似值*x 的相对误差为* **** r e x x e x x δ-=== 而ln x 的误差为()1ln *ln *ln **e x x x e x =-≈ 进而有(ln *)x εδ≈ 2．设x 的相对误差为2%，求n x 的相对误差。 解：设()n f x x =，则函数的条件数为'()||() p xf x C f x = 又1'()n f x nx -= , 1 ||n p x nx C n n -?∴== 又((*))(*)r p r x n C x εε≈? 且(*)r e x 为2 ((*))0.02n r x n ε∴≈ 3．下列各数都是经过四舍五入得到的近似数，即误差限不超过最后一位的半个单位，试指 出它们是几位有效数字：*1 1.1021x =,*20.031x =, *3385.6x =, *456.430x =,*5 7 1.0.x =? 解：*1 1.1021x =是五位有效数字； *20.031x =是二位有效数字； *3385.6x =是四位有效数字； *456.430x =是五位有效数字； *57 1.0.x =?是二位有效数字。 4．利用公式(2.3)求下列各近似值的误差限：(1) *** 124x x x ++,(2) ***123x x x ,(3) **24 /x x . 其中****1234,,,x x x x 均为第3题所给的数。 解：
*4 1*3 2*13*3 4*1 51 ()102 1()102 1()102 1()102 1()102x x x x x εεεεε-----=?=?=?=?=? ***124***1244333 (1)() ()()() 111101010222 1.0510x x x x x x εεεε----++=++=?+?+?=? ***123*********123231132143 (2)() ()()() 1111.10210.031100.031385.610 1.1021385.610222 0.215 x x x x x x x x x x x x εεεε---=++=???+???+???≈ **24****24422 *4 33 5 (3)(/)()() 110.0311056.430102256.43056.430 10x x x x x x x εεε---+≈??+??=?= 5计算球体积要使相对误差限为1，问度量半径R 时允许的相对误差限是多少？ 解：球体体积为343V R π= 则何种函数的条件数为 2 3 '4343 p R V R R C V R ππ=== (*)(*)3(*)r p r r V C R R εεε∴≈= 又(*)1r V ε=
实验4 常微分方程的数值解 【实验目的】 1．掌握用MATLAB软件求微分方程初值问题数值解的方法； 2．通过实例用微分方程模型解决简化的实际问题； 3．了解欧拉方法和龙格-库塔方法的基本思想和计算公式，及稳定性等概念。 【实验内容】 题3 小型火箭初始重量为1400kg，其中包括1080kg燃料。火箭竖直向上发射时燃料燃烧率为18kg/s，由此产生32000N的推力，火箭引擎在燃料用尽时关闭。设火箭上升时空气阻力正比于速度的平方，比例系数为0.4kg/m，求引擎关闭瞬间火箭的高度、速度、加速度，及火箭到达最高点的时的高度和加速度，并画出高度、速度、加速度随时间变化的图形。 模型及其求解 火箭在上升的过程可分为两个阶段，在全过程中假设重力加速度始终保持不变，g=9.8m/s2。 在第一个过程中，火箭通过燃烧燃料产生向上的推力，同时它还受到自身重力（包括自重和该时刻剩余燃料的重量）以及与速度平方成正比的空气阻力的作用，根据牛顿第二定律，三个力的合力产生加速度，方向竖直向上。因此有如下二式： a=dv/dt=(F-mg-0.4v2)/m=(32000-0.4v2)/(1400-18t)-9.8 dh/dt=v 又知初始时刻t=0，v=0，h=0。记x(1)=h，x(2)=v，根据MATLAB 可以求出0到60秒内火箭的速度、高度、加速度随时间的变化情况。程序如下： function [ dx ] = rocket( t,x ) a=[(32000-0.4*x(2)^2)/(1400-18*t)]-9.8; dx=[x(2);a]; end ts=0:1:60;
x0=[0,0]; [t,x]=ode45(@rocket,ts,x0); h=x(:,1); v=x(:,2); a=[(32000-0.4*(v.^2))./(1400-18*t)]-9.8; [t,h,v,a]; 数据如下： t h v a 0 0 0 13.06 1.00 6.57 13.19 13.30 2.00 26.44 26.58 1 3.45 3.00 59.76 40.06 13.50 4.00 106.57 53.54 13.43 5.00 16 6.79 66.89 13.26 6.00 240.27 80.02 12.99 7.00 326.72 92.83 12.61 8.00 425.79 105.22 12.15 9.00 536.99 117.11 11.62 10.00 659.80 128.43 11.02 11.00 793.63 139.14 10.38 12.00 937.85 149.18 9.71 13.00 1091.79 158.55 9.02 14.00 1254.71 167.23 8.33 15.00 1425.93 175.22 7.65 16.00 1604.83 182.55 6.99 17.00 1790.78 189.22 6.36 18.00 1983.13 195.27 5.76 19.00 2181.24 200.75 5.21 20.00 2384.47 205.70 4.69 21.00 2592.36 210.18 4.22 22.00 2804.52 214.19 3.79 23.00 3020.56 217.79 3.41 24.00 3240.08 221.01 3.07 25.00 3462.65 223.92 2.77 26.00 3687.88 226.56 2.50 27.00 3915.58 228.97 2.27
第一章 绪论（12） 1、设0>x ，x 的相对误差为δ，求x ln 的误差。 [解]设0*>x 为x 的近似值，则有相对误差为δε=)(*x r ，绝对误差为**)(x x δε=，从而x ln 的误差为δδεε=='=* ****1)()(ln )(ln x x x x x ， 相对误差为* * ** ln ln ) (ln )(ln x x x x r δ εε= = 。 2、设x 的相对误差为2%，求n x 的相对误差。 [解]设*x 为x 的近似值，则有相对误差为%2)(*=x r ε，绝对误差为**%2)(x x =ε，从而n x 的误差为n n x x n x n x x n x x x ** 1 *** %2%2) ()()()(ln * ?=='=-=εε， 相对误差为%2) () (ln )(ln *** n x x x n r == εε。 3、下列各数都是经过四舍五入得到的近似数，即误差不超过最后一位的半个单位，试指出它们是几位有效数字： 1021.1*1=x ，031.0*2=x ，6.385*3=x ，430.56*4=x ，0.17*5 ?=x 。 [解]1021.1*1 =x 有5位有效数字；0031.0* 2=x 有2位有效数字；6.385*3=x 有4位有效数字；430.56* 4 =x 有5位有效数字；0.17*5?=x 有2位有效数字。 4、利用公式（3.3）求下列各近似值的误差限，其中* 4*3*2*1,,,x x x x 均为第3题所给 的数。 （1）* 4*2*1x x x ++； [解]3 334* 4*2*11** *4*2*1*1005.1102 1 10211021)()()()()(----=?=?+?+?=++=? ??? ????=++∑x x x x x f x x x e n k k k εεεε； （2）* 3*2 *1x x x ；
第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 αβ''-=<
第一章绪论 1.1 引言 常微分方程是现代数学的一个重要分支，是人们解决各种实际问题的有效工具。微分方程的理论和方法从17世纪末开始发展起来，很快成了研究自然现象的强有力工具，在17到18世纪，在力学、天文、科学技术、物理中，就已借助微分方程取得了巨大的成就。1864年Leverrer根据这个方程预见了海王星的存在，并确定出海王星在天空中的位置。现在，常微分方程在许多方面获得了日新月异的应用。这些应用也为常微分方程的进一步发展提供了新的问题，促使人们对微分方程进行更深入的研究，以便适应科学技术飞速发展的需要。 研究常微分方程常用数值解是数学工作者的一项基本的且重要的工作。在国内外众多数学家的不懈努力，使此学科基本上形成了一套完美的体系。微分方程的首要问题是如何求一个给定方程的通解或特解。到目前为止，人们已经对许多微分方程得出了求解的一般方法。由于在生产实际和科学研究中所遇到的微分方程问题比较复杂，使这些问题的解即使能求出解析表达式，也往往因计算量太大而难于求出，而对于一些典型的微分方程则可以运用基本方法求出其解析解，并可以根据初值问题的条件把其中的任意常数确定下来。 由于求通解存在许多困难，人们就开始研究带某种定解条件的特解。首先是Cauchy对微分方程初始解的存在惟一性进行了研究。目前解的存在惟一性、延拓性、大范围的存在性以及解对初始解和参数的延续性和可微性等理论问题都已发展成熟。与此同时，人们开始采取各种近似方法来求微分方程的特解，例如求微分方程数值解的Euler折线法、Runge-Kutta法等，可以求得若干个点上微分方程的近似解。最后，由于当代高科技的发展为数学的广泛应用和深入研究提供了更好的手段。用计算机结合Matlab软件求方程的精确解、近似解，对解的性态进行图示和定性、稳定性研究都十分方便有效。 本章先介绍常微分的一般概念、导出微分方程的一些典型例子及求解微分方程的思路分析。从而得到常微分方程的常用数值解法。
第一章 绪论 p19 2．设x 的相对误差为2%，求n x 的相对误差。 解：设()n f x x =，则函数的条件数为'() | |() p xf x C f x = 又 1 '()n f x nx -=, 1 ||n p x nx C n n -?∴== 又 ((*))(*)r p r x n C x εε≈? 且(*)r e x 为2% ((*))0.02n r x n ε∴≈ 5计算球体积要使相对误差限为1，问度量半径R 时允许的相对误差限是多少？ 解：球体体积为343 V R π= 则何种函数的条件数为 2 3'4343 p R V R R C V R ππ=== (*)(*)3(*)r p r r V C R R εεε∴≈= 又 (*)1r V ε= 故度量半径R 时允许的相对误差限为1 (*)10.333 r R ε= ?≈ 7．求方程2 5610x x -+=的两个根，使它至少具有427.982 =）。 解：2 5610x x -+= ， 故方程的根应为1,228x =故 128 2827.98255.982x = ≈+= 1x ∴具有5位有效数字 211 280.0178632827.98255.982 x =-= ≈ =≈+ 2x 具有5位有效数字
9．正方形的边长大约为了100cm ，应怎样测量才能使其面积误差不超过2 1cm ？ 解：正方形的面积函数为2 ()A x x = p7 当*100x =时，若(*)1A ε≤, 则21 (*)102 x ε-≤ ? 故测量中边长误差限不超过0.005cm 时，才能使其面积误差不超过2 1cm 第二章 插值法p48 1．当1,1,2 x =-时，()0,3,4f x =-, 分别用单项式基底、拉格朗日基底、牛顿基底求() f x 的二次插值多项式。 解： 0120121200102021101201220211,1,2, ()0,()3,()4;()()1 ()(1)(2)()()2()()1 ()(1)(2) ()()6 ()()1 ()(1)(1) ()()3 x x x f x f x f x x x x x l x x x x x x x x x x x l x x x x x x x x x x x l x x x x x x x ==-===-=--==-+-----==------= =-+-- 则二次拉格朗日插值多项式为 2 20 ()()k k k L x y l x ==∑ 0223()4() 14 (1)(2)(1)(1)23537623 l x l x x x x x x x =-+=---+-+=+- 2．给出()ln f x x =的数值表 用线性插值及二次插值计算的近似值。 解：由表格知，
放大电路及其分析方法 2.1 放大电路的基本概念 三极管具有电流放大作用，如何使用三极管构成一个电路，实现对输入信号的放大？本节就来讨论这一问题。 基本放大电路是放大电路中最基本的结构，是构成复杂放大电路的基本单元。它利用双极型半导体三极管输入电流控制输出电流的特性，或场效应半导体三极管输入电压控制输出电流的特性，实现信号的放大。本章基本放大电路的知识是进一步学习电子技术的重要基础。本书中双极型半导体三极管简称三极管，场效应半导体三极管简称场效应管。 2.1.1 放大的概念 基本放大电路一般是指由一个三极管或场效应管组成的放大电路。从电路的角度来看，可以将基本放大电路看成一个双端口网络。放大的作用体现在如下方面： 1．放大电路主要利用三极管或场效应管的控制作用放大微弱信号，输出信号在电压或电流的幅度上得到了放大，输出信号的能量得到了加强。 2．输出信号的能量实际上是由直流电源提供的，只是经过三极管的控制，使之转换成信号能量，提供给负载。放大电路的结构示意图见图2-1-1。 图2-1-1 放大电路结构示意图 2.1.2 基本放大电路的组成及工作原理 一、共射组态基本放大电路的组成 共射组态基本放大电路如图2-1-2所示。在该电路中，输入信号加在加在基极和发射极之间，耦合电容器C1和Ce视为对交流信号短路。输出信号从集电极对地取出，经耦合电容器C2隔除直流量，仅将交流信号加到负载电阻RL之上。放大电路的共射组态实际上是指放大电路中的三极管是共射组态。 图2-1-2 共射组态交流基本放大电路 二、放大原理 在输入信号为零时，直流电源通过各偏置电阻为三极管提供直流的基极电流和直流集电极电流，并在三极管的三个极间形成一定的直流电压。由于耦合电容的隔直流作用，直流电压无法到达放大电路的输入端和输出端。
第四章 数值积分与数值微分 1.确定下列求积公式中的特定参数，使其代数精度尽量高，并指明所构造出的求积公式所具有的代数精度： 101210121 12120 (1)()()(0)(); (2)()()(0)(); (3)()[(1)2()3()]/3; (4)()[(0)()]/2[(0)()]; h h h h h f x dx A f h A f A f h f x dx A f h A f A f h f x dx f f x f x f x dx h f f h ah f f h -----≈-++≈-++≈-++''≈++-?? ?? 解： 求解求积公式的代数精度时，应根据代数精度的定义，即求积公式对于次数不超过m 的多项式均能准确地成立，但对于m+1次多项式就不准确成立，进行验证性求解。 （1）若101(1) ()()(0)()h h f x dx A f h A f A f h --≈-++? 令()1f x =，则 1012h A A A -=++ 令()f x x =，则 110A h A h -=-+ 令2 ()f x x =，则 3 221123 h h A h A -=+ 从而解得 01 1431313A h A h A h -?=?? ?=?? ?=?? 令3 ()f x x =，则 3()0h h h h f x dx x dx --==? ? 101()(0)()0A f h A f A f h --++= 故 101()()(0)()h h f x dx A f h A f A f h --=-++? 成立。 令4 ()f x x =，则
毕业论文 题目抛物型方程的差分解法学院数学科学学院 专业信息与计算科学 班级计算0802 学生王丹丹 学号20080901045 指导教师王宣欣 二〇一二年五月二十五日
摘要 偏微分方程的数值解法在数值分析中占有重要的地位，很多科学技术问题的数值计算包括了偏微分方程的数值解问题【1】。近三十多年来，数值解法的理论和方法都有了很大的发展，而且在各个科学技术的领域中应用也愈来愈广泛。本文的研究主要集中在依赖于时间的问题，借助于简单的常系数扩散方程，介绍抛物型方程的差分解法。本文以基本概念和基本方法为主，同时结合算例实现算法。 第一部分介绍偏微分方程及差分解法的基本概念，引入本文的研究对象——常系 数扩散方程： 2 2 ,,0 u u a x R t t x ?? =∈>?? 第二部分介绍上述方程的几种差分格式及每种格式的相容性、收敛性与稳定性。 第三部分通过算例检验每种差分格式的可行性。 关键词：偏微分方程；抛物型；差分格式；收敛性；稳定性；算例
ABSTRACT The numerical solution of partial differential equation holds an important role in numerical analysis .Many problems of compution in the field of science and techology include the numerical solution of partial differential equation. For more than 30 years, the theory and method of the numerical computation made a great development and its applications in various fields of science and technology are more and more widely. This paper focuses on the problems based on time. I will use object-constant diffusion equation to introduces the finite difference method of parabolic equation. This paper mainly focus on the basic concept ,basic method and simple numerical example. The first part of this paper introduces partial differential equations and basic concepts of finite difference method.I will introduce the object-constant diffusion equation for the first time. 2 2 ,,0 u u a x R t t x ?? =∈>?? The second part of this paper introduces several difference schemes of the above equation and their compatibility ,convergence and stability. The third part tests the accuracy of each scheme. Key words:partial differential equation；parabolic；difference scheme；convergence；stability；application
双曲型方程的有限差分法 线性双曲型方程定解问题： （a ）一阶线性双曲型方程 ()0=??+??x u x a t u （b ）一阶常系数线性双曲型方程组 0=??+??x t u A u 其中A ，s 阶常数方程方阵，u 为未知向量函数。 （c ）二阶线性双曲型方程（波动方程） ()022=?? ? ??????-??x u x a x t u ()x a 为非负函数 （d ）二维，三维空间变量的波动方程 0222222=???? ????+??-??y u x u t u 022222222=???? ????+??+??-??z u y u x u t u §1 波动方程的差分逼近 1.1 波动方程及其特征 线性双曲型偏微方程的最简单模型是一维波动方程： （1.1） 22 222x u a t u ??=?? 其中0>a 是常数。 （1.1）可表示为：022 222=??-??x u a t u ，进一步有
0=??? ????+?????? ????-?? u x a t x a t 由于 x a t ?? ±??当a dt dx ±=时为()t x u ,的全导数 （=dt du dt dx x u t u ???+??x u a t u ??±??=）,故由此定出两个方向 （1.3） a dx dt 1±= 解常微分方程（1.3）得到两族直线 （1．4） 1C t a x =?+ 和 2C t a x =?- 称其为特征。 特征在研究波动方程的各种定解问题时，起着非常重要的作用。 比如，我们可通过特征给出（1.1）的通解。（行波法、特征线法） 将（1.4）视为),(t x 与),(21C C 之间的变量替换。由复合函数的微分法则 2 12211C u C u x C C u x C C u x u ??+ ??=?????+?????=?? x C C u C u C x C C u C u C x u ????? ? ????+????+?????? ????+????=??2 212121122 2 22122212212C u C C u C C u C u ??+???+???+??= 22 22122122C u C C u C u ??+???+??= 同理可得 a t t a t C -=??-=??1，a t C =??2 ???? ????-??=?????+?????=??21 2211C u C u a t C C u t C C u t u
放大电路的基本分析方法（20分钟） 一、参考教材 第二章2.1.4 放大电路的基本分析方法 《模拟电子技术简明教程》张国平、曾高荣主编，电子工业出版社出版 二、教学内容 1.放大电路的直流通路和交流通路 2.估算法确定静态工作点 3.图解法确定静态工作点 三、教学目的 1.掌握放大电路的直流通路与交流通路的画法； 2.掌握估算法确定静态工作点 3.掌握图解法确定静态工作点 四、教学重点、难点 1.放大电路的直流通路与交流通路的画法 2.估算法和图解法确定静态工作点 3.分析静态工作点的意义 五、教学方法 采用课堂讲授加PPT展示的方法，通过例题讲解加深学生对教学内容的理解。 六、教学过程设计 1.旧课复习（3分钟），回顾上一节的知识点，如组成放大电路的基本原则、特 点、主要性能指标等。 2.新课内容（17分钟） 1）首先引入静态和动态两个概念，使学生理解放大电路的分析实际上为直流 通路和交流通路分析的叠加；并且在分析中要采用先静态后动态的分析顺序； 引出静态工作点的概念。 2）放大电路的直流通路和交流通路：详细介绍直流通路和交流通路的画法， 并通过实例分析来加深印象。可以让学生自己进行随堂练习以确保对这一知识 点的领会和掌握。在进行实例分析时，简单介绍放大电路的基本分类（共射、共基、共集）。 3) 通过电路实例分析，介绍如何通过估算法获得静态工作点。 4）图解法是放大电路常用的分析方法之一，简单介绍图解法与微变等效电路 分析法的区别，及适用范围。通过分步解析的方式，详细介绍图解法确定静态 工作点。 七、作业 复习题二（5）；三（3）；习题2.3, 2.4 八、教学后记