MATLAB求解RC电路：从微分方程到批量仿真

RC 电路是电路分析里最基础的一阶动态电路很多人在课程作业或者实验报告里都要用 MATLAB 画出它充放电的曲线。网上的代码其实不少但大多数是“把解析公式抄进去直接画图”换一组 R、C 参数之后就很难复用。这篇文章换个思路从 RC 电路本身的微分方程出发按照“解析解 → 数值解 → 批量参数扫描 → 函数封装”的顺序把 RC 电路在 MATLAB 里完整求解一遍。标题里写“求解 3”可以理解成三条路径第一是用解析公式直接画充放电曲线适合推导验证第二是用 ode45 做微分方程数值求解适合处理初值变化和复杂激励第三是用传递函数做阶跃响应分析适合和《信号与系统》《自动控制原理》衔接。最后再把这三条路径整合成一个可复用的 MATLAB 函数方便做批量仿真。整个过程只依赖 MATLAB 的基础函数不需要额外安装工具箱普通办公电脑就能跑也不涉及 GPU 和显存。如果你正在学《电路》《信号与系统》或者要做 MATLAB 课程设计这篇内容可以直接照着敲。1. 核心能力速览能力项说明项目类型一阶 RC 电路建模与 MATLAB 数值求解核心方法解析解、ode45 数值解、传递函数 step 响应、批量参数扫描输入参数电阻 R、电容 C、电源电压 Us、初始电容电压 U0、仿真时长输出内容电容电压/电流随时间变化曲线、时间常数分析、参数对比图运行环境MATLAB 基础版即可无需额外工具箱批量能力支持for 循环批量修改 R、C结果可以自动保存接口能力可封装为函数也可通过 MATLAB Engine 被 Python 调用硬件门槛普通 CPU 电脑即可无 GPU 要求扩展方向Simulink 建模、控制系统传递函数分析、曲线拟合与参数辨识2. 适用场景与使用边界这套 MATLAB 求解流程适合下面几类场景课程作业和实验报告需要画出 RC 充电、放电曲线并验证时间常数 τ RC 对响应速度的影响。课设和竞赛预研用脚本批量修改 R、C观察不同参数组合下的动态响应辅助选型。信号与系统课程衔接通过传递函数1/(RCs1)和一阶惯性环节理解极点位置与时间常数的关系。控制工程入门RC 电路是最典型的一阶惯性对象后续可以扩展为 PID 控制对象模型。但也有不适合的地方。RC 集总参数模型忽略了导线电感、电容寄生参数和高频效应所以不能用于高频电路、分布参数电路或者开关电源大信号瞬态分析。这类问题需要 SPICE、ADS 或专用电磁场仿真工具。MATLAB 脚本适合验证原理但不能替代专业电路仿真器做完整设计。另外凡是涉及测量数据、实验报告的场合仿真结果需要用真实电路实测或者厂商器件手册的数据做对照不能只拿理论曲线当结论。3. 环境准备与前置条件3.1 MATLAB 版本与工具箱本文用到的函数包括linspace、plot、ode45、tf、step、writetable、saveas这些都属于 MATLAB 基础模块不需要额外安装 Simulink 或 Control System Toolbox。实际上从较早版本开始这些函数就是标配所以无论你电脑上装的是 R2018b 还是更新的 R2024a、R2025a基本都能直接运行。如果要用传递函数tf和step有些低版本需要确认 Control System Toolbox 是否授权。更稳妥的做法是打开 MATLAB 输入以下命令检查。ver % 查看已安装工具箱 which tf % 查看 tf 函数路径如果which tf返回空或者提示未找到说明当前版本没有 Control System Toolbox。这时候可以直接跳过传递函数部分或者干脆使用lsim之前的经典方法自己写exp(-t/tau)画曲线。本质上传递函数也只是一种表达形式不影响 RC 电路物理规律的理解。3.2 工作目录与脚本文件编码建议为每个仿真项目单独建一个目录例如rc_circuit_solver。在 MATLAB 当前文件夹窗口中切换到该目录或者用cd命令cd(D:\work\matlab\rc_circuit_solver)如果脚本中有中文注释在 Windows 上偶尔会遇到编码乱码问题。新版 MATLAB 默认使用 UTF-8 编码基本正常老版本可能按本地编码GBK读取。遇到注释乱码时可以统一把脚本保存为 UTF-8 编码或者在纯计算脚本里直接使用英文注释减少编码导致的显示问题。3.3 常用脚本结构创建一个 MATLAB 脚本文件rc_main.m建议开头用clear; close all; clc;清理工作区避免上一次运行留下的变量影响本次计算。这套清理习惯在写复杂脚本时尤其重要。4. RC 电路模型与求解思路4.1 一阶 RC 电路的微分方程以典型的串联 RC 电路为例电压源 Us 接到电阻 R 和电容 C 组成的串联回路输出取电容电压 Uc。根据基尔霍夫电压定律Us i * R Uc电容电流与电压的关系i C * dUc/dt代入之后得到一阶常系数线性微分方程RC * dUc/dt Uc Us其中时间常数 τ R * C单位是秒。这个 τ 决定了电路响应快慢是 RC 电路最核心的参数。4.2 解析解当 Us 是阶跃电压时方程可以求得解析解。零状态响应初始 Uc0为Uc(t) Us * (1 - exp(-t/tau))当 t τ 时Uc Us * (1 - 1/e) ≈ 0.632 * Us。这是判断曲线是否正确的一个重要参考点。零输入响应Us0初始 UcU0为Uc(t) U0 * exp(-t/tau)当 t τ 时Uc ≈ 0.368 * U0。如果初始电压不为零且电源同时作用可以写成统一形式Uc(t) Us (U0 - Us) * exp(-t/tau)这个公式反映了电容电压从初始值 U0 开始按指数规律趋近稳态值 Us。4.3 数值求解与传递函数思路解析解只适合简单的直流阶跃输入。如果输入是方波、正弦波、分段信号或者电阻值随电压变化解析公式就会失效。这时候用ode45直接对微分方程做数值积分是更通用的方法。另外RC 电路还经常写成传递函数形式G(s) 1 / (RC * s 1)这是一阶惯性环节。在 MATLAB 中可以用tf构造系统再用step观察阶跃响应。这种视角和控制理论中的时间常数分析完全一致适合后续扩展。5. MATLAB 实现解析解与数值解5.1 解析解画充电曲线先写一个简单的解析解脚本。以 R 1 kΩ、C 100 μF 为例时间常数 τ 0.1 s电源 Us 5 V初始电压 U0 0 V。画充电曲线。% 解析求解 RC 充电曲线 clear; close all; clc; R 1000; % 电阻单位 Ω C 100e-6; % 电容单位 F Us 5; % 电源电压单位 V U0 0; % 初始电容电压单位 V tau R * C; % 时间常数单位 s % 从 0 到 5 倍时间常数取 500 个点 t linspace(0, 5 * tau, 500); uc Us (U0 - Us) * exp(-t / tau); figure; plot(t, uc, LineWidth, 1.5); xlabel(时间 (s)); ylabel(电容电压 Uc (V)); title([RC 充电曲线, \tau , num2str(tau), s]); grid on; % 标记时间常数点 hold on; plot(tau, Us * (1 - exp(-1)), ro, MarkerSize, 8, MarkerFaceColor, r); text(tau, Us * (1 - exp(-1)), t \tau, VerticalAlignment, bottom); hold off;运行之后应该看到一条从 0 V 开始指数趋近 5 V 的曲线。在 t 0.1 s 时电压约为 3.16 V也就是 0.632 * Us。5.2 用 ode45 求解微分方程接着用ode45求解同一个方程验证数值解与解析解是否一致。% 用 ode45 数值求解 RC 充电过程 clear; close all; clc; R 1000; C 100e-6; Us 5; U0 0; tau R * C; % 将微分方程整理为 dUc/dt (Us - Uc) / (R*C) rhs (t, Uc) (Us - Uc) / tau; tspan [0, 5 * tau]; % 仿真时间范围 Uc0 U0; % 初始电容电压 [t, Uc] ode45(rhs, tspan, Uc0); % 解析解做对照 Uc_analytic Us (U0 - Us) * exp(-t / tau); figure; plot(t, Uc, b-, LineWidth, 1.5); hold on; plot(t, Uc_analytic, r--, LineWidth, 1.2); hold off; xlabel(时间 (s)); ylabel(电容电压 Uc (V)); legend(ode45 数值解, 解析解, Location, best); title(RC 充电过程ode45 与解析解对比); grid on; % 输出最大误差 max_err max(abs(Uc - Uc_analytic)); fprintf(数值解与解析解最大误差: %e\n, max_err);ode45会自动选择积分步长对大多数非刚性一阶电路问题都能满足精度要求。输出结果中两条曲线应该重合。误差通常在 1e-5 量级以下具体数值取决于RelTol和AbsTol默认设置。5.3 用传递函数做阶跃响应如果安装了 Control System Toolbox还可以用传递函数方式求解。RC 电路的传递函数是1/(R*C*s 1)阶跃响应和解析解本质上是一样的。% 传递函数阶跃响应 clear; close all; clc; R 1000; C 100e-6; tau R * C; sys tf(1, [tau, 1]); figure; step(sys); title(RC 电路单位阶跃响应1/(R*C*s1)); grid on;step默认对单位阶跃输入计算纵坐标最终趋近 1。这里注意因为传递函数默认是单位阶跃不是 5 V 阶跃所以稳态值是 1。如果想要对应 5 V 电源阶跃需要乘以 5figure; step(5 * sys); title(RC 电路 5V 阶跃响应); grid on;6. 批量参数扫描与结果输出6.1 批量扫描电容值实际工程里经常要回答一个问题电容选大一点还是小一点这个问题的答案不需要反复改参数运行脚本可以直接在循环里批量扫描。% 批量扫描不同电容值下的充电曲线 clear; close all; clc; R 1000; Us 5; C_list [10e-6, 100e-6, 470e-6]; % 三种电容值10uF、100uF、470uF colors lines(length(C_list)); figure; hold on; for i 1:length(C_list) C C_list(i); tau R * C; t linspace(0, 5 * tau, 300); uc Us * (1 - exp(-t / tau)); plot(t, uc, Color, colors(i, :), LineWidth, 1.5, ... DisplayName, [C , num2str(C * 1e6), uF]); end hold off; legend(Location, best); xlabel(时间 (s)); ylabel(电容电压 Uc (V)); title(不同电容值下的 RC 充电曲线); grid on;运行后会得到一组曲线电容越大曲线上升越慢达到稳态需要的时间越久。这也符合直观感受时间常数 τ R*C电容越大充电时间越长。6.2 批量保存数据表格与图片批量扫描之后通常还要把结果保存下来方便写报告或者重新绘图。可以用writetable导出 CSV用saveas或exportgraphics保存图片。% 批量计算并保存时间常数和到达稳态时间 R 1000; Us 5; C_list [10e-6, 100e-6, 470e-6]; results table(); for i 1:length(C_list) C C_list(i); tau R * C; % 取电压达到 0.99*Us 的时间作为近似稳态时间 t_settle 5 * tau; results [results; table(R, C, tau, t_settle)]; end disp(results); writetable(results, rc_scan_results.csv); % 保存图片 saveas(gcf, rc_scan_curves.png);如果使用 R2020a 之后的版本更推荐exportgraphics(gcf, rc_scan_curves.png, Resolution, 300)画质更好适合直接放进论文或实验报告。6.3 封装成可复用函数为了后续批量和接口调用把求解逻辑写成一个函数。这个函数输入 R、C、Us、U0、仿真时长和采样点数输出时间向量和电压向量。function [t, uc] rc_solver(R, C, Us, U0, tEnd, npts) % RC_SOLVER 求解一阶 RC 电路时域响应 % 输入: % R 电阻值单位 Ω % C 电容值单位 F % Us 电源电压单位 V % U0 电容初始电压单位 V % tEnd 仿真结束时间单位 s % npts 输出点数默认 500 % 输出: % t 时间向量 % uc 电容电压向量 if nargin 6 || isempty(npts) npts 500; end tau R * C; t linspace(0, tEnd, npts); uc Us (U0 - Us) * exp(-t / tau); end保存为rc_solver.m之后在其他脚本里可以直接调用C_list [10e-6, 100e-6, 470e-6]; trend zeros(length(C_list), 1); for i 1:length(C_list) [t, uc] rc_solver(1000, C_list(i), 5, 0, 0.5, 500); trend(i) uc(end); fprintf(C %.1f uF, t0.5s 时电压 %.3f V\n, C_list(i)*1e6, trend(i)); end7. 功能测试与效果验证写完脚本后不要急着交报告先用几个测试用例把结果验一遍。下面给出一套简单的验证流程。7.1 测试用例 1零状态充电参数R 1 kΩC 100 μFUs 5 VU0 0 V。此时 τ 0.1 s。理论预期t τ 时Uc ≈ 5 * 0.632 3.16 V。t 5τ 时Uc ≈ 5 * (1 - exp(-5)) ≈ 4.966 V基本接近 5 V。用数值结果验证R 1000; C 100e-6; Us 5; U0 0; tau R * C; t_check tau; uc_check Us (U0 - Us) * exp(-t_check / tau); fprintf(t tau 时电压: %.4f V\n, uc_check);如果输出结果在 3.16 V 附近说明公式和代码方向正确。7.2 测试用例 2零输入放电参数R 1 kΩC 100 μFUs 0 VU0 5 V。此时 τ 0.1 s。理论预期t τ 时Uc ≈ 5 * exp(-1) ≈ 1.84 V。t 5τ 时Uc ≈ 5 * exp(-5) ≈ 0.034 V。用rc_solver测试[t, uc] rc_solver(1000, 100e-6, 0, 5, 0.5, 500); idx find(t 0.1, 1, first); fprintf(t0.1s 时放电电压: %.4f V\n, uc(idx));结果应当接近 1.84 V。7.3 测试用例 3ode45 与解析解误差在这个场景里ode45 的解应该无限接近解析解。用最大绝对误差作为指标max_err max(abs(Uc - Uc_analytic)); fprintf(最大绝对误差: %e\n, max_err);如果误差在 1e-6 到 1e-4 之间说明数值积分设置合理。如果误差很大需要检查微分方程是否写错比如符号是否反了时间常数是否算错。完成这三个用例后基本可以确定脚本没有低级错误。8. 批量仿真与 MATLAB 接口调用8.1 脚本内批量仿真批量仿真在上一节已经给出核心就是用循环遍历参数数组。需要注意参数组合多时提前用zeros或cell预分配结果变量避免循环内动态扩展数组。每次循环如果需要绘图用figure新建图窗或者在循环结束后统一绘图。批量结果建议统一写入表格方便事后筛选。8.2 用 Python 调用 MATLAB 函数如果已经有 Python 环境而且希望把 MATLAB 求解能力接到自动测试脚本中可以安装 MATLAB Engine API for Python。安装方式不是 MATLAB 默认自带功能需要在 MATLAB 安装目录中执行安装命令具体路径因版本而异。安装完成后Python 里可以这样调用刚才封装的rc_solver# Python 调用 MATLAB Engine 示例 import matlab.engine eng matlab.engine.start_matlab() eng.cd(rD:\work\matlab\rc_circuit_solver) t, uc eng.rc_solver(1000.0, 100e-6, 5.0, 0.0, 0.5, 500, nargout2) print(最后一个电压值:, uc[-1]) eng.quit()需要提醒一点eng返回的数组是 MATLAB 的 double 数组转成 numpy 数组时要用numpy.array(uc)但注意 MATLAB 默认返回的是列向量维度处理要看实际输出。这种调用方式适合在 Python 侧做参数优化、批量归档、报告生成把 MATLAB 当成计算内核使用。8.3 批量任务目录与失败重试如果是几百上千组参数扫描建议在脚本里增加日志记录。最简单的做法% 批量任务日志 diary(rc_batch_log.txt); for i 1:length(C_list) try [t, uc] rc_solver(R, C_list(i), Us, U0, tEnd, 500); fprintf(参数 %d: C%.3e 成功\n, i, C_list(i)); % 保存结果 T table(t(:), uc(:), VariableNames, {t, uc}); writetable(T, sprintf(result_%d.csv, i)); catch ME fprintf(参数 %d: 失败原因: %s\n, i, ME.message); end end diary off;这样即使某个参数组合有问题也不会中断整个批量任务日志里能看清哪些失败、失败原因是什么。9. 资源占用与性能观察MATLAB 求解 RC 电路属于轻量计算不需要 GPU也不会出现深度学习那种显存占用问题。但批量扫描或者把时间采样点数设得很高时还是要注意内存和执行时间。内存单次linspace(0, 5*tau, 500)和ode45返回结果占用的内存通常只有几 KB 到几十 KB普通电脑可以忽略。但如果把采样点提高到 1e7或者在一个循环里反复拼接矩阵内存会明显上涨。CPUode45会自动改变步长计算量不大。上千组参数扫描时总时间会线性增长可以用tic/toc统计。tic; for i 1:length(C_list) [t, uc] rc_solver(R, C_list(i), Us, U0, tEnd, 500); end elapsed toc; fprintf(批量仿真耗时: %.3f s\n, elapsed);如果发现 ode45 求解很慢可以调整误差容限opts odeset(RelTol, 1e-4, AbsTol, 1e-6); [t, uc] ode45(rhs, tspan, Uc0, opts);误差容限放宽后求解速度会变快但精度会下降。对 RC 电路这类简单问题默认设置通常已经够用。10. 常见问题与排查方法表格形式比较直观问题现象可能原因排查方式解决方案中文注释乱码脚本编码与 MATLAB 默认编码不一致检查文件编码保存为 UTF-8或改用英文注释ode45 报错 Not enough input arguments匿名函数句柄变量未定义检查 rhs 中是否用了工作区变量确保 R、C、Us 在定义函数句柄前已赋值运行后曲线是一条水平直线时间常数太小或仿真时间太长电压变化集中在前段打印 tau 值检查时间范围缩小 tspan 上界或改用对数时间轴曲线快速冲到稳态看不出过渡过程采样点数不足检查 t 的间距增加采样点缩小仿真时间解析解与 ode45 数值解误差很大微分方程写错符号不对打印中间值手推一遍方程把 dUc/dt (Us - Uc)/tau 代入验证无法调用 tf 或 step没有 Control System Toolbox执行which tf改用解析解方式或安装工具箱批量循环运行很慢循环内动态扩展数组重新检查代码用 zeros 预分配结果先存数组后统一写表格输出图片中文乱码或字体异常系统没有对应中文字体检查 figure 字体在绘图前设置字体set(groot,defaultAxesFontName,SimHei)11. 最佳实践与使用建议整个 RC 电路求解看起来简单但工程使用中有几个细节值得记住单位必须统一。电阻用 Ω电容用 F时间用 s。很多人写出错误曲线就是因为电容用了 μF代入公式时忘记换算成 1e-6 单位。先用解析解验证数值解。对于阶跃输入RC 电路解析解很好写先用解析公式算几个特征点再去跑 ode45能快速发现方程写错的问题。参数用结构体管理。当参数很多时建议封装成结构体减少脚本里的全局变量。params.R 1000; params.C 100e-6; params.Us 5; params.U0 0; tau params.R * params.C;批量仿真要留日志。不要只靠“最后一张图”判断结果尽量把每次计算的参数和输出指标写入 CSV。涉及 MATLAB 课程设计提交时代码结构要清晰注释不要只写“这句是画图”而要说明变量含义和公式来源这样才方便复查和答辩。如果脚本要长期使用建议把求解部分和画图部分分开求解函数只负责计算画图脚本负责结果呈现方便测试和复用。12. 总结与下一步RC 电路在 MATLAB 里的求解价值不在于“能画出指数曲线”而在于帮你建立从物理方程到代码实现的完整链路。通过解析解可以快速判断结果方向通过 ode45 可以把方法迁移到任意激励和更复杂的电路方程通过参数扫描可以直观看到 R、C 选型对动态过程的影响。我的建议是先跑通最简单的解析解脚本然后用同样参数跑 ode45确认两者一致再接批量扫描。这样每一步出错都能快速定位。如果你已经完成了这个基础版本下一步可以试着把输入换成方波或正弦波用ode45求稳态响应也可以把电阻换成非线性电阻或者把一阶 RC 扩展到二阶 RLC 电路观察欠阻尼、临界阻尼和过阻尼现象。到时候用类似的分析流程只是微分方程从一阶变成二阶函数封装和批量扫描的思路完全一样。