利用计算机模拟RC电路试验-青岛大学.PDF

维普资讯 第 14卷 第 2期 青 岛大学学报 V()L 14NO 2 JuIv． 20t)】 2001兰 7 月 文章编号 ：1006-1037 (2001)02-0077-079 利用计算机模拟 RC电路实验 孙光辉 (青岛大学物理 系．青岛 266071) 摘要 ：奉文对 电路实验进行了计算机模拟求解，这对现有的Rc电路实验是一个很好的 补充，同时它也表 明用计算机解决物理学 问题是一个很好的途径。 关键词：RC 电路 计算机模拟 ：数值解 中图分类号：04-33 文献标识码 ：^ 引言 RC 电路实验属普通物理 电磁学实验 ．它的主要内容是研究 RC 电路在放 电与充 电过 程 中电容器上的电量与时间的关系 在这个实验 p描述 RC 电路的微 分方程存在精确解， 我们采用计算机模拟方法给出其数值解 ．一方面对 有的RC电路实验是一个很好的补充， 另一方面揭示了物理学与计算机的密切联系，鼠某种程度上甚至可 以起到代替该实验课的 作用 。 1 Rc电路中电容器的放 电 如果把一个 已充电的电容器与一个 电阻串联连接，则 一旦开关闭合，电容器就马上开 始放 电。因为电阻两端的电压在任意时刻部与电容器 琦『极 的电压值相等，所 以可列出下 式 VR+_vr：1【 1‘) 又 VR=TR=Rdq／di V ~q．／C (2) 将 (2)式代入 (1)式，得到： Rdq tq／C-~ (3) 此一阶微分方程整理后变为： dq／dr=一q，『 4‘) 式中 T=RC，其精确解为 q=q}一 “．其中 q．j表示 1=0时电容器 上的电量，T是当减少到 其初始值的 1／e~-．o．368倍时所需要的时 2 RC电路中电容器的充电 在 RC电路中，电动势为E的电池正在对 电容器充电．如果我们把基尔霍夫回路方程 用于此 电路，则 IR+q／C=E (5) 利用 I=dq／dt上式整理后得到：dq／dt=E~ -qR／C (6) 该式的精确解为 q=q(1--C—vRC)，其中q是经过 穷长叫间后 电容器中的电量，它的值 为q-=CEo，假定初始条件为当时间 t=0时电容器一 由电量为零。 3 计算机模拟求解 我们采用 FORTRAN语言，利用四阶 l~unge—Kutta算法编写 了 一个计算程序对上述 问题 ’收稿 日期；2CO14364)4 维普资讯 青岛大 学报 第 14卷 进行了计算。其求解结果存在 rk．o11t文件 中．初始条件存在 rk．inp文件里，该程序可计 算 Rc电路中电容器的充 电过程和放 电过程 中其 电量的数值解 。该程序清单如下： subro|ltiuerk(tO，tm，x0，h，f) real k1，k2，k3，k4 wrire(2．10) 10 format(5x，lht，14x，1hx) Write(220) tO．x0 20 format(1x．f5 I{fl5 5)