SVPWM算法理解、Simulink仿真及程序编写

SVPWM算法理解、Simulink仿真及程序编写

1. 空间矢量理解

三相逆变器中,三相电压若分开描述则较繁琐,故将三相电压量用一个合成量表示(该合成量能保持三相信息的完整性),则三相问题简化为单相(复数)问题。相似于坐标变换。

在这里插入图片描述
三相电压 a = k sin x ; b = k sin (x-120) ; c = k sin (x+120) 。
空间矢量电压 V = Re+ j Im = 1.5k e^j {x-90} ; 其轨迹为圆。
对于任一三相电压状态,只有唯一空间矢量电压点 与之对应;对于任一空间矢量电压点,亦只有唯一三相电压状态 与之对应。两者相互转换。

举一个常见的例子,在两电平三相逆变电路中,如下。
在这里插入图片描述
开关管导通状态有 000~111 ,共8种。其表示意义为 “abc”3桥臂 ,“1”上半桥臂导通。
由开关管导通状态 得到 三相电压状态 如表格中所示,将 三相电压状态 转换为 空间电压矢量 表示如下所示。
在空间电压矢量 中很清晰得到开关管导通顺序:a(110001…); b ( 011100… ) ; c ( 000111 )
在这里插入图片描述

联系第一段的理论基础,6个空间矢量组成的 空间电压矢量 是个不规整的圆,故其得到的三相电流 亦十分失真。然而能够 通过使相邻两个矢量(和适当的零矢量)短时间内 交替作用,似于 合成新的空间矢量,使得 圆规整些、三相电流趋于正弦。这种合成“新”矢量 的方法便是 SVPWM。

2. SVPWM算法

主要包括3个部分:判断参考电压矢量Uref所在扇区;计算相邻空间矢量的作用时间;根据作用时间合成三相PWM信号。

参考电压矢量Uref 所在扇区判断

空间扇区定义如下所示;图中Ud为逆变器主电路的直流母线电压,矢量幅值为2Ud/3。根据恒幅值clack变换推得。
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

计算相邻两电压空间矢量的作用时间

在这里插入图片描述
在这里插入图片描述

根据作用时间合成三相PWM信号

在这里插入图片描述
在这里插入图片描述

3. Simulink仿真

在这里插入图片描述
采用恒幅值变换将 三相变换为两相 作为svpwm1程序的输入;余下2输入分别为 采样时间和直流侧母线电压;
在这里插入图片描述
在这里插入图片描述
svpwm1输出的调制信号
在这里插入图片描述
经三角载波调制后输出的PWM波
在这里插入图片描述
B点电势
在这里插入图片描述
B相负载上的电压、电流波形。
在这里插入图片描述

4. svpwm程序

#define S_FUNCTION_NAME  svpwm1
#define S_FUNCTION_LEVEL 2

#include "simstruc.h"
#include <math.h>

static void mdlInitializeSizes(SimStruct *S)
{
    /* 函数用来设置输入、输出和参数 */
    ssSetNumSFcnParams(S, 0);  /* 设置参数个数,这里为0 */
    if (ssGetNumSFcnParams(S) != ssGetSFcnParamsCount(S)) 
    {
        return;
    }
    ssSetNumContStates(S, 0);  /* 设置连续状态的个数,缺省为0 */
    ssSetNumDiscStates(S, 0);  /* 设置离散状态的个数,缺省为0 */

    if (!ssSetNumInputPorts(S, 1)) return;  /* 设置输入变量的个数,这里为1 */
    ssSetInputPortWidth(S, 0, 4);  /* 设置输入变量0的维数为4 */
    ssSetInputPortRequiredContiguous(S, 0, true); /* 设置input0的访问方式,ture就是邻近访问 */
    ssSetInputPortDirectFeedThrough(S, 0, 1);/* 设置输入端口的信号是否mdloutputs函数中使用,这儿设置为ture */

    if (!ssSetNumOutputPorts(S, 1)) return;   /* 设置输出变量的个数,这里为1 */
    ssSetOutputPortWidth(S, 0, 3);  /* 设置输出变量0的维数为3 */

    ssSetNumSampleTimes(S, 2e-4);  /* 设置采样时间,这里为2e-4S */
    ssSetNumRWork(S, 0);  /* 不管 */
    ssSetNumIWork(S, 0);  /* 不管 */
    ssSetNumPWork(S, 0);  /* 不管 */  
    ssSetNumModes(S, 0);  /* 不管 */
    ssSetNumNonsampledZCs(S, 0);  /* 不管 */

    ssSetOptions(S, SS_OPTION_EXCEPTION_FREE_CODE);  
}

static void mdlInitializeSampleTimes(SimStruct *S)
{
    ssSetSampleTime(S, 0, 2e-4);
    ssSetOffsetTime(S, 0, 0.0);
}

#define MDL_INITIALIZE_CONDITIONS   /* Change to #undef to remove function */
#if defined(MDL_INITIALIZE_CONDITIONS)
  static void mdlInitializeConditions(SimStruct *S)
  {
  }
#endif 

#define MDL_START  /* Change to #undef to remove function */
#if defined(MDL_START) 
  static void mdlStart(SimStruct *S)
  {
  }
#endif /*  MDL_START */

/*
//svpwm算法程序
//svpwm算法程序
//svpwm算法程序
*/
static void mdlOutputs(SimStruct *S, int_T tid)
{
	const double *u = (const double*) ssGetInputPortSignal(S,0);
	double *y = ssGetOutputPortSignal(S,0);
	double Valfa_ref=u[0];
	double Vbeta_ref=u[1];
	double Tpwm=u[2];
	double Vdc=u[3];
	double Vdcinvt=Tpwm/Vdc;
	double Vref1,Vref2,Vref3;
	double X,Y,Z;
	double t1,t2,t1sat,t2sat;
	double taon,tbon,tcon;
	double CMPR1,CMPR2,CMPR3;
    double Tpwm_1=-Tpwm;
	int sector;
	int A,B,C;
	//计算Vref1,Vref2,Vref3,X,Y,Z
	Vref1=Vbeta_ref;
	Vref2=sqrt(3.0)*Valfa_ref-Vbeta_ref;
	Vref3=-sqrt(3.0)*Valfa_ref-Vbeta_ref;
	X=sqrt(3.0)*Vdcinvt*Vbeta_ref;
	Y=sqrt(3.0)*Vdcinvt*Vbeta_ref/2.0+1.5*Vdcinvt*Valfa_ref;
	Z=sqrt(3.0)*Vdcinvt*Vbeta_ref/2.0-1.5*Vdcinvt*Valfa_ref;
	//判断扇区
	if(Vref1>0)
		A=1;
	else
		A=0;
	if(Vref2>0)
		B=1;
	else
		B=0;
	if(Vref3>0)
		C=1;
	else
		C=0;
	sector=A+2*B+4*C;
	switch(sector)
	{
	case 1:
		t1=Z;
		t2=Y;
		break;
	case 2:
		t1=Y;
		t2=-X;
		break;
	case 3:
		t1=-Z;
		t2=X;
		break;
	case 4:
		t1=-X;
		t2=Z;
		break;
	case 5:
		t1=X;
		t2=-Y;
		break;
	default:
		t1=-Y;
		t2=-Z;
	}
	//判断是否进入饱合,如果进入饱合则按相位不变的原则进行处理
	if((t1+t2)>(Tpwm))
	{
		t1sat=t1*Tpwm/(t1+t2);
		t2sat=t2*Tpwm/(t1+t2);
		t1=t1sat;
		t2=t2sat;
	}
	//给taon,tbon,tcon赋值
	taon=(Tpwm-t1-t2)/4.0;
	tbon=taon+t1/2.0;
	tcon=tbon+t2/2.0;
	//给CMPR1,CMPR2,CMPR3赋值
	switch(sector)
	{
	case 1:
		CMPR1=tbon;
		CMPR2=taon;
		CMPR3=tcon;
		break;
	case 2:
		CMPR1=taon;
		CMPR2=tcon;
		CMPR3=tbon;
		break;
	case 3:
		CMPR1=taon;
		CMPR2=tbon;
		CMPR3=tcon;
		break;
	case 4:
		CMPR1=tcon;
		CMPR2=tbon;
		CMPR3=taon;
		break;
	case 5:
		CMPR1=tcon;
		CMPR2=taon;
		CMPR3=tbon;
		break;
	default:
		CMPR1=tbon;
		CMPR2=tcon;
		CMPR3=taon;
	}
  
//输出
	y[0]=CMPR1;
	y[1]=CMPR2;
	y[2]=CMPR3; 
}


#define MDL_UPDATE  /* Change to #undef to remove function */
#if defined(MDL_UPDATE)
  static void mdlUpdate(SimStruct *S, int_T tid)
  {
  }
#endif /* MDL_UPDATE */

#define MDL_DERIVATIVES  /* Change to #undef to remove function */
#if defined(MDL_DERIVATIVES)
  static void mdlDerivatives(SimStruct *S)
  {
  }
#endif /* MDL_DERIVATIVES */

static void mdlTerminate(SimStruct *S)
{
}
#ifdef  MATLAB_MEX_FILE    /* Is this file being compiled as a MEX-file? */
#include "simulink.c"      /* MEX-file interface mechanism */
#else
#include "cg_sfun.h"       /* Code generation registration function */
#endif
©️2020 CSDN 皮肤主题: 大白 设计师: CSDN官方博客 返回首页
实付0元
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、C币套餐、付费专栏及课程。

余额充值