Simulación y cálculo de parámetros electromagnéticos de placas dieléctricas
Domine el cálculo de simulación del coeficiente de reflexión y el coeficiente de transmisión de la onda electromagnética armónica del plano uniforme incidente en la placa dieléctrica; compile el programa matlab para calcular los parámetros electromagnéticos de la placa dieléctrica utilizando el coeficiente de reflexión y el coeficiente de transmisión.
Utilice el software Eastwave para diseñar primero el modelo experimental correspondiente, luego diseñe las condiciones de contorno, la cuadrícula, la fuente de excitación, etc., y ejecute la distribución del campo de observación. El modelo diseñado en este experimento (el objeto utiliza un material dieléctrico, la permitividad relativa es 12,56) y la cuadrícula son:
resultados experimentales
Utilice MATLAB para programar de acuerdo con el principio para obtener los parámetros electromagnéticos correspondientes. el código se muestra a continuación:
gdata.DELTA_T=1;
gdata.UT=1.66782e-17;
gata.UL=1e-8;
gdata.STEP_MAX=1000000;
Ein=load('Electromagnetic_parameter_measurement.ewp.data\Ein.ed');
Er=load('Electromagnetic_parameter_measurement.ewp.data\Erz.ed');
Et=load('Electromagnetic_parameter_measurement.ewp.data\Etz.ed');
Er(1:390)=0;
fmax=1/gdata.DELTA_T/gdata.UT;
freq=linspace(0,1,gdata.STEP_MAX)*fmax;
FEin=fft(Ein,gdata.STEP_MAX);
FEr=fft(Er,gdata.STEP_MAX);
FEt=fft(Et,gdata.STEP_MAX);
R=FEr./FEin;
T=FEt./FEin;
k0=freq*2*pi/3/10^8;
S11=R.*exp(1i*k0*(210*gata.UL));
S21=T.*exp(1i*k0*(210*gata.UL));
Z=sqrt((((1+S11).^2-S21.^2)./((1-S11).^2-S21.^2)));
n=1./(k0*10*gata.UL).*acos((1-S11.^2+S21.^2)./(2.*S21));
subplot(131);plot(freq,real(Z.*n),'g');
legend('相对磁导率');title('相对磁导率');xlim([2e14 10e14]);
subplot(132);plot(freq,real(n./Z),'r');
legend('相对介电常数');title('相对介电常数');xlim([2e14 10e14]);
subplot(133);plot(freq,abs(R),'g');hold on;plot(freq,abs(T),'r');
legend('透射率','反射率');title('反射率与透射率');xlim([2e14 10e14]);
figure;
n0=1./(k0*10*gata.UL).*(2*pi-acos((1-S11.^2+S21.^2)./(2.*S21)));
subplot(131);plot(freq,real(Z.*n0),'g');
legend('相对磁导率');title('相对磁导率');xlim([2e14 10e14]);
subplot(132);plot(freq,real(n0./Z),'r');
legend('相对介电常数');title('相对介电常数');xlim([2e14 10e14]);
subplot(133);plot(freq,abs(R),'g');hold on;plot(freq,abs(T),'r');
legend('透射率','反射率');title('反射率与透射率');xlim([2e14 10e14]);
El resultado es:
0—> pi (la primera sección de resonancia de espesor, la sección izquierda):
Cuando pi—> 2pi (la segunda sección de resonancia de espesor, la sección central):
Se encuentra que el error no es grande. Si desea continuar midiendo el segmento posterior, puede continuar modificando la función del índice de refracción de acuerdo con la forma de onda sinusoidal y combinar las reglas de selección para resolverlo.
El archivo de código del proyecto se detalla en: https://download.csdn.net/download/hyl1181/12646798