半对数y图（`semilogy`）

半对数y图（semilogy）

01_二维线图　5 张　MATLAB R2023b　100 mm 宽 @ 300 dpi

这个方法适合什么

纵轴取对数、横轴线性，把指数衰减或跨数量级的目标值压进一屏。适合：优化收敛曲线、损耗随距离衰减、可靠性随时间的下降。

电气场景

1　潮流计算收敛曲线

2　电池内阻增长曲线

3　谐波幅值衰减

4　电晕损耗电压特性

5　载波信号衰减

代码

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

if any(SEL==3)
    % ---- 3.1 潮流计算收敛曲线 ----
    it = 1:8;
    e1 = [8e-1 6e-2 2e-3 8e-6 1e-9 4e-13 4e-13 4e-13];
    e2 = [9e-1 4e-1 1.2e-1 2.5e-2 4e-3 5e-4 6e-5 7e-6];
    f  = ml_fig(PP,100,80);
    ax = axes('Parent',f,'Position',[0.175 0.190 0.770 0.755]);
    h1 = semilogy(ax,it,e1,'-o','Color',C(1,:),'LineWidth',1.1,'MarkerSize',4, ...
                  'MarkerFaceColor',C(1,:)); hold(ax,'on');
    h2 = semilogy(ax,it,e2,'-s','Color',C(2,:),'LineWidth',1.1,'MarkerSize',4, ...
                  'MarkerFaceColor',C(2,:));
    xlabel(ml_tx('迭代次数 <i>k</i>'));
    ylabel(ml_tx('功率不平衡量 <i>\epsilon</i> / pu'));
    xlim(ax,[1 8]); ylim(ax,[1e-14 1e2]);
    lg = legend(ax,[h1 h2],ml_tx({'牛顿-拉夫逊法','高斯-赛德尔法'}),'Location','northeast');
    ml_style_legend(lg,PP);
    save_png(f,FIGD,'潮流计算收敛曲线');

    % ---- 3.2 电池内阻增长 ----
    N  = 0:200:6000;
    R1 = 18 + 0.0045*N + 1.2e-6*N.^2;
    R2 = 18 + 0.0110*N + 6.0e-6*N.^2;
    f  = ml_fig(PP,100,80);
    ax = axes('Parent',f,'Position',[0.175 0.190 0.770 0.755]);
    h1 = semilogy(ax,N,R1,'Color',C(1,:),'LineWidth',1.2); hold(ax,'on');
    h2 = semilogy(ax,N,R2,'Color',C(4,:),'LineWidth',1.2);
    xlabel(ml_tx('循环次数 <i>N</i>'));
    ylabel(ml_tx('直流内阻 <i>R</i> / mΩ'));
    xlim(ax,[0 6000]); ylim(ax,[10 900]);
    lg = legend(ax,[h1 h2],ml_tx({'25 °C','45 °C'}),'Location','northwest');
    ml_style_legend(lg,PP);
    save_png(f,FIGD,'电池内阻增长曲线');

    % ---- 3.3 谐波幅值衰减 ----
    hh = 1:2:31;
    I1 = 120./hh.^1.3;
    I2 = 45./hh.^2.1;
    f  = ml_fig(PP,100,80);
    ax = axes('Parent',f,'Position',[0.175 0.190 0.770 0.755]);
    h1 = semilogy(ax,hh,I1,'-o','Color',C(2,:),'LineWidth',1.1,'MarkerSize',3.5); hold(ax,'on');
    h2 = semilogy(ax,hh,I2,'-s','Color',C(1,:),'LineWidth',1.1,'MarkerSize',3.5);
    xlabel(ml_tx('谐波次数 <i>h</i>'));
    ylabel(ml_tx('谐波电流 <i>I_h</i> / A'));
    xlim(ax,[1 31]); ylim(ax,[1e-3 3e2]);
    lg = legend(ax,[h1 h2],ml_tx({'六脉整流负载','变频调速负载'}),'Location','northeast');
    ml_style_legend(lg,PP);
    save_png(f,FIGD,'谐波幅值衰减');

    % ---- 3.4 电晕损耗-电压特性 ----
    U  = 200:10:525;
    P1 = 0.05*exp((U-350)/30);
    P2 = 0.08*exp((U-350)/26);
    f  = ml_fig(PP,100,80);
    ax = axes('Parent',f,'Position',[0.175 0.190 0.770 0.755]);
    h1 = semilogy(ax,U,P1,'Color',C(1,:),'LineWidth',1.2); hold(ax,'on');
    h2 = semilogy(ax,U,P2,'Color',C(2,:),'LineWidth',1.2);
    xlabel(ml_tx('线路电压 <i>U</i> / kV'));
    ylabel(ml_tx('电晕损耗 <i>P</i> / (kW/km)'));
    xlim(ax,[200 525]); ylim(ax,[1e-5 3e2]);
    lg = legend(ax,[h1 h2],ml_tx({'LGJ-400 四分裂','LGJ-630 双分裂'}),'Location','northwest');
    ml_style_legend(lg,PP);
    save_png(f,FIGD,'电晕损耗电压特性');

    % ---- 3.5 电力线载波信号衰减 ----
    d  = logspace(-1,1.3,40);
    U1 = 1200*10.^(-0.90*d/20);      % 0.90 dB/km
    U2 = 800 *10.^(-2.20*d/20);      % 2.20 dB/km
    f  = ml_fig(PP,100,80);
    ax = axes('Parent',f,'Position',[0.175 0.190 0.770 0.755]);
    h1 = semilogy(ax,d,U1,'Color',C(1,:),'LineWidth',1.2); hold(ax,'on');
    h2 = semilogy(ax,d,U2,'Color',C(3,:),'LineWidth',1.2);
    xlabel(ml_tx('传输距离 <i>l</i> / km'));
    ylabel(ml_tx('信号电平 <i>U</i> / mV'));
    xlim(ax,[0.1 20]); ylim(ax,[1.2 4e3]);
    lg = legend(ax,[h1 h2],ml_tx({'100 kHz 载波','500 kHz 载波'}),'Location','northeast');
    ml_style_legend(lg,PP);
    save_png(f,FIGD,'载波信号衰减');
end