光子晶体光纤由于其区别于传统光纤而具有的无截止单模传输、可调节色散、高双折射、偏振控制、大的有效面积单模运转和小的有效面积高非线性等特性及其广泛的应用前景,成为当前国内外研究的一个热门课题。为了精确地分析、预测光子晶体光纤的传输特性,人们发展了许多理论分析方法,这些方法是研究光子晶体光纤的基本工具,在光子晶体光纤的研究领域占有很重要的地位。国内外在研究光子晶体方面的文章也非常之多,但是大都注重研究的结论,在众多的文献中能得到好多种光纤的特性信息,但却几乎找不到一种研究方法可以拿来直接用而不用经过和原作者一样的各种知识的繁杂学习的,基于提供一种通用而简单的研究光子晶体光纤的方法,作者通过自己对时域有限差分法(FDTD)和有限元方法的实践探索,总结出利用有限元软件COMSOL Multiphysics实现光子晶体光纤数值模拟的一系列简单可行步骤及后处理过程的MATLAB程序,使一般的研究者只要根据本文给出步骤就可以进行各种光子晶体光纤特性的数值模拟。本论文在系统介绍光子晶体光纤基础知识及各种理论研究方法,并对这些方法优缺点作简单比较的基础上,重点介绍利用有限元软件COMSOL Multiphysics实现光子晶体光纤特性数值模拟的具体方法步骤,并应用该方法计算了空气孔三角形排列、空气孔正方形排列光子晶体光纤的特性参数,并对二者做了简单的比较。最后用该方法计算了类矩形芯光子晶体光纤的偏振特性。即验证了该方法的可行性,又为研究光子晶体光纤研究提供了一定的理论数据。