双对数图（`loglog`）

双对数图（loglog）

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

这个方法适合什么

两轴都取对数，幂律关系的斜率直接可读。适合：阻抗频谱、分形噪声、局部放电幅值分布、电缆衰减特性。

电气场景

1　设备阻抗频谱

2　分形噪声功率谱

3　变压器励磁阻抗

4　电缆特性阻抗

5　电容器ESR频谱

代码

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

if any(SEL==4)
    % ---- 4.1 阻抗频谱 ----
    fr = logspace(1,6,400);
    Z1 = hypot(0.5,2*pi*fr*1.2e-3);
    Z2 = hypot(12 ,2*pi*fr*15e-3);
    Z3 = hypot(2  ,2*pi*fr*45e-3);
    f  = ml_fig(PP,100,80);
    ax = axes('Parent',f,'Position',[0.175 0.190 0.770 0.755]);
    h1 = loglog(ax,fr,Z1,'Color',C(1,:),'LineWidth',1.1); hold(ax,'on');
    h2 = loglog(ax,fr,Z2,'Color',C(2,:),'LineWidth',1.1);
    h3 = loglog(ax,fr,Z3,'Color',C(3,:),'LineWidth',1.1);
    xlabel(ml_tx('频率 <i>f</i> / Hz'));
    ylabel(ml_tx('阻抗模 |<i>Z</i>| / Ω'));
    xlim(ax,[10 1e6]); ylim(ax,[0.5 2e6]);
    lg = legend(ax,[h1 h2 h3],ml_tx({'配电线路','配电变压器','串联电抗器'}),'Location','northwest');
    ml_style_legend(lg,PP);
    save_png(f,FIGD,'设备阻抗频谱');

    % ---- 4.2 分形噪声功率谱 ----
    fr = logspace(-2,3,400);
    S1 = 1e-6./fr;
    S2 = 3e-7./fr.^1.5;
    S3 = 2e-8./fr.^2;
    f  = ml_fig(PP,100,80);
    ax = axes('Parent',f,'Position',[0.175 0.190 0.770 0.755]);
    h1 = loglog(ax,fr,S1,'Color',C(1,:),'LineWidth',1.1); hold(ax,'on');
    h2 = loglog(ax,fr,S2,'Color',C(2,:),'LineWidth',1.1);
    h3 = loglog(ax,fr,S3,'Color',C(4,:),'LineWidth',1.1);
    xlabel(ml_tx('频率 <i>f</i> / Hz'));
    ylabel(ml_tx('功率谱密度 / (A^2/Hz)'));
    xlim(ax,[0.01 1e3]); ylim(ax,[5e-15 1e-2]);
    lg = legend(ax,[h1 h2 h3],ml_tx({'1/f 闪烁噪声','1/f^{1.5} 分形噪声','1/f^2 随机游走'}), ...
                'Location','southwest');
    ml_style_legend(lg,PP);
    save_png(f,FIGD,'分形噪声功率谱');

    % ---- 4.3 变压器励磁阻抗 ----
    fr = logspace(2,4,300);
    Z1 = 1.1e3*(fr/50).^0.85;
    Z2 = 2.4e3*(fr/50).^0.80;
    Z3 = 5.0e3*(fr/50).^0.75;
    f  = ml_fig(PP,100,80);
    ax = axes('Parent',f,'Position',[0.175 0.190 0.770 0.755]);
    h1 = loglog(ax,fr,Z1,'Color',C(1,:),'LineWidth',1.1); hold(ax,'on');
    h2 = loglog(ax,fr,Z2,'Color',C(2,:),'LineWidth',1.1);
    h3 = loglog(ax,fr,Z3,'Color',C(3,:),'LineWidth',1.1);
    xlabel(ml_tx('频率 <i>f</i> / Hz'));
    ylabel(ml_tx('励磁阻抗 <i>Z_m</i> / Ω'));
    xlim(ax,[100 1e4]); ylim(ax,[1e3 8e6]);
    lg = legend(ax,[h1 h2 h3],ml_tx({'630 kVA','1250 kVA','2500 kVA'}),'Location','northwest');
    ml_style_legend(lg,PP);
    save_png(f,FIGD,'变压器励磁阻抗');

    % ---- 4.4 电缆特性阻抗 ----
    fr = logspace(4,8,400);
    Zc = 55*(fr/1e4).^-0.12 + 8;
    f  = ml_fig(PP,100,80);
    ax = axes('Parent',f,'Position',[0.175 0.190 0.770 0.755]);
    h1 = loglog(ax,fr,Zc,'Color',C(1,:),'LineWidth',1.2);
    xlabel(ml_tx('频率 <i>f</i> / Hz'));
    ylabel(ml_tx('特性阻抗 <i>Z_c</i> / Ω'));
    xlim(ax,[1e4 1e8]); ylim(ax,[8 130]);
    lg = legend(ax,h1,ml_tx({'YJV-8.7/15'}),'Location','southwest');
    ml_style_legend(lg,PP);
    save_png(f,FIGD,'电缆特性阻抗');

    % ---- 4.5 电容器等效串联电阻 ----
    fr = logspace(1,6,400);
    E1 = 45*(fr/100).^-0.35 + 3;
    E2 = 30*(fr/100).^-0.30 + 2;
    E3 = 12*(fr/100).^-0.22 + 1;
    f  = ml_fig(PP,100,80);
    ax = axes('Parent',f,'Position',[0.175 0.190 0.770 0.755]);
    h1 = loglog(ax,fr,E1,'Color',C(1,:),'LineWidth',1.1); hold(ax,'on');
    h2 = loglog(ax,fr,E2,'Color',C(2,:),'LineWidth',1.1);
    h3 = loglog(ax,fr,E3,'Color',C(3,:),'LineWidth',1.1);
    xlabel(ml_tx('频率 <i>f</i> / Hz'));
    ylabel(ml_tx('等效串联电阻 / mΩ'));
    xlim(ax,[10 1e6]); ylim(ax,[0.8 200]);
    lg = legend(ax,[h1 h2 h3],ml_tx({'铝电解','金属化薄膜','陶瓷'}),'Location','southwest');
    ml_style_legend(lg,PP);
    save_png(f,FIGD,'电容器ESR频谱');
end