matlab mpopt,MATLAB程序解读。急急急急急!!!!!!!!!!!

该MATLAB程序实现了一个电力系统的潮流计算。通过输入包括发电机、负荷、线路等数据,使用不同的功率流算法(牛顿法、快速解耦法、高斯赛德尔法)求解节点电压和功率注入。程序中涉及到了数据转换、矩阵构建、解算过程及结果输出,对于理解和模拟电力系统运行状态有帮助。
摘要由CSDN通过智能技术生成

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

function [scores bus ] = runpdwga(pop, baseMVA, bus,Main, M, gen, NN, branch, mpopt)

%% bus data

% bus_i type Pd Qd Gs Bs area Vm Va baseKV zone Vmax Vmin

baseMVA=1

bus = [

1 3 0 0 0 0 1 1 0 10 1 1.1 0.9;

2 1 0.0831 0.0304 0 0 1 1 0 10 1 1.1 0.9;

3 1 0 0 0 0 1 1 0 10 1 1.1 0.9;

4 1 1.7451 0.6392 0 0 1 1 0 10 1 1.1 0.9;

5 1 0 0 0 0 1 1 0 10 1 1.1 0.9;

6 1 2.9294 1.0729 0 0 1 1 0 10 1 1.1 0.9;

];

%% branch data

% fbus tbus r x b rateA rateB rateC ratio angle status

Main=[

1 2 0.00084 0.001155 0 250 250 250 0 0 1;

1 6 0.00392 0.00539 0 250 250 250 0 0 1;

2 3 0.0021 0.002888 0 250 250 250 0 0 1;

]

M = [

3 4 0.0007 0.0009625 0 150 150 150 0 0 1;

4 5 0.00042 0.0005775 0 150 150 150 0 0 1;

6 5 0.00028 0.000385 0

以下是一个简单的33节点配电网潮流分区的MATLAB程序,仅供参考: ```matlab % 33节点配电网潮流分区 % 数据输入 % 节点数据 busdata = [ 1 1.06 0 0 0 0 0 0; 2 1.045 0 0 0 0 0 0; 3 1.01 0 0 0 0 0 0; 4 1.07 0 0 0 0 0 0; 5 1.09 0 0 0 0 0 0; 6 1.05 0 0 0 0 0 0; 7 1.08 0 0 0 0 0 0; 8 1.07 0 0 0 0 0 0; 9 1.09 0 0 0 0 0 0; 10 1.03 0 0 0 0 0 0; 11 1.02 0 0 0 0 0 0; 12 1.01 0 0 0 0 0 0; 13 1.06 0 0 0 0 0 0; 14 1.09 0 0 0 0 0 0; 15 1.07 0 0 0 0 0 0; 16 1.05 0 0 0 0 0 0; 17 1.03 0 0 0 0 0 0; 18 1.02 0 0 0 0 0 0; 19 1.01 0 0 0 0 0 0; 20 1.06 0 0 0 0 0 0; 21 1.09 0 0 0 0 0 0; 22 1.07 0 0 0 0 0 0; 23 1.05 0 0 0 0 0 0; 24 1.03 0 0 0 0 0 0; 25 1.02 0 0 0 0 0 0; 26 1.01 0 0 0 0 0 0; 27 1.06 0 0 0 0 0 0; 28 1.09 0 0 0 0 0 0; 29 1.07 0 0 0 0 0 0; 30 1.05 0 0 0 0 0 0; 31 1.03 0 0 0 0 0 0; 32 1.02 0 0 0 0 0 0; 33 1.01 0 0 0 0 0 0; ]; % 支路数据 branchdata = [ 1 2 0.032 0.16 0 0 1 1 0 0 0 0; 1 3 0.01 0.08 0 0 1 1 0 0 0 0; 2 4 0.039 0.17 0 0 1 1 0 0 0 0; 2 5 0.017 0.08 0 0 1 1 0 0 0 0; 2 6 0.032 0.16 0 0 1 1 0 0 0 0; 3 7 0.0085 0.072 0 0 1 1 0 0 0 0; 3 8 0.029 0.13 0 0 1 1 0 0 0 0; 4 9 0.064 0.27 0 0 1 1 0 0 0 0; 5 9 0.045 0.20 0 0 1 1 0 0 0 0; 6 10 0.045 0.20 0 0 1 1 0 0 0 0; 6 7 0.02 0.10 0 0 1 1 0 0 0 0; 7 11 0.024 0.10 0 0 1 1 0 0 0 0; 8 11 0.048 0.22 0 0 1 1 0 0 0 0; 9 11 0.045 0.20 0 0 1 1 0 0 0 0; 9 10 0.046 0.20 0 0 1 1 0 0 0 0; 4 12 0.020 0.10 0 0 1 1 0 0 0 0; 12 13 0.023 0.10 0 0 1 1 0 0 0 0; 12 14 0.015 0.08 0 0 1 1 0 0 0 0; 12 15 0.017 0.10 0 0 1 1 0 0 0 0; 12 16 0.018 0.10 0 0 1 1 0 0 0 0; 14 15 0.017 0.10 0 0 1 1 0 0 0 0; 16 17 0.039 0.17 0 0 1 1 0 0 0 0; 15 18 0.016 0.08 0 0 1 1 0 0 0 0; 18 19 0.045 0.20 0 0 1 1 0 0 0 0; 19 20 0.060 0.26 0 0 1 1 0 0 0 0; 10 20 0.026 0.13 0 0 1 1 0 0 0 0; 10 17 0.032 0.16 0 0 1 1 0 0 0 0; 10 21 0.045 0.20 0 0 1 1 0 0 0 0; 10 22 0.048 0.23 0 0 1 1 0 0 0 0; 21 22 0.011 0.08 0 0 1 1 0 0 0 0; 15 23 0.10 0.50 0.25 0.25 1 1 0 0 0 0; 22 24 0.10 0.50 0.25 0.25 1 1 0 0 0 0; 23 24 0.05 0.40 0.25 0.25 1 1 0 0 0 0; 24 25 0.038 0.20 0 0 1 1 0 0 0 0; 25 26 0.025 0.12 0 0 1 1 0 0 0 0; 25 27 0.022 0.10 0 0 1 1 0 0 0 0; 28 27 0.022 0.10 0 0 1 1 0 0 0 0; 27 29 0.048 0.22 0 0 1 1 0 0 0 0; 27 30 0.035 0.16 0 0 1 1 0 0 0 0; 29 30 0.012 0.08 0 0 1 1 0 0 0 0; 8 28 0.063 0.28 0 0 1 1 0 0 0 0; 6 28 0.016 0.08 0 0 1 1 0 0 0 0; 4 28 0.017 0.08 0 0 1 1 0 0 0 0; ]; % 潮流分析 mpc.bus = busdata; mpc.branch = branchdata; mpc.gen = zeros(0,8); mpopt = mpoption('PF_ALG', 1, 'VERBOSE', 0); results = rundcpf(mpc, mpopt); % 分区划分 area1 = [1 2 3 6 7 10 11 12 13 14 15 16 18 19]; area2 = [4 5 8 9 17 20 21 22 23 24 25 26 27 28 29 30 31 32 33]; % 输出结果 fprintf('Area 1:\n'); disp(area1); fprintf('Area 2:\n'); disp(area2); ``` 这个程序使用了MATPOWER工具箱中的`rundcpf`函数进行潮流分析,并根据分析结果将节点划分为两个区域。请注意,这是一个非常简单的示例程序,可以根据实际情况进行修改和扩展。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值