企业官网建设流程全解析

后台经常有人问我光伏电池模型怎么做，正好我最近整理了一套自己一直在用的 MATLAB/Simulink 光伏电池模型，带说明文档，还有配套的仿真绘图代码。这套东西不复杂，适合在校学生、刚接触光伏仿真的工程师，以及要做 MPPT 算法验证、并网逆变器前级仿真的人直接拿去用。它不是让你在 Simscape 里拖个现成模块然后一脸懵，而是从光伏电池的等效电路出发，把方程、Simulink 框图、批量仿真脚本、绘图模板串成一条完整链路，保证你能理解每个环节在干什么，也能照着复现。

我在实际项目中用过 Simscape 自带的光伏模块，也自己搭过模型。自建模型的优势在于：参数可控、逻辑透明、不依赖额外工具箱授权，并且非常方便改成“组件级、阵列级”或者“带温度的实时变化”模型。下面这套就是我平常教学和自己做仿真时最常用的一套版本，我把细节都展开讲一遍，包括建模思路、参数计算、Simulink 搭建方法、批量绘图代码、文档组织方式，以及在仿真中很容易踩的坑。

1. 建模前的准备：光伏电池的核心方程与参数概念

1.1 单二极管等效电路到底在讲什么

光伏电池的电气特性，工程上最常用的是“单二极管五参数模型”。它的等效电路可以这样理解：光照在 PN 结上会产生光生电流 Iph，这个电流可以看成是一个恒定电流源；PN 结本身又是一个二极管，存在反向饱和电流 I0；半导体材料和电极接触会有接触电阻，对应串联电阻 Rs；电池边缘漏电和内部缺陷会造成漏电流路径，对应并联电阻 Rsh。

根据基尔霍夫电流定律，输出电流 I 与输出电压 V 的关系式为：

I = Iph - I0 * [exp((V + IRs) / (a * Vt)) - 1] - (V + IRs) / Rsh

其中 Vt = Ns * k * T / q，是热电压的组件级表达。Ns 是单体电池串联数量，k 是玻尔兹曼常数，T 是电池温度（开尔文），q 是电子电荷，a 是二极管理想因子。

很多新手看到这个方程第一反应是“好复杂”，但它拆开看并不难：第一项 Iph 是光照产生的电流，第二项是二极管分流掉的电流，第三项是并联电阻漏掉的电流。你从外部能测到的输出电流，就是光生电流减去这两部分损耗。

我在实际教学中经常打一个比方：光生电流是水管里流进来的水，二极管和并联电阻是两个漏水口，串联电阻是管道里的摩擦阻力。你希望水龙头的出水尽量大，就得让两个漏水口尽量小、管道阻力尽量低。这个类比能帮你快速理解后面调参的方向：Rs 大，填充因子差；Rsh 小，曲线在最大功率点附近塌陷严重。

1.2 从 datasheet 拿到哪五个关键参数

搭建模型的第一步不是写代码，而是明确你仿真对象是哪种组件。拿一块常见 260W 多晶硅光伏组件来说，datasheet 上标称参数如下：

这里要强调，短路电流、开路电压、最大功率点电流和电压，都是标准测试条件 STC 下的值，也就是辐照度 1000 W/m²、电池温度 25°C、AM1.5 光谱。你可能会问，我们建模型不是直接需要 Rs、Rsh、I0 这些参数吗？datasheet 为什么不直接给？

原因很简单：厂商不会在铭牌上标这些内部参数，它们是“隐含参数”，只能通过标准参数反推。反推的思路是让模型在三个关键点与 datasheet 完全吻合：短路点（V=0, I=Isc）、开路点（I=0, V=Voc）、最大功率点（V=Vmpp, I=Impp）。只要这三个点能对上，整个 I-V 曲线在中段尤其是最大功率点附近就不会差太多，后续做 MPPT 仿真的可靠性才有保障。

温度系数 Ki 和 Kv 也不能忽略。它们决定了辐照度和温度变化时 I-V 曲线的平移趋势。辐照度主要影响电流，温度主要影响电压，这两个规律在后面的绘图阶段你会在曲线上直观看到。

1.3 工程简化模型和隐含数解法

在 Simulink 里直接解那个超越方程，最关键的是 Iph、I0、Rs、Rsh 这四个内部参数的确定。工程上常用如下简化估计：

Iph = (Isc + Ki * (T - Tref)) * G / Gref

I0 = (Isc + Ki * (T - Tref)) / (exp((Voc + Kv * (T - Tref)) / Vth) - 1)

其中 Vth = Ns * a * k * (T + 273.15) / q。至于 Rs 和 Rsh，很多人直接取 Rs = (Voc - Vmpp) / Impp，Rsh 则取一个较大的值比如 300~1000 欧姆。这种取法在最大功率点附近精度是可以接受的，但如果你要做高精度曲线拟合，建议用最小二乘或者粒子群去辨识。

这里有个数学细节必须注意：方程中 V 和 I 都在指数项和分母项里，属于隐式方程，没办法直接写成 I = f(V) 的显式形式。所以在 MATLAB Function 里我习惯用 fzero 求根，把问题转成：

求 I，使得 Iph - I0 * (exp((V + IRs) / Vth) - 1) - (V + IRs) / Rsh - I = 0

初始猜测值设为 0.9 * Iph。求根区间 [0, Iph] 也可以，因为光伏电池正常工作时输出电流不会超过光生电流。后文代码会直接展示这个做法。你也可以用牛顿迭代自己写几行，效果差不多，但 fzero 在数值稳定性上更省心。

2. Simulink 模型搭建：从方程到仿真框图

2.1 建模方式的取舍：MATLAB Function 模块 vs Simscape 自带模型

聊到 Simulink 光伏建模，很多人第一反应是 Simscape Electrical 里的 PV Array 模块。确实，官方模块在并网逆变器、微电网这种需要和电力电子拓扑实时交互的场景里很方便，它内部已经帮你算好了所有参数，还能直接设置组串并联数、MPPT 接口等。

但我还是建议你掌握“从零自建”的方式，至少在研究和教学场景里，它的优势更明显：

我的建议是：如果你是做光伏电池特性研究、MPPT 算法验证，或者给学生上课讲原理，自建模型是首选；如果你在做完整的并网系统、需要把模型放在一个大的 Simscape 物理域网络里跑，那就老老实实用官方模块，别拿自建模型去和电感电容直接连物理端口。两种方式各有各的用武之地。

2.2 用 MATLAB Function 模块写光伏电池方程

实操步骤很简单，先新建一个 Simulink 模型，我这里命名为 pv_curve.slx。然后从常用模块库里拖入几个关键模块：

• MATLAB Function（用于写光伏电池方程）
• 斜坡信号 Ramp 或 Signal Builder（用于产生电压扫描信号）
• Constant（用于给定辐照度 G 和温度 T）
• Product（用于计算功率 P = V * I）
• To Workspace（用于导出 V、I、P 到工作区）
• Scope（用于快速查看曲线）

MATLAB Function 里核心代码如下：

function I = pv_cell_engine(V, G, T) % 组件参数（260W 多晶硅组件示例） Isc = 8.21; % 短路电流 A Voc = 32.9; % 开路电压 V Impp = 7.61; % 最大功率点电流 A Vmpp = 26.3; % 最大功率点电压 V Ns = 54; % 串联电池片数 Ki = 0.0032; % 短路电流温度系数 A/°C Kv = -0.123; % 开路电压温度系数 V/°C a = 1.3; % 二极管理想因子 Tref = 25; % 参考温度 °C Gref = 1000; % 参考辐照度 W/m² k = 1.3806503e-23; q = 1.60217646e-19; Tk = T + 273.15; % 转开尔文 Vth = Ns * a * k * Tk / q; % 估算内部参数 Iph = (Isc + Ki * (T - Tref)) * G / Gref; I0 = (Isc + Ki * (T - Tref)) / (exp((Voc + Kv * (T - Tref)) / Vth) - 1); Rs = (Voc - Vmpp) / Impp; Rsh = Vmpp / (Iph - I0 * (exp((Vmpp + Impp * Rs) / Vth) - 1) - Impp); % 隐式方程求根 fun = @(i) Iph - I0 * (exp((V + i * Rs) / Vth) - 1) - (V + i * Rs) / Rsh - i; I = fzero(fun, 0.9 * Iph); end

这里有几个细节要交代：输入 V 是由外部斜坡信号产生的 0 到 35V 电压扫描，G 和 T 是两个常量。输出 I 后，用 Product 模块把 V 和 I 相乘得到 P，然后把三者一起送进 To Workspace 模块。

我实际测试过，用 ode45 默认步长跑一次完整扫描只需要零点几秒，效率完全够用。不过注意 fzero 在每个时间步都会被调用一次，如果你最后要把仿真时长拉到几十秒，计算量会线性增加。你在批量仿真之前先把 StopTime 设置合理，比如电压扫描速率是 10V/s，扫描 0~35V 需要 3.5 秒，留点余量设成 4 秒就够了。

2.3 代数环问题：仿真中必踩的一个坑

刚开始搭这个模型的时候，最容易碰到的问题就是代数环。出现代数环的原因是：如果你用“可变负载电阻”来连接模型，输出电压 V = I * Rload，而 I 又依赖 V，模型里就出现了瞬时信号循环。Simulink 会在诊断窗口报一个警告，提示存在 Algebraic Loop，求解时会增加迭代，严重的会导致仿真缓慢甚至报错。

解决方式有三条路，我按推荐程度排序：

第一种，也是最推荐的方式：用“电压源扫描法”而不是“可变电阻法”。也就是说，直接从外部给定电压 V，然后模型计算对应的 I。这样信号流向是单向的，从电压到电流，没有循环依赖，根本不会产生代数环。我前面搭的框图就是这种方式。

第二种，如果你的实际场景是光伏电池接电阻负载或者 DC/DC 变换器，必须让电压和电流双向耦合，那就在反馈回路上加一个 Memory 模块或者单位延迟 Unit Delay，把代数环“打断”。本质上是用上一个仿真步长的输出量代替当前步长的瞬时值，代价是会引入一个步长的延迟，但只要步长足够小，对结果影响可忽略。

第三种，把整个方程整合进一个更大的状态空间，用 Simulink 的求解器统一处理。这个方法理论最优，但对使用者要求高，实际很少人用它来搭光伏模型。

我的经验是：做 I-V 曲线扫描时一律用电压源扫描法；做 MPPT 或并网仿真时，在反馈回路里放一个 Memory。这两种方式我都验证过，仿真结果一致，稳定性很好。

2.4 封装子系统与参数接口设计

模型调通以后，不要嫌麻烦，把它封装成子系统，加上 Mask。这一步很多人会跳过，但等你要做温度扫描、辐照度扫描，或者换一块电池板参数时，封装的价值就体现出来了。

封装操作方法：选中所有模块右键 Create Subsystem，然后右键子系统 Mask > Create Mask。在 Mask Editor 里添加参数，比如：

• 短路电流 Isc
• 开路电压 Voc
• 最大功率点电流 Impp
• 最大功率点电压 Vmpp
• 串联电池数 Ns

这样你每次只需要双击子系统，在弹出的对话框里换参数，而不需要打开 MATLAB Function 去改代码。尤其当你有多个相同结构的子阵列时，封装后的效果就是“参数填进去，直接出结果”。

我习惯在封装说明里也写清楚：输出电压范围不超过 Voc 的 1.2 倍，辐照度单位是 W/m²，温度单位是 °C。这能避免同事或学生后续使用时因为单位问题把曲线跑飞。

3. 仿真绘图代码：把仿真结果变成易读的曲线

3.1 使用 sim 函数批量扫描不同工况

单条 I-V 曲线只是在验证模型，真正有说服力的图是“不同辐照度下的曲线族”和“不同温度下的曲线族”。手动在 Simulink 里改常量再点运行当然可以，但那样效率太低，而且不容易统一管理数据。正确做法是写一个仿真脚本，在循环里修改模型参数，调用 sim 函数批量跑。

下面这段脚本是我常用模板。假设模型里辐照度常量模块名为 G，温度常量模块名为 T，To Workspace 输出变量名分别为 Vout、Iout、Pout。

mdl = 'pv_curve'; load_system(mdl); G_list = [1000 800 600 400 200]; results = struct(); for i = 1:length(G_list) set_param([mdl '/G'], 'Value', num2str(G_list(i))); set_param([mdl '/T'], 'Value', '25'); simOut = sim(mdl, 'StopTime', '4'); results(i).G = G_list(i); results(i).V = simOut.Vout.Data; results(i).I = simOut.Iout.Data; results(i).P = simOut.Pout.Data; end

这段代码有几个容易忽略的细节。set_param 修改的 Value 必须是字符串，所以要用 num2str 转换。sim 函数的返回值建议用 Simulink.SimulationOutput 对象接收，这样数据格式统一，后面绘图不容易出错。建议在命令行先运行 load_system，避免每次循环都在加载模型上浪费时间。

另外一个实战经验是：批量仿真之前把模型的“快速重启 Fast Restart”功能打开，可以在 Tools > Run > Fast Restart 里设置。这样循环里的多次仿真会共享编译结果，仿真速度能提高不少。尤其是你后面要做几十组参数扫描时，这个设置能节省大量时间。

3.2 绘制 I-V 曲线与 P-V 曲线的完整代码

数据拿到手以后，下一步就是画图。我习惯把 I-V 曲线和 P-V 曲线放在同一张图的两个子图里，或者用双 y 轴放在同一张图。双 y 轴的好处是能直观看到最大功率点对应的位置，但需要注意坐标轴刻度范围设定，不然曲线会挤在一起。

figure('Color', 'w', 'Position', [100 100 900 600]); colors = lines(length(results)); hold on; yyaxis left; for i = 1:length(results) plot(results(i).V, results(i).I, 'Color', colors(i, :), 'LineWidth', 1.6); end ylabel('电流 I (A)'); xlabel('电压 V (V)'); ylim([0 9]); xlim([0 35]); yyaxis right; for i = 1:length(results) plot(results(i).V, results(i).P, '--', 'Color', colors(i, :), 'LineWidth', 1.2); end ylabel('功率 P (W)'); ylim([0 280]); legend(arrayfun(@(x) sprintf('G=%d W/m²', x.G), results, 'UniformOutput', false), 'Location', 'northeast'); grid on; set(gca, 'FontName', 'Helvetica', 'FontSize', 11);

在 MATLAB 中文环境下，如果图例或坐标轴要显示中文，你需要在画图前设置系统的默认字体，否则中文会显示成方块。常用的设置命令是：

set(0, 'defaultAxesFontName', 'SimHei'); set(0, 'defaultTextFontName', 'SimHei');

如果是在英文版 MATLAB 上运行，建议直接用英文标签，省得字体问题影响阅读。我平时交付给别人的脚本里会保留中英双语变量注释，但图上的文字统一用中文，因为国内团队看得更直观。

图片导出建议用 print 命令而不是直接截图，这样输出分辨率可控。比如：

print(gcf, 'pv_curve_family.png', '-dpng', '-r600');

600 DPI 足够期刊插图使用。如果你要把图放进 Word 报告，300 DPI 就够了，文件体积还小。

3.3 不同温度/辐照度对比曲线的一次成型技巧

做温度扫描的逻辑完全一样，只需要把常量 T 设为循环变量，G 固定为 1000。温度范围我一般取 0°C、25°C、50°C、75°C 这四档，覆盖大部分户外运行场景。

有个技巧是：把辐照度扫描和温度扫描的绘图函数封装成独立函数 draw_iv_curve(results, title_str)，这样两种工况的绘图脚本共用一个函数，代码量小，风格也统一。我就是这么组织的：

function draw_iv_curve(results, title_str) ... end

另外，在画 P-V 曲线时，我习惯同时标注出每个工况的最大功率点。计算方法很简单，用 max 找到最大功率和对应的电压电流索引：

[Pmax, idx] = max(results(i).P); Vmpp_calc = results(i).V(idx); Impp_calc = results(i).I(idx);

然后在图上用 scatter 画一个实心点，再用 text 标注数值。这样审稿人和领导一眼就能看出你的模型最大功率和 datasheet 标称值的偏差。

4. 说明文档的组织方法：让模型可以被别人复现

4.1 一个能复现模型的说明文档应该包含哪些章节

项目标题里特意强调“附带说明文档”，这一点很多人忽略，但它恰恰是这套模型价值能不能发挥出来的关键。我自己收过很多没有文档的模型文件，拿到手根本不知道该改哪里、参数在哪里设置、输出怎么处理。所以我在交付模型时，说明文档固定包含六个部分：

1. 项目说明：这套模型是什么、解决什么问题、适合用在哪些场景。
2. 文件清单：模型文件、函数文件、脚本、数据文件的目录树，以及每个文件的作用。
3. 快速开始：从打开 MATLAB 到复现一张 I-V 曲线的完整操作步骤。
4. 参数表：所有参数的名称、符号、初始值、单位、来源。
5. 仿真结果示例：标况下的 I-V/P-V 曲线图、不同工况对比图。
6. 注意事项与常见问题：包括版本要求、求解器设置、常见报错处理。

你可以直接用 Word 或者 Markdown 写，也可以做成项目里的 README.md。我的习惯是既放一份 Markdown 在项目根目录，再导出一份 PDF 放进 doc 文件夹，方便给不看代码的人快速了解。

说明文档里最关键的是“快速开始”。不要假设接手你模型的人很熟悉 MATLAB。你要写清：第一步打开 pv_curve.slx，第二步点击运行，第三步在命令行执行 plot_curves.m，第四步看输出图片。每多写一个步骤，用户自己摸索的时间就少十分钟。

4.2 命名规范和变量管理

模型里涉及的参数很多，如果命名随意，后期维护非常痛苦。我的命名规范是这样的：

• 物理量直接用英文字母组合，例如 Isc、Voc、Impp、Vmpp、Ns、Ki、Kv。
• 多个参数可以用结构体打包，例如 params.Isc、params.Voc，这样在 MATLAB Function 里只需要传一个 params 进去。
• Simulink 模型名称和文件名称统一，不用中文和空格。
• 脚本名称用动词开头，例如 plot_curves.m、run_sweep.m、export_figures.m。

另外要注意，Simulink 模型内部的信号名也很重要。双击信号线可以输入名称，比如 V_pv、I_pv、P_pv，这样 To Workspace 导出以后的工作区变量名会带上信号名，方便识别。很多人忽略这个细节，最后导出数据后还要手动猜哪个变量是电压哪个是电流，非常浪费时间。

4.3 关于版本和工具箱的注意事项

我搭建这套模型用的是 R2021a 版本，核心功能在 R2020b 以上都能正常运行。需要特别注意的是，MATLAB Function 模块在早期版本叫 Embedded MATLAB Function，界面略有不同，但代码写法基本一致。如果你的 MATLAB 版本是 2016a 之前的老版本，建议先升级，否则有些模块接口会不兼容。

工具箱方面，这套自建模型只需要 MATLAB 基础包和 Simulink，不需要额外安装 Simscape、Simscape Electrical 等授权工具箱，这对很多没有购买全套工具箱的学生来说非常友好。但如果你要用官方 PV Array 模块对比验证，那你需要确认许可证里包含 Simscape Electrical。

版本兼容性还会影响绘图导出的字体，例如不同操作系统下 SimHei 字体的支持情况不同。为了避免交付后图片字体错乱，我通常在脚本开头统一检查是否有对应字体，没有就自动切换成系统的默认字体。

5. 常见问题与排查技巧实录

5.1 问题1：MATLAB Function 模块报错"Not enough input arguments"

这个问题在我刚把代码拷贝进 MATLAB Function 模块时经常遇到。原因通常是模块的输入引脚数量和函数声明里的输入参数数量不一致。比如你在函数声明里写了 function I = pv_cell_engine(V, G, T)，但模块只有两个输入引脚，那就必然报错。

解决办法：双击 MATLAB Function 模块，在编辑器上方的“编辑数据”里检查输入名称，确保与函数参数名一致，并且数量对得上。另外，输入的数据类型也要留意，如果 G 是 uint8，而函数内部做浮点运算，可能会触发类型转换错误，建议在 Constant 模块的输出数据类型里直接设成 double。

5.2 问题2：仿真曲线出现毛刺或负电流

出现毛刺，大概率是求解器步长过大导致的。光伏电池的 I-V 曲线在最大功率点附近转折很剧烈，如果变步长求解器没有捕捉到这个转折，就会出现数值振荡。解决办法是打开 Configuration Parameters，把 Solver 类型保持 Variable-step，算法选 ode45，同时把相对误差 Relative tolerance 改小，例如 1e-5，最大步长 Max step size 设置成 0.01 或更小。

如果出现负电流，一般是电压扫描超过了开路电压太多，导致二极管反向偏置，模型在某些数值条件下算出很小的负电流。解决办法很简单，把扫描电压上限设成 Voc 的 1.05 倍就够了，不用扫到 40V、50V。同时 fzero 求根区间最好限制在 [0, Iph]，防止在非物理区间迭代发散。

5.3 问题3：最大功率点偏差过大

这是最需要重点排查的问题。模型在短路点、开路点一般都很准，但最大功率点如果不准，通常原因有三个：第一个是 Rs 和 Rsh 估计过于粗糙，Rs 增大最大功率点功率显著下降，Rsh 减小曲线在最大功率点左侧塌陷；第二个是二极管理想因子 a 取值不合适，常见晶硅组件 a 在 1~1.5 之间，但不同厂商工艺有差异；第三个是温度系数 Ki、Kv 设置错误，导致温度变化时最大功率点偏移。

我的经验是先用 datasheet 给的最大功率点反推一遍 Rs 和 Rsh，然后在仿真里扫描 a 的取值，看哪个值让最大功率点功率误差最小。如果误差还是超过 3%，那就得用优化算法辨识参数。别怕这个步骤麻烦，参数辨识一次完成后，模型在宽范围内都稳定。

5.4 问题4：批量仿真数据覆盖或无法保存

批量循环里如果每次都用同一个变量名接收 sim 输出，下一次循环会把上一次的结果覆盖掉。所以必须用结果数组或元胞数组保存。我常用 results(i).V = simOut.Vout.Data 这种方式，每一组数据单独存放。还有就是 To Workspace 模块的输出格式，建议选“Timeseries”，不要用“Array”或“Structure With Time”，因为 Timeseries 自带时间和数据，后续提取更方便。

如果 sim 函数最后一次运行报错，检查模型编译阶段是否有问题，可以在循环外加 try-catch，把出错工况打印出来。这样即使某个辐照度下模型收敛失败，也不会中断整个批量扫描。

5.5 常见排查速查表

6. 一眼定睛：接线信号方向与 To Workspace 的数据类型

最后再分享两个我踩过不少坑的小细节。

第一个是关于信号方向。在 Simulink 里从电压源到 MATLAB Function 模块再到 Product，信号方向一定是从左到右单向流动。如果你不小心把 G 的 Constant 也接到函数输出端，或者连线交叉导致信号回流，仿真结果会非常诡异。检查方法是选中连线，看是否有红色虚线，那是未连接或连接错误的提示。

第二个是 To Workspace 模块的数据类型。在设置界面里有个 Save format 下拉菜单，我推荐选择 Timeseries。很多人默认会选 Structure，等到工作区里发现 simOut.Vout.Data 取不到数据才回头改。这个小坑浪费了我不少时间，所以提前写在这里。

我自己使用这套模型时，已经形成了固定的项目文件结构：模型文件放根目录，脚本放 scripts 文件夹，生成的图片放 figures 文件夹，说明文档放 doc 文件夹。每次做 MP PT 算法对比时，只需替换外围控制器部分，光伏电池模型本身基本不用动。你现在拿这套代码去复现，只要按顺序把步骤走完，半小时内就能跑出完整的 I-V/P-V 曲线族。如果后续要扩展，加上一个 Boost 变换器和扰动观察法，就能直接开始做最大功率点跟踪的仿真实验了。