用matlab模拟一个一维光子晶体(中间插入kerr非线性介质),模拟出透过率图。其他的我做出来了,就是非线性介质要用逆向传输矩阵,这点我不知道如何编写,我可以给参数和公式。希望能有大佬指点。若用COMSOL仿真也可以