马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
终于把MPPT调成功了,还是C语言编程能力不足啊,找了几天才找出问题,不过现在也不是很确定是完全正确还是碰巧波形能够跟需要结果相似。在这里呢就还是先介绍一下,然后劳请各位大侠们指正。首先随便搭了一个Boost电路的测试模型,用的是前面介绍过的PSIM9.0自带的光伏板模型,以下就是模型电路图
- r t( a7 O( K# d
MPPT测试模型
" w# Q7 O- ?3 {- Y! h: j% E
! D: F9 C4 m" w' i6 @# n
其中MPPT是一个DLL模块,进行MPPT计算并产生PWM波的,最下面一个DLL模块就拿来测试一下效率,是输出功率与最大功率的比。
, z9 p% R8 s/ A* D p6 B 对于DLL模块的使用,非常简单,首先打开VC++6.0,然后新建一个Win32 Dynamic-Link Library工程,再新建一个C文件,就可以编程了。编好之后在组建里面选择全部重建,完了之后就可以在你所建的工程文件夹Debug里面找到相应的dll文件了,最后在DLL模块里面选择用该dll文件就可以了。
" r/ ?$ T7 s! s9 [% _7 f7 f 对于编程来讲,首先在c文件里面包含如下语句: ) m# _0 v& ]: G7 i3 N1 _. x
#include <stdio.h>. E9 ?* N8 a/ Z) X+ R7 a
#include <math.h>
) t$ u+ s* Z8 {1 { | __declspec(dllexport) void simuser (double t, double delt, double *in, double *out)2 q) J% B' w2 l/ @/ |( r
{
% _; f; D4 |+ ~. D# o4 w# f4 o }
, l4 C# R7 W0 W$ t5 U! |* W6 } 那么在括号内就可以根据自己的需要编写程序了,t是指时间,delt是仿真步长,这两个都在你的仿真电路的时钟控制里面设置,编程不用管。# m) D. c" c& t# |) n
in和out分别对应输入输出,用的时候直接用用数组形式赋值,比如a=in[0],out[0]=b;就行了,DLL模块有小圆圈的一头从上往下依次是in[0],[1],[2]……out[0],[1],[2]……0 l7 V- X! a% X- A1 R. b
那么我编写的MPPT程序如下:: i5 k; |. b4 N% T" a( E
2 A' f( x" |- r' h- w+ S6 ] Vi=in[0];: K5 o" D2 E+ E! H6 q/ `
% _7 M6 P7 q" M% z! O* K Ii=in[1];6 ^8 U% r) Q0 e8 g+ X7 K
, u9 U4 N, Q/ a1 `% [ Pi=Vi*Ii;( h! ~3 R2 x% S$ O) ]* T
6 Y- S R# p7 @: y% L% j dv=Vi-Vl; //采样电压电流,获得电压,功率增量% ^% M5 T8 v6 }' i6 r
1 z. o& A; o' d6 R
dp=Pi-Pl;$ }5 Z' L6 C7 R5 C
1 N9 s1 K% u. ?$ n" ? if (dp==0)) a& [+ n$ C" g0 v3 w7 h0 ?
1 D/ @5 P% ~5 ?, `/ ?9 e; ] i
{. E1 g2 r) o# }3 W" _ c A
# i; ?/ ~6 n$ z0 b. I
err=0;
/ \4 W) P1 r6 B8 ~) k1 D7 ]1 H
9 ] H- S7 T2 j8 N3 M2 T- {# E3 G" V } z, d5 O0 H0 r+ L% B3 K0 l4 y
% d8 Z1 W, Y3 t
else2 \! |! R7 }. ?: f) ]* [- C$ V
. M/ I" ?( o5 Y$ t% I. f9 V/ C8 m5 ?
{/ G' [% y+ i. w ]
& F- y3 C* b( X$ S* |7 }0 B
if (dp>0)& j* A+ Z# m* g+ t; u7 M% `
3 I% L) y- _( _0 _
{
% O! L, B$ N6 w4 u4 s3 L
3 i* [$ M3 d1 K: _3 Z& o6 } if (dv>0)6 ?7 X: B2 U- K. U0 U# }/ H
6 J* ~; {% m4 G$ n
{err=0.1;}2 x/ Q) A6 E/ f
|* G/ i) I% ]1 l else //MPPT+ @. ^. }. S1 M
- V* ^7 H& I4 @: c. q
{err=-0.1;}
5 `8 j! }* V2 K; x/ N : `+ D. T' \4 k% ^
}
, L$ j7 t0 S: O, S) i/ P
' n$ k' s' x# [9 ` else
/ t3 J) C0 l/ S" g7 K
; k: ~% O0 q3 X: c3 X {
9 ~4 h# @5 O% l+ v$ o
! U# ?, y# c' L3 C* P if (dv>0)
, C% N8 _* C4 Z* f8 B3 [/ _& {* F. Q
" }) P( v& w1 P5 w4 W {err=-0.1;}& g0 ]* |' R! @( X
0 a& `! ]1 a7 i* f( u/ T X: V else
4 T- H' D: H q/ t% R# {2 G * F8 {" D7 V# V' m3 @) A
{err=0.1;}" l: D9 n$ \7 Y9 s8 p) h/ T
) N' F, j1 i7 I. a9 B5 b3 o
}
4 P( n$ e h" J* E. q& }% v( x% s
+ g' R' I: g! K/ A; J0 n }
% p4 U8 f+ W5 p% h1 d# q4 w . l( o4 |) c; G; Q- C
我用的是最简单的PO法,然后设每次Vref=Vin+0.1或者-0.1,所以上面err就是一次计算完了之后的误差值,送入下面的增量型PI控制% w+ d% a$ U" K. K
dVout=20*(err-errf)+err;
4 E' Y9 J" m0 S* J8 L
. W% Z! F! B# X- }) i if (fabs(dVout)>4.1)
) l8 {, d" \' d1 G/ j5 p 7 b4 K; p% D9 E! A$ U
{
9 S* A- N/ H# o* l& ] 1 t2 V0 l, E) J1 T. Q( z
if ((dVout)>0)# a3 d4 t' l$ Z4 @' H
9 Q+ p7 O" i8 |' t! j9 t) C9 L {dVout=4.1;}0 } n" [" R ]: F. K7 b4 K# h
! ^: ]6 m1 {4 j5 ^+ Y' Y
else; Y; ~( z0 |' M6 y" u
" P+ @3 V' Q: z5 } {dVout=-4.1;}6 ]: K8 a a% n3 x& [
8 V# Z) _, x' y! M
}$ m _3 J$ ~: }2 o
$ @% R) @$ g! s4 @5 y G1 ?+ G, C Vtemp=Vout+dVout; //PI0 H0 H* }; d; K+ C" I4 W d
s" W) r5 W5 p% y, T0 f: A
" _' S: u" f+ [$ p
. W9 D" h, |# f4 X r% S6 {; O1 i if (Vtemp>10)
# m; |( G3 a. [. Z6 e
0 `' a9 G; I a8 o' c, D: d( L {Vtemp=10;}
( l5 @- H9 ]' A# V ( M, J6 H: p$ A" C1 p& b2 e
else
- P. t2 I6 B% h+ M& D & R8 Y9 D. s1 ~
{0 z8 u; F f( D* F& v5 a! ]: m
3 x$ b4 _+ [ Y! s0 q
if (Vtemp<0) S: D: a4 N; n( H3 P- k
) A1 L" d& R9 i6 ^ {Vtemp=0;}' N: _( F6 y) w) E! ?
! ?4 h' h% L( B
else! o4 l' L7 |" H; o1 r
5 u* m9 s1 V1 {; O; q- P- @ {Vtemp=Vtemp;}
: V# l. o% @5 r
/ U' N) ^0 L% V, b( L) N8 L }
# Y, i" x d$ J1 s, W8 }* k9 l, l! G 2 b6 T, `7 H% o" m
Vout=Vtemp;: N& G. `5 C4 P$ ~+ v4 x1 r q
9 m1 j3 ~3 B! O& `/ Y errf=err;
9 w2 S8 j7 k- T: w |$ Y
2 B, ^# ^% g( }1 t! U9 e' \: E 用来得到比较值的变化,最后跟锯齿波比较之后输出占空比。. b0 E0 g, }4 w! r' p- X! m
以上程序跑完之后的仿真结果如下所示:
; @& O: c% `" m) H
; P( x% |1 X: w1 Z
0 N, s7 `) G6 k9 t
感觉还是蛮精确的,效率算下来大概在99.87%,但是我总觉得哪儿有问题,不是太对,上面误差输出的值就只有0.1,-0.1,0三个值
0 }$ a J$ p5 j1 T$ d 不过不知道问题到底原因是什么,仿真效果还蛮好嘛,呵呵
$ J- ~9 g2 E/ R; K) [+ ]) O 大家给我指导一下呗,多谢!~7 [% m7 \/ o% U/ j$ b( B' L" b
楼主热帖