图一是正图
图二是错在什么地方，怎么数据差那么大。有没有大佬求解
% ============= 1. 论文指定参数（与原文完全一致） =============
c33 = 208e9; % 弹性刚度 (Pa)
e33 = 0.18; % 压电常数 (C/m²)
eps33 = 8.5 * 8.854e-12; % 介电常数 (F/m)
p3 = 9.4e-6; % 热电常数 (C/(m·K))
lambda33 = 3.017e-6; % 热膨胀系数 (1/K)
n0 = 1e21; % 初始电子浓度 (m⁻³)
q = 1.602e-19; % 元电荷 (C)
kB = 1.380649e-23; % 玻尔兹曼常数 (J/K)
Theta0 = 300; % 参考温度 (K)
L = 0.6e-6; % 纤维半长 (m，对应x3范围±0.6μm)
theta = [lbk]0.05, 0.10, 0.20[rbk]; % 温度变化 (K，论文中三组值)
% ============= 2. 计算论文公式中的系数 =============
eps33_hat = eps33 - (e33^2)/c33; % 有效介电常数
kappa = sqrt( (n0 * q^2) / (kB * Theta0 * eps33_hat) ); % 特征系数κ
G0 = (e33 * lambda33 + p3 * c33) / (eps33 * c33 + e33^2); % 电场系数G0
% ============= 3. 生成x3坐标并计算E3（严格使用论文公式） =============
x3 = linspace(-L, L, 500); % 加密坐标点（保证曲线平滑）
E3 = -G0 * theta' * (cosh(kappa * x3) ./ cosh(kappa * L)); % 论文E3公式
E3 = abs(E3) / 1e4; % 取绝对值（匹配图中正值）+ 缩放至×10⁴ V/m
% ============= 4. 绘图（1:1匹配论文图2第二张图样式） =============
figure('Position', [lbk]100, 100, 700, 500[rbk]);
plot(x3*1e6, E3(1,:), 'b--', 'LineWidth', 1.2); hold on; % θ=0.05K（蓝虚线）
plot(x3*1e6, E3(2,:), 'm--', 'LineWidth', 1.2); % θ=0.10K（品红虚线）
plot(x3*1e6, E3(3,:), 'g-', 'LineWidth', 1.5); % θ=0.20K（绿实线）
% 坐标轴与图例（完全对应论文图）
xlabel('$x_3$ ($\mu$m)', 'FontSize', 12, 'Interpreter', 'Latex');
ylabel('$E_3$ (V/m) $\times 10^4$', 'FontSize', 12, 'Interpreter', 'Latex');
legend({'$\theta=0.05$K', '$\theta=0.10$K', '$\theta=0.20$K'}, ...
'Location', 'northwest', 'FontSize', 10, 'Interpreter', 'Latex');
xlim([lbk]-0.6, 0.6[rbk]); % x3范围匹配论文图
ylim([lbk]0, 1.5[rbk]); % E3范围匹配论文图
grid on; % 显示网格（与论文图一致）