这人写错了,分母应该是减号。w02-w^2 Real and imaginary parts of the permittivity for (a) STO under different temperatures
%% figure; w = linspace(0*pi,20*pi,5000); %%for T = 250:50:400 T=400; temp = 1; w_0 = sqrt(31.2*(T-42.5)); gamma = -3.3+0.094*T; epsilon_inf =9.6; f = 2.3e6; epsilon_w =epsilon_inf +f./(w_0^2-w.^2-(1i*w*gamma)); subplot(2,1,temp); temp = temp + 1; plot(w./(2*pi),real(epsilon_w),'-'); hold on subplot(2,1,temp); plot(w./(2*pi),imag(epsilon_w),'-'); hold on a=real(epsilon_w); b=imag(epsilon_w);