气泡图（`bubblechart`）

气泡图（bubblechart）

02_柱与散点　5 张　MATLAB R2023b　100 mm 宽 @ 300 dpi

这个方法适合什么

点的面积编码第三个变量，一张图读三个量。适合：节点电压—负荷—短路容量、机组出力—煤耗—排放的三维比较。

电气场景

1　负载率效率容量气泡图

2　电压偏差谐波畸变负荷气泡图

3　年份线损率供电量气泡图

4　供电半径电压降负荷矩气泡图

5　功率温升运行年限气泡图

代码

脚本 p2_bar_scatter.m 中对应 bubblechart 的一段，SEL 用于选择只跑某几个方法。

if any(SEL==6)
    rng(21);
    % ---- 6.1 负载率 - 效率 - 容量 ----
    b1 = [22 35 48 62 74 86];
    b2 = [14 26 38 52 66 80];
    b3 = [18 30 44 58 70 82];
    e1 = 94.2 + 4.6*(1-exp(-(b1-10)/26)) - 0.0012*max(b1-80,0).^2;
    e2 = 95.4 + 3.4*(1-exp(-(b2-10)/32)) - 0.0009*max(b2-84,0).^2;
    e3 = 93.4 + 5.4*(1-exp(-(b3-10)/20)) - 0.0016*max(b3-76,0).^2;
    s1 = [8 12 18 24 30 36]; s2 = [20 30 42 55 68 80]; s3 = [4 7 11 15 20 26];
    f  = ml_fig(PP,100,80);
    ax = axes('Parent',f,'Position',[0.175 0.190 0.770 0.755]);
    h1 = bubblechart(ax,b1,e1,s1/4,repmat(C(1,:),6,1)); hold(ax,'on');
    h2 = bubblechart(ax,b2,e2,s2/4,repmat(C(2,:),6,1));
    h3 = bubblechart(ax,b3,e3,s3/4,repmat(C(3,:),6,1));
    set([h1 h2 h3],'MarkerFaceAlpha',0.55,'MarkerEdgeColor','none');
    xlabel(ml_tx('负载率 <i>\beta</i> / %'));
    ylabel(ml_tx('效率 <i>\epsilon</i> / %'));
    xlim(ax,[5 100]); ylim(ax,[92.5 104]);
    lg = legend(ax,[h1 h2 h3],ml_tx({'10 kV 配变','35 kV 主变','380 V 逆变器'}), ...
                'Location','northeast');
    ml_style_legend(lg,PP);
    save_png(f,FIGD,'负载率效率容量气泡图');

    % ---- 6.2 电压偏差 - 谐波畸变 - 负荷 ----
    v1 = [-4.2 -3.1 -2.4 -1.6 -0.8 -0.2];
    v2 = [-3.6 -2.7 -1.8 -1.1 -0.5  0.3];
    v3 = [-2.8 -1.9 -1.2 -0.6 -0.1  0.8];
    t1 = [6.8 5.9 5.1 4.4 3.8 3.2];
    t2 = [4.6 4.1 3.5 3.0 2.6 2.1];
    t3 = [8.4 7.6 6.8 6.0 5.4 4.8];
    q1 = [0.6 1.1 1.8 2.4 3.1 3.8]; q2 = [4 7 11 16 21 26]; q3 = [1.2 2.0 2.8 3.6 4.5 5.4];
    f  = ml_fig(PP,100,80);
    ax = axes('Parent',f,'Position',[0.175 0.190 0.770 0.755]);
    h1 = bubblechart(ax,v1,t1,q1/4,repmat(C(4,:),6,1)); hold(ax,'on');
    h2 = bubblechart(ax,v2,t2,q2/4,repmat(C(1,:),6,1));
    h3 = bubblechart(ax,v3,t3,q3/4,repmat(C(6,:),6,1));
    set([h1 h2 h3],'MarkerFaceAlpha',0.55,'MarkerEdgeColor','none');
    xlabel(ml_tx('母线电压偏差 <i>\delta U</i> / %'));
    ylabel(ml_tx('电压总谐波畸变率 <i>U</i>_{THD} / %'));
    xlim(ax,[-5.6 1.9]); ylim(ax,[0.2 11.2]);
    lg = legend(ax,[h1 h2 h3],ml_tx({'居民馈线','工业馈线','商业馈线'}), ...
                'Location','northeast');
    ml_style_legend(lg,PP);
    save_png(f,FIGD,'电压偏差谐波畸变负荷气泡图');

    % ---- 6.3 年份 - 线损率 - 供电量 ----
    yy = [2017 2018 2019 2020 2021 2022 2023 2024];
    l1 = [6.8 6.5 6.1 5.8 5.4 5.1 4.8 4.5];
    l2 = [5.2 5.0 4.7 4.4 4.1 3.8 3.5 3.2];
    l3 = [7.6 7.3 6.9 6.5 6.0 5.6 5.2 4.9];
    g1 = [86 94 102 110 118 127 136 145]; g2 = [42 46 50 55 60 65 70 76]; ...
    g3 = [16 18 20 22 24 26 29 32];
    f  = ml_fig(PP,100,80);
    ax = axes('Parent',f,'Position',[0.175 0.190 0.770 0.755]);
    h1 = bubblechart(ax,yy,l1,g1/4,repmat(C(1,:),8,1)); hold(ax,'on');
    h2 = bubblechart(ax,yy,l2,g2/4,repmat(C(3,:),8,1));
    h3 = bubblechart(ax,yy,l3,g3/4,repmat(C(2,:),8,1));
    set([h1 h2 h3],'MarkerFaceAlpha',0.55,'MarkerEdgeColor','none');
    xlabel(ml_tx('年份'));
    ylabel(ml_tx('线损率 <i>\epsilon</i> / %'));
    xlim(ax,[2016.4 2024.6]); ylim(ax,[2.6 8.4]);
    lg = legend(ax,[h1 h2 h3],ml_tx({'城网','农网','工业园'}), ...
                'Location','northeast','Orientation','horizontal');
    ml_style_legend(lg,PP);
    save_png(f,FIGD,'年份线损率供电量气泡图');

    % ---- 6.4 供电半径 - 电压降 - 负荷矩 ----
    r1 = [0.8 1.6 2.4 3.2 4.0 4.8 5.6];
    r2 = [1.2 2.0 2.9 3.8 4.7 5.6 6.5];
    r3 = [2.0 3.0 4.1 5.2 6.3 7.4 8.5];
    d1 = [0.6 1.3 2.1 3.0 4.0 5.1 6.3];
    d2 = [0.9 1.6 2.4 3.4 4.5 5.7 7.0];
    d3 = [1.4 2.2 3.2 4.4 5.8 7.4 9.2];
    m1 = [4 8 12 17 22 27 33]; m2 = [6 10 15 21 27 34 41]; m3 = [10 16 23 31 40 50 61];
    f  = ml_fig(PP,100,80);
    ax = axes('Parent',f,'Position',[0.175 0.190 0.770 0.755]);
    h1 = bubblechart(ax,r1,d1,m1/4,repmat(C(5,:),7,1)); hold(ax,'on');
    h2 = bubblechart(ax,r2,d2,m2/4,repmat(C(6,:),7,1));
    h3 = bubblechart(ax,r3,d3,m3/4,repmat(C(4,:),7,1));
    set([h1 h2 h3],'MarkerFaceAlpha',0.55,'MarkerEdgeColor','none');
    xlabel(ml_tx('供电半径 <i>r</i> / km'));
    ylabel(ml_tx('线路电压降 <i>\Delta U</i> / %'));
    xlim(ax,[0 10.8]); ylim(ax,[0 12.6]);
    lg = legend(ax,[h1 h2 h3],ml_tx({'LGJ-95','LGJ-150','YJV-240'}), ...
                'Location','northwest');
    ml_style_legend(lg,PP);
    save_png(f,FIGD,'供电半径电压降负荷矩气泡图');

    % ---- 6.5 功率 - 温升 - 运行年限 ----
    p1 = [12 26 40 54 68 82 96];
    p2 = [16 30 45 60 74 88 100];
    p3 = [8 22 36 50 64 78 92];
    w1 = [14 20 26 33 41 50 60];
    w2 = [18 24 31 38 46 55 65];
    w3 = [10 16 22 29 36 44 53];
    a1 = [3 5 7 10 13 16 20]; a2 = [1 2 4 6 9 12 15]; a3 = [8 11 14 17 20 24 28];
    f  = ml_fig(PP,100,80);
    ax = axes('Parent',f,'Position',[0.175 0.190 0.770 0.755]);
    h1 = bubblechart(ax,p1,w1,a1/4,repmat(C(1,:),7,1)); hold(ax,'on');
    h2 = bubblechart(ax,p2,w2,a2/4,repmat(C(3,:),7,1));
    h3 = bubblechart(ax,p3,w3,a3/4,repmat(C(2,:),7,1));
    set([h1 h2 h3],'MarkerFaceAlpha',0.55,'MarkerEdgeColor','none');
    xlabel(ml_tx('有功功率 <i>P</i> / kW'));
    ylabel(ml_tx('绕组温升 <i>\Delta T</i> / K'));
    xlim(ax,[-2 114]); ylim(ax,[2 82]);
    lg = legend(ax,[h1 h2 h3],ml_tx({'新投运','运行 8 年','运行 15 年'}), ...
                'Location','northwest');
    ml_style_legend(lg,PP);
    save_png(f,FIGD,'功率温升运行年限气泡图');
end