Portfolio optimization with insider s initial information and counterparty risk

Size: px
Start display at page:

Download "Portfolio optimization with insider s initial information and counterparty risk"

Transcription

1 Portfoio optimization with insider s initia information and counterparty risk Caroine HILLAIRET Ying JIAO Juy 11, 212 Abstract We study the gain of an insider having private information which concerns the defaut risk of a counterparty. More precisey, the defaut time τ is modeed as the first time a stochastic process hits a random barrier L. The insider knows this barrier (as it can be the case for exampe for the manager of the counterparty), whereas standard investors ony observe its vaue at the defaut time. A investors aim to maximize the expected utiity from termina weath, on a financia market where the risky asset price is exposed to a sudden oss at the defaut time of the counterparty. In this framework, the insider s information is modeed by using an initia enargement of fitration and τ is a stopping time with respect to this enarged fitration. We prove that the reguator must impose short seing constraints for the insider, in order to excude the vaue process to reach infinity. We then sove the optimization probem and we study the gain of the insider, theoreticay and numericay. In genera, the insider achieves a arger vaue of expected utiity than the standard investor. But in extreme situations for the defaut and oss risks, a standard investor may in average outperform the insider, by taking advantage of an aggressive short seing position which is not aowed for the insider, but at the risk of big osses if the defaut finay occurs after the maturity. Keywords : asymmetric information, enargement of fitrations, counterparty risk, optima investment, duaity, dynamic programming. MSC 21 : 6H3 91B28 91G4 93E2 This research is part of a project of Europace Institute of Finance. We thank Laurent Denis, Nicoe E Karoui, Monique Jeanbanc, Huyên Pham, Abass Sagna, Nizar Touzi and Lioudmia Vostrikova for discussions. CMAP Ecoe Poytechnique, Emai: caroine.hiairet@poytechnique.edu. Financia support by Chair Financia Risks of the Risk Foundation, Chair Derivatives of the Future sponsored by the Fédération Bancaire Française, Chair Finance and Sustainabe Deveopment sponsored by EDF and Cayon. LPMA Université Paris Diderot and Peking University, Emai: jiao@math.univ-paris-diderot.fr. Financia support by Ama Recherche. 1

2 1 Introduction The insider s optima investment is a cassica probem where an investor possessing some extra fow of information aims to maximize the expected utiity on the fina vaue of her portfoio. As the insider has more information, she has access to a arger set of avaiabe trading strategies, eading to a higher expected utiity from termina weath. In the iterature, an interesting question has been studied: what is the cost of the extra information? From an indifference point of view, we search for the vaue at which the investor accepts to buy the information at the initia time, that is, the amount of money she is ready to pay such that this cost is offset by the increase of the maxima expected utiity. This is the approach adopted by Amendinger et a. 1, where the authors study the vaue of an initia information in the setting of a compete defaut free market. The extra information they consider is a termina information distorted by an independent noise, for exampe, a noisy signa of a functiona of the fina vaue of the assets. We adopt a more direct manner: we are interested in the gain of the insider from her investment strategy on the portfoio compared to other investors not having access to the extra information. The originaity of our paper is to study this probem in the context of credit risks: the insider s information concerns the defaut risk of a counterparty firm. During the financia crisis, the counterparty defaut has become an important source of risk we need to take into account. Jiao and Pham 14 have considered an optima investment probem where the risky asset in the portfoio is subjected to the defaut risk of a counterparty firm and its vaue may suffer a sudden oss at the counterparty defaut time τ. This paper is a good benchmark in our study in order to quantify the vaue of the extra information. The accessibe information for a standard investor is described as in the cassica credit risk modeing by Bieecki and Rutkowski 4, using the progressive enargement of a reference defaut-free fitration F = (F t ) t by the defaut τ. To anayze the impact of defaut, the defaut density framework deveoped in E Karoui et a. 5 has been adopted. This current paper concentrates on an insider in comparison with a standard investor. Both agents can invest in the same risk-free asset and risky one and they observe the same market price for each asset. However, the insider possesses more information on the risky asset since it is infuenced by the counterparty defaut on which the insider has additiona knowedge. Due to the extra information, the insider may gain arger profit. The insider s information is modeed by using an initia enargement of fitration as in 1 and in Grorud and Pontier 7. More precisey, in the credit risk context, we mode the defaut time τ as the first time that a stochastic process hits a random barrier L. The insider knows the barrier from the initia time and the other investors ony see its vaue at the defaut time. This extra information is caed the insider s information, or the fu information in Hiairet and Jiao 9. We sha consider the insider s optimization probem in parae with the one studied in 14. The canonica decomposition of processes adapted to the enarged fitration induces to specify the investment strategies on the two sets: before-defaut one {t < τ} and after-defaut one {t τ}, which is a simiar point to 14. However, due to the extra knowedge on the defaut barrier L, the insider s strategy depends on L before the counterparty defaut, which is not the case for the standard investor. If the defaut occurs, the insider s strategy wi depend on the defaut time τ. From the methodoogy point of view, the main difference here is that for the insider, the defaut time is modeed as in the cassica structura approach mode since the random barrier L is known, so that τ becomes a stopping time w.r.t. the reference fitration F. 2

3 Therefore, the defaut density hypothesis, which is crucia in 14, fais to hod for the insider and we can no onger adopt the conditiona density approach in this situation. We appy the theory of initia enargement of fitration, assuming that the conditiona aw of L given F t is equivaent to the aw of L. The corresponding Radon Nikodym derivative process, (p t (.), t ) wi pay a key roe in our methodoogy. The main observation of our study is that, if the short-seing is not reguated, then the insider can obtain unbounded termina weath. This justifies the necessity to consider the optimization probem with portfoio strategies where the short-seing is imited to a given eve. We decompose the optimization probem as an after-defaut probem and a goba before-defaut one, that we sove respectivey by using the dua and the dynamic programming methods. To make comparison with the standard investor in 14, we choose to consider CRRA utiity function. The paper is organized as foows. In section 2, we introduce the mode for the counterparty defaut, and we define and compare the informationa structure of an insider with respect to a standard investor. In Section 3 we present the insider s investment probem and we decompose it into an after defaut and goba before defaut ones, using the Radon Nikodym derivative process. We aso prove the necessity of imposing short seing constraint for the insider to excude the vaue process to reach infinity. In Section 4 we sove the two optimization probems: the after defaut one through duaity methods in a defaut free compete market, and the goba before defaut one through dynamic programming approach. We perform theorica comparison of the vaue process of the after defaut optimization probem, and for the goba before defaut optimization probem, the comparison is done in Section 5 through numerica iustrations. 2 Counterparty defaut mode and information We first introduce the mode for the counterparty defaut, which is a genera and standard mode in the credit risk anaysis. Let us fix a probabiity space (Ω, A, P) equipped with a reference fitration F = (F t ) t satisfying the usua conditions, which represents the defautfree information. Let τ be a positive random time denoting the defaut time of the counterparty, which is not necessariy an F stopping time. The defaut mode We consider the defaut risk of the counterparty in a genera barrier mode. Let (λ t, t ) be a positive F-adapted process representing the defaut intensity process of the counterparty. Denote by Λ t = t λ sds. It is an increasing process. We mode the defaut time as the first passage time of the process Λ to a positive random barrier L, i.e. (2.1) τ = inf{t : Λ t L} where the defaut threshod L is a positive A-measurabe random variabe. In the particuar but widey used case of Cox process mode, L is independent of F and foows the uni-exponentia aw. In the case where L is constant or deterministic, τ is an F-stopping time as in the cassica structura defaut modes. Information of the insider We suppose that, besides the information on the defaut-free market, the insider has compete information on L: this is the case for exampe of the counterparty firm s managers who 3

4 determine the defaut threshod. This fu information is modeed as the initia enargement of the fitration F by L and denoted by G M = (G M t ) t, G M t = F t σ(l). Without oss of generaity, we assume that a the fitrations we dea with in the foowing satisfy the usua conditions. We suppose that a standard investor on the market observes whether the defaut has occurred or not and if so, the defaut time τ, together with the information contained in the fitration F. Mathematicay, this information is represented by the progressive enargement of fitration F by τ, or more precisey, by the fitration G = (G t ) t where G t = F t D t, D t = σ(1 τ s, s t). This is the standard credit risk modeing for an market investor as in 4. The investor s information is incuded in the insider s information fow. We have G t Gt M for any t. In fact, before the defaut τ, i.e., on the set {t < τ}, the insider has additiona information on L, so her information Gt M is in genera stricty arger than G t. After the defaut occurs, both of them observe the defaut event and subsequenty the vaue of L so that they have equa information fow. We reca the canonica decomposition of G M -adapted (respectivey G M -predictabe) processes (see Jeuin 13 Lemma 3.13 and 4.4). Lemma For t, any Gt M -measurabe random variabe can be written in the form Y t = 1 τ>t Yt (L) + 1 τ t Yt 1 (τ) where Yt ( ) and Yt 1 ( ) are F t B(R + )-measurabe. 2. Any G M -adapted process Y admits the decomposition form Y t = 1 τ>t Y t (L) + 1 τ t Y 1 t (τ) where Y ( ) and Y 1 ( ) are F B(R + )-adapted Any G M -predictabe process Y admits the decomposition form Y t = 1 τ t Y t (L)+1 τ<t Y 1 t (τ) where Y ( ) and Y 1 ( ) are P(F) B(R + )-measurabe, P(F) being the predictabe σ-agebra associated with the fitration F. Remark 2.2 To compare with the case of a standard investor, we reca that any G t -measurabe random variabe Z t can be written as Z t = 1 τ>t Zt + 1 τ t Zt 1 (τ) where Zt and Zt 1 ( ) are respectivey F t -measurabe and F t B(R + )-measurabe. 3 Insider s optimization probem 3.1 Portfoio investment strategy and weath process From now on, a finite horizon T is fixed and a investment strategies take pace from time to time T. The insider has access to the same financia market as the standard investor, more precisey, she can invest in two types of financia assets. The first one is a risk-free bond with stricty positive vaues. We choose it as the numéraire and assume, without oss of generaity that the vaue of this bond equa to 1. The other asset is a risky one which is affected by defaut risk of the counterparty firm on which the insider has extra information. The price of this risky asset is observabe by a investors on market at any time t, T. Since it is subject to the counterparty defaut risk, the price process is modeed by a G-adapted 1 Namey for any t, the function Y i t ( ) is F t B(R +)-measurabe. 4

5 process S, which admits the decomposition form (3.1) S t = S t 1 t<τ + S 1 t (τ)1 t τ, t T where S is F-adapted and S 1 ( ) is F B(R + )-adapted. We suppose that the asset suffers a contagious oss at the defaut time of the counterparty, that is, S 1 θ (θ) = S θ (1 γ θ), and that < γ θ < 1 for any θ so that the asset price remains stricty positive. The process γ is F-adapted and represents the proportiona oss at defaut. We consider the trading strategy of the insider, who chooses to adjust the portfoio of assets according to information accessibiity. Therefore, the investment strategy process is characterized by a G M -predictabe process π which represents the proportion of weath invested in the risky asset and is of the form π t = 1 t τ π t (L) + 1 t>τ π 1 t (τ), where π ( ) and π 1 ( ) are P(F) B(R + )-measurabe processes. Starting from an initia weath X R +, the tota weath of the insider s portfoio is then a G M -adapted process given by (3.2) X t = 1 t<τ X t (L) + 1 t τ X 1 t (τ) where the before-defaut weath process satisfies the sef-financing equation (3.3) dx t (L) = X t (L)π t (L) ds t S t, t T and after the defaut τ, the weath process has a change of regime in its dynamics and satisfies (3.4) dxt 1 (τ) = Xt 1 (τ)πt 1 (τ) ds1 t (τ) St 1, t τ, T. (τ) At the defaut time, the weath has a jump in its vaue. Therefore, at time τ, the initia vaue of the after-defaut weath process is (3.5) X 1 τ (τ) = X τ (L) ( 1 π τ (L)γ τ ). We suppose that π τ (L)γ τ < 1, so that the weath remains stricty positive after the jump due to the counterparty defaut. We consider the foowing dynamics for the asset price S on before-defaut set {t < τ} for S and on after-defaut set {t τ} for S 1 : ds t = S t (µ t dt + σ t dw t ), t T ds 1 t (θ) = S 1 t (θ)(µ 1 t (θ)dt + σ 1 t (θ)dw t ), θ t T where the coefficients µ and σ are F-adapted processes, µ 1 (θ) and σ 1 (θ) are F B(R + )-adapted processes, and W is an F-Brownian motion. In addition, we suppose the integrabiity condition T µ t σt 2 dt + T θ µ1 t (θ) 2 dt + σt 1(θ) T 5 σ t 2 dt + T θ σ 1 t (θ) 2 dt <.

6 So the vaues of the before-defaut and after-defaut weath satisfy the dynamics (3.6) (3.7) dx t (L) = X t (L)π t (L)(µ t dt + σ t dw t ), t T dx 1 t (τ) = X 1 t (τ)π 1 t (τ)(µ 1 t (τ)dt + σ 1 t (τ)dw t ), t τ, T and the jump at defaut of the weath process is given by the equaity (3.5). Finay, we define the admissibe trading strategy famiy A L as the set of pairs (π ( ), π 1 ( )), where π ( ) and π 1 ( ) are P(F) B(R + )-measurabe processes such that ( τ T (3.8) >, T ) πt ()σt 2 dt + πt 1 ( )σt 1 ( ) 2 dt < and πτ ()γ τ < 1, a.s., T where is the F-stopping time defined by := inf{t : Λ t }. Remark 3.1 Let A denote the set of a G M -predictabe processes π such that T π tσ t 2 dt < and π τ γ τ < 1. If (π ( ), π 1 ( )) is an eement in A L, then (π t = π t (L)1 τ t +π 1 t (τ L )1 τ<t, t ) is a processus in the set A. Conversy, given a process π A, there exists a pair (π ( ), π 1 ( )) A L such that π t = π t (L)1 τ t + π 1 t (τ L )1 τ<t for any t, thanks to Lemma The optimization probem The insider has the objective to maximize her expected utiity function on the termina weath of the portfoio. Let U be a utiity function defined on (, + ), stricty increasing, stricty concave and of cass C 1 on (, + ), and satisfying im x + U (x) = + and im x U (x) =. We sha consider the probem (3.9) V = sup π A L EU(X T ) and search for the optima strategy π for the insider. A simiar probem has been studied in 14 for a standard investor with G-predictabe strategy π = (π, π 1 ( )) and G-adapted weath X = (X, X 1 ( )). The admissibe strategy set A consists of pairs π = (π, π 1 ( )) where π and π 1 ( ) are respectivey F-predictabe and P(F) B(R + )-measurabe, and such that θ R +, T T π t σt 2 dt + π 1 t (θ)σt 1 (θ) 2 dt < and π θ γ θ < 1, a.s. θ In this paper, we concentrate on the insider s optimization probem (3.9) and we are interested in the information fow impact on the trading strategies. Intuitivey, the insider shoud have a arger gain of investment due to the extra information. Indeed, if the investor and the insider have the same utiity function and if they can invest in the same financia assets, then the ony difference between them reies on their avaiabe information in the sense that the corresponding fitrations satisfy G G M, which impies the same incusion for the sets of admissibe strategies: A A L. Thus the corresponding vaue functions satisfy (3.1) V := sup EU(X T ) V. π A 6

7 since the supremum in V is taken on a smaer set than the one in V. Remark here that on a given sampe path, it may happen that U(X T ) > U(X T ) if both investor and insider foow their optima strategies π and π respectivey, but one aways have in expectation for optima strategies EU(X π T ) EU(X π T ). as The expectation EU(X T ) can be written, by using the weath decomposition formua (3.2), (3.11) EU(X T ) = E1 T <τ U(X T (L)) + 1 T τ U(X 1 T (τ)). The aim of the above formuation, simiar as in 14, is to reduce the initia optimization probem in an incompete market into two probems : the after-defaut and before-defaut ones. Nevertheess, the approach we adopt here is different since, as mentioned previousy, the random time τ is not a totay inaccessibe random time for the insider and we can no onger use the conditiona density approach to sove the probem. Our approach wi use the theory of initia enargement of fitration (aso caed the strong information modeing in 2) by using the vaue of the random defaut barrier L known to the insider. More precisey, we introduce a famiy of F-stopping times = inf{t : Λ t } for a > which are possibe reaizations of L and we work under an equivaent probabiity measure P L under which L is independent to F T. Thus in our framework, we sha need the Radon-Nikodym derivative process p t (L) which is the density of the historica probabiity measure P with respect to this equivaent probabiity measure P L and it wi pay a simiar roe as the defaut density process in 14. This probabiity density hypothesis is given beow. It is a standard hypothesis in the theory of initia enargement of fitration due to Jacod 11, 12. Assumption 3.1 We assume that L is an A-measurabe random variabe with vaues in, +, which satisfies the assumption : P(L F t )(ω) P(L ), t, T, P a.s.. We denote by P L t (ω, dx) a reguar version of the conditiona aw of L given F t and by P L the aw of L (under the probabiity P). According to 12, there exists a measurabe version of the conditiona density (3.12) p t (x)(ω) = dp t L (ω, x) dp L which is an positive (F, P)-martingae and hence can be written as (3.13) p t (x) = p (x) + t p s (x)ρ s (x)dw s, x, +, t, T for some F-predictabe process (ρ t (x)) t,t. It is proved in 7 that Assumption 3.1 is satisfied if and ony if there exists a probabiity measure equivaent to P and under which F T and σ(l) are independent. Among such equivaent probabiity measures, the probabiity P L defined by the Radon-Nikodym derivative process dp E P L G M dp L t = p t (L) 7

8 is the ony one that is identica to P on F. For exampes of L and expicit computations of corresponding p t (L), interest reader may refer to 1. We reca that we consider the optimization probem (3.9) V = sup π AL EU(X T ). Since the initia σ-fied is non-trivia, it is usefu to consider the conditiona optimization probem (3.14) ess sup π A L EU(X T ) G M, where G M = σ(l). The ink between those two optimization probems (3.9) and (3.14) is given by 1 : if the supremum in (3.14) is attained by some strategy in A L, then the ω-wise optimum is aso a soution to (3.9). Athough the supremum is not necessariy attained in our probem, we wi see in Proposition 4.12 that there exists a sequence of admissibe strategies π n such that EU(X πn T ) GM converges in L1 to (3.14) and we can prove that for the same sequence, EU(X πn T ) converges to V. Thus, we wi first sove the optimization probem (3.14), and then deduce the soution of (3.9) by taking the expectation, as expained in the foowing Proposition: Proposition 3.2 Under the Assumption 3.1, one has EU(X T ) G M = E p T () ( 1 T <τ U(XT ()) + 1 T τ U(XT 1 ( )) ) =L EU(X T ) = E p T () ( 1 T <τ U(XT ()) + 1 T τ U(XT 1 ( )) ) P L (d) R + where for >, := inf{t : Λ t }. Proof: We wi use the change of probabiity to P L in order to reduce to the case where L and F T are independent. Firsty, by (3.11), EU(X T ) G M = E 1 T <τ U(XT (L)) + 1 T τ U(XT 1 (τ)) L = E P L pt () ( 1 T <τ U(XT ()) + 1 T τ U(XT 1 ( )) ) =L = E p T () ( 1 T <τ U(X T ()) + 1 T τ U(X 1 T ( )) ) =L where the ast two equaities foow respectivey from the facts that F and σ(l) are independent under P L and that P L is identica to P on F T. Thus EU(X T ) = E E 1 T <τ U(XT (L)) + 1 T τ U(XT 1 (τ)) L = E p T () ( 1 T <τ U(XT ()) + 1 T τ U(XT 1 ( )) ) P L (d). R + This motivates to introduce, for any >, the set A of pairs π = (π, π 1 ( )), where π and π 1 ( ) are respectivey F-predictabe and P(F) B(R + )-measurabe processes, such that (3.15) τ T π t σ t 2 dt + T T and consider the foowing optimization probem π 1 t ( )σ 1 t ( ) 2 dt < and π γ τ < 1, a.s. (3.16) V () = sup π A E p T () ( 1 T <τ U(X T ) + 1 T τ U(X 1 T ( )) ), where is the F-stopping time inf{t : Λ t }. 8

9 3.3 The necessity of imiting short seing for the insider before defaut In this subsection, we show that if the reguators do not impose any constraint on short seing for the insider before the defaut of the counterparty, then the insider can achieve a termina weath that is not bounded in L 1. Proposition 3.3 We suppose that the foowing conditions are satisfied: (1) the process Λ is a.s. stricty increasing on, T, (2) for any in the support of the distribution of the aw of L, one has P(Λ T ) >. Then we have ess sup EX T G M = + a.s. π A L In addition, for any utiity function U such that im U(x) = +, x + ess sup EU(X T ) G M = + a.s. π A L Proof: Let ϕ :, +, + be an increasing function such that ϕ() < for any, +. Let ψ > be a constant. For each, +, we define a strategy π() = (π (), π 1 ( )) A as foows π t () = ψ1 τϕ() <t, π 1 ( ). Note that (π ( ), π 1 ( )) is an admissibe strategy in A L. The vaue at the time of the corresponding weath process X ϕ,ψ () is equa to (1 + γ τ ψ)x ϕ,ψ. By the dynamics of the weath τ process (3.6) and (3.7), on { T }, we have ( ) Moreover, X ϕ,ψ = X (1 + γ τ ψ) exp τ τ ϕ() ( µ t ψ (σ t ) 2 ψ 2) dt τ τ ϕ() σ t ψdw t EX T G M = E p T () ( 1 T <τ XT () + 1 T τ XT 1 ( ) ) =L E 1 T τ p T ()Xτ ϕ,ψ Now fix an increasing sequence (ϕ n ) n 1 of functions such that ϕ n () < for, and that im n + ϕ n () =. By the condition (1), one obtains that τ ϕn() converges a.s. to when n +. The sequence of random variabes T τ ϕn() T σ t ψ dw t, n 1 converges a.s. to. Then by Fatou s emma im inf E1 T τ p T ()X ϕn,ψ n + E1 T τ p T ()X (1 + γ τ ψ) which impies the first assertion since P( T ) > by condition (2). argument and the assumption on U to obtain im ψ + E1 T τ p T ()U(X (1 + γ τ ψ)) =L =.. =L We use the simiar 9

10 Remark 3.4 The strategies mentionned in this proof are not arbitrage strategies because for any fixed function ϕ as in the proof, P(T τ ϕ(l), τ ) > and on this event, the strategy of the insider that consists of betting on the defaut before maturity T turns out to be a wrong bet. Thus, on a non nu probabiity set, the strategy of a standard investor outperforms the one of the insiders, athough the converse inequaity hods on expectation. Proposition 3.3 justifies to consider, instead of A L defined in (3.8), the foowing admissibe trading strategy sets for the insider : for a δ, (3.17) A δ L = {(π ( ), π 1 ( )) A L such that π δ} and to quantify the impact of the ower bound δ. Simiar to (3.15), we define A δ = {(π, π 1 ( )) A such that π δ}. Before considering the main optimization probem (3.18) ess sup EU(X T ) G M π A δ L we first consider an aternative famiy of optimization probems depending on the parameter, +, (3.19) V δ () := sup E p T () ( 1 T <τ U(XT ()) + 1 T τ U(XT 1 ( )) ). π A δ The foowing theorem shows that the optima vaue of the optimization probem (3.18) is actuay equa to V δ (L). Theorem 3.5 With the above notation, we have V δ (L) = ess sup EU(X T ) G M π A δ L a.s. Proof: 3.2 we obtain that Assume that (π ( ), π 1 ( )) is an eement in A δ L, then (π (), π 1 ( )) A δ. By Proposition ess sup EU(X T ) G M V δ (L). π A δ L For the converse inequaity, we sha use measurabe seection theorem. For any ε > and any,, et F ε () be the set of strategies (π, π 1 ( )) A δ which are ε-optima with respect to the probem (3.19), namey such that E p T () ( 1 T <τ U(X T ()) + 1 T τ U(X 1 T ( )) ) { V δ 1/ε, () ε, if V δ () < +, if V δ () = +. By a measurabe seection theorem (cf. Benes 3, Lemma 1), there exists a measurabe (with respect to ) famiy {(π (), π 1 (, ))} R+ such that (π (), π 1 (, )) F ε () for any >. Finay et π ( ) := π ( ), π 1 t (x) := 1 t>x π 1 t (x, Λ x ). 1

11 We have ( π ( ), π 1 ( )) A L and π 1 t ( ) = π 1 t (, ) for any > on { < t}. Therefore, by Proposition 3.2, EU( X T ) G M = E p T () ( 1 T <τ U(XT ()) + 1 T τ U(XT 1 ( )) ) =L { V δ (L) ε, if V δ (L) < +, 1/ε, if V δ (L) = +. Since ε is arbitrary, we obtain EU( X T ) G M V δ (L). 4 Soving the optimization probems In this section, we concentrate on soving the optimization probem (3.19) sup E p T () ( 1 T <τ U(XT ()) + 1 T τ U(XT 1 ( )) ) π A δ for any fixed. We reca that the before-defaut and after-defaut weath processes X and X 1 are governed by two contro processes π and π 1 respectivey, so we need to search for a coupe of optima contros ˆπ = (ˆπ, ˆπ 1 ). In Theorem 4.1 we expain how to decompose the optimization probem into two probems each depending ony on π and on π 1 respectivey. The after-defaut optimization probem wi be soved firsty, using the after-defaut fitration F 1 F 1 := (F τ t) t,t. Remark that the initia σ-fied of the fitration F 1 is not trivia : F 1 = F τ = {A A : A { t} F t, t, T } and is F τ -measurabe. A the F 1 -adapted processes invoved in the after-defaut optimization probem are indexed on the right-upper side by the symbo 1. In particuar, we denote by X 1,x ( ) the soution of the SDE (3.7) defined on the stochastic interva, T starting from the F-stopping time with F τ -measurabe initia vaue x. We define by A 1 the admissibe predictabe strategy set (πt 1 ( ), t, T ) such that T πt 1 ( )σt 1 ( ) 2 dt < a.s.. The goba before defaut optimization probem invoves the soution of the after-defaut optimization probem and wi be soved in a second step, using the stopped fitration F F := (F τ t) t,t. A the F -adapted processes invoved in the after-defaut optimization probem are indexed on the right-upper side by the symbo. The admissibe predictabe strategy set A,δ is (πt (), t, T ) such that T πt ()σt 2 dt <, πt () δ and 1 > πτ ()γ τ, a.s.. 11

12 Theorem 4.1 For any >, we denote by V 1 (x ) the optima vaue of the after-defaut optimization probem (4.1) Vτ 1 (x ) = ess sup Ep T ()U(X 1,x T ( )) F τ, π 1 ( ) A 1 then the goba optima vaue V δ () defined in (3.19) can be written as the soution of a beforedefaut optimization probem as (4.2) V δ () = sup E 1 T <τ p T ()U(XT ()) + 1 T τ V 1 ( X τ ()(1 πτ ()γ τ ) ). π A,δ Proof: Consider firsty an arbitrary admissibe strategy (π, π 1 ( )) A δ for a fixed >. By definition, (πt, t, T ) A,δ and (πt 1 ( ), t, T ) A 1. Taking the conditiona expectation with respect to F τ eads to the foowing inequaities E p T () ( 1 T <τ U(XT ()) + 1 T τ U(XT 1 ( )) ) = E 1 T <τ p T ()U(X T ()) + 1 T τ E ( p T ()U(X 1 T ( )) F τ ) E 1 T <τ p T ()U(X T ()) + 1 T τ V 1 (X ()(1 π ()γ τ )) V δ (). For the converse inequaity, et us assume for the moment that the esssup in the definition (4.1) is achieved for a given ˆπ 1 ( ) (see section 4.1 for the proof). Then for any (π t, t, ) in A,δ, by a measurabe seection theorem, there exists a P(F) B(R + )-measurabe process π 1 ( ) such that π 1 ( ) = ˆπ 1 ( ) on, T and that (π, π 1 ( )) A δ, where we have extended π to an F-predictabe process on R +. Thus V δ () E 1 T <τ p T ()U(X T ()) + 1 T τ V 1 (X ()(1 π ()γ τ )) By taking the supremum over a (π t (), t, ) A,δ, we obtain the desired inequaity. Remark 4.2 The supremum in V δ () can be approached by a sequence of admissibe strategies in A,δ (see Proposition 4.12), which induces a sequence of strategies in A δ L such that the corresponding vaue functions converge to V δ := sup π A δ EU(X T ). Thus V δ = L R + V δ()p L (d). Remark 4.3 The process (p t (), t, T ) is essentia in our approach of initia information and pays a simiar roe to the defaut density process in 14 (α t (θ), t, T ) defined as α t (θ)dθ = P(τ dθ F t ). Thus, to quantify the gain of the insider (as in the forthcoming Proposition 4.5), it is interesting to compare those two processes. In the particuar case where the F t -conditiona aw of L admits a density g t () w.r.t. the Lebesgue measure, the defaut density can be competey deduced (see 6, Proposition 3) in this framework as α t (θ) = λ θ g t (Λ θ ), t θ and α t (θ) = Eα θ (θ) F t for t θ where λ is the process given in Section 2. In the genera case, the aw of L can have atoms, then the defaut density does not exist and the approach in 14 is no onger vaid, whereas the insider s optimization probem can be soved with the process p. (). 12

13 4.1 The after-defaut optimization In this section, we focus on the after defaut optimization probem (4.1) in the fitration F 1 = (F τ t) t,t Vτ 1 (x ) = ess sup Ep T ()U(X 1,x T ( )) F τ π 1 ( ) A 1 where is an F stopping time and the initia weath x is F τ -measurabe. This probem is simiar to a standard optimization probem, we wi extend the resuts in our framework where the initia time is a random time (and is an F-stopping time). We define the process ( τ t µ 1 Z t ( ) = exp u( ) σu(τ 1 ) dw u 1 2 τ t µ 1 u( ) σu(τ 1 ) 2 ) du, t, T. This process is an F 1 -oca martingae (cf. 15, page 2), we assume that the coefficients µ 1 ( ) and σ 1 ( ) satisfy a Novikov criterion (see Theorem 4.4 beow) so that (Z t ( )) t,t is an F 1 - martingae. Theorem 4.4 We assume that for any, the coefficients µ 1 u( ) and σu(τ 1 ) satisfy the Novikov criterion ( 1 τ T E exp µ 1 2 ) u( ) 2 σu(τ 1 ) du <. Then the vaue function process to probem (4.1) is a.s. finite and is given by ( 1 (4.3) ˆV τ (x ) = E p T ()U I ( ŷ τ (x ) Z T ( )) ) F 1 p T () where I = (U ) 1 and the Lagrange mutipier ŷ τ ( ) is the unique F τ B(R + )-measurabe soution of the equation 1 Z t ( ) E Z T ( )I ( ŷ τ (x ) Z T ( )) F 1 = x. p T () The corresponding optima weath is equa to (4.4) ˆX1,x t ( ) = 1 Z t ( ) E Z T ( )I ( ŷ τ (x ) Z T ( )) Ft 1 p T (), t, T. Proof: Note that after the defaut, the market is compete. The process (Z t ( )X 1,x t ( )) t,t is a positive oca F 1 -martingae, and thus a supermartingae, eading to the foowing budget constraint : ( ) (4.5) E Z T ( )X 1,x T ( ) F 1 x. Conversey, the martingae representation theorem on the brownian fitration impies that for any F T B(R + )-measurabe X T ( ), there exists a P(F) B(R + )- measurabe process φ( ) such that τ T X T ( )1 τ <T = E(X T ( )1 τ <T F τ ) + φ u ( )du 13

14 Therefore the after defaut optimization probem is soved by the mean of the Lagrange mutipier ( Vτ 1 (x ) = E p T ()U I ( ŷ τ (x ) Z T ( )) ) F 1 p T () and the optima weath is given by ˆX 1,x t ( ) = 1 Z t ( ) E Z T ( )I ( ŷ τ (x ) Z T ( ) p T () ) F 1 t where I = (U ) 1 and the Lagrange mutipier ŷ τ (x ) is F τ B(R + )-measurabe and satisfies ˆX 1,x ( ) = x. The existence, uniqueness and measurabiity of the Lagrange mutipier ŷ τ (x ) in the case of a non trivia initia σ-fied is proved in Proposition 4.5 of Hiairet 8. We can aready give a first anaytica comparison of the vaue function process of the after defaut optimization probem between the initia information and the progressive information. We reca (see 14, Theorem 4.1) that for the progressive information, this vaue function process at the defaut time θ = is ( V τ 1 (x) = E α T ( )U I ( ȳ τ (x) Z T ( )) ) (4.6) F τ α T ( ) ZT ( ) with y τ such that E Z t ( ) I( ȳ τ (x) Z T ( )) F τ = x α T ( ) which has to be compared to the vaue function (4.3) for an initia information. To do this, we assume that the F t -conditiona aw of L admits a density with respect to the Lebesgue measure, denoted as g t () (see Remark 4.3). Using Assumption 3.1, it is sufficient to assume that the aw of L admits a density with respect to the Lebesgue measure. Proposition 4.5 The vaue function processes of the after defaut optimization probem for the initia information ˆV τ 1 (x) and for the progressive information V τ 1 (x) satisfy V 1 (x) ˆV 1 (x) = g ()λ τ where g (.) denotes the density of the aw of L with respect to the Lebesgue measure and is supposed to be stricty positive. In addition, the ratio α T ( ) p T () is F τ -measurabe and is aso equa to V τ 1 (x) ˆV. τ 1 (x) Proof: On the stochastic interva, T, the density of probabiity Z t ( ) is the same for the initia and the progressive information, but not necessariy the Lagrange mutipiers, athough they satisfy the same type of equation. More precisey, ȳ τ (x) is the unique F τ B(R + )- measurabe soution of ZT ( ) (4.7) E Z t ( ) I( ȳ τ (x) Z T ( )) Fτ = x α T ( ) 14

15 and ŷ τ (x) is the unique F τ B(R + )-measurabe soution of ZT ( ) (4.8) E Z t ( ) I( ŷ τ (x) Z T ( )) Fτ = x. p T () We reca that g T () is the density of the F T -conditiona aw of L with respect to the Lebesgue measure. According to (3.12) and Remark 4.3, p T () = g T () g () and α T ( ) = λ τ g T (Λ τ ) = λ τ g T (). Thus the ratio α T ( ) p T () = λ τ g () is F τ -measurabe. Furthermore, (4.7) is equivaent to ZT ( ) E Z t ( ) I( ȳ τ (x) p T () α T ( ) Z T ( )) Fτ = x p T () so the unicity of the F τ -measurabe soution of the equation (4.8) impies that the ȳ τ (x) p T () α T (τ) = ŷ τ (x). Therefore, ȳ (x) ŷ τ (x) = λ g () and U ( I ( ȳ τ (x) Z T ( ))) ( α T ( ) = U I (ŷτ (x) Z T ( ))) p T (). We concude the proof by using again that the ratio α T ( ) p T () = λ τ g () is F τ -measurabe. 1 Remark The constant g () is a threshod to compare with λ. For a given scenario, if λ τ is smaer than 1 g (), which is often the case in practice, then the absoute vaue of the vaue function for the initia information is greater than the one for the progressive information. In this case, the insider puts a higher weight p T () on the after defautoptimization probem (4.1), compared to a standard investor whose weight is α T ( ). This means that p t () conveys more information than α t ( ). Moreover, the smaer g () is, the greater is the gain of an insider. The interpretation is as foows: sma vaue of g () means that P(L d) is sma, thus the insider, who knows the rea vaue of L, has very reevant information. 2. Concerning the optima weath ˆX 1 ( ) in (4.4), we can prove, using the same argument as in the proof of Proposition 4.5 that starting from a same weath x at the defaut time, the optima weath process of the after-defaut optimization probem is the same for the initia and the progressive information, i.e. ˆX1,x ( ) = X 1,x ( ). This resut is not surprising, since after the defaut, the two information fows coincide. But naturay the input weath of the after-defaut optimization probem wi not be the same for the two information fows since they are not the same before. We wi quantify numericay the gain of an insider compared to an standard investor for the goba optimization probem. We now consider, as in 14, Constant Reative Risk Aversion (CRRA) utiity functions U(x) = xp p, < p < 1, x > and I(x) = x 1 p 1. Direct computations from the previous theorem yied the optima weath ˆX 1,x t ( ) = x Z t ( ) E E p T () p T () ( ) p ZT ( ) p 1 p T () ( ZT ( ) p T () ) p p 1 Ft 1, t, T F τ 15

16 and the optima vaue function (4.9) ˆV 1 τ (x ) = xp p ( E p T () ( ) p ZT ( ) p 1 p T () F ) 1 p =: xp p K ony de- ( ( ) p where the F τ -measurabe random variabe K τ = E p T () ZT ( ) p 1 p T () pends on the stopping time and on market parameters. ) 1 p F 4.2 The before-defaut optimization We now consider the optimization probem (4.2) with CRRA utiity functions. Using (4.9), we have to sove : (4.1) V δ () = sup E 1 T <τ p T ()U(X T ()) + 1 T τ K τ U ( X τ ()(1 π τ ()γ τ ) ) π A,δ where the F τ -measurabe random variabe K τ does not depend on the contro process π A,δ. We wi use a dynamic programming approach. Reca that F = (F τ t) t,t is the stopped fitration at the defaut time. Since 1 t<τ is F τ t-measurabe, we have E 1 T <τ p T ()U(XT ()) + 1 T τ K τ U ( Xτ ()(1 πτ ()γ τ ) ) F τ t = 1 t τ K τ U ( Xτ ()(1 πτ ()γ τ ) ) + E 1 T <τ p T ()U(XT ()) + 1 t<τ T K τ U ( Xτ ()(1 πτ ()γ τ ) ) F τ t For any ν A,δ, we introduce the famiy of F -adapted processes (4.11) X t (ν) := ess sup E π A,δ (t,ν) 1 T <τ p T ()U(X T ()) + 1 t<τ T K τ U ( X ()(1 π ()γ τ ) ) F τ t where A,δ (t, ν) is the set of contros coinciding with ν unti time t: for any t, T, ν A,δ, A (t, ν) = {π A,δ : π. t = ν. t }, X ν, denotes the weath process derived from the contro ν A,δ. We have V δ () = X (ν) for any ν A,δ. In the foowing resut, we show that the short seing constraint pays a crucia roe in the optimization. In fact, it is the optima strategy at the defaut time. Proposition 4.7 For any π A,δ, there exists a sequence of strategies (πn A,δ ) n N such that πn,τ = δ and im X (π n + n) X (π ). 16

17 Proof: Let (τ n ) n N where τ n < be a increasing sequence of F stopping times that converge to. Starting from a strategy π A,δ, we define another strategy πn = 1,τnπ 1 τn,τ δ that remains in A,δ. We denotes as X () and Xn() the corresponding weath before defaut, and as X (π ) and X (πn) the corresponding vaue function for those strategies of the before defaut goba optimization probem. On the one hand, by dominated convergence theorem, it is easy to check that On the other hand, on the event {T } X n, ()(1 π n, ()γ τ ) X ()(1 π ()γ τ ) π A,δ im E 1 T <τ p T () U(X n + T ()) U(Xn,T ()) =. = (1 + δγ ( ) τ (1 πτ ()γ τ ) exp ( (πs + δ)µ s + 1 τ ) τ n 2 (σ s) 2 ((πs) 2 δ 2 ))ds (πs + δ)σsdw s τ n impies that π () δ and (1+δγ τ ) (1 π ()γ τ ) 1 (because γ > ). τ n impies that the exponentia term tends to 1 a.s. Thus im n + X n, ()(1 π n, ()γ τ ) X ()(1 π ()γ τ ) 1 and im E 1 T τ K τ U(X n + n, ()(1 πn,τ ()γ τ )) E 1 T τ K τ U(Xτ ()(1 πτ ()γ τ )). Consequenty, im X (π n + n) X (π ). We now characterize the optima strategy process. From the dynamic programming principe, the foowing resut hods: Lemma 4.8 For any ν A,δ, the process ξ ν t := X t (µ) + 1 t τ K τ U(X ν, ()(1 ν τ ()γ τ )), t T is an F -supermartingae. Furthermore, the optima strategy π is characterized by the martingae property : (ξ π t ) t T is an F -martingae. Proof: Let s, t be two times such that s t T. E X t (µ) + 1 t τ K τ U(Xτ ν, ()(1 ν τ ()γ τ )) F τ s = E X t (µ) + 1 s<τ tk τ U(Xτ ν, ()(1 ν τ ()γ τ )) F τ s + 1 s τ K τ U(Xτ ν, ()(1 ν τ ()γ τ )) 17

18 We make expicit the conditiona expectation : E X t (µ) + 1 s<τ tk τ U(Xτ ν, ()(1 ν τ ()γ τ )) F τ s = E ess sup (4.12) π A,δ (t,ν) ess sup π A,δ (s,ν) E 1 T <τ p T ()U(XT ()) + 1 t<τ T K τ U(Xτ ()(1 πτ ()γ τ )) F τ t + 1 s<τ tk τ U(Xτ ν, ()(1 ν τ ()γ τ )) F τ s E 1 T <τ p T ()U(XT ()) + 1 s<τ T K τ U(Xτ ()(1 πτ ()γ τ )) F τ s the ast inequaity foowing from the fact that in the ast esssup, the optima contro is taken from the date s t. Thus E X t (µ) + 1 t τ K τ U(Xτ ν, ()(1 ν τ ()γ τ )) F τ s X s (µ) + 1 s τ K τ U(Xτ ν, ()(1 ν τ ()γ τ )) and (ξt ν ) t T is an F -supermartingae. It is an F -martingae if and ony if the inequaity (4.12) is an equaity for a t, T, meaning that ν is the optima contro on, t, for a t T. This characterizes the optima strategy. Remark that the F -adapted process (4.13) Y t := X t(µ) U(X ν, t ()) ( X = ess sup E 1 T <τ p T () T () ) p ( X τ () ) p(1 + X ν, 1t<τ T K τ π t () X ν, τ t () ()γ ) F τ t π A,δ (s,ν) t T does not depend on ν A,δ, and is constant after. We wi give a characterization of the process Y in terms of a backward stochastic differentia equation (BSDE) and of the optima strategy. Before this, we give a characterization of F -martingae. Lemma 4.9 Let (M t ) t,t be an F -martingae. Then there exists an F-predictabe process φ in L 2 oc (W ) such that M t = M + t φ s1 s τ dw s, t, T. Proof: We first prove that (M t ) t,t is aso an F-martingae. Indeed, for s t T M s = E(M t F τ s) = E(E(M t F τ ) F s ) = E(M t F s ) because M t is F τ -measurabe. Thus, by representation theorem for the F-martingae, and since (M t ) t,t is stopped at time, there exists φ an F-predictabe process such that M t = M + t φ s1 s τ dw s. We are now ready to characterize the optima strategy. Remark that Y t = on, (and zero after ) thus Y L + (F ) where Xt(µ) U(X ν, t ()) is positive L + (F ) := {Ỹ : F -adapted processes s.t. Ỹ t > for t, and Ỹt = for t, }. 18

19 Theorem 4.1 The process Y defined in (4.13) is the smaest soution in L + (F ) to the BSDE (4.14) (1 + δγ τ ) p T τ T τ Y t = 1 T <τ p T () + 1 t<τ T K τ + f(θ, Y θ, φ θ )dθ φ θ dw θ, t, T, p t t for some φ L 2 oc (W ), and where (4.15) f(s, Y s, φ s ) = p ess sup ν A,δ,s.t. ν τ = δ (µ s Y s + σsφ s )ν s 1 p 2 Y s ν s σs 2 The optima vaue function is attained as the imit of a sequence of admissibe strategies (ˆπ t ) t,t τ converging to the supremum in (4.15). Remark 4.11 As in Theorem 4.2 in 14, the optima strategy before defaut is characterized through the optimization of the driver of a BSDE. However, the main difference reies in the fact that in our case, the driver has a jump at the defaut time. Nevertheess, since the jump occurs (if it occurs) ony at the termina date of the BSDE, standard theory on BSDE sti appy. Proof: By Lemma 4.8, for any ν A,δ ξt ν = U(X ν, t ())Y t + 1 t τ K τ U(X ν, ()(1 ν τ ()γ τ ) t T is an F -supermartingae. In particuar, by taking ν =, we see that the processes (Y t + K τ 1 t τ ) t T, and thus (Y t ) t T are F -supermartingaes. By the Doob-Meyer decomposition and Lemma 4.9, there exists φ L 2 oc (W ), and a finite variation increasing F -predictabe process A such that: (4.16) dy t = φ t dw t da t, t, T. From Itô s formua, we deduce that the finite variation process in the decomposition of the F -supermartingae ξ ν, ν A,δ, is given by A ν with da ν t = (X ν, t ()) p{ 1 p da t (µ t Y t + σt φ t 1 t τ )ν t dt 1 p 2 Y t ν t σt 2 (1 ν t γ t ) p dt K t d1 t τ }. p A ν is nondecreasing and the martingae property of ξˆπ and impies that da t = p (µ t Y t + σt φ t 1 t τ ))ˆπ t dt 1 p 2 Y t ˆπ t σt 2 (1 ˆπ t γ t ) p dt + K t d1 t τ p A t = A + p ess sup ν A,δ t ( (µ s Y s + σsφ s )ν s 1 p 2 Y s ν s σs 2) (1 ν τ γ τ ) p ds + 1 t τ K τ. p Maximizing at eads to ν τ = δ (see Proposition 4.7) and A t = A + p ess sup ν A,δ s.t. ν τ = δ t ( µ s Y s + σsφ s )ν s 1 p 2 Y s ν s σs 2) (1 + δγ τ ) p ds + 1 t τ K τ. p 19

20 Furthermore, Y T = 1 T <τ p T () and (Y t ) t T is constant after, thus (Y, φ) soves the BSDE (4.14). Note that Y is not a continuous process, it may jump at time. We now prove that Y is smaest soution in the L + (F ) to the BSDE (4.14). Let Ỹ L+ (F ) be another soution, and we define the famiy of nonnegative F -adapted processes ξ ν (Ỹ ), ν A,δ, as ξt ν (Ỹ ) = U(Xν, t ())Ỹt + 1 t τ K τ U(Xτ ν, ()(1 ν τ ()γ τ )), t, T. By simiar cacuations as above, dξt ν (Ỹ ) = d M t ν dãν t, where Ãν is a nondecreasing F -adapted process, and M ν is a F - oca martingae. By Fatou s emma, this impies that the process ξ ν (Ỹ ) is a F -supermartingae, for any ν A,δ. Since ỸT = 1 T <τ p T (), we deduce that for a ν A,δ E U(X ν, T )1 T <τ p T () + 1 t τ K τ U(Xτ ν, ()(1 ν τ ()γ τ ) Ft U(X ν, t )Ỹt, t, T. Since p >, U(X ν, t ) is positive. By dividing the above inequaities by U(X ν, t ), we deduce by definition of Y (see (4.13)), and arbitrariness of ν A,δ, that Y t Ỹt, t T. This shows that Y is the smaest soution to the BSDE (4.14). For optimizing (4.13) via the BSDE (4.14), a naive approach wi consist in optimizing π at time, eading to an π = δ, and then optimizing for s < the driver f (s, Ys, φ s) = ess sup p ( µ sys + σsφ s)ν s 1 p 2 Y s ν s σs 2 where Y is soution to the BSDE Y t ν s δ = 1 T <τ p T ()+1 t<τ T K τ (1 + δγ τ ) p p T τ T τ + f (θ, Y θ, φ θ )dθ φ θ dw θ, t, T, t t eading to the optima portfoio ˆπ s. Thus, the natura candidate to be the optima strategy before defaut is (4.17) π np := 1,τ ˆπ δ1 τ, but unfortunatey π np is not a predictabe process. Nevertheess, we wi prove the existence of a sequence of predictabe strategies in A,δ such that the corresponding vaue function tends to the vaue function reative to this non predictabe strategy. To do this, for any strategy π A,δ, we reca the corresponding vaue function of the before defaut goba optimization probem (4.18) X (π ) = E 1 T <τ p T ()U(XT π ()) + 1 T τ K τ U(Xτ π ()(1 πτ ()γ τ )). Note that (4.18) can aso be defined for a strategy π that is predictabe ony on, (and not necessary on, ). Using Proposition 4.7, we have the foowing resut: Proposition 4.12 The strategies (πn = 1,τnˆπ 1 τn,τ δ) where ˆπ is the optima process for the driver of the BSDE Yt (1 + δγ τ ) p T τ T τ = 1 T <τ p T ()+1 t<τ T K τ + f (θ, Y θ, φ θ )dθ φ θ p dw θ, t, T, t t 2

21 are in A,δ and satisfy im X (π n + n) = V () = E (µ f (s, y, φ) = p sup s y + σsφ)ν 1 p ν δ 2 y νσ s 2 1 T <τ p T ()U(XT ˆπ ()) + 1 T τ K τ U(Xτ ˆπ ()(1 + δγ τ )). Proof: Let (τ n ) n N be an increasing sequence of F-predictabe stopping times that converge to. For any integer n 1, the strategy πn := 1,τnˆπ 1 τn,τ δ is in A,δ and πn converges to the non-predictabe optima strategy π np defined in (4.17) when n. Moreover, for any n N, X (πn) X (π np ) and by Proposition 4.7 X (π np ) im X (π n + n) V (ˆπ ). But the proof of Proposition 4.7 sti hods if we change the vaue at time of the portfoio π, thus the converse inequaity X (π np ) im n + X (πn) hods and E 1 T <τ p T ()U(XT ˆπ ()) + 1 T τ K τ U(Xτ ˆπ ()(1 + δγ τ )) = X (π np ) = X (1,τ ˆπ δ1 τ ) = im n + X (1,τnˆπ 1 τn, δ) = im n + X (π n) 5 Numerica iustrations We now iustrate our previous resuts by expicit modes and we aim to compare the optimization resuts for an insider and a standard investor. We reca that a investors start with an initia weath X. For the purpose of comparison, we choose a simiar mode with the one studied in 14. More precisey, we et the parameters µ, σ, γ to be constant, and µ 1 (θ), σ 1 (θ) are deterministic functions of θ given by (5.1) µ 1 (θ) = µ θ T, σ1 (θ) = σ (2 θ ), θ, T, T which means that the ratio of the after-defaut and before-defaut for the return rate of the asset is smaer than 1 and for the voatiity is arger than 1. Moreover, these ratios increases or decreases ineary with the defaut time respectivey: the after-defaut rate of return drops to zero, when the defaut time occurs near the initia date, and converges to the before-defaut rate of return, when the defaut time occurs near the finite investment horizon. For the voatiity, this ratio converges to the doube (resp. initia) vaue of the before-defaut voatiity, when the defaut time goes to the initia (resp. termina horizon) time. Moreover, in order to satisfy the hypothesis in the simuation part of 14, we have to assume that the defaut barrier L has no atoms (to ensure the density hypothesis, see Remark 4.3) and that L is independent of the fitration F (so that the defaut density is a deterministic function). In this case, p T (L) = 1. 21

22 Consider the CRRA utiity U(x) = xp p, < p < 1, the after-defaut vaue function is given from (4.9) by Vτ 1 (x) = K τ U(x) where K τ = ) ( (EZ T ( ) p 1 p p 1 1 = exp 2 p 1 p ( µ 1 ( ) ) 2(τ σ 1 T )) ( ) Furthermore, the soution of the before-defaut optimization probem is given by V () = Y U(X ) where Y is the soution of the BSDE (4.14) when etting φ =, i.e. (1 + δγ) p (5.2) Y t = 1 T <τ + 1 T τ K τ + p where f(t, y) = p ess sup ν A,δ,ν τ = δ τ T t f(θ, Y θ )dθ {µ ν t 1 p 2 (ν tσ ) 2 }y. We notice that in the case where the defaut time occurs after the maturity T, the optima strategy coincides with the cassica Merton strategy with constraint π δ, 1 γ. In the case where occurs before T, the process Y is stopped at, with the termina vaue depending on the quantity K τ, and the short-seing strategy δ at. We use an iterative agorithm to sove the equation (5.2). The foowing resuts are based on the mode parameters described beow: µ =.3, σ =.2, T = 1, the risk aversion parameter p =.8. For the standard investor, we use the deterministic mode as in 14 etting P(τ > t) = e λt, λ > so that the density function is α(θ) = λe λθ for a θ. Figure 1 compares the optima vaue function for insider, investor and Merton strategy. We fix the short-seing constraint δ =.5, the oss given defaut γ =.2 and the defaut intensity λ =.3. This corresponds to a reativey high risk of defaut. At the defaut time, the vaue function suffers a bruta oss for a the three strategies. The insider outperforms the other two strategies before and after the defaut occurs. Before the defaut, the vaue function for the standard investor is smaer than the Merton one because the atter does not consider at a the potentia defaut risk. However, when the defaut occurs, the investor outperforms the Merton strategy since the defaut risk is taken into account from the beginning. In Figure 2, we fix a smaer short-seing constraint δ =.1 and we study the impact of the oss given defaut γ. We observe that the vaue functions for both insider and investor are increasing with respect to the oss vaue γ. It is interesting to emphasize this phenomenon as a consequence of the short-seing where both insider and the investor wi profit the defaut event and obtain a arger vaue function. More precisey, in the eft-hand figure, we consider a reativey ow defaut risk (with the defaut intensity λ =.1) and we observe that the gain of the insider with respect to the investor remains stabe in time and aso for different vaues of γ. Whereas in the right-hand figure with a higher defaut risk (λ =.3), the impact of γ is more important for the investor: before the defaut the vaue function of the investor increases more significanty as γ increases, and there is no onger a drop at the defaut time. This can be expained by the fact that the investor has no imit for the short-seing strategy. 22

Appendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model

Appendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model Appendix of the Paper The Roe of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Mode Caio Ameida cameida@fgv.br José Vicente jose.vaentim@bcb.gov.br June 008 1 Introduction In this

More information

CS229 Lecture notes. Andrew Ng

CS229 Lecture notes. Andrew Ng CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view

More information

Homework 5 Solutions

Homework 5 Solutions Stat 310B/Math 230B Theory of Probabiity Homework 5 Soutions Andrea Montanari Due on 2/19/2014 Exercise [5.3.20] 1. We caim that n 2 [ E[h F n ] = 2 n i=1 A i,n h(u)du ] I Ai,n (t). (1) Indeed, integrabiity

More information

u(x) s.t. px w x 0 Denote the solution to this problem by ˆx(p, x). In order to obtain ˆx we may simply solve the standard problem max x 0

u(x) s.t. px w x 0 Denote the solution to this problem by ˆx(p, x). In order to obtain ˆx we may simply solve the standard problem max x 0 Bocconi University PhD in Economics - Microeconomics I Prof M Messner Probem Set 4 - Soution Probem : If an individua has an endowment instead of a monetary income his weath depends on price eves In particuar,

More information

Provisions estimation for portfolio of CDO in Gaussian financial environment

Provisions estimation for portfolio of CDO in Gaussian financial environment Technica report, IDE1123, October 27, 2011 Provisions estimation for portfoio of CDO in Gaussian financia environment Master s Thesis in Financia Mathematics Oeg Maximchuk and Yury Vokov Schoo of Information

More information

Mat 1501 lecture notes, penultimate installment

Mat 1501 lecture notes, penultimate installment Mat 1501 ecture notes, penutimate instament 1. bounded variation: functions of a singe variabe optiona) I beieve that we wi not actuay use the materia in this section the point is mainy to motivate the

More information

Coupling of LWR and phase transition models at boundary

Coupling of LWR and phase transition models at boundary Couping of LW and phase transition modes at boundary Mauro Garaveo Dipartimento di Matematica e Appicazioni, Università di Miano Bicocca, via. Cozzi 53, 20125 Miano Itay. Benedetto Piccoi Department of

More information

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is

More information

NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS

NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS NOISE-INDUCED STABILIZATION OF STOCHASTIC DIFFERENTIAL EQUATIONS TONY ALLEN, EMILY GEBHARDT, AND ADAM KLUBALL 3 ADVISOR: DR. TIFFANY KOLBA 4 Abstract. The phenomenon of noise-induced stabiization occurs

More information

Lecture Note 3: Stationary Iterative Methods

Lecture Note 3: Stationary Iterative Methods MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or

More information

Bourgain s Theorem. Computational and Metric Geometry. Instructor: Yury Makarychev. d(s 1, s 2 ).

Bourgain s Theorem. Computational and Metric Geometry. Instructor: Yury Makarychev. d(s 1, s 2 ). Bourgain s Theorem Computationa and Metric Geometry Instructor: Yury Makarychev 1 Notation Given a metric space (X, d) and S X, the distance from x X to S equas d(x, S) = inf d(x, s). s S The distance

More information

Manipulation in Financial Markets and the Implications for Debt Financing

Manipulation in Financial Markets and the Implications for Debt Financing Manipuation in Financia Markets and the Impications for Debt Financing Leonid Spesivtsev This paper studies the situation when the firm is in financia distress and faces bankruptcy or debt restructuring.

More information

A Brief Introduction to Markov Chains and Hidden Markov Models

A Brief Introduction to Markov Chains and Hidden Markov Models A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,

More information

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents

MARKOV CHAINS AND MARKOV DECISION THEORY. Contents MARKOV CHAINS AND MARKOV DECISION THEORY ARINDRIMA DATTA Abstract. In this paper, we begin with a forma introduction to probabiity and expain the concept of random variabes and stochastic processes. After

More information

XSAT of linear CNF formulas

XSAT of linear CNF formulas XSAT of inear CN formuas Bernd R. Schuh Dr. Bernd Schuh, D-50968 Kön, Germany; bernd.schuh@netcoogne.de eywords: compexity, XSAT, exact inear formua, -reguarity, -uniformity, NPcompeteness Abstract. Open

More information

4 Separation of Variables

4 Separation of Variables 4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE

More information

A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC

A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC (January 8, 2003) A NOTE ON QUASI-STATIONARY DISTRIBUTIONS OF BIRTH-DEATH PROCESSES AND THE SIS LOGISTIC EPIDEMIC DAMIAN CLANCY, University of Liverpoo PHILIP K. POLLETT, University of Queensand Abstract

More information

14 Separation of Variables Method

14 Separation of Variables Method 14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt

More information

4 1-D Boundary Value Problems Heat Equation

4 1-D Boundary Value Problems Heat Equation 4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x

More information

A. Distribution of the test statistic

A. Distribution of the test statistic A. Distribution of the test statistic In the sequentia test, we first compute the test statistic from a mini-batch of size m. If a decision cannot be made with this statistic, we keep increasing the mini-batch

More information

The distribution of the number of nodes in the relative interior of the typical I-segment in homogeneous planar anisotropic STIT Tessellations

The distribution of the number of nodes in the relative interior of the typical I-segment in homogeneous planar anisotropic STIT Tessellations Comment.Math.Univ.Caroin. 51,3(21) 53 512 53 The distribution of the number of nodes in the reative interior of the typica I-segment in homogeneous panar anisotropic STIT Tesseations Christoph Thäe Abstract.

More information

THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES

THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES THE REACHABILITY CONES OF ESSENTIALLY NONNEGATIVE MATRICES by Michae Neumann Department of Mathematics, University of Connecticut, Storrs, CT 06269 3009 and Ronad J. Stern Department of Mathematics, Concordia

More information

6 Wave Equation on an Interval: Separation of Variables

6 Wave Equation on an Interval: Separation of Variables 6 Wave Equation on an Interva: Separation of Variabes 6.1 Dirichet Boundary Conditions Ref: Strauss, Chapter 4 We now use the separation of variabes technique to study the wave equation on a finite interva.

More information

Separation of Variables and a Spherical Shell with Surface Charge

Separation of Variables and a Spherical Shell with Surface Charge Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation

More information

Explicit overall risk minimization transductive bound

Explicit overall risk minimization transductive bound 1 Expicit overa risk minimization transductive bound Sergio Decherchi, Paoo Gastado, Sandro Ridea, Rodofo Zunino Dept. of Biophysica and Eectronic Engineering (DIBE), Genoa University Via Opera Pia 11a,

More information

Another Class of Admissible Perturbations of Special Expressions

Another Class of Admissible Perturbations of Special Expressions Int. Journa of Math. Anaysis, Vo. 8, 014, no. 1, 1-8 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.1988/ijma.014.31187 Another Cass of Admissibe Perturbations of Specia Expressions Jerico B. Bacani

More information

Density models for credit risk

Density models for credit risk Density models for credit risk Nicole El Karoui, Ecole Polytechnique, France Monique Jeanblanc, Université d Évry; Institut Europlace de Finance Ying Jiao, Université Paris VII Workshop on stochastic calculus

More information

An approximate method for solving the inverse scattering problem with fixed-energy data

An approximate method for solving the inverse scattering problem with fixed-energy data J. Inv. I-Posed Probems, Vo. 7, No. 6, pp. 561 571 (1999) c VSP 1999 An approximate method for soving the inverse scattering probem with fixed-energy data A. G. Ramm and W. Scheid Received May 12, 1999

More information

(1 ) = 1 for some 2 (0; 1); (1 + ) = 0 for some > 0:

(1 ) = 1 for some 2 (0; 1); (1 + ) = 0 for some > 0: Answers, na. Economics 4 Fa, 2009. Christiano.. The typica househod can engage in two types of activities producing current output and studying at home. Athough time spent on studying at home sacrices

More information

On the evaluation of saving-consumption plans

On the evaluation of saving-consumption plans On the evauation of saving-consumption pans Steven Vanduffe Jan Dhaene Marc Goovaerts Juy 13, 2004 Abstract Knowedge of the distribution function of the stochasticay compounded vaue of a series of future

More information

Proceedings of the 2012 Winter Simulation Conference C. Laroque, J. Himmelspach, R. Pasupathy, O. Rose, and A. M. Uhrmacher, eds.

Proceedings of the 2012 Winter Simulation Conference C. Laroque, J. Himmelspach, R. Pasupathy, O. Rose, and A. M. Uhrmacher, eds. Proceedings of the 2012 Winter Simuation Conference C. aroque, J. Himmespach, R. Pasupathy, O. Rose, and A. M. Uhrmacher, eds. TIGHT BOUNDS FOR AMERICAN OPTIONS VIA MUTIEVE MONTE CARO Denis Beomestny Duisburg-Essen

More information

Indirect Optimal Control of Dynamical Systems

Indirect Optimal Control of Dynamical Systems Computationa Mathematics and Mathematica Physics, Vo. 44, No. 3, 24, pp. 48 439. Transated from Zhurna Vychisite noi Matematiki i Matematicheskoi Fiziki, Vo. 44, No. 3, 24, pp. 444 466. Origina Russian

More information

$, (2.1) n="# #. (2.2)

$, (2.1) n=# #. (2.2) Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier

More information

Physics 235 Chapter 8. Chapter 8 Central-Force Motion

Physics 235 Chapter 8. Chapter 8 Central-Force Motion Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion

More information

SVM: Terminology 1(6) SVM: Terminology 2(6)

SVM: Terminology 1(6) SVM: Terminology 2(6) Andrew Kusiak Inteigent Systems Laboratory 39 Seamans Center he University of Iowa Iowa City, IA 54-57 SVM he maxima margin cassifier is simiar to the perceptron: It aso assumes that the data points are

More information

Module 22: Simple Harmonic Oscillation and Torque

Module 22: Simple Harmonic Oscillation and Torque Modue : Simpe Harmonic Osciation and Torque.1 Introduction We have aready used Newton s Second Law or Conservation of Energy to anayze systems ike the boc-spring system that osciate. We sha now use torque

More information

More Scattering: the Partial Wave Expansion

More Scattering: the Partial Wave Expansion More Scattering: the Partia Wave Expansion Michae Fower /7/8 Pane Waves and Partia Waves We are considering the soution to Schrödinger s equation for scattering of an incoming pane wave in the z-direction

More information

Cryptanalysis of PKP: A New Approach

Cryptanalysis of PKP: A New Approach Cryptanaysis of PKP: A New Approach Éiane Jaumes and Antoine Joux DCSSI 18, rue du Dr. Zamenhoff F-92131 Issy-es-Mx Cedex France eiane.jaumes@wanadoo.fr Antoine.Joux@ens.fr Abstract. Quite recenty, in

More information

B. Brown, M. Griebel, F.Y. Kuo and I.H. Sloan

B. Brown, M. Griebel, F.Y. Kuo and I.H. Sloan Wegeerstraße 6 53115 Bonn Germany phone +49 8 73-347 fax +49 8 73-757 www.ins.uni-bonn.de B. Brown, M. Griebe, F.Y. Kuo and I.H. Soan On the expected uniform error of geometric Brownian motion approximated

More information

Primal and dual active-set methods for convex quadratic programming

Primal and dual active-set methods for convex quadratic programming Math. Program., Ser. A 216) 159:469 58 DOI 1.17/s117-15-966-2 FULL LENGTH PAPER Prima and dua active-set methods for convex quadratic programming Anders Forsgren 1 Phiip E. Gi 2 Eizabeth Wong 2 Received:

More information

Automobile Prices in Market Equilibrium. Berry, Pakes and Levinsohn

Automobile Prices in Market Equilibrium. Berry, Pakes and Levinsohn Automobie Prices in Market Equiibrium Berry, Pakes and Levinsohn Empirica Anaysis of demand and suppy in a differentiated products market: equiibrium in the U.S. automobie market. Oigopoistic Differentiated

More information

Reichenbachian Common Cause Systems

Reichenbachian Common Cause Systems Reichenbachian Common Cause Systems G. Hofer-Szabó Department of Phiosophy Technica University of Budapest e-mai: gszabo@hps.ete.hu Mikós Rédei Department of History and Phiosophy of Science Eötvös University,

More information

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l Probem 7. Simpe Penduum SEMINAR. PENDULUMS A simpe penduum means a mass m suspended by a string weightess rigid rod of ength so that it can swing in a pane. The y-axis is directed down, x-axis is directed

More information

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes

12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes Maima and Minima 1. Introduction In this Section we anayse curves in the oca neighbourhood of a stationary point and, from this anaysis, deduce necessary conditions satisfied by oca maima and oca minima.

More information

Disturbance decoupling by measurement feedback

Disturbance decoupling by measurement feedback Preprints of the 19th Word Congress The Internationa Federation of Automatic Contro Disturbance decouping by measurement feedback Arvo Kadmäe, Üe Kotta Institute of Cybernetics at TUT, Akadeemia tee 21,

More information

Age of Information: The Gamma Awakening

Age of Information: The Gamma Awakening Age of Information: The Gamma Awakening Eie Najm and Rajai Nasser LTHI, EPFL, Lausanne, Switzerand Emai: {eie.najm, rajai.nasser}@epf.ch arxiv:604.086v [cs.it] 5 Apr 06 Abstract Status update systems is

More information

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s

More information

Volume 13, MAIN ARTICLES

Volume 13, MAIN ARTICLES Voume 13, 2009 1 MAIN ARTICLES THE BASIC BVPs OF THE THEORY OF ELASTIC BINARY MIXTURES FOR A HALF-PLANE WITH CURVILINEAR CUTS Bitsadze L. I. Vekua Institute of Appied Mathematics of Iv. Javakhishvii Tbiisi

More information

An Algorithm for Pruning Redundant Modules in Min-Max Modular Network

An Algorithm for Pruning Redundant Modules in Min-Max Modular Network An Agorithm for Pruning Redundant Modues in Min-Max Moduar Network Hui-Cheng Lian and Bao-Liang Lu Department of Computer Science and Engineering, Shanghai Jiao Tong University 1954 Hua Shan Rd., Shanghai

More information

STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION

STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION Journa of Sound and Vibration (996) 98(5), 643 65 STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM G. ERDOS AND T. SINGH Department of Mechanica and Aerospace Engineering, SUNY at Buffao,

More information

Analysis of the problem of intervention control in the economy on the basis of solving the problem of tuning

Analysis of the problem of intervention control in the economy on the basis of solving the problem of tuning Anaysis of the probem of intervention contro in the economy on the basis of soving the probem of tuning arxiv:1811.10993v1 q-fin.gn] 9 Nov 2018 1 st eter Shnurov Nationa Research University Higher Schoo

More information

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische

More information

Lecture 6: Moderately Large Deflection Theory of Beams

Lecture 6: Moderately Large Deflection Theory of Beams Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey

More information

arxiv: v4 [math.pr] 1 Feb 2018

arxiv: v4 [math.pr] 1 Feb 2018 DIFFUSION TRANSFORMATIONS, BLACK-SCHOLES EQUATION AND OPTIMAL STOPPING UMUT ÇETIN arxiv:7.85v4 [math.pr] Feb 28 Abstract. We deveop a new cass of path transformations for one-dimensiona diffusions that

More information

Math 124B January 17, 2012

Math 124B January 17, 2012 Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia

More information

Integrating Factor Methods as Exponential Integrators

Integrating Factor Methods as Exponential Integrators Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been

More information

Asymptotic Properties of a Generalized Cross Entropy Optimization Algorithm

Asymptotic Properties of a Generalized Cross Entropy Optimization Algorithm 1 Asymptotic Properties of a Generaized Cross Entropy Optimization Agorithm Zijun Wu, Michae Koonko, Institute for Appied Stochastics and Operations Research, Caustha Technica University Abstract The discrete

More information

Convergence Property of the Iri-Imai Algorithm for Some Smooth Convex Programming Problems

Convergence Property of the Iri-Imai Algorithm for Some Smooth Convex Programming Problems Convergence Property of the Iri-Imai Agorithm for Some Smooth Convex Programming Probems S. Zhang Communicated by Z.Q. Luo Assistant Professor, Department of Econometrics, University of Groningen, Groningen,

More information

Problem set 6 The Perron Frobenius theorem.

Problem set 6 The Perron Frobenius theorem. Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator

More information

Symbolic models for nonlinear control systems using approximate bisimulation

Symbolic models for nonlinear control systems using approximate bisimulation Symboic modes for noninear contro systems using approximate bisimuation Giordano Poa, Antoine Girard and Pauo Tabuada Abstract Contro systems are usuay modeed by differentia equations describing how physica

More information

Expectation-Maximization for Estimating Parameters for a Mixture of Poissons

Expectation-Maximization for Estimating Parameters for a Mixture of Poissons Expectation-Maximization for Estimating Parameters for a Mixture of Poissons Brandon Maone Department of Computer Science University of Hesini February 18, 2014 Abstract This document derives, in excrutiating

More information

Uniprocessor Feasibility of Sporadic Tasks with Constrained Deadlines is Strongly conp-complete

Uniprocessor Feasibility of Sporadic Tasks with Constrained Deadlines is Strongly conp-complete Uniprocessor Feasibiity of Sporadic Tasks with Constrained Deadines is Strongy conp-compete Pontus Ekberg and Wang Yi Uppsaa University, Sweden Emai: {pontus.ekberg yi}@it.uu.se Abstract Deciding the feasibiity

More information

Alberto Maydeu Olivares Instituto de Empresa Marketing Dept. C/Maria de Molina Madrid Spain

Alberto Maydeu Olivares Instituto de Empresa Marketing Dept. C/Maria de Molina Madrid Spain CORRECTIONS TO CLASSICAL PROCEDURES FOR ESTIMATING THURSTONE S CASE V MODEL FOR RANKING DATA Aberto Maydeu Oivares Instituto de Empresa Marketing Dept. C/Maria de Moina -5 28006 Madrid Spain Aberto.Maydeu@ie.edu

More information

On Some Basic Properties of Geometric Real Sequences

On Some Basic Properties of Geometric Real Sequences On Some Basic Properties of eometric Rea Sequences Khirod Boruah Research Schoar, Department of Mathematics, Rajiv andhi University Rono His, Doimukh-791112, Arunacha Pradesh, India Abstract Objective

More information

Some Measures for Asymmetry of Distributions

Some Measures for Asymmetry of Distributions Some Measures for Asymmetry of Distributions Georgi N. Boshnakov First version: 31 January 2006 Research Report No. 5, 2006, Probabiity and Statistics Group Schoo of Mathematics, The University of Manchester

More information

King Fahd University of Petroleum & Minerals

King Fahd University of Petroleum & Minerals King Fahd University of Petroeum & Mineras DEPARTMENT OF MATHEMATICAL SCIENCES Technica Report Series TR 369 December 6 Genera decay of soutions of a viscoeastic equation Saim A. Messaoudi DHAHRAN 3161

More information

BASIC NOTIONS AND RESULTS IN TOPOLOGY. 1. Metric spaces. Sets with finite diameter are called bounded sets. For x X and r > 0 the set

BASIC NOTIONS AND RESULTS IN TOPOLOGY. 1. Metric spaces. Sets with finite diameter are called bounded sets. For x X and r > 0 the set BASIC NOTIONS AND RESULTS IN TOPOLOGY 1. Metric spaces A metric on a set X is a map d : X X R + with the properties: d(x, y) 0 and d(x, y) = 0 x = y, d(x, y) = d(y, x), d(x, y) d(x, z) + d(z, y), for a

More information

Componentwise Determination of the Interval Hull Solution for Linear Interval Parameter Systems

Componentwise Determination of the Interval Hull Solution for Linear Interval Parameter Systems Componentwise Determination of the Interva Hu Soution for Linear Interva Parameter Systems L. V. Koev Dept. of Theoretica Eectrotechnics, Facuty of Automatics, Technica University of Sofia, 1000 Sofia,

More information

Competitive Diffusion in Social Networks: Quality or Seeding?

Competitive Diffusion in Social Networks: Quality or Seeding? Competitive Diffusion in Socia Networks: Quaity or Seeding? Arastoo Fazei Amir Ajorou Ai Jadbabaie arxiv:1503.01220v1 [cs.gt] 4 Mar 2015 Abstract In this paper, we study a strategic mode of marketing and

More information

Input-to-state stability for a class of Lurie systems

Input-to-state stability for a class of Lurie systems Automatica 38 (2002) 945 949 www.esevier.com/ocate/automatica Brief Paper Input-to-state stabiity for a cass of Lurie systems Murat Arcak a;, Andrew Tee b a Department of Eectrica, Computer and Systems

More information

BALANCING REGULAR MATRIX PENCILS

BALANCING REGULAR MATRIX PENCILS BALANCING REGULAR MATRIX PENCILS DAMIEN LEMONNIER AND PAUL VAN DOOREN Abstract. In this paper we present a new diagona baancing technique for reguar matrix pencis λb A, which aims at reducing the sensitivity

More information

Equilibrium of Heterogeneous Congestion Control Protocols

Equilibrium of Heterogeneous Congestion Control Protocols Equiibrium of Heterogeneous Congestion Contro Protocos Ao Tang Jiantao Wang Steven H. Low EAS Division, Caifornia Institute of Technoogy Mung Chiang EE Department, Princeton University Abstract When heterogeneous

More information

Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.

Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n. Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier

More information

Pricing Multiple Products with the Multinomial Logit and Nested Logit Models: Concavity and Implications

Pricing Multiple Products with the Multinomial Logit and Nested Logit Models: Concavity and Implications Pricing Mutipe Products with the Mutinomia Logit and Nested Logit Modes: Concavity and Impications Hongmin Li Woonghee Tim Huh WP Carey Schoo of Business Arizona State University Tempe Arizona 85287 USA

More information

arxiv:math/ v2 [math.pr] 6 Mar 2005

arxiv:math/ v2 [math.pr] 6 Mar 2005 ASYMPTOTIC BEHAVIOR OF RANDOM HEAPS arxiv:math/0407286v2 [math.pr] 6 Mar 2005 J. BEN HOUGH Abstract. We consider a random wa W n on the ocay free group or equivaenty a signed random heap) with m generators

More information

Restricted weak type on maximal linear and multilinear integral maps.

Restricted weak type on maximal linear and multilinear integral maps. Restricted weak type on maxima inear and mutiinear integra maps. Oscar Basco Abstract It is shown that mutiinear operators of the form T (f 1,..., f k )(x) = R K(x, y n 1,..., y k )f 1 (y 1 )...f k (y

More information

An explicit Jordan Decomposition of Companion matrices

An explicit Jordan Decomposition of Companion matrices An expicit Jordan Decomposition of Companion matrices Fermín S V Bazán Departamento de Matemática CFM UFSC 88040-900 Forianópois SC E-mai: fermin@mtmufscbr S Gratton CERFACS 42 Av Gaspard Coriois 31057

More information

Week 6 Lectures, Math 6451, Tanveer

Week 6 Lectures, Math 6451, Tanveer Fourier Series Week 6 Lectures, Math 645, Tanveer In the context of separation of variabe to find soutions of PDEs, we encountered or and in other cases f(x = f(x = a 0 + f(x = a 0 + b n sin nπx { a n

More information

Haar Decomposition and Reconstruction Algorithms

Haar Decomposition and Reconstruction Algorithms Jim Lambers MAT 773 Fa Semester 018-19 Lecture 15 and 16 Notes These notes correspond to Sections 4.3 and 4.4 in the text. Haar Decomposition and Reconstruction Agorithms Decomposition Suppose we approximate

More information

The Construction of a Pfaff System with Arbitrary Piecewise Continuous Characteristic Power-Law Functions

The Construction of a Pfaff System with Arbitrary Piecewise Continuous Characteristic Power-Law Functions Differentia Equations, Vo. 41, No. 2, 2005, pp. 184 194. Transated from Differentsia nye Uravneniya, Vo. 41, No. 2, 2005, pp. 177 185. Origina Russian Text Copyright c 2005 by Izobov, Krupchik. ORDINARY

More information

Lecture Notes 4: Fourier Series and PDE s

Lecture Notes 4: Fourier Series and PDE s Lecture Notes 4: Fourier Series and PDE s 1. Periodic Functions A function fx defined on R is caed a periodic function if there exists a number T > such that fx + T = fx, x R. 1.1 The smaest number T for

More information

17 Lecture 17: Recombination and Dark Matter Production

17 Lecture 17: Recombination and Dark Matter Production PYS 652: Astrophysics 88 17 Lecture 17: Recombination and Dark Matter Production New ideas pass through three periods: It can t be done. It probaby can be done, but it s not worth doing. I knew it was

More information

Higher dimensional PDEs and multidimensional eigenvalue problems

Higher dimensional PDEs and multidimensional eigenvalue problems Higher dimensiona PEs and mutidimensiona eigenvaue probems 1 Probems with three independent variabes Consider the prototypica equations u t = u (iffusion) u tt = u (W ave) u zz = u (Lapace) where u = u

More information

Iterative Decoding Performance Bounds for LDPC Codes on Noisy Channels

Iterative Decoding Performance Bounds for LDPC Codes on Noisy Channels Iterative Decoding Performance Bounds for LDPC Codes on Noisy Channes arxiv:cs/060700v1 [cs.it] 6 Ju 006 Chun-Hao Hsu and Achieas Anastasopouos Eectrica Engineering and Computer Science Department University

More information

The EM Algorithm applied to determining new limit points of Mahler measures

The EM Algorithm applied to determining new limit points of Mahler measures Contro and Cybernetics vo. 39 (2010) No. 4 The EM Agorithm appied to determining new imit points of Maher measures by Souad E Otmani, Georges Rhin and Jean-Marc Sac-Épée Université Pau Veraine-Metz, LMAM,

More information

David Eigen. MA112 Final Paper. May 10, 2002

David Eigen. MA112 Final Paper. May 10, 2002 David Eigen MA112 Fina Paper May 1, 22 The Schrodinger equation describes the position of an eectron as a wave. The wave function Ψ(t, x is interpreted as a probabiity density for the position of the eectron.

More information

Unconditional security of differential phase shift quantum key distribution

Unconditional security of differential phase shift quantum key distribution Unconditiona security of differentia phase shift quantum key distribution Kai Wen, Yoshihisa Yamamoto Ginzton Lab and Dept of Eectrica Engineering Stanford University Basic idea of DPS-QKD Protoco. Aice

More information

Asynchronous Control for Coupled Markov Decision Systems

Asynchronous Control for Coupled Markov Decision Systems INFORMATION THEORY WORKSHOP (ITW) 22 Asynchronous Contro for Couped Marov Decision Systems Michae J. Neey University of Southern Caifornia Abstract This paper considers optima contro for a coection of

More information

arxiv: v1 [math.co] 17 Dec 2018

arxiv: v1 [math.co] 17 Dec 2018 On the Extrema Maximum Agreement Subtree Probem arxiv:1812.06951v1 [math.o] 17 Dec 2018 Aexey Markin Department of omputer Science, Iowa State University, USA amarkin@iastate.edu Abstract Given two phyogenetic

More information

Stochastic Automata Networks (SAN) - Modelling. and Evaluation. Paulo Fernandes 1. Brigitte Plateau 2. May 29, 1997

Stochastic Automata Networks (SAN) - Modelling. and Evaluation. Paulo Fernandes 1. Brigitte Plateau 2. May 29, 1997 Stochastic utomata etworks (S) - Modeing and Evauation Pauo Fernandes rigitte Pateau 2 May 29, 997 Institut ationa Poytechnique de Grenobe { IPG Ecoe ationae Superieure d'informatique et de Mathematiques

More information

Rate-Distortion Theory of Finite Point Processes

Rate-Distortion Theory of Finite Point Processes Rate-Distortion Theory of Finite Point Processes Günther Koiander, Dominic Schuhmacher, and Franz Hawatsch, Feow, IEEE Abstract We study the compression of data in the case where the usefu information

More information

Many-Help-One Problem for Gaussian Sources with a Tree Structure on their Correlation

Many-Help-One Problem for Gaussian Sources with a Tree Structure on their Correlation Many-Hep-One Probem for Gaussian Sources with a Tree Structure on their Correation Yasutada Oohama arxiv:090988v csit 4 Jan 009 Abstract In this paper we consider the separate coding probem for + correated

More information

1 Heat Equation Dirichlet Boundary Conditions

1 Heat Equation Dirichlet Boundary Conditions Chapter 3 Heat Exampes in Rectanges Heat Equation Dirichet Boundary Conditions u t (x, t) = ku xx (x, t), < x (.) u(, t) =, u(, t) = u(x, ) = f(x). Separate Variabes Look for simpe soutions in the

More information

FRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA)

FRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA) 1 FRST 531 -- Mutivariate Statistics Mutivariate Discriminant Anaysis (MDA) Purpose: 1. To predict which group (Y) an observation beongs to based on the characteristics of p predictor (X) variabes, using

More information

Statistical Inference, Econometric Analysis and Matrix Algebra

Statistical Inference, Econometric Analysis and Matrix Algebra Statistica Inference, Econometric Anaysis and Matrix Agebra Bernhard Schipp Water Krämer Editors Statistica Inference, Econometric Anaysis and Matrix Agebra Festschrift in Honour of Götz Trenker Physica-Verag

More information

Lecture 6 Povh Krane Enge Williams Properties of 2-nucleon potential

Lecture 6 Povh Krane Enge Williams Properties of 2-nucleon potential Lecture 6 Povh Krane Enge Wiiams Properties of -nuceon potentia 16.1 4.4 3.6 9.9 Meson Theory of Nucear potentia 4.5 3.11 9.10 I recommend Eisberg and Resnik notes as distributed Probems, Lecture 6 1 Consider

More information

Akaike Information Criterion for ANOVA Model with a Simple Order Restriction

Akaike Information Criterion for ANOVA Model with a Simple Order Restriction Akaike Information Criterion for ANOVA Mode with a Simpe Order Restriction Yu Inatsu * Department of Mathematics, Graduate Schoo of Science, Hiroshima University ABSTRACT In this paper, we consider Akaike

More information

Schedulability Analysis of Deferrable Scheduling Algorithms for Maintaining Real-Time Data Freshness

Schedulability Analysis of Deferrable Scheduling Algorithms for Maintaining Real-Time Data Freshness 1 Scheduabiity Anaysis of Deferrabe Scheduing Agorithms for Maintaining Rea-Time Data Freshness Song Han, Deji Chen, Ming Xiong, Kam-yiu Lam, Aoysius K. Mok, Krithi Ramamritham UT Austin, Emerson Process

More information

ON THE POSITIVITY OF SOLUTIONS OF SYSTEMS OF STOCHASTIC PDES

ON THE POSITIVITY OF SOLUTIONS OF SYSTEMS OF STOCHASTIC PDES ON THE POSITIVITY OF SOLUTIONS OF SYSTEMS OF STOCHASTIC PDES JACKY CRESSON 1,2, MESSOUD EFENDIEV 3, AND STEFANIE SONNER 3,4 On the occasion of the 75 th birthday of Prof. Dr. Dr.h.c. Wofgang L. Wendand

More information

Numerical solution of one dimensional contaminant transport equation with variable coefficient (temporal) by using Haar wavelet

Numerical solution of one dimensional contaminant transport equation with variable coefficient (temporal) by using Haar wavelet Goba Journa of Pure and Appied Mathematics. ISSN 973-1768 Voume 1, Number (16), pp. 183-19 Research India Pubications http://www.ripubication.com Numerica soution of one dimensiona contaminant transport

More information