24小时热门版块排行榜    

查看: 512  |  回复: 6
当前主题已经存档。
【悬赏金币】回答本帖问题,作者哈佳将赠送您 100 个金币

哈佳

金虫 (小有名气)

[求助] 如何用matlab确定VG model上的三个系数(parameter)呢?

如何用matlab确定VG model上的三个系数(parameter)呢?

我已经想到怎样倒出VG model,但怎么用市场上实际的call option data 去推出来呢?是用fminsearch吗?怎么用呢?向各位请教了。

function S=myfun8(nst,Snul,r,q,T,n,vg)
dt=T/n;
tt=[0:dt:T];
omega=log(1-nst(2)^2*nst(1)/2-nst(3)*nst(1))/nst(1);
for s = 1:n
g= gamrnd(dt/nst(1),nst(1));
vg(s+1) = vg(s) + nst(3)*g+ nst(2)*normrnd(0,sqrt(nst(3)));
%S(s+1) = Snul*exp(((r-q)+log(1-nst(2)^2*nst(1)/2-nst(3)*nst(1))/nst(1)).*tt(s+1)+vg(s+1));
end;

S1 = Snul*exp((r-q+omega).*tt+vg);

S=mean(S1);
这个是我写的VG model代码,但是怎么用数据,例如( S&P index data),  去得到优化的三个系数呢?
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nucleus01

至尊木虫 (著名写手)

小木虫北庭节度使

你给出的信息不明确,VG模型的形式,哪个是自变量,哪几个是拟合的参数,最好给出数据。这样才能帮你试。
By the way, VG模型是van Genuchten关于土壤水分的模型吗?
爽是要付出代价的!
2楼2008-04-08 22:30:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxh8372

荣誉版主 (文坛精英)

小木虫微笑研究僧

看着头大

帮顶一下
____/\~~~/\,----(..)/\____//|(\|(^\/___\/\||__||__|-"
3楼2008-04-08 23:25:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

哈佳

金虫 (小有名气)

是的, 说的有点不清楚,抱歉
我这个模型是经济模型,就是
wirte a pricer for european call anf puts under the variance gamma model using the characteristic function approach,
上面是我写的VG model的代码,nst(1)是VG model的nu,nst(2) 是sigma,nst(3)是theta
4楼2008-04-09 02:22:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

哈佳

金虫 (小有名气)

我现在是没有头绪,就是不明白这个characteristic function approach,所以问问大家
5楼2008-04-09 02:24:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

johnliu1983

至尊木虫 (著名写手)

你可以试一下这个函数:lsqcurvefit,怎么用,你查一下help。
6楼2008-04-09 09:48:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

哈佳

金虫 (小有名气)

用fminsearch 求这三个系数的话,问问大家,我订的初始值变化的话,为什么系数值也跟着变化呢?
7楼2008-04-11 02:27:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 哈佳 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 ttH3f1oA6ace 2026-11-04 8/400 2026-11-05 13:44 by My20PsvIExYO
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-05 13:33 by My20PsvIExYO
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 kvewwGtKktqB 2026-11-04 6/300 2026-11-05 13:22 by My20PsvIExYO
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 9/450 2026-11-05 13:22 by My20PsvIExYO
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-05 13:13 by My20PsvIExYO
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 ncTVBM6AOJHF 2026-11-04 9/450 2026-11-05 13:11 by My20PsvIExYO
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ncTVBM6AOJHF 2026-11-04 9/450 2026-11-05 13:11 by My20PsvIExYO
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 13:11 by My20PsvIExYO
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 ncTVBM6AOJHF 2026-11-04 7/350 2026-11-05 13:00 by My20PsvIExYO
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-05 12:51 by My20PsvIExYO
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 6/300 2026-11-05 12:16 by My20PsvIExYO
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 ncTVBM6AOJHF 2026-11-04 7/350 2026-11-05 11:12 by PYZmcHr5kSaq
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 ncTVBM6AOJHF 2026-11-04 8/400 2026-11-05 11:01 by PYZmcHr5kSaq
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ncTVBM6AOJHF 2026-11-04 8/400 2026-11-05 10:59 by PYZmcHr5kSaq
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-05 10:48 by PYZmcHr5kSaq
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 10:48 by PYZmcHr5kSaq
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-05 10:39 by PYZmcHr5kSaq
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 hNpp4Y2CGvSC 2026-11-04 9/450 2026-11-05 09:20 by PYZmcHr5kSaq
[基金申请] 各位基础老师们,看完中央《从严治党》的决定后,对学术界更加有信心了吧 +5 iaeyuan 2026-11-01 5/250 2026-11-04 14:43 by jinchengzi
信息提示
请填处理意见