太阳能光伏系统的Matlab Simulink仿真搭建与探索-程序员充电站

Matlab simulink仿真太阳能光伏系统搭建包含有功、无功模块、3-2变换等。通过改变id、iq数值可得出不同的波形，图2为纯有功，a相电压于相电流重合；图3为纯无功时，a相电压与a相电流相差90度，可通过改变id、iq大小控制功率因数角的大小。图4、5、6为在某一时间段突然增加负载，直流母线电压、a相电压与相电流的变化、功率变化。

在太阳能光伏系统的研究与设计中，Matlab Simulink是一个强大的工具。今天就来聊聊如何搭建一个包含有功、无功模块以及3 - 2变换等的太阳能光伏系统仿真模型。

系统搭建

有功、无功模块

在Simulink中搭建有功和无功模块是关键部分。有功功率（$P$）和无功功率（$Q$）在电力系统中有着重要意义。从理论上来说，$P = VI\cos\varphi$，$Q = VI\sin\varphi$，其中$V$是电压，$I$是电流，$\varphi$是功率因数角。

在Simulink里，我们可以通过控制电流的直轴分量$id$和交轴分量$iq$来实现对有功和无功功率的调节。以一个简单的电流控制模块为例（假设使用S函数实现）：

function [sys,x0,str,ts] = space_vectors(t,x,u,flag) switch flag, case 0, [sys,x0,str,ts]=mdlInitializeSizes; case 3, sys=mdlOutputs(t,x,u); case {2,4,9} sys = []; otherwise error(['Unhandled flag = ',num2str(flag)]); end function [sys,x0,str,ts]=mdlInitializeSizes sizes = simsizes; sizes.NumContStates = 0; sizes.NumDiscStates = 0; sizes.NumOutputs = 3; sizes.NumInputs = 2; sizes.DirFeedthrough = 1; sizes.NumSampleTimes = 1; sys = simsizes(sizes); x0 = []; str = []; ts = [0 0]; function sys=mdlOutputs(t,x,u) id = u(1); iq = u(2); % 这里省略复杂的空间矢量计算 % 简单假设直接输出三相电流分量 ia = id; ib = -0.5*id + sqrt(3)/2*iq; ic = -0.5*id - sqrt(3)/2*iq; sys = [ia;ib;ic];

这段代码接收$id$和$iq$作为输入，然后通过简单的计算输出三相电流分量。实际应用中，会涉及到更复杂的空间矢量调制等算法，但基本原理类似。通过改变输入的$id$和$iq$数值，我们就能灵活调整有功和无功功率。

3 - 2变换

3 - 2变换在三相系统分析中非常重要，它能将三相静止坐标系下的量转换到两相旋转坐标系下，方便分析和控制。其变换矩阵如下：

\[ C_{3s/2r} = \sqrt{\frac{2}{3}}\begin{bmatrix} \cos\theta & \cos(\theta - \frac{2\pi}{3}) & \cos(\theta + \frac{2\pi}{3}) \\ -\sin\theta & -\sin(\theta - \frac{2\pi}{3}) & -\sin(\theta + \frac{2\pi}{3}) \end{bmatrix} \]

在Simulink中，可以通过一些现成的模块或者自定义模块来实现这个变换。比如利用Matlab Function模块：

function [id,iq] = abc_to_dq(ia,ib,ic, theta) alpha = (2/3)*(ia - 0.5*ib - 0.5*ic); beta = (2/3)*(sqrt(3)/2*ib - sqrt(3)/2*ic); id = alpha*cos(theta) + beta*sin(theta); iq = -alpha*sin(theta) + beta*cos(theta);

这段代码实现了从三相电流$ia$，$ib$，$ic$到直轴电流$id$和交轴电流$iq$的变换。通过这个变换，我们能更好地控制有功和无功功率，因为在两相旋转坐标系下，$id$主要控制有功功率，$i_q$主要控制无功功率。

仿真结果分析

改变$i_d$，$i_q$得出不同波形

当我们改变$id$和$iq$的数值时，系统输出的波形会发生明显变化。图2展示的是纯有功的情况，此时a相电压与相电流重合。这是因为在这种情况下，功率因数角$\varphi = 0$，$Q = 0$，只有有功功率在起作用。从前面提到的$P = VI\cos\varphi$公式也能理解，$\cos\varphi = 1$时，就是纯有功状态。

而图3呈现的是纯无功的情况，a相电压与a相电流相差90度。此时$\varphi = 90^{\circ}$，$\cos\varphi = 0$，$P = 0$，只有无功功率存在。我们可以通过灵活改变$id$和$iq$的大小，轻松控制功率因数角$\varphi$的大小，从而实现对系统功率特性的精确调整。