Simulation and calculation of electromagnetic parameters of dielectric plates
Master the simulation calculation of the reflection coefficient and transmission coefficient of the harmonic uniform plane electromagnetic wave incident on the dielectric plate; compile the matlab program to calculate the electromagnetic parameters of the dielectric plate by using the reflection coefficient and transmission coefficient.
Use Eastwave software to design the corresponding experimental model first, then design the boundary conditions, grid, excitation source, etc., and run the observation field distribution. The model designed in this experiment (the object uses a dielectric material, the relative permittivity is 12.56) and the grid are:
experimental results
Use MATLAB to program according to the principle to obtain the corresponding electromagnetic parameters. code show as below:
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]);
The result is:
0—>pi (the first section of thickness resonance, the left section):
When pi—>2pi (the second section of thickness resonance, the middle section):
It is found that the error is not large. If you want to continue to measure the back segment, you can continue to modify the refractive index function according to the sinusoidal waveform and combine the selection rules to solve it.
The project code file is detailed at: https://download.csdn.net/download/hyl1181/12646798