求助B样条计算利率期限结构……我快崩溃了

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

《基于B 样条对国债利率期限结构的实证研究》
用的和这篇文献几乎一模一样的数据，但运行出来的结果前三个基函数系数很相近，但后三个差距很大，翻来覆去地检查了很多遍代码，问题可能是出现在基函数上，但我实在是不知道自己错在哪，完全是照着文献写的啊
我用的matlab编写的程序，具体如下：
function g=bsplinenew(T,t,k)
g=zeros(length(T)-k-1,length(t));%第i行代表第i个样条基函数在不同的时间所得到的样条函数值
for i=1:length(T)-k-1
    for j=1:length(t)
    g(i,j)=1/(T(i+1)-T(i))/(T(i+2)-T(i))/(T(i+3)-T(i))/(T(i+4)-T(i))*max([(t(j)-T(i)),0])^3+...
        1/(T(i)-T(i+1))/(T(i+2)-T(i+1))/(T(i+3)-T(i+1))/(T(i+4)-T(i+1))*max([(t(j)-T(i+1)),0])^3+...
        1/(T(i)-T(i+2))/(T(i+1)-T(i+2))/(T(i+3)-T(i+2))/(T(i+4)-T(i+2))*max([(t(j)-T(i+2)),0])^3+...
        1/(T(i)-T(i+3))/(T(i+1)-T(i+3))/(T(i+2)-T(i+3))/(T(i+4)-T(i+3))*max([(t(j)-T(i+3)),0])^3+...
        1/(T(i)-T(i+4))/(T(i+1)-T(i+4))/(T(i+2)-T(i+4))/(T(i+3)-T(i+4))*max([(t(j)-T(i+4)),0])^3;
    end
end

上面是基函数的计算，三次B样条，k就是3，T长度为10，所以样条函数个数为6，T是债券时段划分，四个节点，前三后三，分别是-3，-2，-1，0,5,8,20,25,30,35，t是债券现金流时间的序列，比如0.5，1，1.5，2
基函数矩阵就是g，第一行为第一个基函数在所有现金流时间下的基函数取值

哪位大侠能帮帮小弟{:soso_e109:}论文很快就要交了，现在真是欲哭无泪

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏4 回帖

关键词：利率期限结构 期限结构 利率期限 function bspline 计算 结构