主程序FEMV_main.m
当前程序处于计算非均匀波导的状态。
如需测试均匀波导,修改:
第7、8行:注释掉
R2=[3,4,0.8,1.2,1.2,0.8,0.2,0.2,0.3,0.3]';
R3=[3,4,0,2,2,0,0,0,0.2,0.2]';
第9-14行:换用注释掉的部分
%gm=[R1];
gm=[R1,R2,R3];
%sf='R1';
sf='R1+R2+R3';
%ns=char('R1');
ns=char('R1','R2','R3');
观察图1,获取边界编号。
第30-32行:提供非零狄利克雷边界值及其编号以及有源边界(指有沿边界的线电流源)编号,如没有置零
Nonzero_Dedge=0;
E=0;
Source_Dedge=3;
第40-41行:修改各子域介电常数
%eps_domain=[1];
eps_domain=[1,1.2,3.8];
第43-44行:设定k0与电流源大小
k0=pi*5;
J=[1,0]';