PSIM9.0学习笔记2——光伏MPPT测试和DLL模块的使用

scott_sl

终于把MPPT调成功了，还是C语言编程能力不足啊，找了几天才找出问题，不过现在也不是很确定是完全正确还是碰巧波形能够跟需要结果相似。在这里呢就还是先介绍一下，然后劳请各位大侠们指正。首先随便搭了一个Boost电路的测试模型，用的是前面介绍过的PSIM9.0自带的光伏板模型，以下就是模型电路图
' N$ D! W5 }2 `& e0 U

MPPT测试模型

& e4 J! i/ v  v5 {& ?$ Z( R
其中MPPT是一个DLL模块，进行MPPT计算并产生PWM波的，最下面一个DLL模块就拿来测试一下效率，是输出功率与最大功率的比。
对于DLL模块的使用，非常简单，首先打开VC++6.0，然后新建一个Win32 Dynamic-Link Library工程，再新建一个C文件，就可以编程了。编好之后在组建里面选择全部重建，完了之后就可以在你所建的工程文件夹Debug里面找到相应的dll文件了，最后在DLL模块里面选择用该dll文件就可以了。
对于编程来讲，首先在c文件里面包含如下语句：
" }8 C4 N; B9 J1 }$ g6 z* p#include <stdio.h>
#include <math.h>
* F6 K4 H% L3 {+ n7 ?__declspec(dllexport) void simuser (double t, double delt, double *in, double *out)
{
}
; L$ i- `$ z' h* P& ^那么在括号内就可以根据自己的需要编写程序了，t是指时间，delt是仿真步长，这两个都在你的仿真电路的时钟控制里面设置，编程不用管。
in和out分别对应输入输出，用的时候直接用用数组形式赋值，比如a=in[0],out[0]=b;就行了，DLL模块有小圆圈的一头从上往下依次是in[0],[1],[2]……out[0],[1],[2]……
那么我编写的MPPT程序如下：
. r3 P. Y. t5 m3 m6 w8 N* l' k
/ h6 O! d5 P4 X% n9 F        Vi=in[0];

3 T1 y1 J: Y. F' ^9 tIi=in[1];

Pi=Vi*Ii;

/ z3 _+ a1 o5 L% \dv=Vi-Vl;                                          //采样电压电流，获得电压，功率增量
6 X+ @! J9 h- d4 L: o
) X$ j: v- X9 u9 Rdp=Pi-Pl;
" z3 A- S  p: ?) Y4 d5 o
if (dp==0)

{
4 a6 i; R% c2 Y' e4 n
2 Y5 @  U  g0 {% Verr=0;
/ O" t0 c7 K" `0 u' v
}

4 ~! O8 X0 |: Y' T$ e7 [else
  a. V. b# b- T; M
{

, Y, O: W$ v2 o4 mif (dp>0)

{
% N4 L2 U% V0 ]1 S$ L1 P) ~7 R
if (dv>0)
5 x" k  G  Q- c) D
5 U; b/ d6 z4 N! s! D{err=0.1;}
  d3 @( c# Q/ B7 i* O4 P+ W1 L) D  j
else                                   //MPPT

{err=-0.1;}

}

9 \5 U- h" `% r4 [( R6 J/ Z6 b6 ?else

{

' M, f; k9 y3 t3 F9 t) s; Aif (dv>0)

{err=-0.1;}

* ^( n7 n+ a5 L) q/ C) `else

{err=0.1;}
0 o- d9 _6 ^+ E% V+ Q
1 ^9 I+ O5 k9 _7 |+ v; y( g}
$ ?& F5 }5 \5 _6 V4 E( _: z& J0 j
- p5 k2 M. C; X0 x+ d) J}

我用的是最简单的PO法，然后设每次Vref=Vin+0.1或者-0.1，所以上面err就是一次计算完了之后的误差值，送入下面的增量型PI控制
2 Q: o% {8 Q, P5 e7 \' T. n  h                dVout=20*(err-errf)+err;                 
0 v$ u$ X/ ]  H# n+ j( _
- I- }& |8 g5 G' j: ^' U3 V% vif (fabs(dVout)>4.1)
+ P) @* T0 W/ m3 \- _2 v  |& S' r1 H
5 \6 ]( j2 Q" T; Q{

if ((dVout)>0)
6 f2 \$ {$ Y( f% t8 h( a, p6 _( v
- A: U7 z3 ~6 f3 |# V: g3 l. a{dVout=4.1;}
9 F0 L: _# C4 t3 G. c' n6 q
else
: r( L$ e  w. I; |; N& Z! @6 N
0 C- J$ X' _: Z{dVout=-4.1;}
) }  b: D5 F  N( j. u
}

7 K3 d! i, @& [# G3 f0 |Vtemp=Vout+dVout;                          //PI
" y2 D& N; N. @' _


if (Vtemp>10)                                

) K! z& T1 R/ i' O9 H& ]8 p{Vtemp=10;}
6 J- a0 d- V! ^
$ c3 N# c( P6 ielse

{

2 N. A! w8 c7 \' U% S; ^if (Vtemp<0)
. b. A" s# L: ^0 K/ l3 \5 B# u, E
{Vtemp=0;}

else
! S+ d( I( m& B( b
+ b, K/ F# o# C: R3 d/ ?{Vtemp=Vtemp;}

6 _9 ?4 ]# p" ?8 [  G! W. l3 G}
- X2 I0 d, a6 H1 Q9 d' b& Q
Vout=Vtemp;

errf=err;

  p  ^& D1 x" Y! Z& ^: ]! t用来得到比较值的变化，最后跟锯齿波比较之后输出占空比。
以上程序跑完之后的仿真结果如下所示：

: q: W; j7 ]+ T6 G. b2 s; g
感觉还是蛮精确的，效率算下来大概在99.87%，但是我总觉得哪儿有问题，不是太对，上面误差输出的值就只有0.1，-0.1,0三个值
5 c+ i9 ]# n0 k: y4 I5 b不过不知道问题到底原因是什么，仿真效果还蛮好嘛，呵呵
大家给我指导一下呗，多谢！~
, p" P$ i, g6 E

白露为霜

楼主上面的各个变量具体的定义是什么呢？还有就是为什么要送进增量型PI控制中去呢？

丝雨如愁

不错不错 一起学习

mingyu

个人兴趣关注一下，希望有这方面的高人指点！

kongxincailei

啥时候上线呀，交流下可以吗？

天锁斩月

希望经常交流

天锁斩月

你说的那个误差值应该是3点震荡引起的吧，还有一种2点震荡的。

eden7

不错不错~

mainbody

好帖！学习下强大的操作！

redplum

好详细的资料啊

liyuwen

你好，我想请问一下最后一步“最后在DLL模块里面选择用该dll文件就可以了”，是怎样选择啊？

tonyemperor

感谢兄弟分享