用息票剥离法编程求债券即期利率太复杂了，本人现有一方法，请牛人帮忙看看是否可行：
利用债券原始交易数据，直接用MATLAB 的nlsn（）函数估计NS模型，得到折现因子D（t）数据，再根据公式R（t）=(1/D(t))^(1/t)-1换算成即期利率。