花椒一点都不好
关注数: 0 粉丝数: 1 发帖数: 4 关注贴吧数: 2
图一是正图 图二是错在什么地方,怎么数据差那么大。有没有大佬求解 % ============= 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; % 显示网格(与论文图一致)
1 下一页