切片流线图（`streamslice`）

切片流线图（streamslice）

05_场与矢量　5 张　MATLAB R2023b　100 mm 宽 @ 300 dpi

这个方法适合什么

在指定切片上画带箭头的流线。适合：二维流场、油流分布、气流组织，能同时看出流向与疏密。

电气场景

1　变压器油流切片流线

2　变电站风速切片流线

3　电缆沟自然对流切片流线

4　开关柜气流切片流线

5　电除尘器气流切片流线

代码

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

if any(SEL==6)
    % ---- 6.1 变压器油流切片流线 ----
    try
        [X,Y,U,V] = gen_oil2(161,121);

        f  = ml_fig(PP,100,90);
        ax = ax2(f);
        h = streamslice(ax,X,Y,U,V,2);
        set(h,'Color',C(1,:),'LineWidth',0.55);
        xlabel(ml_tx('油道长度 <i>x</i> / m'));
        ylabel(ml_tx('油道高度 <i>y</i> / m'));
        xlim(ax,[0 4]); ylim(ax,[0 2.4]);
        save_png(f,FIGD,'变压器油流切片流线');
    catch ME
        fprintf('[FAIL] 6-1 %s\n',ME.message); close all
    end

    % ---- 6.2 变电站风速切片流线 ----
    try
        [X,Y,U,V] = gen_wind2(161,121);

        f  = ml_fig(PP,100,90);
        ax = ax2(f);
        h = streamslice(ax,X,Y,U,V,2);
        set(h,'Color',C(4,:),'LineWidth',0.55);
        xlabel(ml_tx('站区长度 <i>x</i> / m'));
        ylabel(ml_tx('站区宽度 <i>y</i> / m'));
        xlim(ax,[0 40]); ylim(ax,[0 30]);
        save_png(f,FIGD,'变电站风速切片流线');
    catch ME
        fprintf('[FAIL] 6-2 %s\n',ME.message); close all
    end

    % ---- 6.3 电缆沟自然对流切片流线 ----
    try
        [X,Y,U,V] = gen_duct2(121,101);

        f  = ml_fig(PP,100,90);
        ax = ax2(f);
        h = streamslice(ax,X,Y,U,V,2.2);
        set(h,'Color',C(3,:),'LineWidth',0.55);
        xlabel(ml_tx('电缆沟宽度 <i>x</i> / m'));
        ylabel(ml_tx('沟深 <i>d</i> / m'));
        xlim(ax,[0 1.2]); ylim(ax,[0 1.0]);
        save_png(f,FIGD,'电缆沟自然对流切片流线');
    catch ME
        fprintf('[FAIL] 6-3 %s\n',ME.message); close all
    end

    % ---- 6.4 开关柜气流切片流线 ----
    try
        [X,Y,U,V] = gen_cab2(101,141);

        f  = ml_fig(PP,100,90);
        ax = ax2(f);
        h = streamslice(ax,X,Y,U,V,2);
        set(h,'Color',C(2,:),'LineWidth',0.55);
        xlabel(ml_tx('柜宽 <i>x</i> / m'));
        ylabel(ml_tx('柜高 <i>z</i> / m'));
        xlim(ax,[0 1.2]); ylim(ax,[0 1.8]);
        save_png(f,FIGD,'开关柜气流切片流线');
    catch ME
        fprintf('[FAIL] 6-4 %s\n',ME.message); close all
    end

    % ---- 6.5 电除尘器气流切片流线 ----
    try
        [X,Y,U,V] = gen_esp2(161,101);

        f  = ml_fig(PP,100,90);
        ax = ax2(f);
        h = streamslice(ax,X,Y,U,V,2);
        set(h,'Color',C(5,:),'LineWidth',0.55);
        xlabel(ml_tx('电场长度 <i>x</i> / m'));
        ylabel(ml_tx('电场高度 <i>y</i> / m'));
        xlim(ax,[0 6]); ylim(ax,[0 2.4]);
        save_png(f,FIGD,'电除尘器气流切片流线');
    catch ME
        fprintf('[FAIL] 6-5 %s\n',ME.message); close all
    end
end