%===========================================================
%电磁场与电磁波实验课：
%对导电煤质分界面平面的垂直入射的Matlab动态仿真图
%学生：杰作
%学号：xx 改：201011-20-13-33
%===========================================================
clc; %清除命令窗口代码
clear all; %清除系统所用变量
close all; %关掉系统所用临时窗口
%==========================================================
%系统变量的定义部分
%==========================================================
% w: 极化频率
% t: 时间
% k: 波数 即在2π空间距离内包含的波长个数
% z: 沿传播方向的位置坐标
% faix：x方向上的初始相位
% faiy：y方向上的初始相位
% Exm：x方向上的电场振幅
% Eym：y方向上的电场振幅
% Ex：电场在x方向上的瞬时分量
% Ey：电场在y方向上的瞬时分量
% epsilon:介电常数ε
% mu:磁导率μ
% sigma:电导率σ
% gamma:传播常数γ
% eta:本征阻抗η
syms w t %波特征
w=300e6;
t=[0:1e-9:1e-7];
z=[0:0.00001:0.002];
%=======================================================
%介质1中的电参数
syms gamma1 mu1 epsilon1 sigma1 eta1%介质1中的电参数的定义
mu1=1.0000004;%介质1磁导率
epsilon1=56;%介质1的相对介电常数
sigma1=100;%空气的电导率为
%介质2中的电参数
syms gamma2 mu2 epsilon2 sigma2 eta2%介质2中的电参数的定义
mu2=0.99999;%介质2磁导率
epsilon2=100;%介质2的相对介电常数
sigma2=5.8e10;%介质2的电导率为
%=======================================================
Eim=0.5;
gamma1=j*w*sqrt(mu1*epsilon1*(1-j*sigma1/(w*epsilon1)))
gamma2=j*w*sqrt(mu2*epsilon2*(1-j*sigma2/(w*epsilon2)))
% real(gamma1); %求虚数的实部
% real(gamma2); %求虚数的实部
