{"categories":["Stochastic Calculus"],"contentHtml":"<h5>Content Review</h5><p>Problem Set 5</p><p><math><mi>𝑌</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>=</mo><msup><mi>𝑒</mi><mrow><mi>𝑡</mi><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></mrow></msup></math></p><p>Find the differential</p><p><math><mi>𝑑</mi><mi>𝑌</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></math></p><p>Recognize that this is some function that looks like <math><mi>𝑓</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑎</mi><mo>,</mo><mi>𝑏</mi><mo stretchy=\"false\">)</mo></math>, like</p><p><math><mi>𝑓</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑎</mi><mo>,</mo><mi>𝑏</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>=</mo><msup><mi>𝑒</mi><mrow><mi>𝑎</mi><mi>𝑏</mi></mrow></msup></math></p><p>where <math><mi>𝑎</mi></math> is a deterministic component and <math><mi>𝑏</mi></math> is a stochastic component.</p><p>Then solve your partials</p><p><math><msub><mi>𝑓</mi><mi>𝑎</mi></msub><mo>=</mo><mi>𝑏</mi><msup><mi>𝑒</mi><mrow><mi>𝑎</mi><mi>𝑏</mi></mrow></msup></math> <math><msub><mi>𝑓</mi><mi>𝑏</mi></msub><mo>=</mo><mi>𝑎</mi><msup><mi>𝑒</mi><mrow><mi>𝑎</mi><mi>𝑏</mi></mrow></msup></math></p><p><math><msub><mi>𝑓</mi><mrow><mi>𝑏</mi><mi>𝑏</mi></mrow></msub><mo>=</mo><msup><mi>𝑎</mi><mn>2</mn></msup><msup><mi>𝑒</mi><mrow><mi>𝑎</mi><mi>𝑏</mi></mrow></msup></math></p><p>then using the Ito formula, this is</p><p><math><mi>𝑑</mi><mi>𝑌</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>=</mo><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><msup><mi>𝑒</mi><mrow><mi>𝑡</mi><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></mrow></msup><mi>𝑑</mi><mi>𝑡</mi><mo>+</mo><mi>𝑡</mi><msup><mi>𝑒</mi><mrow><mi>𝑡</mi><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></mrow></msup><mi>𝑑</mi><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>𝑡</mi><mn>2</mn></msup><msup><mi>𝑒</mi><mrow><mi>𝑡</mi><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></mrow></msup><mi>𝑑</mi><mi>𝑡</mi></math></p><p>Question on why we’re able to get away with this? <math><mi>𝔼</mi><mo form=\"prefix\" stretchy=\"false\">[</mo><msup><mi>𝑒</mi><mrow><mi>𝜎</mi><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></mrow></msup><mo form=\"postfix\" stretchy=\"false\">]</mo><mo>=</mo><msup><mi>𝑒</mi><mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>𝜎</mi><mn>2</mn></msup><mi>𝑡</mi></mrow></msup></math></p><p>This is because of a result, confirmed in the result, that</p><p><math><mi>𝑋</mi><mo>∼</mo><mi>𝑁</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑎</mi><mo>,</mo><msup><mi>𝑏</mi><mn>2</mn></msup><mo stretchy=\"false\">)</mo></math> implies <math><mi>𝔼</mi><mo form=\"prefix\" stretchy=\"false\">[</mo><msup><mi>𝑒</mi><mrow><mi>𝑢</mi><mi>𝑋</mi></mrow></msup><mo stretchy=\"false\">]</mo></math></p><p>Which is the moment generating function with respect to the dummy variable <math><mi>𝑢</mi></math>, is</p><p><math><mo form=\"infix\" lspace=\"0em\">=</mo><msup><mi>𝑒</mi><mrow><mi>𝑎</mi><mi>𝑢</mi><mo lspace=\"0em\" rspace=\"0em\">+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>𝑏</mi><mn>2</mn></msup><msup><mi>𝑢</mi><mn>2</mn></msup></mrow></msup></math></p><p>We’re using this off-label. So if we plug it into this we can do this. What if <math><mi>𝜎</mi></math> is adapted to time? Can we use this result? Not necessarily. Being an adapted process means that you’re measurable with respect to filtration at that time. But if <math><mi mathvariant=\"normal\">Δ</mi><mi>𝑡</mi></math> is a random process, then this is a totally different problem.</p><p><math><mfrac><mrow><mi>𝑑</mi><mi>𝑆</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></mrow><mrow><mi>𝑆</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></mrow></mfrac><mo>=</mo><mi>𝜇</mi><mi>𝑑</mi><mi>𝑡</mi><mo>+</mo><mi>𝜎</mi><mi>𝑑</mi><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></math></p><p>Can you take the integral from both sides? You can, but it causes some issues.</p><p><math><msubsup><mo lspace=\"0em\">∫</mo><mn>0</mn><mi>𝑡</mi></msubsup><mfrac><mrow><mi>𝑑</mi><mi>𝑆</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo stretchy=\"false\">)</mo></mrow><mrow><mi>𝑆</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo stretchy=\"false\">)</mo></mrow></mfrac><mo>=</mo><mi>𝛼</mi><mi>𝑡</mi><mo>+</mo><mi>𝜎</mi><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></math></p><p>Doesn’t really give you a closed form approximation.</p><p>From the derivation of Black Scholes</p><p>Expectation of the discount factor, <math><msub><mi>𝑆</mi><mi>𝑡</mi></msub><mo>−</mo><mi>𝐾</mi></math>, then you have an indicator function.</p><p>Indicator is transformed, where does that came from</p><p><math><mi>𝔼</mi><mrow><mo>[</mo><mrow><msup><mi>𝑒</mi><mrow><mo>−</mo><mi>𝑟</mi><mi>𝑇</mi></mrow></msup><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑆</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑇</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>−</mo><mi>𝐾</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><msub><mi>𝕀</mi><mrow><mo>{</mo><mrow><mi>𝑑</mi><mo lspace=\"0em\" rspace=\"0em\">∗</mo><mo lspace=\"0em\" rspace=\"0em\">&lt;</mo><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑇</mi><mo stretchy=\"false\">)</mo></mrow><mo>}</mo></mrow></msub></mrow><mo>]</mo></mrow></math></p><p>We can substitute</p><p>What is random? <math><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></math>. What it’s distribution? normal over real line.</p><p>Then take the infinite integral.</p><p>Question from sample midterm: <math><mi>𝐼</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>=</mo><msubsup><mo form=\"prefix\" lspace=\"0em\">∫</mo><mn>0</mn><mi>𝑡</mi></msubsup><mi mathvariant=\"normal\">Δ</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mi>𝑑</mi><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo stretchy=\"false\">)</mo></math></p><p>Find <math><mi>𝔼</mi><mo form=\"prefix\" stretchy=\"false\">[</mo><msup><mi>𝐼</mi><mn>4</mn></msup><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo stretchy=\"false\">]</mo></math></p><p>And <math><mi mathvariant=\"normal\">Δ</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo stretchy=\"false\">)</mo></math> is <math><mn>3</mn><mi>𝑢</mi></math></p><p>There are a couple ways to do this! There’s the way with the hint, and other ways. Any time you’re doing an expected value of a non-well-defined process, like an Ito process raised to a power, look at the Ito decomposition instead of the process as a whole.</p><p><math><mi>𝑑</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><msup><mi>𝐼</mi><mn>4</mn></msup><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>=</mo><mn>4</mn><msup><mi>𝐼</mi><mn>3</mn></msup><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mi>𝑑</mi><mi>𝐼</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>+</mo><mn>6</mn><msup><mi>𝐼</mi><mn>2</mn></msup><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑑</mi><mi>𝐼</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><msup><mo stretchy=\"false\">)</mo><mn>2</mn></msup></math></p><p>And then</p><p><math><mi>𝑑</mi><mi>𝐼</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>=</mo><mn>3</mn><mi>𝑡</mi><mi>𝑑</mi><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></math></p><p>based on what we already know…</p><p><math><mo stretchy=\"false\">(</mo><mi>𝑑</mi><mi>𝐼</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><msup><mo form=\"postfix\" stretchy=\"false\">)</mo><mn>2</mn></msup><mo>=</mo><mn>9</mn><msup><mi>𝑡</mi><mn>2</mn></msup><mi>𝑑</mi><mi>𝑡</mi></math> based on what we know</p><p><math><msup><mi>𝐼</mi><mn>4</mn></msup><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>=</mo><msubsup><mo form=\"prefix\" lspace=\"0em\">∫</mo><mn>0</mn><mi>𝑡</mi></msubsup><mn>4</mn><msup><mi>𝐼</mi><mn>3</mn></msup><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mi mathvariant=\"normal\">Δ</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mi>𝑑</mi><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>+</mo><msubsup><mo form=\"prefix\" lspace=\"0em\">∫</mo><mn>0</mn><mi>𝑡</mi></msubsup><mn>6</mn><msup><mi>𝐼</mi><mn>2</mn></msup><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo form=\"prefix\" stretchy=\"false\">(</mo><mn>9</mn><msup><mi>𝑢</mi><mn>2</mn></msup><mo form=\"postfix\" stretchy=\"false\">)</mo><mi>𝑑</mi><mi>𝑢</mi></math></p><p>When I take the expectation, I can see that an Ito integral with respect to a martingale is a martingale, so the first term is gone.</p><p><math><mi>𝔼</mi><mo form=\"prefix\" stretchy=\"false\">[</mo><msup><mi>𝐼</mi><mn>4</mn></msup><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo form=\"postfix\" stretchy=\"false\">]</mo><mo>=</mo><msubsup><mo form=\"prefix\" lspace=\"0em\">∫</mo><mn>0</mn><mi>𝑡</mi></msubsup><mn>54</mn><msup><mi>𝑢</mi><mn>2</mn></msup><mi>𝔼</mi><mo form=\"prefix\" stretchy=\"false\">[</mo><msup><mi>𝐼</mi><mn>2</mn></msup><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo form=\"postfix\" stretchy=\"false\">]</mo><mi>𝑑</mi><mi>𝑢</mi></math></p><p>Could you be more clever about this?</p><p>This is an Ito integral of a deterministic integrand, so it has to be normally distributed. The fourth moment of any normal variable is three times the variance.</p><p>This would give you the answer. But if it was like in the homework problem where you want to do <math><msup><mi>𝐼</mi><mn>6</mn></msup></math>, you would have to do it in the Ito decomposition way.</p><p>An Ito decomposition results in A nonrandom initial constant, a Riemann integral <math><mi>𝑑</mi><mi>𝑢</mi></math>, and an Ito integral <math><mi>𝑑</mi><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></math>.</p><p>Understand that this formula is</p><p>Generalized geometric Brownian motion.</p><p>If you recognize that if you have a process that can be rewritten as e^Ito process, do we need to do our partials? If you can reliably get the correct answer without doing that, you can, but it’s better to do the partials.</p><p>Itos’ formula:</p><p><math><mi>𝑓</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo>,</mo><mi>𝑋</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>=</mo><mi>𝑓</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mn>0</mn><mo>,</mo><mi>𝑋</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mn>0</mn><mo form=\"postfix\" stretchy=\"false\">)</mo><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>+</mo><msubsup><mo form=\"prefix\" lspace=\"0em\">∫</mo><mn>0</mn><mi>𝑡</mi></msubsup><msub><mi>𝑓</mi><mi>𝑎</mi></msub><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo>,</mo><mi>𝑋</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo form=\"postfix\" stretchy=\"false\">)</mo><mi>𝑑</mi><mi>𝑢</mi><mo>+</mo><msubsup><mo form=\"prefix\" lspace=\"0em\">∫</mo><mn>0</mn><mi>𝑡</mi></msubsup><msub><mi>𝑓</mi><mi>𝑏</mi></msub><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo>,</mo><mi>𝑋</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo form=\"postfix\" stretchy=\"false\">)</mo><mi>𝑑</mi><mi>𝑋</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msubsup><mo form=\"prefix\">∫</mo><mn>0</mn><mi>𝑡</mi></msubsup><msub><mi>𝑓</mi><mrow><mi>𝑏</mi><mi>𝑏</mi></mrow></msub><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo>,</mo><mi>𝑋</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo form=\"postfix\" stretchy=\"false\">)</mo><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑑</mi><mi>𝑋</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><msup><mo stretchy=\"false\">)</mo><mn>2</mn></msup></math></p><p><math><mi>𝑋</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></math> is an Ito process</p><p><math><mi>𝑢</mi></math> is your index</p><p>Brownian motion is definitely an Ito process. In order to be an Ito process, you need to be able to be expressed as the nonrandom constant + Riemann + stochastic.</p><p>Quadratic variation of a process</p><p><math><mi>𝑋</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>=</mo><mn>3</mn><mo>−</mo><msubsup><mo form=\"prefix\" lspace=\"0em\">∫</mo><mn>0</mn><mi>𝑡</mi></msubsup><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mi>𝑑</mi><mi>𝑢</mi><mo>+</mo><msubsup><mo form=\"prefix\" lspace=\"0em\">∫</mo><mn>0</mn><mi>𝑡</mi></msubsup><mi>𝑢</mi><mi>𝑑</mi><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo stretchy=\"false\">)</mo></math></p><p>Quickest way to find this is recognizing that this is an Ito process.</p><p>Quadratic variation of an Ito process is</p><p><math><msubsup><mo lspace=\"0em\">∫</mo><mn>0</mn><mi>𝑡</mi></msubsup><msup><mi mathvariant=\"normal\">Δ</mi><mn>2</mn></msup><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mi>𝑑</mi><mi>𝑢</mi></math> REMEMBER THIS</p><p>How do we get the expected value of an Ito integral squared? The Ito isometry.</p><p>If I have the expression</p><p>If I tell you <math><mi>𝑢</mi><mo>=</mo><mn>2</mn></math> and <math><mi>𝑡</mi><mo>=</mo><mn>4</mn></math>,</p><p><math><msubsup><mo lspace=\"0em\">∫</mo><mn>0</mn><mi>𝑡</mi></msubsup><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><msup><mi>𝑢</mi><mn>2</mn></msup><mspace width=\"0.1667em\"></mspace><mi>𝑑</mi><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo stretchy=\"false\">)</mo></math></p><p>THAT IS NOT AN ADAPTED PROCESS.</p><h3>Week 9</h3><h4>General Probability Theory</h4><p><math><mi mathvariant=\"normal\">Ω</mi></math> is the <strong>sample space</strong> of all possible outcomes <math><mi>𝜔</mi></math>.</p><p>A <strong><math><mi>𝜎</mi></math>-algebra</strong> is the set of all sets generated from a starting set of outcomes, their complements, and the unions of all complements. Conceptually, it represents an addressable subset in a potentially uncountably infinite sample space. The <math><mi>𝜎</mi></math>-algebra is represented by <math><mi>ℱ︀</mi></math>.</p><p>The <strong>probability measure</strong> <math><mi>ℙ</mi></math> distributes values of probability on every set in the <math><mi>𝜎</mi></math>-algebra <math><mi>ℱ︀</mi></math>, which must add up to 1 across <math><mi mathvariant=\"normal\">Ω</mi></math>.</p><p>The <strong>Borel set</strong> is the <math><mi>𝜎</mi></math>-algebra generated by the set of all closed intervals on the real number line. Conceptually, it represents an addressable subset of real numbers.</p><p>A <strong>random variable</strong> is a mapping from outcomes <math><mi>𝜔</mi></math> to a Borel subset. This could be a single number or an interval.</p><p>The <strong>distribution measure</strong> of a random variable <math><mi>𝑋</mi></math> is the probability measure that <math><mi>𝑋</mi></math> is in each Borel subset. This is computed by finding the <em>preimage</em> of a Borel subset, which is an event in the sample space, and finding the probability measure of that event.</p><p>Two different random variables could have the same distribution measure, and the same random variable could have two different distribution measures under two different probability measures.</p><p>The <strong>expectation</strong> of a random variable is the average result of the random variable. This can be computed in many ways. For a discrete random variable, you can directly compute the result of each outcome mapped through the random variable, multiply them by the probability of the outcome, and sum them. For a continuous random variable, you must use an integral.</p><p><math><mi>𝔼</mi><mo form=\"prefix\" stretchy=\"false\">[</mo><mi>𝑋</mi><mo form=\"postfix\" stretchy=\"false\">]</mo><mo>=</mo><msub><mo form=\"prefix\" lspace=\"0em\">∫</mo><mi mathvariant=\"normal\">Ω</mi></msub><mi>𝑋</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝜔</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mi>𝑑</mi><mi>ℙ</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝜔</mi><mo stretchy=\"false\">)</mo></math></p><p>(Definition 1.3.3)</p><p>Expectations have the property of linearity.</p><h4>Information and Conditioning</h4><p>We can model the information we have regarding an uncertain outcome with <math><mi>𝜎</mi></math>-algebras.</p><p>If we know which set of outcomes that the outcome must belong to based on the current information, we can generate a <math><mi>𝜎</mi></math>-algebra from that set of outcomes.</p><p>A <strong>filtration</strong> is a sequence of such <math><mi>𝜎</mi></math>-algebras that represent the information available over time until time <math><mi>𝑡</mi></math>.</p><p>Filtrations always increase in information over time.</p><p>The filtration at time <math><mi>𝑡</mi></math> represents all the information at time <math><mi>𝑡</mi></math>, and it is used for conditioning. We can say that a random variable is <strong>measurable</strong> with respect to a filtration if the information in the filtration can be used to determine the value of the random variable.</p><p>For any two sets, <strong>independence</strong> is generally defined by</p><p><math><mi>ℙ</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝐴</mi><mo>∩</mo><mi>𝐵</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>=</mo><mi>ℙ</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝐴</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo lspace=\"0.2222222222222222em\" rspace=\"0.2222222222222222em\">⋅</mo><mi>ℙ</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝐵</mi><mo stretchy=\"false\">)</mo></math></p><p>Two <math><mi>𝜎</mi></math>-algebras are independent if this is true for all sets in each <math><mi>𝜎</mi></math>-algebra.</p><p>Two random variables are independent if the <math><mi>𝜎</mi></math>-algebras generated from them are independent.</p><p>Any functions on independent random variables will also result in independent random variables.</p><p>Not gonna bother with joint densities and crap.</p><p>A <strong>conditional expectation</strong> allows us to use information from a filtration to estimate the expectation of a random variable without it being completely measurable. In shorthand, we might use a random variable as a condition to indicate that the filtration is the <math><mi>𝜎</mi></math>-algebra generated by the random variable.</p><p>Conditional expectations have the following properties (Theorem 2.3.2)</p><ul><li>Linearity</li><li>Taking out what is known: if a random variable <math><mi>𝑋</mi></math> is measurable with respect to the condition, it can be taken out. <math><mi>𝔼</mi><mo form=\"prefix\" stretchy=\"false\">[</mo><mi>𝑋</mi><mi>𝑌</mi><mo lspace=\"0em\" rspace=\"0em\">|</mo><mi>𝒢︀</mi><mo form=\"postfix\" stretchy=\"false\">]</mo><mo>=</mo><mi>𝑋</mi><mi>𝔼</mi><mo form=\"prefix\" stretchy=\"false\">[</mo><mi>𝑌</mi><mo lspace=\"0em\" rspace=\"0em\">|</mo><mi>𝒢︀</mi><mo stretchy=\"false\">]</mo></math></li><li>Iterated conditioning: if we have a <math><mi>𝜎</mi></math>-algebra <math><mi>𝐻</mi></math> with <em>less</em> information, we can subsume the <math><mi>𝜎</mi></math>-algebra <math><mi>𝐺</mi></math> with <em>more</em> information into it. <math><mi>𝔼</mi><mo form=\"prefix\" stretchy=\"false\">[</mo><mi>𝔼</mi><mo form=\"prefix\" stretchy=\"false\">[</mo><mi>𝑋</mi><mo lspace=\"0em\" rspace=\"0em\">|</mo><mi>𝒢︀</mi><mo form=\"postfix\" stretchy=\"false\">]</mo><mo lspace=\"0em\" rspace=\"0em\">|</mo><mi>ℋ︀</mi><mo form=\"postfix\" stretchy=\"false\">]</mo><mo>=</mo><mi>𝔼</mi><mo form=\"prefix\" stretchy=\"false\">[</mo><mi>𝑋</mi><mo lspace=\"0em\" rspace=\"0em\">|</mo><mi>ℋ︀</mi><mo stretchy=\"false\">]</mo></math></li><li>Independence: if a random variable <math><mi>𝑋</mi></math> is independent of the condition, the condition is irrelevant. <math><mi>𝔼</mi><mo form=\"prefix\" stretchy=\"false\">[</mo><mi>𝑋</mi><mo lspace=\"0em\" rspace=\"0em\">|</mo><mi>𝒢︀</mi><mo form=\"postfix\" stretchy=\"false\">]</mo><mo>=</mo><mi>𝔼</mi><mo form=\"prefix\" stretchy=\"false\">[</mo><mi>𝑋</mi><mo stretchy=\"false\">]</mo></math></li></ul><p>A <strong>stochastic process</strong> is a collection of random variables indexed at time <math><mi>𝑡</mi></math>.</p><p>An <strong>adapted stochastic process</strong> is a stochastic process where each random variable is measurable by the filtration at time <math><mi>𝑡</mi></math>. Conceptually, it is a process which only depends on information available at each time <math><mi>𝑡</mi></math>.</p><p>A <strong>martingale</strong> is an adapted stochastic process which is expected to remain the same value at all times. Formally, (Definition 2.3.5)</p><p><math><mi>𝔼</mi><mo form=\"prefix\" stretchy=\"false\">[</mo><mi>𝑀</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo lspace=\"0em\" rspace=\"0em\">|</mo><mi>ℱ︀</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑠</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo form=\"postfix\" stretchy=\"false\">]</mo><mo>=</mo><mi>𝑀</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑠</mi><mo stretchy=\"false\">)</mo></math> for all <math><mn>0</mn><mo>≤</mo><mi>𝑠</mi><mo>≤</mo><mi>𝑡</mi><mo>≤</mo><mi>𝑇</mi></math></p><p>Submartingales and supermartingales observe this property for <math><mo lspace=\"0em\" rspace=\"0em\">≥</mo></math> and <math><mo lspace=\"0em\" rspace=\"0em\">≤</mo></math>, respectively.</p><p>A <strong>Markov</strong> is an adapted stochastic process where we can determine the expectation of the process using only the previous value. In other words, all of the information (filtration) up to time <math><mi>𝑠</mi></math> is encoded in the value. Formally (2.3.6)</p><p><math><mi>𝔼</mi><mo form=\"prefix\" stretchy=\"false\">[</mo><mi>𝑓</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑋</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo form=\"postfix\" stretchy=\"false\">)</mo><mo lspace=\"0em\" rspace=\"0em\">|</mo><mi>ℱ︀</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑠</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo form=\"postfix\" stretchy=\"false\">]</mo><mo>=</mo><mi>𝑔</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑋</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑠</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo stretchy=\"false\">)</mo></math></p><h4>Brownian Motion</h4><p>The <strong>symmetric random walk</strong> is an adapted stochastic process which moves up and down by 1 by the flip of a coin at each time step. It has the property that each of its increments are independent. The expectation of the increment is 0, and its variance is <math><mi mathvariant=\"normal\">Δ</mi><mi>𝑡</mi></math>. The walk has the martingale property, which can be shown by splitting the current value into an increment, which is independent, and the previous value, which is measurable.</p><p>The <strong>scaled symmetric random walk</strong> approximates the Brownian motion.</p><p><math><msup><mi>𝑊</mi><mrow><mo stretchy=\"false\">(</mo><mi>𝑛</mi><mo stretchy=\"false\">)</mo></mrow></msup><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>=</mo><mfrac><mn>1</mn><msqrt><mi>𝑛</mi></msqrt></mfrac><msub><mi>𝑀</mi><mrow><mi>𝑛</mi><mi>𝑡</mi></mrow></msub></math></p><p>as <math><mi>𝑛</mi></math> approaches <math><mi>∞</mi></math>.</p><p>I’m going to skip a lot of stuff relating to random walks, assuming that it won’t be relevant.</p><p><strong>Brownian motion</strong> <math><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></math> has the following properties:</p><ul><li>Independent increments</li><li>Normally distributed increments</li><li>Expectation of increments is 0</li><li>Variance of increments is <math><mi mathvariant=\"normal\">Δ</mi><mi>𝑡</mi></math></li><li>Martingale</li></ul><p>The <strong>first-order variation</strong> of a function is the definite integral from time 0 to <math><mi>𝑇</mi></math> of the absolute value of its derivative. It can also be computed by partitioning time <math><mi>𝑇</mi></math> into partition <math><mi mathvariant=\"normal\">Π</mi></math>. As <math><mo form=\"infix\" lspace=\"0em\">∥</mo><mi mathvariant=\"normal\">Π</mi><mo form=\"infix\" rspace=\"0em\">∥</mo></math> (the norm of <math><mi mathvariant=\"normal\">Π</mi></math>, aka size of its largest increment) approaches 0, we’re interested in summing the absolute value of each of these increments. Essentially we’re measuring how much change there is in the function.</p><p>The <strong>quadratic variation</strong> is computed in the same way, but with the increments being squared. Formally (3.4.1)</p><p><math><mo stretchy=\"false\">[</mo><mi>𝑓</mi><mo>,</mo><mi>𝑓</mi><mo form=\"postfix\" stretchy=\"false\">]</mo><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑇</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>=</mo><msub><mo form=\"prefix\" lspace=\"0em\" rspace=\"0.16666666666666666em\">lim</mo><mrow><mo form=\"infix\" lspace=\"0em\" rspace=\"0em\">∥</mo><mi mathvariant=\"normal\">Π</mi><mo lspace=\"0em\" rspace=\"0em\">∥</mo><mo lspace=\"0em\" rspace=\"0em\" stretchy=\"false\">→</mo><mn>0</mn></mrow></msub><msubsup><mo form=\"prefix\" lspace=\"0em\" rspace=\"0em\">∑</mo><mrow><mi>𝑗</mi><mo lspace=\"0em\" rspace=\"0em\">=</mo><mn>0</mn></mrow><mrow><mi>𝑛</mi><mo lspace=\"0em\" rspace=\"0em\">−</mo><mn>1</mn></mrow></msubsup><mo form=\"prefix\" stretchy=\"false\">[</mo><mi>𝑓</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><msub><mi>𝑡</mi><mrow><mi>𝑗</mi><mo lspace=\"0em\" rspace=\"0em\">+</mo><mn>1</mn></mrow></msub><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>−</mo><mi>𝑓</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><msub><mi>𝑡</mi><mi>𝑗</mi></msub><mo form=\"postfix\" stretchy=\"false\">)</mo><msup><mo stretchy=\"false\">]</mo><mn>2</mn></msup></math></p><p>Typically, quadratic variation is 0 for any differentiable function. The quadratic variation of Brownian motion <math><mo stretchy=\"false\">[</mo><mi>𝑊</mi><mo>,</mo><mi>𝑊</mi><mo form=\"postfix\" stretchy=\"false\">]</mo><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑇</mi><mo stretchy=\"false\">)</mo></math> is <math><mi>𝑇</mi></math>. Conceptually, we can understand this as Brownian motion accumulating quadratic variation at rate one per unit time.</p><p>The differential of time <math><mi>𝑑</mi><mi>𝑡</mi></math> multiplied by any other differential is 0. The differential of Brownian motion <math><mi>𝑑</mi><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></math> multiplied by itself is <math><mi>𝑑</mi><mi>𝑡</mi></math>. This relates to quadratic variation because</p><p><math><mi>𝑑</mi><mo form=\"prefix\" stretchy=\"false\">[</mo><mi>𝑊</mi><mo>,</mo><mi>𝑊</mi><mo form=\"postfix\" stretchy=\"false\">]</mo><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>=</mo><mi>𝑑</mi><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mi>𝑑</mi><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>=</mo><mi>𝑑</mi><mi>𝑡</mi></math></p><p><strong>Geometric Brownian motion</strong> is an application of Brownian motion to stock prices. The process is described by</p><p><math><mi>𝑆</mi><mrow><mo>(</mo><mi>𝑡</mi><mo>)</mo></mrow><mo>=</mo><mi>𝑆</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo form=\"prefix\" lspace=\"0.16666666666666666em\" rspace=\"0em\">exp</mo><mrow><mo>{</mo><mrow><mi>𝜎</mi><mi>𝑊</mi><mrow><mo>(</mo><mi>𝑡</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mrow><mi>𝛼</mi><mo>−</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>𝜎</mi><mn>2</mn></msup></mrow><mo>)</mo></mrow><mi>𝑡</mi></mrow><mo>}</mo></mrow></math></p><p>And log returns in this process are described by</p><p><math><mo lspace=\"0em\" rspace=\"0.16666666666666666em\">log</mo><mfrac><mrow><mi>𝑆</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><msub><mi>𝑡</mi><mrow><mi>𝑗</mi><mo lspace=\"0em\" rspace=\"0em\">+</mo><mn>1</mn></mrow></msub><mo stretchy=\"false\">)</mo></mrow><mrow><mi>𝑆</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><msub><mi>𝑡</mi><mi>𝑗</mi></msub><mo stretchy=\"false\">)</mo></mrow></mfrac><mo>=</mo><mi>𝜎</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><msub><mi>𝑡</mi><mrow><mi>𝑗</mi><mo lspace=\"0em\" rspace=\"0em\">+</mo><mn>1</mn></mrow></msub><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>−</mo><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><msub><mi>𝑡</mi><mi>𝑗</mi></msub><mo form=\"postfix\" stretchy=\"false\">)</mo><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>+</mo><mrow><mo>(</mo><mrow><mi>𝛼</mi><mo>−</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>𝜎</mi><mn>2</mn></msup></mrow><mo>)</mo></mrow><mo form=\"prefix\" stretchy=\"false\">(</mo><msub><mi>𝑡</mi><mrow><mi>𝑗</mi><mo lspace=\"0em\" rspace=\"0em\">+</mo><mn>1</mn></mrow></msub><mo>−</mo><msub><mi>𝑡</mi><mi>𝑗</mi></msub><mo stretchy=\"false\">)</mo></math></p><p>Skipping first passage time distribution, reflection principle, maximum to date</p><h4>Stochastic Calculus</h4><p>The purpose of stochastic calculus is to integrate an adapted stochastic process on a Brownian motion.</p><p><math><msubsup><mo lspace=\"0em\">∫</mo><mn>0</mn><mi>𝑇</mi></msubsup><mi mathvariant=\"normal\">Δ</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mi>𝑑</mi><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></math> This <math><mi mathvariant=\"normal\">Δ</mi></math> describes our position in an asset at time <math><mi>𝑡</mi></math>, which is dependent on the change in price. This integral is the value of the portfolio at time <math><mi>𝑇</mi></math>. The change in price is a stochastic process based on Brownian motion, which means that we cannot integrate with a Riemann/Lebesgue integral.</p><p>In order to construct the <strong>Itô integral</strong>, we add up increments of <math><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></math> and multiply them by the stock position <math><mi mathvariant=\"normal\">Δ</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></math> at that time. The discrete analogy is that of <math><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo stretchy=\"false\">)</mo></math> as stock price and the time intervals being trading dates. The value of the portfolio is therefore the sum of the product of: the change in stock price each day, and the position held in stock on that day. Taking this discrete analogy continuously creates the integral.</p><p>The Itô integral is described as follows:</p><p><math><mi>𝐼</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑡</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mo>=</mo><msubsup><mo form=\"prefix\" lspace=\"0em\">∫</mo><mn>0</mn><mi>𝑡</mi></msubsup><mi mathvariant=\"normal\">Δ</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo form=\"postfix\" stretchy=\"false\">)</mo><mi>𝑑</mi><mi>𝑊</mi><mo form=\"prefix\" stretchy=\"false\">(</mo><mi>𝑢</mi><mo stretchy=\"false\">)</mo></math></p><p>The Itô integral is a martingale. This is proved here.</p><h3>Final Methods</h3>","contentSourceFormat":"typst","contentTypst":"==== Content Review\n<content-review>\nProblem Set 5\n\n$Y\\(t\\)= e^(t W\\(t\\))$\n\nFind the differential\n\n$d Y\\(t\\)$\n\nRecognize that this is some function that looks like $f\\(a\\,b\\)$, like\n\n$f\\(a\\,b\\)= e^(a b)$\n\nwhere $a$ is a deterministic component and $b$ is a stochastic component.\n\nThen solve your partials\n\n$f_a = b e^(a b)$ $f_b = a e^(a b)$\n\n$f_(b b) = a^2 e^(a b)$\n\nthen using the Ito formula, this is\n\n$d Y\\(t\\)= W\\(t\\)e^(t W\\(t\\)) d t + t e^(t W\\(t\\)) d W\\(t\\)+ 1 / 2 t^2 e^(t W\\(t\\)) d t$\n\nQuestion on why we're able to get away with this? $bb(E)\\[e^(sigma W\\(t\\))\\]= e^(1 / 2 sigma^2 t)$\n\nThis is because of a result, confirmed in the result, that\n\n$X tilde.op N\\(a\\,b^2\\)$ implies $bb(E)\\[e^(u X)\\]$\n\nWhich is the moment generating function with respect to the dummy variable $u$, is\n\n$= e^(a u + 1 / 2 b^2 u^2)$\n\nWe're using this off-label. So if we plug it into this we can do this. What if $sigma$ is adapted to time? Can we use this result? Not necessarily. Being an adapted process means that you're measurable with respect to filtration at that time. But if $Delta t$ is a random process, then this is a totally different problem.\n\n$frac(d S\\(t\\), S\\(t\\)) = mu d t + sigma d W\\(t\\)$\n\nCan you take the integral from both sides? You can, but it causes some issues.\n\n$integral_0^t frac(d S\\(u\\), S\\(u\\)) = alpha t + sigma W\\(t\\)$\n\nDoesn't really give you a closed form approximation.\n\nFrom the derivation of Black Scholes\n\nExpectation of the discount factor, $S_t - K$, then you have an indicator function.\n\nIndicator is transformed, where does that came from\n\n$bb(E) [e^(- r T) \\( S \\( T \\) - K \\) bb(I)_({ d * < W\\(T\\)})]$\n\nWe can substitute\n\nWhat is random? $W\\(t\\)$. What it's distribution? normal over real line.\n\nThen take the infinite integral.\n\nQuestion from sample midterm: $I\\(t\\)= integral_0^t Delta\\(u\\)d W\\(u\\)$\n\nFind $bb(E)\\[I^4\\(t\\)\\]$\n\nAnd $Delta\\(u\\)$ is $3 u$\n\nThere are a couple ways to do this! There's the way with the hint, and other ways. Any time you're doing an expected value of a non-well-defined process, like an Ito process raised to a power, look at the Ito decomposition instead of the process as a whole.\n\n$d\\(I^4\\(t\\)\\)= 4 I^3\\(t\\)d I\\(t\\)+ 6 I^2\\(t\\)\\(d I\\(t\\)\\)^2$\n\nAnd then\n\n$d I\\(t\\)= 3 t d W\\(t\\)$\n\nbased on what we already know…\n\n$\\(d I\\(t\\)\\)^2= 9 t^2 d t$ based on what we know\n\n$I^4\\(t\\)= integral_0^t 4 I^3\\(u\\)Delta\\(u\\)d W\\(u\\)+ integral_0^t 6 I^2\\(u\\)\\(9 u^2\\)d u$\n\nWhen I take the expectation, I can see that an Ito integral with respect to a martingale is a martingale, so the first term is gone.\n\n$bb(E)\\[I^4\\(t\\)\\]= integral_0^t 54 u^2 bb(E)\\[I^2\\(u\\)\\]d u$\n\nCould you be more clever about this?\n\nThis is an Ito integral of a deterministic integrand, so it has to be normally distributed. The fourth moment of any normal variable is three times the variance.\n\nThis would give you the answer. But if it was like in the homework problem where you want to do $I^6$, you would have to do it in the Ito decomposition way.\n\nAn Ito decomposition results in A nonrandom initial constant, a Riemann integral $d u$, and an Ito integral $d W\\(t\\)$.\n\nUnderstand that this formula is\n\nGeneralized geometric Brownian motion.\n\nIf you recognize that if you have a process that can be rewritten as e^Ito process, do we need to do our partials? If you can reliably get the correct answer without doing that, you can, but it's better to do the partials.\n\nItos' formula:\n\n$f\\(t\\,X\\(t\\)\\)= f\\(0\\,X\\(0\\)\\)+ integral_0^t f_a\\(u\\,X\\(u\\)\\)d u + integral_0^t f_b\\(u\\,X\\(u\\)\\)d X\\(u\\)+ 1 / 2 integral_0^t f_(b b)\\(u\\,X\\(u\\)\\)\\(d X\\(u\\)\\)^2$\n\n$X\\(t\\)$ is an Ito process\n\n$u$ is your index\n\nBrownian motion is definitely an Ito process. In order to be an Ito process, you need to be able to be expressed as the nonrandom constant + Riemann + stochastic.\n\nQuadratic variation of a process\n\n$X\\(t\\)= 3 - integral_0^t W\\(u\\)d u + integral_0^t u d W\\(u\\)$\n\nQuickest way to find this is recognizing that this is an Ito process.\n\nQuadratic variation of an Ito process is\n\n$integral_0^t Delta^2\\(u\\)d u$ REMEMBER THIS\n\nHow do we get the expected value of an Ito integral squared? The Ito isometry.\n\nIf I have the expression\n\nIf I tell you $u = 2$ and $t = 4$,\n\n$integral_0^t W\\(t\\)u^2 thin d W\\(u\\)$\n\nTHAT IS NOT AN ADAPTED PROCESS.\n\n== Week 9\n<week-9>\n=== General Probability Theory\n<general-probability-theory>\n$Omega$ is the #strong[sample space] of all possible outcomes $omega$.\n\nA #strong[$sigma$-algebra] is the set of all sets generated from a starting set of outcomes, their complements, and the unions of all complements. Conceptually, it represents an addressable subset in a potentially uncountably infinite sample space. The $sigma$-algebra is represented by $cal(F)$.\n\nThe #strong[probability measure] $bb(P)$ distributes values of probability on every set in the $sigma$-algebra $cal(F)$, which must add up to 1 across $Omega$.\n\nThe #strong[Borel set] is the $sigma$-algebra generated by the set of all closed intervals on the real number line. Conceptually, it represents an addressable subset of real numbers.\n\nA #strong[random variable] is a mapping from outcomes $omega$ to a Borel subset. This could be a single number or an interval.\n\nThe #strong[distribution measure] of a random variable $X$ is the probability measure that $X$ is in each Borel subset. This is computed by finding the #emph[preimage] of a Borel subset, which is an event in the sample space, and finding the probability measure of that event.\n\nTwo different random variables could have the same distribution measure, and the same random variable could have two different distribution measures under two different probability measures.\n\nThe #strong[expectation] of a random variable is the average result of the random variable. This can be computed in many ways. For a discrete random variable, you can directly compute the result of each outcome mapped through the random variable, multiply them by the probability of the outcome, and sum them. For a continuous random variable, you must use an integral.\n\n$bb(E)\\[X\\]= integral_Omega X\\(omega\\)d bb(P)\\(omega\\)$\n\n\\(Definition 1.3.3)\n\nExpectations have the property of linearity.\n\n=== Information and Conditioning\n<information-and-conditioning>\nWe can model the information we have regarding an uncertain outcome with $sigma$-algebras.\n\nIf we know which set of outcomes that the outcome must belong to based on the current information, we can generate a $sigma$-algebra from that set of outcomes.\n\nA #strong[filtration] is a sequence of such $sigma$-algebras that represent the information available over time until time $t$.\n\nFiltrations always increase in information over time.\n\nThe filtration at time $t$ represents all the information at time $t$, and it is used for conditioning. We can say that a random variable is #strong[measurable] with respect to a filtration if the information in the filtration can be used to determine the value of the random variable.\n\nFor any two sets, #strong[independence] is generally defined by\n\n$bb(P)\\(A inter B\\)= bb(P)\\(A\\)dot.op bb(P)\\(B\\)$\n\nTwo $sigma$-algebras are independent if this is true for all sets in each $sigma$-algebra.\n\nTwo random variables are independent if the $sigma$-algebras generated from them are independent.\n\nAny functions on independent random variables will also result in independent random variables.\n\nNot gonna bother with joint densities and crap.\n\nA #strong[conditional expectation] allows us to use information from a filtration to estimate the expectation of a random variable without it being completely measurable. In shorthand, we might use a random variable as a condition to indicate that the filtration is the $sigma$-algebra generated by the random variable.\n\nConditional expectations have the following properties (Theorem 2.3.2)\n\n- Linearity\n- Taking out what is known: if a random variable $X$ is measurable with respect to the condition, it can be taken out. $bb(E)\\[X Y\\|cal(G)\\]= X bb(E)\\[Y\\|cal(G)\\]$\n- Iterated conditioning: if we have a $sigma$-algebra $H$ with #emph[less] information, we can subsume the $sigma$-algebra $G$ with #emph[more] information into it. $bb(E)\\[bb(E)\\[X\\|cal(G)\\]\\|cal(H)\\]= bb(E)\\[X\\|cal(H)\\]$\n- Independence: if a random variable $X$ is independent of the condition, the condition is irrelevant. $bb(E)\\[X\\|cal(G)\\]= bb(E)\\[X\\]$\n\nA #strong[stochastic process] is a collection of random variables indexed at time $t$.\n\nAn #strong[adapted stochastic process] is a stochastic process where each random variable is measurable by the filtration at time $t$. Conceptually, it is a process which only depends on information available at each time $t$.\n\nA #strong[martingale] is an adapted stochastic process which is expected to remain the same value at all times. Formally, (Definition 2.3.5)\n\n$bb(E)\\[M\\(t\\)\\|cal(F)\\(s\\)\\]= M\\(s\\)$ for all $0 lt.eq s lt.eq t lt.eq T$\n\nSubmartingales and supermartingales observe this property for $gt.eq$ and $lt.eq$, respectively.\n\nA #strong[Markov] is an adapted stochastic process where we can determine the expectation of the process using only the previous value. In other words, all of the information (filtration) up to time $s$ is encoded in the value. Formally (2.3.6)\n\n$bb(E)\\[f\\(X\\(t\\)\\)\\|cal(F)\\(s\\)\\]= g\\(X\\(s\\)\\)$\n\n=== Brownian Motion\n<brownian-motion>\nThe #strong[symmetric random walk] is an adapted stochastic process which moves up and down by 1 by the flip of a coin at each time step. It has the property that each of its increments are independent. The expectation of the increment is 0, and its variance is $Delta t$. The walk has the martingale property, which can be shown by splitting the current value into an increment, which is independent, and the previous value, which is measurable.\n\nThe #strong[scaled symmetric random walk] approximates the Brownian motion.\n\n$W^(\\(n\\))\\(t\\)= 1 / sqrt(n) M_(n t)$\n\nas $n$ approaches $oo$.\n\nI'm going to skip a lot of stuff relating to random walks, assuming that it won't be relevant.\n\n#strong[Brownian motion] $W\\(t\\)$ has the following properties:\n\n- Independent increments\n- Normally distributed increments\n- Expectation of increments is 0\n- Variance of increments is $Delta t$\n- Martingale\n\nThe #strong[first-order variation] of a function is the definite integral from time 0 to $T$ of the absolute value of its derivative. It can also be computed by partitioning time $T$ into partition $Pi$. As $parallel Pi parallel$ (the norm of $Pi$, aka size of its largest increment) approaches 0, we're interested in summing the absolute value of each of these increments. Essentially we're measuring how much change there is in the function.\n\nThe #strong[quadratic variation] is computed in the same way, but with the increments being squared. Formally (3.4.1)\n\n$\\[f\\,f\\]\\(T\\)= lim_(parallel Pi parallel arrow.r 0) sum_(j = 0)^(n - 1)\\[f\\(t_(j + 1)\\)- f\\(t_j\\)\\]^2$\n\nTypically, quadratic variation is 0 for any differentiable function. The quadratic variation of Brownian motion $\\[W\\,W\\]\\(T\\)$ is $T$. Conceptually, we can understand this as Brownian motion accumulating quadratic variation at rate one per unit time.\n\nThe differential of time $d t$ multiplied by any other differential is 0. The differential of Brownian motion $d W\\(t\\)$ multiplied by itself is $d t$. This relates to quadratic variation because\n\n$d\\[W\\,W\\]\\(t\\)= d W\\(t\\)d W\\(t\\)= d t$\n\n#strong[Geometric Brownian motion] is an application of Brownian motion to stock prices. The process is described by\n\n$S(t) = S(0) exp{sigma W(t) + (alpha - 1 / 2 sigma^2) t}$\n\nAnd log returns in this process are described by\n\n$log frac(S\\(t_(j + 1)\\), S\\(t_j\\)) = sigma\\(W\\(t_(j + 1)\\)- W\\(t_j\\)\\)+ (alpha - 1 / 2 sigma^2)\\(t_(j + 1) - t_j\\)$\n\nSkipping first passage time distribution, reflection principle, maximum to date\n\n=== Stochastic Calculus\n<stochastic-calculus>\nThe purpose of stochastic calculus is to integrate an adapted stochastic process on a Brownian motion.\n\n$integral_0^T Delta\\(t\\)d W\\(t\\)$ This $Delta$ describes our position in an asset at time $t$, which is dependent on the change in price. This integral is the value of the portfolio at time $T$. The change in price is a stochastic process based on Brownian motion, which means that we cannot integrate with a Riemann/Lebesgue integral.\n\nIn order to construct the #strong[Itô integral], we add up increments of $W\\(t\\)$ and multiply them by the stock position $Delta\\(t\\)$ at that time. The discrete analogy is that of $W\\(t\\)$ as stock price and the time intervals being trading dates. The value of the portfolio is therefore the sum of the product of: the change in stock price each day, and the position held in stock on that day. Taking this discrete analogy continuously creates the integral.\n\nThe Itô integral is described as follows:\n\n$I\\(t\\)= integral_0^t Delta\\(u\\)d W\\(u\\)$\n\nThe Itô integral is a martingale. This is proved here.\n\n== Final Methods\n<final-methods>","dataUrl":"https://sharifhsn.dev/api/posts/stochastic-calculus-review.json","date":"2024-10-31","datePublished":"2024-10-31","description":"Recognize that this is some function that looks like 𝑓 ( 𝑎 , 𝑏 ) , like where 𝑎 is a deterministic component and 𝑏 is a stochastic component. Then solve your partials","site":"https://sharifhsn.dev","slug":"stochastic-calculus-review","source":"FE-610 | Stochastic Calculus","sourceUrl":null,"tags":["Stochastic Calculus","Review","Pricing"],"title":"Stochastic Calculus Review","url":"https://sharifhsn.dev/blog/stochastic-calculus-review/","version":"1","wordCount":2592}