AFD主动反孤岛

    这里，我先对AFD主动反孤岛进行算法仿真。

0.相关理论知识；

首先，这个电力系统，按照之前你所给出的要求，我们做如下的设置，即：

230V/50Hz，单相（single phase or 1 phase）。

    然后，仿真的具体要求为：

负载：小型电网 single phase 1kW grid-connected PV systems 

根据 IEEE 1547 standard 对反孤岛的要求

·Voltage： 240V ± 10%， 

·frequency：50Hz ± 1%

·Quality Factor, Qf:  Qf < 2.5 (Qf = 0 , 0.02,0.04, … , 2.4)

·Required Islanding detection time: tdetect<0.2s  and isolated to the load

·THD: <5%

关于AFD的相关理论，见具体的参考文献。

1.Simulink模型；

    下面首先建立一个单相光伏发电系统（single phase PV system），然后结合AFD的基本算法进行仿真。

注意：

之前你提供的那个AFD的S-function，我仔细看了，觉得你写的应该是的对的（然后补充了相应的代码注释，在最后给你的改进方案中，我也会对AFD的算法进行改进，改为自适应的AFD，在这里，先为传统的AFD），但是simulink模型里面没有内容，所以我基于你的这个S-funtion建立一个仿真模型。

    对S函数和simulink做了修改之后，最后得到的结构如下所示：

2.仿真波形的数据保存；

    Simulink最后的输出波形，其显示的时候没法进行label的修改，不太方便，这里我们将simulink仿真的结果保存到workspace中，然后有MATLAB调用，显示波形，保存的结果如下所示：

3. MATLAB调用数据进行波形的完整显示；

    运行文件夹中的run_figure.m就能见simulink仿真的结果进行显示，具体的仿真结果如下所示：

因为采用电压和频率边界检测，即UVP,OVP，UFP,OFP这个算法的时候，其对于边界具有一定的模糊性，我们举个例子，比如IEEE1547规定的49hz和51hz的范围之内，当发生孤岛的时候，其频率可能也是在51hz，但是这个是理论上的值，在实际中，我们会取一个极小值作为判别依据。

比如实际频率为F，F-51的绝对值或者F-49的绝对值小于0.05。或者别的一个值，我们就认为检测到了，然后就给出断开指令。

对于电压的情况，实际也是这么处理的。

所以修改之后，最后的结构图如下所示：

2.仿真波形的数据保存；

    同taskA，具体看workspace中的数据。

3. MATLAB调用数据进行波形的完整显示；