手把手教你学simulink（58.3）--并网逆变器场景示例：基于Simulink的孤岛检测的并网逆变器控制

 function setup_IslandDetectionInverter()

    
     % 打开Simulink模型（如果模型已存在，则加载；否则新建）
    
     if exist('IslandDetectionInverter', 'file') == 2
    
     open_system('IslandDetectionInverter');
    
     else
    
     new_system('IslandDetectionInverter');
    
     end
    
  
    
     % 添加必要的模块
    
     add_block('simscape/Electrical/Specialized Power Systems/Power Electronics/Three-Phase VSI', 'IslandDetectionInverter/Inverter');
    
     add_block('simscape/Electrical/Specialized Power Systems/Additional Components/PLL', 'IslandDetectionInverter/PLL');
    
     add_block('simulink/User-Defined Functions/MATLAB Function', 'IslandDetectionInverter/IslandDetectionController');
    
     add_block('simulink/User-Defined Functions/MATLAB Function', 'IslandDetectionInverter/PowerController');
    
     add_block('simulink/Math Operations/Sum', 'IslandDetectionInverter/CurrentError');
    
     add_block('simulink/Sinks/Scope', 'IslandDetectionInverter/Scope');
    
     add_block('simulink/Sinks/To Workspace', 'IslandDetectionInverter/ToWorkspace');
    
  
    
     % 连接模块
    
     connect_system('IslandDetectionInverter', 'Inverter/AC Voltage', 'PLL/Input Voltage');
    
     connect_system('IslandDetectionInverter', 'PLL/Phase', 'IslandDetectionController/phase');
    
     connect_system('IslandDetectionInverter', 'PLL/Frequency', 'IslandDetectionController/frequency');
    
     connect_system('IslandDetectionInverter', 'IslandDetectionController/island_detection_signal', 'Inverter/Island Detection Signal');
    
     connect_system('IslandDetectionInverter', 'PowerController/power_signal', 'CurrentError/1');
    
     connect_system('IslandDetectionInverter', 'CurrentError/2', 'Inverter/Reference Current');
    
     connect_system('IslandDetectionInverter', 'Scope/+', 'Inverter/Output Current');
    
     connect_system('IslandDetectionInverter', 'ToWorkspace/+', 'Inverter/Output Current');
    
  
    
     % 设置仿真时长
    
     set_param('IslandDetectionInverter', 'StopTime', '60');
    
  
    
     % 运行仿真
    
     sim('IslandDetectionInverter');
    
  
    
     % 获取仿真结果
    
     out = simout;
    
 end

 % 添加直流电源模型

    
 add_block('simscape/Electrical/Specialized Power Systems/Sources/DC Voltage Source', 'IslandDetectionInverter/DCSource');
    
  
    
 % 配置DCSource参数
    
 set_param('IslandDetectionInverter/DCSource', ...
    
     'Amplitude', '400', ... % 直流电压 (V)
    
     'InternalResistance', '0.1' ... % 内阻 (Ω)
    
 );

 % 配置Inverter参数

    
 set_param('IslandDetectionInverter/Inverter', ...
    
     'SwitchingFrequency', '10e3', ... % 开关频率 (Hz)
    
     'ModulationIndex', '0.9', ... % 调制指数
    
     'FilterInductance', '0.1e-3', ... % 滤波电感 (H)
    
     'FilterCapacitance', '10e-6', ... % 滤波电容 (F)
    
     'IslandDetectionSignal', 'island_detection_signal' ... % 孤岛检测信号输入
    
 );

 % 配置PLL参数

    
 set_param('IslandDetectionInverter/PLL', ...
    
     'Kp', '10', ... % 比例增益
    
     'Ki', '100', ... % 积分增益
    
     'Kf', '1', ... % 频率增益
    
     'MaxFrequency', '60', ... % 最大频率 (Hz)
    
     'MinFrequency', '50' ... % 最小频率 (Hz)
    
 );

 % 添加孤岛检测控制器模块

    
 add_block('simulink/User-Defined Functions/MATLAB Function', 'IslandDetectionInverter/IslandDetectionController');
    
  
    
 % 编写孤岛检测控制器代码
    
 set_param('IslandDetectionInverter/IslandDetectionController', 'Function', 'island_detection_controller');
    
  
    
 function island_detection_signal = island_detection_controller(phase, frequency, grid_voltage, inverter_current)
    
     % 定义孤岛检测参数
    
     voltage_threshold = 230;  % 电压阈值 (V)
    
     frequency_threshold = 0.5;  % 频率阈值 (Hz)
    
     phase_jump_threshold = 0.1;  % 相位跳变阈值 (rad)
    
     
    
     % 被动式孤岛检测
    
     passive_detection = false;
    
     if abs(grid_voltage - 230) > voltage_threshold || abs(frequency - 60) > frequency_threshold
    
     passive_detection = true;
    
     end
    
     
    
     % 主动式孤岛检测
    
     active_detection = false;
    
     if abs(diff(phase)) > phase_jump_threshold
    
     active_detection = true;
    
     end
    
     
    
     % 综合孤岛检测信号
    
     island_detection_signal = passive_detection || active_detection;
    
     
    
     % 如果检测到孤岛效应，停止逆变器运行
    
     if island_detection_signal
    
     inverter_current = 0;
    
     end
    
 end

 % 添加功率控制器模块

    
 add_block('simulink/User-Defined Functions/MATLAB Function', 'IslandDetectionInverter/PowerController');
    
  
    
 % 编写功率控制器代码
    
 set_param('IslandDetectionInverter/PowerController', 'Function', 'power_controller');
    
  
    
 function power_signal = power_controller(load_power, grid_voltage, grid_frequency)
    
     % 定义功率控制参数
    
     P_ref = 5000;  % 有功功率参考值 (W)
    
     Q_ref = 0;  % 无功功率参考值 (Var)
    
     
    
     % 根据负载需求和电网情况生成功率控制信号
    
     power_signal = [P_ref, Q_ref];
    
 end

 % 设置限幅器

    
 add_block('simulink/Math Operations/Saturation', 'IslandDetectionInverter/CurrentSaturation');
    
 set_param('IslandDetectionInverter/CurrentSaturation', ...
    
     'LowerLimit', '-100', ... % 下限 (-100A)
    
     'UpperLimit', '100' ... % 上限 (100A)
    
 );
    
  
    
 % 连接连限幅器
    
 connect_system('IslandDetectionInverter', 'PowerController/power_signal', 'CurrentSaturation/in');
    
 connect_system('IslandDetectionInverter', 'CurrentSaturation/out', 'CurrentError/1');

 % 添加过压保护逻辑

    
 add_block('simulink/Logic and Bit Operations/Relational Operator', 'IslandDetectionInverter/OvervoltageProtection');
    
 set_param('IslandDetectionInverter/OvervoltageProtection', ...
    
     'Operator', '>', ... % 比较操作符
    
     'Inputs', '2', ... % 输入数量
    
     'ZeroCross', 'on' ... % 启用零交叉检测
    
 );
    
  
    
 % 设置过压阈值
    
 add_block('simulink/Sources/Constant', 'IslandDetectionInverter/OvervoltageThreshold');
    
 set_param('IslandDetectionInverter/OvervoltageThreshold', ...
    
     'Value', '260' ... % 过压阈值 (V)
    
 );
    
  
    
 % 连接过压保护逻辑
    
 connect_system('IslandDetectionInverter', 'Inverter/AC Voltage', 'OvervoltageProtection/1');
    
 connect_system('IslandDetectionInverter', 'OvervoltageThreshold/+', 'OvervoltageProtection/2');
    
 connect_system('IslandDetectionInverter', 'OvervoltageProtection/+', 'Scope/2');
    
  
    
 % 添加过流保护逻辑
    
 add_block('simulink/Logic and Bit Operations/Relational Operator', 'IslandDetectionInverter/OvercurrentProtection');
    
 set_param('IslandDetectionInverter/OvercurrentProtection', ...
    
     'Operator', '>', ... % 比较操作符
    
     'Inputs', '2', ... % 输入数量
    
     'ZeroCross', 'on' ... % 启用零交叉检测
    
 );
    
  
    
 % 设置过流阈值
    
 add_block('simulink/Sources/Constant', 'IslandDetectionInverter/OvercurrentThreshold');
    
 set_param('IslandDetectionInverter/OvercurrentThreshold', ...
    
     'Value', '100' ... % 过流阈值 (A)
    
 );
    
  
    
 % 连接过流保护逻辑
    
 connect_system('IslandDetectionInverter', 'Inverter/Output Current', 'OvercurrentProtection/1');
    
 connect_system('IslandDetectionInverter', 'OvercurrentThreshold/+', 'OvercurrentProtection/2');
    
 connect_system('IslandDetectionInverter', 'OvercurrentProtection/+', 'Scope/3');

 % 提取时间向量和输出

    
 time = out.time;
    
 grid_voltage = out.signals.values(:, 2); % 假设电网电压是第二个信号
    
 inverter_current = out.signals.values(:, 3); % 假设逆变器输出电流是第三个信号
    
 island_detection_signal = out.signals.values(:, 4); % 假设孤岛检测信号是第四个信号
    
  
    
 % 绘制电网电压随时间变化的图形
    
 figure;
    
 subplot(3,1,1);
    
 plot(time, grid_voltage);
    
 xlabel('Time (s)');
    
 ylabel('Grid Voltage (V)');
    
 title('Grid Voltage Response');
    
 grid on;
    
  
    
 % 绘制逆变器输出电流随时间变化的图形
    
 subplot(3,1,2);
    
 plot(time, inverter_current);
    
 xlabel('Time (s)');
    
 ylabel('Inverter Output Current (A)');
    
 title('Inverter Output Current Response');
    
 grid on;
    
  
    
 % 绘制孤岛检测信号随时间变化的图形
    
 subplot(3,1,3);
    
 plot(time, island_detection_signal);
    
 xlabel('Time (s)');
    
 ylabel('Island Detection Signal');
    
 title('Island Detection Signal Response');
    
 grid on;

 % 选择目标硬件平台

    
 targetPlatform = 'Arduino';
    
  
    
 % 配置硬件连接
    
 if strcmp(targetPlatform, 'Arduino')
    
     arduinoObj = arduino('COM3', 'Uno', 'Libraries', {'Servo'});
    
 else
    
     error('Unsupported platform');
    
 end
    
  
    
 % 下载生成的代码到目标硬件
    
 deploy(arduinoObj, 'IslandDetectionInverter');
    
  
    
 % 开始运行
    
 start(arduinoObj);
    
  
    
 % 监控实时数据并记录日志
    
 logFile = fopen('real_time_data.log', 'w');
    
 if logFile == -1
    
     error('Failed to open log file');
    
 end
    
  
    
 while true
    
     % 获取实时数据
    
     data = getSignal(arduinoObj, 'Scope');
    
     
    
     % 写入日志文件
    
     fprintf(logFile, '%.2f, %.2f, %.2f, %.2f\n', ...
    
     data.Time, data.Grid_Voltage, data.Inverter_Current, data.Island_Detection_Signal);
    
     
    
     % 暂停一段时间以减少数据量
    
     pause(1);
    
 end
    
  
    
 % 关闭日志文件
    
 fclose(logFile);

 % 配置代码生成参数

    
 cfg = coder.config('hdl');
    
 cfg.HDLSubsystem = 'IslandDetectionInverter/Inverter';
    
 cfg.GenerateExampleDesign = true;
    
  
    
 % 生成逆变器代码
    
 codegen -config cfg IslandDetectionInverter
    
  
    
 % 配置孤岛检测控制器代码生成参数
    
 cfg.HDLSubsystem = 'IslandDetectionInverter/IslandDetectionController';
    
  
    
 % 生成孤岛检测控制器代码
    
 codegen -config cfg IslandDetectionInverter
    
  
    
 % 配置PLL代码生成参数
    
 cfg.HDLSubsystem = 'IslandDetectionInverter/PLL';
    
  
    
 % 生成PLL代码
    
 codegen -config cfg IslandDetectionInverter
    
  
    
 % 配置功率控制器代码生成参数
    
 cfg.HDLSubsystem = 'IslandDetectionInverter/PowerController';
    
  
    
 % 生成功率控制器代码
    
 codegen -config cfg IslandDetectionInverter
    
  
    
 % 配置过压保护逻辑代码生成参数
    
 cfg.HDLSubsystem = 'IslandDetectionInverter/OvervoltageProtection';
    
  
    
 % 生成过压保护逻辑代码
    
 codegen -config cfg IslandDetectionInverter
    
  
    
 % 配置过流保护逻辑代码生成参数
    
 cfg.HDLSubsystem = 'IslandDetectionInverter/OvercurrentProtection';
    
  
    
 % 生成过流保护逻辑代码
    
 codegen -config cfg IslandDetectionInverter

 % 提取时间向量和输出

    
 time = out.time;
    
 grid_voltage = out.signals.values(:, 2); % 假设电网电压是第二个信号
    
 inverter_current = out.signals.values(:, 3); % 假设逆变器输出电流是第三个信号
    
 island_detection_signal = out.signals.values(:, 4); % 假设孤岛检测信号是第四个信号
    
  
    
 % 绘制电网电压随时间变化的图形
    
 figure;
    
 subplot(3,1,1);
    
 plot(time, grid_voltage);
    
 xlabel('Time (s)');
    
 ylabel('Grid Voltage (V)');
    
 title('Grid Voltage Response');
    
 grid on;
    
  
    
 % 绘制逆变器输出电流随时间变化的图形
    
 subplot(3,1,2);
    
 plot(time, inverter_current);
    
 xlabel('Time (s)');
    
 ylabel('Inverter Output Current (A)');
    
 title('Inverter Output Current Response');
    
 grid on;
    
  
    
 % 绘制孤岛检测信号随时间变化的图形
    
 subplot(3,1,3);
    
 plot(time, island_detection_signal);
    
 xlabel('Time (s)');
    
 ylabel('Island Detection Signal');
    
 title('Island Detection Signal Response');
    
 grid on;

 % 选择目标硬件平台

    
 targetPlatform = 'Arduino';
    
  
    
 % 配置硬件连接
    
 if strcmp(targetPlatform, 'Arduino')
    
     arduinoObj = arduino('COM3', 'Uno', 'Libraries', {'Servo'});
    
 else
    
     error('Unsupported platform');
    
 end
    
  
    
 % 下载生成的代码到目标硬件
    
 deploy(arduinoObj, 'IslandDetectionInverter');
    
  
    
 % 开始运行
    
 start(arduinoObj);
    
  
    
 % 监控实时数据并记录日志
    
 logFile = fopen('real_time_data.log', 'w');
    
 if logFile == -1
    
     error('Failed to open log file');
    
 end
    
  
    
 while true
    
     % 获取实时数据
    
     data = getSignal(arduinoObj, 'Scope');
    
     
    
     % 写入日志文件
    
     fprintf(logFile, '%.2f, %.2f, %.2f, %.2f\n', ...
    
     data.Time, data.Grid_Voltage, data.Inverter_Current, data.Island_Detection_Signal);
    
     
    
     % 暂停一段时间以减少数据量
    
     pause(1);
    
 end
    
  
    
 % 关闭日志文件
    
 fclose(logFile);