按照论文给的公式，matlab复现，结果跟论文不一样。求大佬看看什么问题。

matlab代码
clear all;
close all;
wavenumber=0:1:50;
T=293;
omega0=sqrt(31.2*(T-42.5));
gamma=-3.3+0.094*T;
f=2.3e6;
epsilon_inf=9.6;
epsilon_omega=epsilon_inf+f./(omega0^2-wavenumber.^2-1i*wavenumber*gamma);
figure;plot(wavenumber,real(epsilon_omega));

 返回小木虫查看更多