<?xml version="1.0" encoding="UTF-8"?>
<rss version="2.0"
	xmlns:content="http://purl.org/rss/1.0/modules/content/"
	xmlns:wfw="http://wellformedweb.org/CommentAPI/"
	xmlns:dc="http://purl.org/dc/elements/1.1/"
	xmlns:atom="http://www.w3.org/2005/Atom"
	xmlns:sy="http://purl.org/rss/1.0/modules/syndication/"
	xmlns:slash="http://purl.org/rss/1.0/modules/slash/"
	>

<channel>
	<title>Stand Alone Complex</title>
	<atom:link href="http://ichentao.com/feed/" rel="self" type="application/rss+xml" />
	<link>http://ichentao.com</link>
	<description>不存在的原创产生原创的拷贝</description>
	<lastBuildDate>Sun, 13 May 2012 22:09:23 +0000</lastBuildDate>
	<language>en</language>
	<sy:updatePeriod>hourly</sy:updatePeriod>
	<sy:updateFrequency>1</sy:updateFrequency>
	<generator>http://wordpress.org/?v=3.2.1</generator>
		<item>
		<title>The Journey to the West S02E05</title>
		<link>http://ichentao.com/blog/archives/2232/</link>
		<comments>http://ichentao.com/blog/archives/2232/#comments</comments>
		<pubDate>Sun, 13 May 2012 22:09:23 +0000</pubDate>
		<dc:creator>Tao</dc:creator>
				<category><![CDATA[日记]]></category>
		<category><![CDATA[CME]]></category>
		<category><![CDATA[Research]]></category>

		<guid isPermaLink="false">http://ichentao.com/?p=2232</guid>
		<description><![CDATA[大概一个月前, 带我实习的法国 boss 开始表达希望我暑假能留下来继续实习的想法, 并且很热心的要帮忙解决工作许可的问题, 于是我一段曲折起伏的心路历程开始了. 在他和我谈话之前, 我的想... ]]></description>
			<content:encoded><![CDATA[<p><img src="https://blufiles.storage.live.com/y1p_h2n76S63g5FhQBsePFCHzEWm1dTkeeNyMEXltMNH61Yc6Zi-xd1BtCzSV-kewLzjz9ery9fXYQ/3199.jpg" alt="CME Group" /></p>

<p>大概一个月前, 带我实习的法国 boss 开始表达希望我暑假能留下来继续实习的想法, 并且很热心的要帮忙解决工作许可的问题, 于是我一段曲折起伏的心路历程开始了.</p>

<p>在他和我谈话之前, 我的想法是暑假专心做研究, 可是法国人偏偏提出了一个让我很难拒绝的提议: 就是暑期工作时可以用50%的时间做数学. 之所以诱人, 是因为这样的安排已经是一个普通 quant 的时间分配了(我顿时体会到 PhD 三个字母在对方心目中的分量), 可另一方面, 实在也担心老板会反对我暑假继续工作. 于是带着犹豫的心理, 回学校找老板谈话. 老板当然是不用听我说完就能明白我的心思, 跟我分析了好久目光长短的问题, 然后表态说可以在不超过某一限度的范围内继续实习. 我像吃了定心丸一样, 开始着手要解决工作许可的问题. 可没想到过两天, 系里一年一度 Menger Day 大会上, 颁奖嘉宾突然叫道我的名字让我上台去领暑期研究的奖金, 之前老板说因为我还没开始做, 所以这次奖金的申请是不指望的, 结果自然是让我大感意外, 同时是真想不通为什么最近好事多的让我无法抉择. 说实话, 钱对现在的我来说是很有用的, 但更重要的是, 数学系没有选择那么些已经在读的博士, 偏偏把奖金给了我, 真的让人很有压力, 再一想导师对我整个博士期间的规划, 千里之行, 始于足下, 为了开一个好头, 我还是放弃了暑假实习的想法.</p>

<p>现在回顾这一段实习经历还有点过早, 因为有很多东西在忙碌的实习期内来不及消化, 特别是数学上的. 非常幸运, 在读博前可以有机会了解真实的市场, 的确学到不少学校内学不到的东西, 而且能加入 Clearing 部门参与 OTC Clearing 的项目, 我记得 Quant 组的 Director 在第一次跟我谈话时就说, OTC Clearing 现在是金融业的一件大事, 不光如此, 未来 OTC 交易转 Exchange 肯定也会成为趋势, 所以, 能有机会参与到最前沿的工作当中, 心里很高兴. 这两个多月里, 我发现大部分 quant 们都是跟普通人不太一样的 math geek, 自己也受到一些影响, 感觉如此过一生, 大概就是最适合自己也感觉最舒服的一条路.</p>
]]></content:encoded>
			<wfw:commentRss>http://ichentao.com/blog/archives/2232/feed/</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>The Black Keys &amp; Arctic Monkeys @ United Center</title>
		<link>http://ichentao.com/blog/archives/2223/</link>
		<comments>http://ichentao.com/blog/archives/2223/#comments</comments>
		<pubDate>Thu, 29 Mar 2012 02:11:33 +0000</pubDate>
		<dc:creator>Tao</dc:creator>
				<category><![CDATA[影音]]></category>
		<category><![CDATA[Arctic Monkeys]]></category>
		<category><![CDATA[The Black Keys]]></category>

		<guid isPermaLink="false">http://ichentao.com/?p=2223</guid>
		<description><![CDATA[其实是上个星期一去的, 后来一直比较忙, 所以拖到今天才写. 冬天和同学去看 Feist 的时候, 在一个破破烂烂的剧场里, 我们坐在二层, 当时我还说坐着听听 Feist 可以, 但要是看摇滚乐队演出就比... ]]></description>
			<content:encoded><![CDATA[<p>其实是上个星期一去的, 后来一直比较忙, 所以拖到今天才写.</p>

<p>冬天和同学去看 Feist 的时候, 在一个破破烂烂的剧场里, 我们坐在二层, 当时我还说坐着听听 Feist 可以, 但要是看摇滚乐队演出就比较奇怪了, 结果没想到这次 The Black Keys 和 Arctic Monkeys 的演出, 我们便是在联合中心球馆的 300 Level 坐着看完的. 之前因为在附近的意大利餐馆吃饭耽误了时间, 到场的时候已经迟了, Arctic Monkeys 开始唱 I Bet You Look Good on the Dance Floor 的时候我们还在上楼梯找座位, 懊悔不已.</p>

<p>想到去年 Lollapalooza 上看 Muse 的经历, 结合网上一些 Placebo 在美国演出的视频, 再加上这次看 Arctic Monkeys, 真的发现英国乐队在美国的不易, Muse 已经是我看到的在美国最成功的英国乐队了. 就拿这次 Arctic Monkeys 来说, 我错过的部分不算, 全场也就 Fluorescent Adolescent 响起时欢呼声稍微大一点, 其它时间观众太安静了, 这还是在猴子们唱了所有他们最流行的歌曲的情况下. 我只能很欣慰的讲, 中场出去抽烟时, 听到有人说他们是来看北极猴的.</p>

<p><a href="http://www.flickr.com/photos/cherntodd/6852875902/" title="DSCN0210 by Todd - 魔群月光, on Flickr"><img src="http://farm8.staticflickr.com/7198/6852875902_05c8d2d3e6.jpg" width="500" height="375" alt="DSCN0210"></a></p>

<p>The Black Keys 受到的待遇明显不同, 从头至尾高分贝欢呼不断, 尽管他们很多首歌听起来很重复, 但可以看出他们绝对是一支踏实认真的乐队, 再加上风格特别合适现场演出, 所以, 我觉得他们还是配得上观众的热情. 之所以认为 The Black Keys 朴实努力, 是因为 1. 主唱的话很少, 不像一些乐队中途要说很多, 他就一个劲的唱; 2. 唱的歌很多, 让观众过足了瘾, 而且返场更炫, 能让人不由自主的起立跟着他们摆动. 同学说, 感觉这次 The Black Keys 把所有的歌都唱了一遍.</p>

<p>我会嘲笑美国人那么粉 Muse, 可我又非常赞同他们对于 The Black Keys 的好感, 事实上, 我还很希望能在今年的 Lolla 上再看他们一次.</p>

<p><a href="http://www.flickr.com/photos/cherntodd/6852904876/" title="DSCN0261 by Todd - 魔群月光, on Flickr"><img src="http://farm7.staticflickr.com/6111/6852904876_e2df8584d9.jpg" width="500" height="375" alt="DSCN0261"></a></p>
]]></content:encoded>
			<wfw:commentRss>http://ichentao.com/blog/archives/2223/feed/</wfw:commentRss>
		<slash:comments>7</slash:comments>
		</item>
		<item>
		<title>The Journey to the West S02E04 -- PhD</title>
		<link>http://ichentao.com/blog/archives/2220/</link>
		<comments>http://ichentao.com/blog/archives/2220/#comments</comments>
		<pubDate>Mon, 12 Mar 2012 23:30:02 +0000</pubDate>
		<dc:creator>Tao</dc:creator>
				<category><![CDATA[日记]]></category>
		<category><![CDATA[PhD]]></category>

		<guid isPermaLink="false">http://ichentao.com/?p=2220</guid>
		<description><![CDATA[今天上班很轻松, 全组和在纽约的部门开了一个视频会议. 剩下的时间我改了几个 forward curve, 于是就下班回学校上课去了. 随机偏微分方程的课上问了老师几个问题, 然后同样很轻松的水过. 不过... ]]></description>
			<content:encoded><![CDATA[<p>今天上班很轻松, 全组和在纽约的部门开了一个视频会议. 剩下的时间我改了几个 forward curve, 于是就下班回学校上课去了. 随机偏微分方程的课上问了老师几个问题, 然后同样很轻松的水过.</p>

<p>不过, 今天正式确认下来, 拿到了数学系的 offer, 接着的三至四年将进行我的 PhD 生活, 这意味着从小学算起到30岁左右毕业时, 我将整整在校园里待25年, 很恐怖的一个数字.</p>

<p>为什么读 PhD, 因为心里有着做 pure quant 的梦想, 一方面此类工作对学历的要求都是博士, 另一方面我也深感自己还不够 quantitative, 所以希望能够得到更多的数学训练. 在美国待了快两年, 发现了和国内不同的地方. 出国前, 听到最多的声音是, 继续读书再不出来工作就和社会脱节了; 提起学数学, 很多人第一反应是"你要当老师啊", 总之我认为国内弥漫着"读书无用论"的风气. 在美国的这段时间里, 参加过一些招聘, 听过一些业内人士的讲座, 在实习的时候和同事有一些交流, 大家总把"PhD"挂在嘴边:"你们不是 PhD, 所以你们只能怎样怎样", "这个问题以前也有 PhD 来问, 尚且花了很多时间才搞懂", 大概意思就是 PhD 在人们的眼中是一群 super smart guys, 这个行业的顶尖人才都是 PhD 出身. 导师在跟我谈的时候说的没错, 有了博士的头衔, 别人对待你的方式就不同了, 比如某某某, 现在在公司只有他能和公司的一号人物在一起研究论文. 美国好吗? 我觉得其实一般, 这里有无形的歧视, 有难吃的食物, 人们的品味普遍出现问题, 但我能确信一点, 美国是最发达的国家, 因为这个国家真正的尊重知识和追求创新. 想好好干一翻工作, 还是这个国家最合适.</p>

<p>还有选择学校的问题, 事实上, 我觉得除了整体排名低一点, 每个月生活费比一些学校少一点以外, 其它各方面都算非常理想了. 我想做的是金融数学, 数遍美国各个学校的数学系, 基本上在这个方向上都不超过两位教授, 最厉害的是卡耐基梅隆, 有两个大牛, 接下来可能是普林斯顿, 哥伦比亚等等, 选择面其实很窄, 除了现在在的学校以外就只剩下顶级牛校. 在这种情况下, 本校也可以接受了, 有两个金融数学方向的老师, 两个人都算我的导师, 真正官方的导师是和 Brigo 这样的大牛齐名的人物, 同时从主观角度来说, 当这样的大牛能够从你的身上发现别人看不到的优点, 能够交给你很难的题目并且说 You have the brain for it 的时候, 我觉得无法离开. 古人说, 知遇之恩当涌泉相报, 我想就是这个道理.</p>

<p>Anyway, 接下来的三至四年会比较辛苦, 因为这一组是号称数学系最苦逼的一组, 导师对我发火是迟早的一件事情, 但是不付出努力, 别的只能免谈了, 在国内的时候荒废了那么多年, 那就接下来补回来好了.</p>

<p>其实, 想到每年都还有 Pitchfork, Lollapalooza 和数不清的演唱会, 还可以在这里等着 Placebo 来, 真是一件非常开心的事情.</p>
]]></content:encoded>
			<wfw:commentRss>http://ichentao.com/blog/archives/2220/feed/</wfw:commentRss>
		<slash:comments>14</slash:comments>
		</item>
		<item>
		<title>The Journey to the West S02E03 -- Internship</title>
		<link>http://ichentao.com/blog/archives/2168/</link>
		<comments>http://ichentao.com/blog/archives/2168/#comments</comments>
		<pubDate>Fri, 17 Feb 2012 03:18:44 +0000</pubDate>
		<dc:creator>Tao</dc:creator>
				<category><![CDATA[日记]]></category>
		<category><![CDATA[CME]]></category>
		<category><![CDATA[工作]]></category>

		<guid isPermaLink="false">http://ichentao.com/?p=2168</guid>
		<description><![CDATA[三月份开始要在 CME Group(芝加哥商业交易所集团) 实习, 职位是 Quantitative Risk Research Intern, 所在小组负责利率衍生品场外交易清算的模型建立, 开发和测试. 本来给自己定的目标是通过前五, 六次... ]]></description>
			<content:encoded><![CDATA[<p><img alt="" src="https://blufiles.storage.live.com/y1p0OlHbkQ1pPwn3c8qYPNUw3r4Tdp23pfbcMSbtfuMOJmDUbYUgOBo-_6fOcsEN-n4jkNUj9KJZ-I/CME-NYMEX.jpg" title="CME" class="alignnone" width="500" height="333" /></p>

<p>三月份开始要在 <a href="http://www.cmegroup.com/">CME Group</a>(<a href="http://www.cmegroup.com/cn-s/">芝加哥商业交易所集团</a>) 实习, 职位是 Quantitative Risk Research Intern, 所在小组负责利率衍生品场外交易清算的模型建立, 开发和测试. 本来给自己定的目标是通过前五, 六次面试进行锻炼, 没想到第三次就成功了.</p>

<p>说实话, 自己在找实习方面运气还是不错的: 第二学期投的简历相对多一点, 得到了美国一家顶尖的自营交易公司的电面, 面试内容是两位数乘两位数的心算和概率题的心算, 结果由于不熟悉面试的流程没发挥好; 上学期只投了一份简历, 得到了面试机会, 但因为没什么统计和大规模数据处理的经验, 再加上紧张, 又失败了, 决定毕业读 PhD 以后, 去 NYU 全美金融工程专业招聘会上打了一圈酱油, 把普华永道的 representative 搞的 pissed off, 我也被他搞的 pissed off; 这学期要不是老师帮忙把简历递到 CME, 估计是一份都不会投的.</p>

<p>因为心里还有别的小九九, 所以去 CME 面试时抱的是无所谓的态度, 心态也前所未有的轻松. 虽然经历了四个人轮番轰炸, 但还算应付自如. 第一个面试官本应该是未来的 boss, 结果他当天有事, 临时换成一个刚从纽约调来的女孩儿, 此人麻省理工数学系代数方向毕业, 我有两个问题没答好, 一是短期和长期 yield curve 中间的 gap 可以通过哪种衍生品价格补全, 我想了半天估摸着说可能是掉期期权, 她说我把问题搞复杂了, 掉期期货就可以; 二是掉期期权的 Black-Scholes 定价公式, 我这种从来不记公式的人当场就傻了, 又是搞了半天在黑板上从头到尾把公式推了一遍, 然后她建议我有空还是要把 John Hull 的那本"圣经"好好看看... 第二位是个有些诡异的美国人, 见面也不和我打招呼, 然后一开始的气氛就有些尴尬, 我们在我的简历上花了好多时间, 后来我发现跟他说话还是挺轻松的, 于是面试就有往聊天方向发展的趋势, 他问我对 Black-Scholes model 的看法, 我扯了一通然后往自己比较擅长的 Heston model 上靠, 末了不忘加一句, "其实 Heston model 还有完善的余地", 他又问我为什么从数学转到金融, 我说我想把数学运用在实际当中, 他说那还有别的方向可以用数学, 我说"你指工程吗, 我感觉也就用用傅里叶变换, 挺无趣的", 他想了半天居然说了一句"我觉得还要用偏微分方程"... 快结束时他问我知不知道 Garch model, 我就完全放肆了, 心想 Heston 我都知道, Garch 不是反而还差些吗, 于是说仅仅听到过这个模型, 还反问他"这个模型好不好呀"... 第三第四个面试官都比较随和, 我答的也不错, 问题大概就是 Heston model 好在哪里, 为什么要用蒙特卡洛方法, 怎么用蒙特卡洛方法, 还有几个期权方面的基本数学问题. 走出办公楼的时候我比较轻松, 心想反正过不了我就继续留在学校看美女. 等了一个星期, 中间听说班上另一个去面试的同学录取了, 我则是收到邮件说周四再去一次, 这次接我上楼的人上次面过我, 对我说不要担心, 今天跟我谈的是我 boss, 随便聊一聊. 我一开始都没反应过来, 后来面着面着发现他开始跟我谈工作时间, 工作前景的问题, 我想这下有戏了, 也的确是果不其然.</p>

<p>走出办公楼, 给家里打了一个电话通报一声, 兴奋了大概一个小时. 然后就是纠结实习和讨论班的冲突问题, 还有半个月, 希望船到桥头自然直.</p>

<p>有实习经验的话, 后面工作也会比较好找, 比如利率衍生品的方向, QRM 这样的公司应该是很有希望的. 不过, 因为心里面还抱着做100%时间用来推公式的 pure quant 的梦想, 所以下定决心要读 PhD, 希望再经历一个三到四年的训练, 能够达到自己期望的标准. 美国的同学曾经问, 什么样的工作可以让我放弃读博. 我说, 至少起薪10万刀吧.</p>

<p>说的有些过分了, "堤高于岸, 浪必催之", 低调, 慎行.</p>
]]></content:encoded>
			<wfw:commentRss>http://ichentao.com/blog/archives/2168/feed/</wfw:commentRss>
		<slash:comments>9</slash:comments>
		</item>
		<item>
		<title>Fixed Income Modeling Review 8</title>
		<link>http://ichentao.com/blog/archives/2161/</link>
		<comments>http://ichentao.com/blog/archives/2161/#comments</comments>
		<pubDate>Sat, 04 Feb 2012 23:19:17 +0000</pubDate>
		<dc:creator>Tao</dc:creator>
				<category><![CDATA[技术]]></category>
		<category><![CDATA[Fixed Income]]></category>
		<category><![CDATA[Mathematical Finance]]></category>

		<guid isPermaLink="false">http://ichentao.com/?p=2161</guid>
		<description><![CDATA[Feynman-Kac formula, named after Richard Feynman and Mark Kac, establishes a link between parabolic partial differential equations (PDEs) and stochastic processes. Suppose follows the stochastic process and if is defined as then satisfies , . The theorem... ]]></description>
			<content:encoded><![CDATA[<p>Feynman-Kac formula, named after Richard Feynman and Mark Kac, establishes a link between parabolic partial differential equations (PDEs) and stochastic processes. Suppose <img src="http://l.wordpress.com/latex.php?latex=S%28t%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="S(t)" style="vertical-align:-20%;" class="tex" alt="S(t)" /> follows the stochastic process</p>

<p><img src="http://l.wordpress.com/latex.php?latex=dS%28t%29%3D%5Cmu%28t%2C%5Comega%29dt%2B%5Csigma%28t%2C%5Comega%29dW%28t%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="dS(t)=\mu(t,\omega)dt+\sigma(t,\omega)dW(t)" style="vertical-align:-20%;" class="tex" alt="dS(t)=\mu(t,\omega)dt+\sigma(t,\omega)dW(t)" /></p>

<p>and if <img src="http://l.wordpress.com/latex.php?latex=V&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="V" style="vertical-align:-20%;" class="tex" alt="V" /> is defined as</p>

<p><img src="http://l.wordpress.com/latex.php?latex=V%28t%2CS%28t%29%29%3DE%5EQ_t%5BH%28T%2CS%28T%29%29%5D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="V(t,S(t))=E^Q_t[H(T,S(T))]" style="vertical-align:-20%;" class="tex" alt="V(t,S(t))=E^Q_t[H(T,S(T))]" /></p>

<p>then <img src="http://l.wordpress.com/latex.php?latex=V&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="V" style="vertical-align:-20%;" class="tex" alt="V" /> satisfies</p>

<p><img src="http://l.wordpress.com/latex.php?latex=%5Cfrac%7B%5Cpartial%20V%7D%7B%5Cpartial%20t%7D%2B%5Cmu%5Cfrac%7B%5Cpartial%20V%7D%7B%5Cpartial%20S%7D%2B%5Cfrac%7B1%7D%7B2%7D%5Csigma%5E2%5Cfrac%7B%5Cpartial%5E2%20V%7D%7B%5Cpartial%20S%5E2%7D%3D0&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\frac{\partial V}{\partial t}+\mu\frac{\partial V}{\partial S}+\frac{1}{2}\sigma^2\frac{\partial^2 V}{\partial S^2}=0" style="vertical-align:-20%;" class="tex" alt="\frac{\partial V}{\partial t}+\mu\frac{\partial V}{\partial S}+\frac{1}{2}\sigma^2\frac{\partial^2 V}{\partial S^2}=0" />, <img src="http://l.wordpress.com/latex.php?latex=V%28T%2CS%28T%29%29%3DH%28T%2CS%28T%29%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="V(T,S(T))=H(T,S(T))" style="vertical-align:-20%;" class="tex" alt="V(T,S(T))=H(T,S(T))" />.</p>

<p>The theorem can be stated the other way around. Thus, Feynman-Kac formula could be applied to interest rate models, for example, under Hull-White model, zero coupon bond price <img src="http://l.wordpress.com/latex.php?latex=P%28t%2CT%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="P(t,T)" style="vertical-align:-20%;" class="tex" alt="P(t,T)" /> satisfies</p>

<p><img src="http://l.wordpress.com/latex.php?latex=%5Cfrac%7B%5Cpartial%20P%7D%7B%5Cpartial%20t%7D%2B%28%5Ctheta%28t%29-ar%28t%29%29%5Cfrac%7B%5Cpartial%20P%7D%7B%5Cpartial%20r%7D%2B%5Cfrac%7B1%7D%7B2%7D%5Csigma%5E2%5Cfrac%7B%5Cpartial%5E2%20P%7D%7B%5Cpartial%20r%5E2%7D%3Dr%28t%29P&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\frac{\partial P}{\partial t}+(\theta(t)-ar(t))\frac{\partial P}{\partial r}+\frac{1}{2}\sigma^2\frac{\partial^2 P}{\partial r^2}=r(t)P" style="vertical-align:-20%;" class="tex" alt="\frac{\partial P}{\partial t}+(\theta(t)-ar(t))\frac{\partial P}{\partial r}+\frac{1}{2}\sigma^2\frac{\partial^2 P}{\partial r^2}=r(t)P" />, <img src="http://l.wordpress.com/latex.php?latex=P%28T%2CT%29%3D1&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="P(T,T)=1" style="vertical-align:-20%;" class="tex" alt="P(T,T)=1" /></p>

<p>we can deduce the analytic solution for <img src="http://l.wordpress.com/latex.php?latex=P%28t%2CT%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="P(t,T)" style="vertical-align:-20%;" class="tex" alt="P(t,T)" /> by using the theorem described above.</p>

<p>Another application of Feynman-Kac formula in finance is discretization: explicit schemes and implicit schemes. For explicit schemes, we use such discretization method:</p>

<p><img src="http://l.wordpress.com/latex.php?latex=%5Cleft%28%5Cfrac%7B%5Cpartial%20S%7D%7B%5Cpartial%20r%7D%5Cright%29_%7Bi%2Cj%7D%5Capprox%5Cfrac%7BS_%7Bi%2Cj%2B1%7D-S_%7Bi%2Cj-1%7D%7D%7B2%5CDelta%20r%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\left(\frac{\partial S}{\partial r}\right)_{i,j}\approx\frac{S_{i,j+1}-S_{i,j-1}}{2\Delta r}" style="vertical-align:-20%;" class="tex" alt="\left(\frac{\partial S}{\partial r}\right)_{i,j}\approx\frac{S_{i,j+1}-S_{i,j-1}}{2\Delta r}" />,<br />
<img src="http://l.wordpress.com/latex.php?latex=%5Cleft%28%5Cfrac%7B%5Cpartial%5E2%20S%7D%7B%5Cpartial%20r%5E2%7D%5Cright%29_%7Bi%2Cj%7D%5Capprox%5Cfrac%7BS_%7Bi%2Cj%2B1%7D-2S_%7Bi%2Cj%7D%2BS_%7Bi%2Cj-1%7D%7D%7B%28%5CDelta%20r%29%5E2%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\left(\frac{\partial^2 S}{\partial r^2}\right)_{i,j}\approx\frac{S_{i,j+1}-2S_{i,j}+S_{i,j-1}}{(\Delta r)^2}" style="vertical-align:-20%;" class="tex" alt="\left(\frac{\partial^2 S}{\partial r^2}\right)_{i,j}\approx\frac{S_{i,j+1}-2S_{i,j}+S_{i,j-1}}{(\Delta r)^2}" />,<br />
<img src="http://l.wordpress.com/latex.php?latex=%5Cleft%28%5Cfrac%7B%5Cpartial%20S%7D%7B%5Cpartial%20t%7D%5Cright%29_%7Bi%2Cj%7D%5Capprox%5Cfrac%7BS_%7Bi%2Cj%7D-S_%7Bi-1%2Cj%7D%7D%7B%5CDelta%20t%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\left(\frac{\partial S}{\partial t}\right)_{i,j}\approx\frac{S_{i,j}-S_{i-1,j}}{\Delta t}" style="vertical-align:-20%;" class="tex" alt="\left(\frac{\partial S}{\partial t}\right)_{i,j}\approx\frac{S_{i,j}-S_{i-1,j}}{\Delta t}" /></p>

<p>substitute these into the PDE then the equation could be written as</p>

<p><img src="http://l.wordpress.com/latex.php?latex=S_%7Bi-1%2Cj%7D%3DA_%7Bi%2Cj%7DS_%7Bi%2Cj%2B1%7D%2BB_%7Bi%2Cj%7DS_%7Bi%2Cj%7D%2BC_%7Bi%2Cj%7DS_%7Bi%2Cj-1%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="S_{i-1,j}=A_{i,j}S_{i,j+1}+B_{i,j}S_{i,j}+C_{i,j}S_{i,j-1}" style="vertical-align:-20%;" class="tex" alt="S_{i-1,j}=A_{i,j}S_{i,j+1}+B_{i,j}S_{i,j}+C_{i,j}S_{i,j-1}" /></p>

<p>This kind of scheme has equivalence to trinomial tree valuation</p>

<p><img src="http://l.wordpress.com/latex.php?latex=S_%7Bi-1%2Cj%7D%3Dp_uS_%7Bi%2Cj%2B1%7D%2Bp_mS_%7Bi%2Cj%7D%2Bp_dS_%7Bi%2Cj-1%7D-r_j%5CDelta%20tS_%7Bi%2Cj%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="S_{i-1,j}=p_uS_{i,j+1}+p_mS_{i,j}+p_dS_{i,j-1}-r_j\Delta tS_{i,j}" style="vertical-align:-20%;" class="tex" alt="S_{i-1,j}=p_uS_{i,j+1}+p_mS_{i,j}+p_dS_{i,j-1}-r_j\Delta tS_{i,j}" /></p>

<p>and setting the size of rate step <img src="http://l.wordpress.com/latex.php?latex=%5CDelta%20r%3D%5Csigma%283%5CDelta%20t%29%5E%7B1%2F2%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\Delta r=\sigma(3\Delta t)^{1/2}" style="vertical-align:-20%;" class="tex" alt="\Delta r=\sigma(3\Delta t)^{1/2}" /> leads to the most efficient approximation, under which the scheme becomes 4th order w.r.t. to <img src="http://l.wordpress.com/latex.php?latex=r&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="r" style="vertical-align:-20%;" class="tex" alt="r" />.</p>

<p>While applying explicit finite difference methods, boundary conditions have very little effect. We need to have</p>

<p><img src="http://l.wordpress.com/latex.php?latex=%5Csigma%5E2_%7Bi%2Cj%7D%5Cfrac%7B%5CDelta%20t%7D%7B%5CDelta%20r%29%5E2%7D%5Cleq%5Cfrac%7B1%7D%7B2%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\sigma^2_{i,j}\frac{\Delta t}{\Delta r)^2}\leq\frac{1}{2}" style="vertical-align:-20%;" class="tex" alt="\sigma^2_{i,j}\frac{\Delta t}{\Delta r)^2}\leq\frac{1}{2}" /></p>

<p>hold cause otherwise small numerical errors can grow uncontrollably.</p>

<p>For implicit scheme, we discretize the PDE like this:</p>

<p><img src="http://l.wordpress.com/latex.php?latex=%5Cleft%28%5Cfrac%7B%5Cpartial%20S%7D%7B%5Cpartial%20r%7D%5Cright%29_%7Bi%2Cj%7D%5Capprox%5Cfrac%7BS_%7Bi-1%2Cj%2B1%7D-S_%7Bi-1%2Cj-1%7D%7D%7B2%5CDelta%20r%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\left(\frac{\partial S}{\partial r}\right)_{i,j}\approx\frac{S_{i-1,j+1}-S_{i-1,j-1}}{2\Delta r}" style="vertical-align:-20%;" class="tex" alt="\left(\frac{\partial S}{\partial r}\right)_{i,j}\approx\frac{S_{i-1,j+1}-S_{i-1,j-1}}{2\Delta r}" />,<br />
<img src="http://l.wordpress.com/latex.php?latex=%5Cleft%28%5Cfrac%7B%5Cpartial%5E2%20S%7D%7B%5Cpartial%20r%5E2%7D%5Cright%29_%7Bi%2Cj%7D%5Capprox%5Cfrac%7BS_%7Bi-1%2Cj%2B1%7D-2S_%7Bi-1%2Cj%7D%2BS_%7Bi-1%2Cj-1%7D%7D%7B%28%5CDelta%20r%29%5E2%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\left(\frac{\partial^2 S}{\partial r^2}\right)_{i,j}\approx\frac{S_{i-1,j+1}-2S_{i-1,j}+S_{i-1,j-1}}{(\Delta r)^2}" style="vertical-align:-20%;" class="tex" alt="\left(\frac{\partial^2 S}{\partial r^2}\right)_{i,j}\approx\frac{S_{i-1,j+1}-2S_{i-1,j}+S_{i-1,j-1}}{(\Delta r)^2}" />,<br />
<img src="http://l.wordpress.com/latex.php?latex=%5Cleft%28%5Cfrac%7B%5Cpartial%20S%7D%7B%5Cpartial%20t%7D%5Cright%29_%7Bi%2Cj%7D%5Capprox%5Cfrac%7BS_%7Bi%2Cj%7D-S_%7Bi-1%2Cj%7D%7D%7B%5CDelta%20t%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\left(\frac{\partial S}{\partial t}\right)_{i,j}\approx\frac{S_{i,j}-S_{i-1,j}}{\Delta t}" style="vertical-align:-20%;" class="tex" alt="\left(\frac{\partial S}{\partial t}\right)_{i,j}\approx\frac{S_{i,j}-S_{i-1,j}}{\Delta t}" /></p>

<p>then Feynman-Kac formula becomes</p>

<p><img src="http://l.wordpress.com/latex.php?latex=A_%7Bi%2Cj%7DS_%7Bi-1%2Cj%2B1%7D%2BB_%7Bi%2Cj%7DS_%7Bi-1%2Cj%7D%2BC_%7Bi%2Cj%7DS_%7Bi-1%2Cj-1%7D%3DS_%7Bi%2Cj%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="A_{i,j}S_{i-1,j+1}+B_{i,j}S_{i-1,j}+C_{i,j}S_{i-1,j-1}=S_{i,j}" style="vertical-align:-20%;" class="tex" alt="A_{i,j}S_{i-1,j+1}+B_{i,j}S_{i-1,j}+C_{i,j}S_{i-1,j-1}=S_{i,j}" /></p>

<p>It has a relationship to trinomial trees:</p>

<p><img src="http://l.wordpress.com/latex.php?latex=%281-r_j%5CDelta%20t%29S_%7Bi%2Cj%7D%5Capprox%20p_uS_%7Bi-1%2Cj%2B1%7D%2Bp_mS_%7Bi-1%2Cj%7D%2Bp_dS_%7Bi-1%2Cj-1%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="(1-r_j\Delta t)S_{i,j}\approx p_uS_{i-1,j+1}+p_mS_{i-1,j}+p_dS_{i-1,j-1}" style="vertical-align:-20%;" class="tex" alt="(1-r_j\Delta t)S_{i,j}\approx p_uS_{i-1,j+1}+p_mS_{i-1,j}+p_dS_{i-1,j-1}" /></p>

<p>Since we use time <img src="http://l.wordpress.com/latex.php?latex=i-1&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="i-1" style="vertical-align:-20%;" class="tex" alt="i-1" /> value to calculate <img src="http://l.wordpress.com/latex.php?latex=S_%7Bi%2Cj%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="S_{i,j}" style="vertical-align:-20%;" class="tex" alt="S_{i,j}" />, then if we want to do pricing backward, we need to solve a system of linear equations. Also, we need to care about the boundary conditions where we usually set <img src="http://l.wordpress.com/latex.php?latex=r_0%3D0&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="r_0=0" style="vertical-align:-20%;" class="tex" alt="r_0=0" /> and <img src="http://l.wordpress.com/latex.php?latex=r_M%3Dr_%5Cinfty&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="r_M=r_\infty" style="vertical-align:-20%;" class="tex" alt="r_M=r_\infty" /> which is large enough. Hence, implicit scheme is slightly more difficult to implement than the explicit one. However, it is always stable and convergent. In practice, one method that is used very often is called Crank-Nicolson method which combines both explicit scheme and implicit scheme:</p>

<p><img src="http://l.wordpress.com/latex.php?latex=%5Cleft%28%5Cfrac%7B%5Cpartial%20S%7D%7B%5Cpartial%20r%7D%5Cright%29_%7Bi%2Cj%7D%5Capprox%281-%5Ctheta%29%5Cfrac%7BS_%7Bi%2Cj%2B1%7D-S_%7Bi%2Cj-1%7D%7D%7B2%5CDelta%20r%7D%2B%5Ctheta%20%5Cfrac%7BS_%7Bi-1%2Cj%2B1%7D-S_%7Bi-1%2Cj-1%7D%7D%7B2%5CDelta%20r%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\left(\frac{\partial S}{\partial r}\right)_{i,j}\approx(1-\theta)\frac{S_{i,j+1}-S_{i,j-1}}{2\Delta r}+\theta \frac{S_{i-1,j+1}-S_{i-1,j-1}}{2\Delta r}" style="vertical-align:-20%;" class="tex" alt="\left(\frac{\partial S}{\partial r}\right)_{i,j}\approx(1-\theta)\frac{S_{i,j+1}-S_{i,j-1}}{2\Delta r}+\theta \frac{S_{i-1,j+1}-S_{i-1,j-1}}{2\Delta r}" />,<br />
<img src="http://l.wordpress.com/latex.php?latex=%5Cleft%28%5Cfrac%7B%5Cpartial%5E2%20S%7D%7B%5Cpartial%20r%5E2%7D%5Cright%29_%7Bi%2Cj%7D%5Capprox%281-%5Ctheta%29%5Cfrac%7BS_%7Bi%2Cj%2B1%7D-2S_%7Bi%2Cj%7D%2BS_%7Bi%2Cj-1%7D%7D%7B%28%5CDelta%20r%29%5E2%7D%2B%5Ctheta%20%5Cfrac%7BS_%7Bi-1%2Cj%2B1%7D-2S_%7Bi-1%2Cj%7D%2BS_%7Bi-1%2Cj-1%7D%7D%7B%28%5CDelta%20r%29%5E2%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\left(\frac{\partial^2 S}{\partial r^2}\right)_{i,j}\approx(1-\theta)\frac{S_{i,j+1}-2S_{i,j}+S_{i,j-1}}{(\Delta r)^2}+\theta \frac{S_{i-1,j+1}-2S_{i-1,j}+S_{i-1,j-1}}{(\Delta r)^2}" style="vertical-align:-20%;" class="tex" alt="\left(\frac{\partial^2 S}{\partial r^2}\right)_{i,j}\approx(1-\theta)\frac{S_{i,j+1}-2S_{i,j}+S_{i,j-1}}{(\Delta r)^2}+\theta \frac{S_{i-1,j+1}-2S_{i-1,j}+S_{i-1,j-1}}{(\Delta r)^2}" />,<br />
<img src="http://l.wordpress.com/latex.php?latex=%5Cleft%28%5Cfrac%7B%5Cpartial%20S%7D%7B%5Cpartial%20t%7D%5Cright%29_%7Bi%2Cj%7D%5Capprox%281-%5Ctheta%29%5Cfrac%7BS_%7Bi%2Cj%7D-S_%7Bi-1%2Cj%7D%7D%7B%5CDelta%20t%7D%2B%5Ctheta%20%5Cfrac%7BS_%7Bi%2Cj%7D-S_%7Bi-1%2Cj%7D%7D%7B%5CDelta%20t%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\left(\frac{\partial S}{\partial t}\right)_{i,j}\approx(1-\theta)\frac{S_{i,j}-S_{i-1,j}}{\Delta t}+\theta \frac{S_{i,j}-S_{i-1,j}}{\Delta t}" style="vertical-align:-20%;" class="tex" alt="\left(\frac{\partial S}{\partial t}\right)_{i,j}\approx(1-\theta)\frac{S_{i,j}-S_{i-1,j}}{\Delta t}+\theta \frac{S_{i,j}-S_{i-1,j}}{\Delta t}" /></p>

<p>Although Crank-Nicolson method is a little bit harder to implement than both the explicit and implicit scheme, it is as stable as the fully implicit method and is the second-order accurate w.r.t. time step which is better than explicit and implicit schemes.</p>
]]></content:encoded>
			<wfw:commentRss>http://ichentao.com/blog/archives/2161/feed/</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>Fixed Income Modeling Review 7</title>
		<link>http://ichentao.com/blog/archives/2150/</link>
		<comments>http://ichentao.com/blog/archives/2150/#comments</comments>
		<pubDate>Sat, 04 Feb 2012 03:48:57 +0000</pubDate>
		<dc:creator>Tao</dc:creator>
				<category><![CDATA[技术]]></category>
		<category><![CDATA[Fixed Income]]></category>
		<category><![CDATA[Mathematical Finance]]></category>

		<guid isPermaLink="false">http://ichentao.com/?p=2150</guid>
		<description><![CDATA[In finance, there are two major applications of the Monte Carlo simulation: -- Generating stochastic paths for interest rates, exchange rates, and stock prices; -- Numerical valuation of derivative instruments; consider the risk-neutral pricing equation ... ]]></description>
			<content:encoded><![CDATA[<p>In finance, there are two major applications of the Monte Carlo simulation:</p>

<p>-- Generating stochastic paths for interest rates, exchange rates, and stock prices;<br />
-- Numerical valuation of derivative instruments;</p>

<p>consider the risk-neutral pricing equation for security S with payoff <img src="http://l.wordpress.com/latex.php?latex=H%28T%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="H(T)" style="vertical-align:-20%;" class="tex" alt="H(T)" /></p>

<p><img src="http://l.wordpress.com/latex.php?latex=S%28t%29%3DE_t%5B%5Cexp%28-%5Cint%5ET_tr%28s%2C%5Comega%29ds%29H%28T%2C%5Comega%29%5D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="S(t)=E_t[\exp(-\int^T_tr(s,\omega)ds)H(T,\omega)]" style="vertical-align:-20%;" class="tex" alt="S(t)=E_t[\exp(-\int^T_tr(s,\omega)ds)H(T,\omega)]" /></p>

<p>we can generate a large number of equally probable sample paths then the security value at time <img src="http://l.wordpress.com/latex.php?latex=t%3D0&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="t=0" style="vertical-align:-20%;" class="tex" alt="t=0" /> can be approximated as</p>

<p><img src="http://l.wordpress.com/latex.php?latex=S%280%29%3D%5Cfrac%7B1%7D%7BN%7D%5Csum%5EN_%7Bn%3D1%7D%5Cexp%28-%5Cint%5ET_tr%28s%2C%5Comega_n%29ds%29H%28T%2C%5Comega_n%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="S(0)=\frac{1}{N}\sum^N_{n=1}\exp(-\int^T_tr(s,\omega_n)ds)H(T,\omega_n)" style="vertical-align:-20%;" class="tex" alt="S(0)=\frac{1}{N}\sum^N_{n=1}\exp(-\int^T_tr(s,\omega_n)ds)H(T,\omega_n)" />.</p>

<p>For rate process <img src="http://l.wordpress.com/latex.php?latex=dr%3D%5Cmu%28r%2Ct%29dt%2B%5Csigma%28r%2Ct%29dW%28t%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="dr=\mu(r,t)dt+\sigma(r,t)dW(t)" style="vertical-align:-20%;" class="tex" alt="dr=\mu(r,t)dt+\sigma(r,t)dW(t)" />, path generation is straightforward:</p>

<p><img src="http://l.wordpress.com/latex.php?latex=r%28t_%7Bn%2B1%7D%29%3Dr%28t_n%29%2B%5Cmu%28r%28t_n%29%2Ct_n%29%5Ccdot%5CDelta%20t%2B%5Csigma%28r%28t_n%29%2Ct_n%29%5Ccdot%5Csqrt%7B%5CDelta%20t%7D%5Ccdot%20%5Cvarepsilon_%7Bn%2B1%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="r(t_{n+1})=r(t_n)+\mu(r(t_n),t_n)\cdot\Delta t+\sigma(r(t_n),t_n)\cdot\sqrt{\Delta t}\cdot \varepsilon_{n+1}" style="vertical-align:-20%;" class="tex" alt="r(t_{n+1})=r(t_n)+\mu(r(t_n),t_n)\cdot\Delta t+\sigma(r(t_n),t_n)\cdot\sqrt{\Delta t}\cdot \varepsilon_{n+1}" /></p>

<p>where <img src="http://l.wordpress.com/latex.php?latex=%5Cvarepsilon_n&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\varepsilon_n" style="vertical-align:-20%;" class="tex" alt="\varepsilon_n" /> is standard normal random number. For stochastic driver like this</p>

<p><img src="http://l.wordpress.com/latex.php?latex=r%28t%29%3DF%28%5Cvarphi%28t%29%2Bu%28t%29%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="r(t)=F(\varphi(t)+u(t))" style="vertical-align:-20%;" class="tex" alt="r(t)=F(\varphi(t)+u(t))" /><br />
<img src="http://l.wordpress.com/latex.php?latex=du%3D%5Cmu%28t%2Cu%29dt%2B%5Csigma%28t%2Cu%29dW%28t%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="du=\mu(t,u)dt+\sigma(t,u)dW(t)" style="vertical-align:-20%;" class="tex" alt="du=\mu(t,u)dt+\sigma(t,u)dW(t)" /></p>

<p>we sample <img src="http://l.wordpress.com/latex.php?latex=u%28t%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="u(t)" style="vertical-align:-20%;" class="tex" alt="u(t)" /> first and set <img src="http://l.wordpress.com/latex.php?latex=r_n%28t%29%3DF%28%5Cvarphi%28t%29%2Bu_n%28t%29%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="r_n(t)=F(\varphi(t)+u_n(t))" style="vertical-align:-20%;" class="tex" alt="r_n(t)=F(\varphi(t)+u_n(t))" />.</p>

<p>When compared to the lattice valuation method, Monte Carlo approach has two important advantages:</p>

<p>-- Monte Carlo method easily handles path dependent instruments<br />
-- Monte Carlo method is well suited to be used with multi-factor models</p>

<p>however, it also has drawbacks:</p>

<p>-- It is ill-suited for pricing interest rate derivatives with embedded exercise rights<br />
-- It converges to true value very slowly<br />
-- Longer rates along short rate paths cannot be implied from the paths<br />
-- It does not give the same result twice since its value is random<br />
-- Generating sample paths for high dimensional problems has considerable practical difficulties, which is an important feature of interest rate derivatives<br />
-- MC has the opposite problem of recombining lattice: the full future is inaccessible</p>

<p>To implement Monte Carlo method, first we generate uniform random numbers then use methods such as Box-Muller or inverse transform to get normal random numbers. Because of the perfect foresight problem, zero coupon bonds prices cannot be computed without the use of either analytic pricing formulas or interest rate trees. Although MC method will not be used to calibrate the model, MC paths should be calibrated on the top of model calibration to make sure they price zero coupon bonds correctly. Instead of the continuous sampling as described above, we could first discretize the underlying continuous short rate process by means of a short rate lattice and then sample the lattice instead, which is call discrete sampling. The advantage of discrete sampling is that longer yields are readily available at each MC node and we also have a complete view of the stochastic future at each node on a path.</p>

<p>We want to decease the error of MC method as much as possible. One way is to increase the number of paths, the other is variance reduction. There are also two ways to do variance reduction: control variates and improving the sampling quality. The latter's advantage is it does not depend on characteristics of each instrument as the first one does. To improve the sampling quality, there is no other way better than improving the quality of Brownian motion sampled. Since we generate uniform random numbers first, then we can apply an approach called stratified sampling to improve their uniformity by</p>

<p><img src="http://l.wordpress.com/latex.php?latex=u%5E%2A_n%3D%5Cfrac%7Bn-1%7D%7BN%7D%2B%5Cfrac%7Bu_n%7D%7BN%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="u^*_n=\frac{n-1}{N}+\frac{u_n}{N}" style="vertical-align:-20%;" class="tex" alt="u^*_n=\frac{n-1}{N}+\frac{u_n}{N}" /></p>

<p>Based on the fact that Brownian motion is perfectly symmetric about zero, we generate Wiener sample paths which are symmetric to each other w.r.t. zero:</p>

<p><img src="http://l.wordpress.com/latex.php?latex=W_%7B-n%7D%28t%29%3D-W_n%28t%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="W_{-n}(t)=-W_n(t)" style="vertical-align:-20%;" class="tex" alt="W_{-n}(t)=-W_n(t)" /><br />
<img src="http://l.wordpress.com/latex.php?latex=V_%7BMC%7D%3D%5Cfrac%7B1%7D%7BN%7D%5Csum%5EN_%7Bn%3D1%7D%5Cfrac%7BV%28W_n%29%2BV%28W_%7B-n%7D%29%7D%7B2%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="V_{MC}=\frac{1}{N}\sum^N_{n=1}\frac{V(W_n)+V(W_{-n})}{2}" style="vertical-align:-20%;" class="tex" alt="V_{MC}=\frac{1}{N}\sum^N_{n=1}\frac{V(W_n)+V(W_{-n})}{2}" /></p>

<p>and this is called antithetic sampling. After the sample paths <img src="http://l.wordpress.com/latex.php?latex=W_n%28t_i%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="W_n(t_i)" style="vertical-align:-20%;" class="tex" alt="W_n(t_i)" /> are generated, one can make a simple adjustment called moment matching to ensure the correct mean and standard deviation:</p>

<p><img src="http://l.wordpress.com/latex.php?latex=W%5E%2A_n%28t_i%29%3D%5Csqrt%7Bt_i%7D%5Ccdot%5Cfrac%7BW_n%28t_i%29-M_i%7D%7BS_i%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="W^*_n(t_i)=\sqrt{t_i}\cdot\frac{W_n(t_i)-M_i}{S_i}" style="vertical-align:-20%;" class="tex" alt="W^*_n(t_i)=\sqrt{t_i}\cdot\frac{W_n(t_i)-M_i}{S_i}" /></p>

<p>where <img src="http://l.wordpress.com/latex.php?latex=M_i&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="M_i" style="vertical-align:-20%;" class="tex" alt="M_i" /> and <img src="http://l.wordpress.com/latex.php?latex=S_i&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="S_i" style="vertical-align:-20%;" class="tex" alt="S_i" /> are the actual mean and standard deviation of the paths <img src="http://l.wordpress.com/latex.php?latex=W_n&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="W_n" style="vertical-align:-20%;" class="tex" alt="W_n" /> at time <img src="http://l.wordpress.com/latex.php?latex=t_i&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="t_i" style="vertical-align:-20%;" class="tex" alt="t_i" />.</p>
]]></content:encoded>
			<wfw:commentRss>http://ichentao.com/blog/archives/2150/feed/</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>Fixed Income Modeling Review 6</title>
		<link>http://ichentao.com/blog/archives/2140/</link>
		<comments>http://ichentao.com/blog/archives/2140/#comments</comments>
		<pubDate>Sat, 04 Feb 2012 01:27:38 +0000</pubDate>
		<dc:creator>Tao</dc:creator>
				<category><![CDATA[技术]]></category>
		<category><![CDATA[Fixed Income]]></category>
		<category><![CDATA[Mathematical Finance]]></category>

		<guid isPermaLink="false">http://ichentao.com/?p=2140</guid>
		<description><![CDATA[To apply the theoretical models, let us first see how to use interest rate trees to price and calibrate. As mentioned in Review 5, markovian property of short rate models is needed to implement recombining lattice. Usually, a short rate model can be writ... ]]></description>
			<content:encoded><![CDATA[<p>To apply the theoretical models, let us first see how to use interest rate trees to price and calibrate. As mentioned in Review 5, markovian property of short rate models is needed to implement recombining lattice.</p>

<p>Usually, a short rate model can be written as</p>

<p><img src="http://l.wordpress.com/latex.php?latex=r%28t%29%3DF%28u%28t%29%2B%5Cvarphi%28t%29%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="r(t)=F(u(t)+\varphi(t))" style="vertical-align:-20%;" class="tex" alt="r(t)=F(u(t)+\varphi(t))" /><br />
<img src="http://l.wordpress.com/latex.php?latex=du%28t%29%3D-au%28t%29dt%2B%5Csigma%20dW%28t%29.&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="du(t)=-au(t)dt+\sigma dW(t)." style="vertical-align:-20%;" class="tex" alt="du(t)=-au(t)dt+\sigma dW(t)." /></p>

<p>Hence, we first discretize the process <img src="http://l.wordpress.com/latex.php?latex=u%28t%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="u(t)" style="vertical-align:-20%;" class="tex" alt="u(t)" />. Typically, people use binomial tree and trinomial tree but for mean reverting model, binomial tree cannot be used. The time step <img src="http://l.wordpress.com/latex.php?latex=%5CDelta%20t&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\Delta t" style="vertical-align:-20%;" class="tex" alt="\Delta t" /> is arbitrary, the state step <img src="http://l.wordpress.com/latex.php?latex=%5CDelta%20u&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\Delta u" style="vertical-align:-20%;" class="tex" alt="\Delta u" /> and branch probabilities <img src="http://l.wordpress.com/latex.php?latex=p_u&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="p_u" style="vertical-align:-20%;" class="tex" alt="p_u" />, <img src="http://l.wordpress.com/latex.php?latex=p_m&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="p_m" style="vertical-align:-20%;" class="tex" alt="p_m" /> and <img src="http://l.wordpress.com/latex.php?latex=p_d&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="p_d" style="vertical-align:-20%;" class="tex" alt="p_d" /> should be chosen in such a way that the discrete dynamics has the first few moments matched. Since we could solve</p>

<p><img src="http://l.wordpress.com/latex.php?latex=u%28t%2B%5CDelta%20t%29%3De%5E%7B-a%5CDelta%20t%7Du%28t%29%2B%5Csigma%5Cint%5E%7Bt%2B%5CDelta%20t%7D_te%5E%7B-a%5CDelta%20t%7DdW%28s%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="u(t+\Delta t)=e^{-a\Delta t}u(t)+\sigma\int^{t+\Delta t}_te^{-a\Delta t}dW(s)" style="vertical-align:-20%;" class="tex" alt="u(t+\Delta t)=e^{-a\Delta t}u(t)+\sigma\int^{t+\Delta t}_te^{-a\Delta t}dW(s)" /></p>

<p>then by moment-matching, we can set</p>

<p><img src="http://l.wordpress.com/latex.php?latex=%5CDelta%20u%3D%5Csqrt%7B3V%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\Delta u=\sqrt{3V}" style="vertical-align:-20%;" class="tex" alt="\Delta u=\sqrt{3V}" />, <img src="http://l.wordpress.com/latex.php?latex=V%3D%5Csigma%5E2%5Cfrac%7B1-%5Cexp%28-2a%5CDelta%20t%29%7D%7B2a%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="V=\sigma^2\frac{1-\exp(-2a\Delta t)}{2a}" style="vertical-align:-20%;" class="tex" alt="V=\sigma^2\frac{1-\exp(-2a\Delta t)}{2a}" /></p>

<p>and we select the middle branching node at next step to be the closest one to the mean of the continuous process:</p>

<p><img src="http://l.wordpress.com/latex.php?latex=k%3D%5Cint%5Bje%5E%7B-a%5CDelta%20t%7D%2B%5Cfrac%7B1%7D%7B2%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="k=\int[je^{-a\Delta t}+\frac{1}{2}" style="vertical-align:-20%;" class="tex" alt="k=\int[je^{-a\Delta t}+\frac{1}{2}" /></p>

<p>Branching probabilities are determined as</p>

<p><img src="http://l.wordpress.com/latex.php?latex=p_u%3D%5Cfrac%7B1%7D%7B6%7D%2B%5Cfrac%7B1%7D%7B2%7D%28%5Cbeta%5E2%2B%5Cbeta%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="p_u=\frac{1}{6}+\frac{1}{2}(\beta^2+\beta)" style="vertical-align:-20%;" class="tex" alt="p_u=\frac{1}{6}+\frac{1}{2}(\beta^2+\beta)" /><br />
<img src="http://l.wordpress.com/latex.php?latex=p_m%3D%5Cfrac%7B2%7D%7B3%7D-%5Cbeta%5E2&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="p_m=\frac{2}{3}-\beta^2" style="vertical-align:-20%;" class="tex" alt="p_m=\frac{2}{3}-\beta^2" /><br />
<img src="http://l.wordpress.com/latex.php?latex=p_d%3D%5Cfrac%7B1%7D%7B6%7D%2B%5Cfrac%7B1%7D%7B2%7D%28%5Cbeta%5E2-%5Cbeta%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="p_d=\frac{1}{6}+\frac{1}{2}(\beta^2-\beta)" style="vertical-align:-20%;" class="tex" alt="p_d=\frac{1}{6}+\frac{1}{2}(\beta^2-\beta)" /><br />
<img src="http://l.wordpress.com/latex.php?latex=beta%3Dje%5E%7B-a%5CDelta%20t%7D-k&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="beta=je^{-a\Delta t}-k" style="vertical-align:-20%;" class="tex" alt="beta=je^{-a\Delta t}-k" /></p>

<p>Once tree for the stochastic driver <img src="http://l.wordpress.com/latex.php?latex=u%28t%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="u(t)" style="vertical-align:-20%;" class="tex" alt="u(t)" /> is built, we need to convert it to the short rate tree according to the functions:</p>

<p>normal: <img src="http://l.wordpress.com/latex.php?latex=r%28i%2Cj%29%3Dr_0%28i%29%2Bu%28i%2Cj%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="r(i,j)=r_0(i)+u(i,j)" style="vertical-align:-20%;" class="tex" alt="r(i,j)=r_0(i)+u(i,j)" /><br />
lognormal: <img src="http://l.wordpress.com/latex.php?latex=r%28i%2Cj%29%3Dr_0%28i%29%5Ccdot%5Cexp%28u%28i%2Cj%29%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="r(i,j)=r_0(i)\cdot\exp(u(i,j))" style="vertical-align:-20%;" class="tex" alt="r(i,j)=r_0(i)\cdot\exp(u(i,j))" /></p>

<p>where mean level vector <img src="http://l.wordpress.com/latex.php?latex=r_0%28i%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="r_0(i)" style="vertical-align:-20%;" class="tex" alt="r_0(i)" /> is calibrated to the term structure of zero coupon bond prices.</p>

<p>Since trees represent discretization of the rate dynamics, then within each time step the evolution of short rate is not described by tree and thus needs to be specified. Usually we assume short rate is constant within each period. To price derivatives, we do backwards</p>

<p><img src="http://l.wordpress.com/latex.php?latex=S%28t%29%3DE%5EB_t%5B%5Cexp%28-%5Cint_t%5E%7Bt%2B%5CDelta%20t%7Dr%28s%29ds%29%28S%28t%2B%5CDelta%20t%29%2BCF%28t%2B%5CDelta%20t%29%29%5D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="S(t)=E^B_t[\exp(-\int_t^{t+\Delta t}r(s)ds)(S(t+\Delta t)+CF(t+\Delta t))]" style="vertical-align:-20%;" class="tex" alt="S(t)=E^B_t[\exp(-\int_t^{t+\Delta t}r(s)ds)(S(t+\Delta t)+CF(t+\Delta t))]" /></p>

<p>where <img src="http://l.wordpress.com/latex.php?latex=CF%28t%2B%5CDelta%20t%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="CF(t+\Delta t)" style="vertical-align:-20%;" class="tex" alt="CF(t+\Delta t)" /> is the cash flow at time <img src="http://l.wordpress.com/latex.php?latex=t%2B%5CDelta%20t&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="t+\Delta t" style="vertical-align:-20%;" class="tex" alt="t+\Delta t" />. And now the remaining problem is how to calibrate to get mean level vector <img src="http://l.wordpress.com/latex.php?latex=r_0%28i%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="r_0(i)" style="vertical-align:-20%;" class="tex" alt="r_0(i)" />. Denote prices of zero coupon bonds with maturities <img src="http://l.wordpress.com/latex.php?latex=t_i%3Di%5CDelta%20t&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="t_i=i\Delta t" style="vertical-align:-20%;" class="tex" alt="t_i=i\Delta t" /> by <img src="http://l.wordpress.com/latex.php?latex=P%28i%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="P(i)" style="vertical-align:-20%;" class="tex" alt="P(i)" />, then apply the following iterative search process forward:</p>

<p>(1) Set <img src="http://l.wordpress.com/latex.php?latex=r_0%280%29%3D-%5Cln%28P%281%29%29%2F%5CDelta%20t&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="r_0(0)=-\ln(P(1))/\Delta t" style="vertical-align:-20%;" class="tex" alt="r_0(0)=-\ln(P(1))/\Delta t" /><br />
(2) Search for the value of <img src="http://l.wordpress.com/latex.php?latex=r_0%281%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="r_0(1)" style="vertical-align:-20%;" class="tex" alt="r_0(1)" /> such that lattice price of the zero coupon bond maturing at <img src="http://l.wordpress.com/latex.php?latex=t_2&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="t_2" style="vertical-align:-20%;" class="tex" alt="t_2" /> equals <img src="http://l.wordpress.com/latex.php?latex=P%282%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="P(2)" style="vertical-align:-20%;" class="tex" alt="P(2)" /><br />
(3) ...<br />
(4) Assuming <img src="http://l.wordpress.com/latex.php?latex=r_0%280%29%2C%5Cldots%2Cr_0%28i-1%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="r_0(0),\ldots,r_0(i-1)" style="vertical-align:-20%;" class="tex" alt="r_0(0),\ldots,r_0(i-1)" /> are done, search for the value of <img src="http://l.wordpress.com/latex.php?latex=r_0%28i%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="r_0(i)" style="vertical-align:-20%;" class="tex" alt="r_0(i)" /> to fit <img src="http://l.wordpress.com/latex.php?latex=P%28i%2B1%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="P(i+1)" style="vertical-align:-20%;" class="tex" alt="P(i+1)" /><br />
(5) ...</p>

<p>We need to keep in mind that we always need to do numerical calibration for <img src="http://l.wordpress.com/latex.php?latex=r_0%28i%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="r_0(i)" style="vertical-align:-20%;" class="tex" alt="r_0(i)" /> even when there are analytic formulas for the no-arbitrage drift because the use of trees alters the continuous dynamics. To be more specific about step (2) and step (4), we introduce Arrow-Debreu Security which pays 1$ once the node (i,j) is reached and pays nothing otherwise. Thus in trees we have</p>

<p><img src="http://l.wordpress.com/latex.php?latex=AD%28i%2B1%2Cj%29%3D%5Csum_k%20AD%28i%2Ck%29%5Cexp%28-r%28i%2Ck%29%5CDelta%20t%29P_%7Bi%2Ck%5Crightarrow%20i%2B1%2Cj%7D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="AD(i+1,j)=\sum_k AD(i,k)\exp(-r(i,k)\Delta t)P_{i,k\rightarrow i+1,j}" style="vertical-align:-20%;" class="tex" alt="AD(i+1,j)=\sum_k AD(i,k)\exp(-r(i,k)\Delta t)P_{i,k\rightarrow i+1,j}" /></p>

<p>One very useful property of Arrow-Debreu prices is that</p>

<p><img src="http://l.wordpress.com/latex.php?latex=P%280%2Ct_i%29%3D%5Csum_j%20AD%28i%2Cj%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="P(0,t_i)=\sum_j AD(i,j)" style="vertical-align:-20%;" class="tex" alt="P(0,t_i)=\sum_j AD(i,j)" /></p>

<p>Recombining interest rate trees have a serious shortcoming that the rate history is totally lost and thus it is incompatible with any path dependent instrument. To overcome this problem, one can use the tower law and record the cash flows at the lattice nodes not when they are paid but rather when they are certain and express them with appropriate discounting.</p>
]]></content:encoded>
			<wfw:commentRss>http://ichentao.com/blog/archives/2140/feed/</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>Fixed Income Modeling Review 5</title>
		<link>http://ichentao.com/blog/archives/2127/</link>
		<comments>http://ichentao.com/blog/archives/2127/#comments</comments>
		<pubDate>Fri, 03 Feb 2012 20:19:21 +0000</pubDate>
		<dc:creator>Tao</dc:creator>
				<category><![CDATA[技术]]></category>
		<category><![CDATA[科学人]]></category>
		<category><![CDATA[Fixed Income]]></category>
		<category><![CDATA[Mathematical Finance]]></category>

		<guid isPermaLink="false">http://ichentao.com/?p=2127</guid>
		<description><![CDATA[This topic is mainly about modeling of short rate which plays a central role in the theory and practice because zero coupon bonds can be priced as and so could discount factor. Assume there is only one source of uncertainty represented by a Brownian moti... ]]></description>
			<content:encoded><![CDATA[<p>This topic is mainly about modeling of short rate <img src="http://l.wordpress.com/latex.php?latex=r%28t%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="r(t)" style="vertical-align:-20%;" class="tex" alt="r(t)" /> which plays a central role in the theory and practice because zero coupon bonds can be priced as</p>

<p><img src="http://l.wordpress.com/latex.php?latex=P%28t%2CT%29%3DE_t%5B%5Cexp%28-%5Cint%5ET_tr%28s%29ds%29%5D&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="P(t,T)=E_t[\exp(-\int^T_tr(s)ds)]" style="vertical-align:-20%;" class="tex" alt="P(t,T)=E_t[\exp(-\int^T_tr(s)ds)]" /></p>

<p>and so could discount factor.</p>

<p>Assume there is only one source of uncertainty represented by a Brownian motion <img src="http://l.wordpress.com/latex.php?latex=W%28t%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="W(t)" style="vertical-align:-20%;" class="tex" alt="W(t)" /> under the risk-neutral measure <img src="http://l.wordpress.com/latex.php?latex=Q&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="Q" style="vertical-align:-20%;" class="tex" alt="Q" />. We want the stochastic process to be markovian which is referred to as path independence and hence allows for the use of recombining lattice. We don't want to use non-recombining lattice cause the number of nodes growing exponentially w.r.t. the number of time steps and thus is not practical. There are several important things we need to think about when choosing a model:</p>

<p>-- Does the dynamics guarantee positive rates?<br />
-- What distribution does the dynamics imply for the short rate?<br />
-- Can the model be calibrated to fit today's structure of zero coupon bond prices, how easy, and does the calibration stable?<br />
-- Are the zero bond prices computable by means of analytic formulas?<br />
-- Are there analytic pricing formulas available for call/put bond option?<br />
-- Is the model mean-reverting?<br />
-- How do the volatility structures implied by the model look like?<br />
-- How suited is the model for building recombining lattices?</p>

<p>If the parameters of a model are constant like <img src="http://l.wordpress.com/latex.php?latex=dr%3Da%28b-r%28t%29%29dt%2B%5Csigma%20dW&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="dr=a(b-r(t))dt+\sigma dW" style="vertical-align:-20%;" class="tex" alt="dr=a(b-r(t))dt+\sigma dW" /> (Vasicek) then it is called equilibrium model. It's simple but could not be calibrated to fit the current term structure of rates. If we make the parameter of the drift term time dependent then we can overcome this problem and make the model no-arbitrage.</p>

<p>Based on empirical observation, people find that interest rates (unlike stock prices) tend to be pulled back to some mean level. So we want our model to be mean reverting; we also want to have the analytical formulas for zero coupon bond, call/put bond options, etc., and hence we require the model be affine which means the drift term is linear to short rate and square of the volatility is also linear to short rate. Regarding the volatility, typically we have two kinds of model for volatility: normal and lognormal. It is very hard to say which one is better, normal assumption results in analytic tractability and lognormal assumption leads to positive rates. Also, normal dynamics suit low rate environments and for lognormal, the other way around.</p>

<p>Several important models:</p>

<p>Ho-Lee: <img src="http://l.wordpress.com/latex.php?latex=dr%3D%5Ctheta%28t%29dt%2B%5Csigma%20dW&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="dr=\theta(t)dt+\sigma dW" style="vertical-align:-20%;" class="tex" alt="dr=\theta(t)dt+\sigma dW" /><br />
Hull-White: <img src="http://l.wordpress.com/latex.php?latex=dr%3Da%28%5Ctheta%28t%29-r%29dt%2B%5Csigma%20dW&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="dr=a(\theta(t)-r)dt+\sigma dW" style="vertical-align:-20%;" class="tex" alt="dr=a(\theta(t)-r)dt+\sigma dW" /><br />
CIR model: <img src="http://l.wordpress.com/latex.php?latex=dr%3Da%28%5Ctheta%28t%29-r%29dt%2B%5Csigma%5Csqrt%7Br%7DdW&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="dr=a(\theta(t)-r)dt+\sigma\sqrt{r}dW" style="vertical-align:-20%;" class="tex" alt="dr=a(\theta(t)-r)dt+\sigma\sqrt{r}dW" /><br />
BDT model: <img src="http://l.wordpress.com/latex.php?latex=d%5Cln%28r%29%3D%28%5Ctheta%28t%29%2B%28%5Csigma%27%28t%29%2F%5Csigma%28t%29%29%5Cln%28r%29%29dt%2B%5Csigma%28t%29dW%28t%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="d\ln(r)=(\theta(t)+(\sigma'(t)/\sigma(t))\ln(r))dt+\sigma(t)dW(t)" style="vertical-align:-20%;" class="tex" alt="d\ln(r)=(\theta(t)+(\sigma'(t)/\sigma(t))\ln(r))dt+\sigma(t)dW(t)" /></p>

<p>Among these models, Hull-White has a positive probability to generate negative rates but we can cut that and it does not affect pricing very much. CIR model is the only model with full analytic tractability which also has positive rates. While doing calibration, <img src="http://l.wordpress.com/latex.php?latex=%5Ctheta%28t%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\theta(t)" style="vertical-align:-20%;" class="tex" alt="\theta(t)" /> term in CIR must be solved for numerically. Under lognormal assumption, the total volatility is proportional to the rate level and thus it grows/declines with the rates, which produces rate path distribution skewed upward. Theoretically, we could deduce such equation under lognormal model:</p>

<p><img src="http://l.wordpress.com/latex.php?latex=E_0%5BB%28%5CDelta%20t%29%5D%3D%5Cinfty&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="E_0[B(\Delta t)]=\infty" style="vertical-align:-20%;" class="tex" alt="E_0[B(\Delta t)]=\infty" /></p>

<p>therefore we can make infinite money in an arbitrary small time from $1 initial investment, however in numerical implementations of a short rate model on a tree this problem does not appear as the number of states is finite.</p>

<p>In order to make our models better, we could add a deterministic shift:</p>

<p><img src="http://l.wordpress.com/latex.php?latex=r%28t%29%3Dr_0%28t%29%2B%5Cvarphi%28t%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="r(t)=r_0(t)+\varphi(t)" style="vertical-align:-20%;" class="tex" alt="r(t)=r_0(t)+\varphi(t)" /><br />
<img src="http://l.wordpress.com/latex.php?latex=dr_0%28t%29%3Da%28%5Ctheta%28t%29-r_0%29dt%2B%5Csigma%5Csqrt%7Br_0%7DdW%28t%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="dr_0(t)=a(\theta(t)-r_0)dt+\sigma\sqrt{r_0}dW(t)" style="vertical-align:-20%;" class="tex" alt="dr_0(t)=a(\theta(t)-r_0)dt+\sigma\sqrt{r_0}dW(t)" />.</p>

<p>This kind of extension could preserve the tractability but not the positiveness of rates. Shifted CIR model could perfectly fit the market price of zero coupon bond and according to this:</p>

<p>when rates are low: <img src="http://l.wordpress.com/latex.php?latex=%5Csigma%28t%29%28r%28t%29%2B%5Calpha%29%5Capprox%5Csigma%28t%29%5Calpha&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\sigma(t)(r(t)+\alpha)\approx\sigma(t)\alpha" style="vertical-align:-20%;" class="tex" alt="\sigma(t)(r(t)+\alpha)\approx\sigma(t)\alpha" /><br />
when rates are high: <img src="http://l.wordpress.com/latex.php?latex=%5Csigma%28t%29%28r%28t%29%2B%5Calpha%29%5Capprox%5Csigma%28t%29r%28t%29&#038;bg=FFFFFF&#038;fg=000000&#038;s=0" title="\sigma(t)(r(t)+\alpha)\approx\sigma(t)r(t)" style="vertical-align:-20%;" class="tex" alt="\sigma(t)(r(t)+\alpha)\approx\sigma(t)r(t)" /></p>

<p>we know that shifted lognormal models capture the phenomenon that rates follow close-to-normal dynamics in low rate environments and close-to-lognormal one in high rate environments.</p>
]]></content:encoded>
			<wfw:commentRss>http://ichentao.com/blog/archives/2127/feed/</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
		<item>
		<title>The Best of 2011 - Uzak</title>
		<link>http://ichentao.com/blog/archives/2113/</link>
		<comments>http://ichentao.com/blog/archives/2113/#comments</comments>
		<pubDate>Sat, 31 Dec 2011 22:21:48 +0000</pubDate>
		<dc:creator>Tao</dc:creator>
				<category><![CDATA[电影评论]]></category>
		<category><![CDATA[Nuri Ceylan]]></category>

		<guid isPermaLink="false">http://ichentao.com/?p=2113</guid>
		<description><![CDATA[首先, 我觉得中文译名不够好, 对于本片, 从涵盖主题的角度来说, "远方"并不如英文译名"Distant"来的精确. 之所以喜欢这部电影, 还是在于导演对"疏离"这个元素的表达和描述. 我一直相信, 一部真... ]]></description>
			<content:encoded><![CDATA[<p><img alt="" src="https://blufiles.storage.live.com/y1psFcedTQQ_y9_Oejqtc8f1qFlrvjH0cdPqr05OeypMyMk_9OfLNtT99n5ccMlFuSjPULE7xTtddU/p461194043.jpg" title="Uzak" class="alignnone" width="300" height="438" /></p>

<p>首先, 我觉得中文译名不够好, 对于本片, 从涵盖主题的角度来说, "远方"并不如英文译名"Distant"来的精确.</p>

<p>之所以喜欢这部电影, 还是在于导演对"疏离"这个元素的表达和描述. 我一直相信, 一部真正的好电影, 一定是重述和探索人们日常生活中的真实的. 就像摄影, 我们都喜欢那些捕捉到有趣, 动人瞬间的片子. 不是说商业剧情片不好, 我认为它们也有趣, 但代表了另一种乐趣, 如同摄影师用 Photoshop 把羊头天衣无缝的接到牛身上, 一种表象的刺激, 我们看到都会笑, 但笑完之后, 几乎没留下任何回味的余地. 回到"疏离感"上来, 这是美国影评家非常喜欢的一个词, 尽管他们的导演们似乎并不太喜欢向观众表现它. 其实这是多么普遍的一种体会和感受呢, 它来自于工业化, 来自于社会的进步, 来自于人内心的骄傲等等一切可以设置障碍的有形或无形的变化. 比如说已故的安东尼奥尼, "疏离感"深入到他电影的骨髓里, 建筑, 色彩, 主角的行为, 心理活动, 一切外在或内在无不表现着这一矛盾.</p>

<p>当代欧洲有很多出色的导演, Almodóvar 关注女性; Nanni Moretti 关注男性的气质危机; Fatih Akin 关注不同文化间的冲击; 唯一深入研究超越种族, 文化的人类的本质性矛盾的, 就我目前了解, 好像只有锡兰 (Nuri Ceylan) 一个人. 在他的作品中, 他总是把这种矛盾投射到各个方面: 伊斯坦布尔与小镇; 男人与女人; "文化"与"文盲"; 努力工作与游手好闲, 在这些差别中他展现和探讨人与人, 人与生活的距离. 我们喜欢欧洲, 不仅因为他们有很长的历史, 很多的雕塑与建筑, 有多元的文化, 有摇滚乐与足球, 还因为他们关注人, 有人文精神. 锡兰正是这种人文关怀的表达者.</p>

<p>锡兰本身还是一位出色的摄影师, "诗意的镜头语言"是塔科夫斯基的标志, 锡兰的影片有同样的特点, 这也是人们称他为"土耳其的塔科夫斯基"的原因, 将电影拍成照片, 又将照片连接成电影, 所以他的影片总在诗意中透出那种"愉悦的忧伤", 这种忧伤符合他的身份: 一个扎根于土耳其的导演, 他的国家由弱到强, 又盛极而衰. 他的作品偶尔提及历史 (在极少的情况下), 而聚焦于当下, 历史变迁, 形成了当代土耳其人的气质, 锡兰在作品中将这一切回归于普通的土耳其人, 做了平淡而又饱满的表达.</p>
]]></content:encoded>
			<wfw:commentRss>http://ichentao.com/blog/archives/2113/feed/</wfw:commentRss>
		<slash:comments>2</slash:comments>
		</item>
		<item>
		<title>The Best of 2011 (候选)</title>
		<link>http://ichentao.com/blog/archives/2084/</link>
		<comments>http://ichentao.com/blog/archives/2084/#comments</comments>
		<pubDate>Wed, 28 Dec 2011 07:41:14 +0000</pubDate>
		<dc:creator>Tao</dc:creator>
				<category><![CDATA[影音]]></category>
		<category><![CDATA[Fatih Akin]]></category>
		<category><![CDATA[Ferzan Özpetek]]></category>
		<category><![CDATA[Nuri Ceylan]]></category>

		<guid isPermaLink="false">http://ichentao.com/?p=2084</guid>
		<description><![CDATA[又到年末, 照例会挑出一部当年最喜欢的电影, "新浪潮"的老板说我每次都在碟子卖完后才去他那问, 不光买碟, 看电影也一样, 今年豆瓣上打5颗星的电影最近也是2010年的, 后知后觉到我自己都不... ]]></description>
			<content:encoded><![CDATA[<p><img alt="" src="https://blufiles.storage.live.com/y1pmcsm2obOHxUH6MBdBnw4wV7FJisjJ5T-iZuZ9onFYWTtaYJH9GhTi8Wj4VtLVGwI7u7ocN8ByW0/Akin_solino14.jpg" title="Fatih Akin" class="alignnone" width="500" height="341" /></p>

<p><img alt="" src="https://blufiles.storage.live.com/y1p3Q0X_5QY8l5NCSsSn6Y8hWhsJ_YB4hL4bnbA6UVgJeR8kZzKDsMnrTINDvRGBpMLqNHg3zj8MVk/Kasim2008-2_1273090093.jpg" title="Ferzan Ozpetek" class="alignnone" width="500" height="331" /></p>

<p><img alt="" src="https://blufiles.storage.live.com/y1p2azYaZ8iegYGf0s7r1O1FlBlywuKNOYAxY-I1mQb2Qo-Qs91FAXdb0UEcgjNj2rxniDRohkviq4/slider-anatolia-nuri.jpg" title="Ceylan" class="alignnone" width="500" height="282" /></p>

<p>又到年末, 照例会挑出一部当年最喜欢的电影, "新浪潮"的老板说我每次都在碟子卖完后才去他那问, 不光买碟, 看电影也一样, 今年豆瓣上打5颗星的电影最近也是2010年的, 后知后觉到我自己都不能忍了.</p>

<p>这年的年初看过几部电影, 学期中途总是很忙无暇看片, 这个寒假算是过了瘾, 差不多是以每天不止一部的速度在看, 其间不停的感慨, 不停的流泪, 但观影速度太快不能静下心来思考, 明知是对好电影的糟蹋, 又忍不住想看下一部是怎样的故事.</p>

<p>缘起 Fatih Akin 的 <a href="http://movie.douban.com/subject/3104216/">Soul Kitchen</a>, 最近这段时间看了相当多的土耳其导演作品, 还有 Ferzan Ozpetek 和 Nuri Ceylan, 打了5颗星的有: <a href="http://movie.douban.com/subject/1442128/">Solino</a>, <a href="http://movie.douban.com/subject/1298331/">Im Juli.</a>, <a href="http://movie.douban.com/subject/1309012/">Head-on</a>, <a href="http://movie.douban.com/subject/1307921/">Uzak</a>, <a href="http://movie.douban.com/subject/1300827/">Mayıs Sıkıntısı</a>, <a href="http://movie.douban.com/subject/3637237/">Mine vaganti</a>, <a href="http://movie.douban.com/subject/1449431/">Hamam</a>. 他们的风格互不相同, 这些不同又是缘于个人出身和经历的差别. Fatih Akin 是出生在汉堡的土耳其后裔, 他的电影始终关注着德国的土耳其移民, 移民在异乡的位置以及他们与故乡的关系; Ferzan Ozpetek 本人是一位同性恋(如果我没有理解错 an openly gay director 这句话的话), 他的作品不时聚焦在同性恋的身上, 描写社会对同性恋的排斥与宽容, 与同性恋角色与生俱来的问题是他们对自身的探索, 种种原因导致 Ozpetek 的电影中总是隐藏着秘密, 带有悬疑的味道; Nuri Ceylan 是很多人追捧的一位导演, 因为他的镜头有着与阿巴斯或者塔可夫斯基类似的风格, 他本人同时是摄影家, 在他的掌控下, 影片或像一张张流动的胶片, 或者就是一张定格的摄影作品, 也由于如此, 他的电影极其擅长通过人物面部的特写去展现心理活动.</p>

<p>除了这些土耳其制作, 09年的法国电影<a href="http://movie.douban.com/subject/3285622/">非法入境 (Welcome)</a> 讲述伊拉克难民试图偷渡到英国的故事, 开放性的结尾画龙点睛, 温暖动人; 法斯宾德的<a href="http://movie.douban.com/subject/1294129/">罗拉 (Lola)</a> 是他"女性三部曲"中我唯一没看过的, 赞叹导演的天才: 通过一个柔弱的女性角色爆发出巨大的能量, 完成对社会的强烈反讽; <a href="http://movie.douban.com/subject/3667102/">乐与路</a>则是一个简单, 带有明显日本风格的主流爱情片, 偏偏人就容易被这样的电影催泪; <a href="http://movie.douban.com/subject/4843480/">Inside Job</a> 是关于金融危机的纪录片, 我们不能指望艺术工作者能够完全客观的看待这次危机: 有些剪辑显得刻意, 有些评论过于激进, 但总体而言, 很清晰明了的解释了金融危机的来龙去脉; <a href="http://movie.douban.com/subject/3040026/">靖大爷和他的老主顾们</a>也是一部纪录片, 完成于02年, 片长不足一小时, 年过八旬的老理发匠四处奔走为年纪相仿的老主顾们理发, 他见证了政权, 城市和人生的变迁, 送走一位位主顾, 也不得不随着北京城的发展而告别一个个传统, 八十多岁, 知天命, 按他自己的话说到了该死的年龄, 每天乐观而平静的活着.</p>

<p>接下来的几天会好好回味一下这些片子, 跨年之前挑出年度最喜爱的作品.</p>
]]></content:encoded>
			<wfw:commentRss>http://ichentao.com/blog/archives/2084/feed/</wfw:commentRss>
		<slash:comments>6</slash:comments>
		</item>
	</channel>
</rss>

