A new combination of techniques for the numerical computation of incompressible flow is presented. The temporal discretization bases on the discontinuous Galerkin-formulation. Both constant (DG(0)) and linear approximation (DG(1)) in time is discussed. In case of DG(1) an iterative method reduces the problem to a sequence of problems each with the dimension of the DG(0) approach. For the semi-discrete problems a Galerkin/least-squares method is applied. Furthermore a non-overlapping domain decomposition method can be used for a parallelized computation. The main advantage of this approach is the low amount of information which must be exchanged between the subdomains. Due to the slight bandwidth a workstation-cluster is a suitable platform. Otherwise this method is efficient only for a small number of subdomains. The interface condition is of the Robin/Robin-type and for the Navier-Stokes equation a formulation introducing a further pressure interface condition is used. Additionally a suggestion for the implementation of the standard k-epsilon turbulence model with special wall function is done in this context. All the features mentioned above are implemented in a code called ParallelNS. Using this code the verification of this approach was done on a large number of examples ranging from simple advection-diffusion problems to turbulent convection in a closed cavity.

許多公司在發行公司債的時候，會給此公司債一個可提前贖回的特性，此種公司債稱為可贖回公司債(Callable Bond)，用來規避利率變動風險的金融商品也與我們熟知的利率交換不同，稱為可取消利率交換(Cancelable Swap)。其實可取消利率交換可以拆解成百慕達利率交換選擇權(Bermudan Swaption)加上利率交換，由於利率交換之評價較簡單也有市場一致的評價方法，因此百慕達利率交換選擇權便成為評價的重點。 評價的部分，由於百慕達式的商品有提前履約的特性，造成其封閉解不存在，因此需要利用其他的近似解或是數值方法來求它的價格。由於本文採用BGM(1997)的市場利率模型(Libor Market Model)，其高維度的性質導致數狀方法與有限差分法使用起來較無效率，因此本文選擇使用蒙地卡羅法做為評價的方法，同時利用Longstaff and Schwartz(2001)的最小平方蒙地卡羅法(Least Squares Monte Carlo Method)來解決提前履約的問題。 最後，本文將採用2種利率波動度假設與2種不同利率間相關係數的假設，共4種組合，在歐式利率交換選擇權的市場波動度下進行校準，使用校準出來的參數進行評價來得到4種價格。再進行商品的敏感度分析(Sensitivity Analysis)和風險值(Value at Risk)的計算。

No description available.

Spacecraft missions such as spacecraft docking and formation flying require high precision relative position and attitude data. Although Global Positioining Systems can provide this capability near the earth, deep space missions require the use of alternative technologies. One such technology is the vision-based navigation (VISNAV) sensor system developed at Texas A&M University. VISNAV comprises an electro-optical sensor combined with light sources or beacons. This patented sensor has an analog detector in the focal plane with a rise time of a few microseconds. Accuracies better than one part in 2000 of the field of view have been obtained. This research presents a new approach involving simultaneous activation of beacons with frequency division multiplexing as part of the VISNAV sensor system. In addition, it discusses the synchronous demodulation process using digital heterodyning and decimating filter banks on a low-power fixed point DSP, which improves the accuracy of the sensor measurements and the reliability of the system. This research also presents an optimal and computationally efficient six-degree-of-freedom estimation algorithm using a new measurement model based on the attitude representation of Modified Rodrigues Parameters.

In many fields of science, engineering, and economics large amounts of data are stored and there is a need to analyze these data in order to extract information for various purposes. Data mining is a general concept involving different tools for performing this kind of analysis. The development of mathematical models and efficient algorithms is of key importance. In this thesis we discuss algorithms for the reduced rank regression problem and algorithms for the computation of the best multilinear rank approximation of tensors. The first two papers deal with the reduced rank regression problem, which is encountered in the field of state-space subspace system identification. More specifically the problem is \[ \min_{\rank(X) = k} \det (B - X A)(B - X A)\tp, \] where $A$ and $B$ are given matrices and we want to find $X$ under a certain rank condition that minimizes the determinant. This problem is not properly stated since it involves implicit assumptions on $A$ and $B$ so that $(B - X A)(B - X A)\tp$ is never singular. This deficiency of the determinant criterion is fixed by generalizing the minimization criterion to rank reduction and volume minimization of the objective matrix. The volume of a matrix is defined as the product of its nonzero singular values. We give an algorithm that solves the generalized problem and identify properties of the input and output signals causing a singular objective matrix. Classification problems occur in many applications. The task is to determine the label or class of an unknown object. The third paper concerns with classification of handwritten digits in the context of tensors or multidimensional data arrays. Tensor and multilinear algebra is an area that attracts more and more attention because of the multidimensional structure of the collected data in various applications. Two classification algorithms are given based on the higher order singular value decomposition (HOSVD). The main algorithm makes a data reduction using HOSVD of 98--99 \% prior the construction of the class models. The models are computed as a set of orthonormal bases spanning the dominant subspaces for the different classes. An unknown digit is expressed as a linear combination of the basis vectors. The resulting algorithm achieves 5\% in classification error with fairly low amount of computations. The remaining two papers discuss computational methods for the best multilinear rank approximation problem \[ \min_{\cB} \| \cA - \cB\| \] where $\cA$ is a given tensor and we seek the best low multilinear rank approximation tensor $\cB$. This is a generalization of the best low rank matrix approximation problem. It is well known that for matrices the solution is given by truncating the singular values in the singular value decomposition (SVD) of the matrix. But for tensors in general the truncated HOSVD does not give an optimal approximation. For example, a third order tensor $\cB \in \RR^{I \x J \x K}$ with rank$(\cB) = (r_1,r_2,r_3)$ can be written as the product \[ \cB = \tml{X,Y,Z}{\cC}, \qquad b_{ijk}=\sum_{\lambda,\mu,\nu} x_{i\lambda} y_{j\mu} z_{k\nu} c_{\lambda\mu\nu}, \] where $\cC \in \RR^{r_1 \x r_2 \x r_3}$ and $X \in \RR^{I \times r_1}$, $Y \in \RR^{J \times r_2}$, and $Z \in \RR^{K \times r_3}$ are matrices of full column rank. Since it is no restriction to assume that $X$, $Y$, and $Z$ have orthonormal columns and due to these constraints, the approximation problem can be considered as a nonlinear optimization problem defined on a product of Grassmann manifolds. We introduce novel techniques for multilinear algebraic manipulations enabling means for theoretical analysis and algorithmic implementation. These techniques are used to solve the approximation problem using Newton and Quasi-Newton methods specifically adapted to operate on products of Grassmann manifolds. The presented algorithms are suited for small, large and sparse problems and, when applied on difficult problems, they clearly outperform alternating least squares methods, which are standard in the field.

本文實證探討黃金期貨報酬率的特性並在標的黃金期貨價格遵循狀態轉換跳躍擴散過程時實現美式選擇權之評價。在這樣的動態過程下，跳躍事件被一個複合普瓦松過程與對數常態跳躍振幅所描述，以及狀態轉換到達強度是由一個其狀態代表經濟狀態的隱藏馬可夫鏈所捕捉。考量不同的跳躍風險假設，我們使用Merton測度與Esscher轉換推導出在一個不完全市場設定下的風險中立黃金期貨價格動態過程。為了達到所需的精確度，最小平方蒙地卡羅法被用來近似美式黃金期貨選擇權的價值。基於實際市場資料，我們提供實證與數值結果來說明這個動態模型的優點。 / This dissertation empirically investigates the characteristics of gold futures returns and achieves the valuation of American-style options when the underlying gold futures price follows a regime-switching jump-diffusion process. Under such dynamics, the jump events are described as a compound Poisson process with a log-normal jump amplitude, and the regime-switching arrival intensity is captured by a hidden Markov chain whose states represent the economic states. Considering the different jump risk assumptions, we use the Merton measure and Esscher transform to derive risk-neutral gold futures price dynamics under an incomplete market setting. To achieve a desired accuracy level, the least-squares Monte Carlo method is used to approximate the values of American gold futures options. Our empirical and numerical results based on actual market data are provided to illustrate the advantages of this dynamic model.

進行調查時，最後的抽樣單位常是從不同的群集取得的，而同一群集內的樣本對象，因背景類似而對於某些問題常會傾向相同或類似的反應，研究者若忽略這種群內相關性，仍以獨立性樣本進行分析時，因其共變異數矩陣通常會與多項模式的共變異數矩陣相差懸殊，而造成所謂的額外變異或不足變異的現象。本文在不同的情況下，提出了Dirichlet-Multinomial模式(簡稱DM模式)、擴展的DM模式、以及兩種平均數-共變異數矩陣模式，以適當的彙整所有的群集資料。並討論DM與EDM模式中相關之參數及格機率之最大概似估計法，且分別對此兩種平均數-共變異數矩陣模式，提出求導廣義最小平方估計的程序。此外，也針對幾種特殊的二維表及三維表結構，探討對應的參數及格機率之估計方法。並提出計算簡易的Score統計檢定量以判斷群內相關(intra-cluster correlation)之存在性，及判斷資料集具有額外或不足變異，而對於不同母體的群內相關同質性檢定亦提出討論。 / This paper presents a modelling method of analyzing categorical data with overdispersion or underdispersion.　In many studies, data are collected from differ clusters, and members within the same cluster behave similary.　Thus, the responses of members within the same cluster are not independent and the multinomial distribution is not the correct distribution for the observed counts.　Therefore, the covariance matrix of the sample proportion vector tends to be much different from that of the multinomial model.　We discuss four different models to fit counts data with overdispersion or underdispersion feature, witch include Dirichlet-Multinomial model (DM model), extended DM model (EDM model), and two mean-covariance models.　Method of maximum-likelihood estimation is discussed for DM and EDM models.　Procedures to derive generalized least squares estimates are proposed for the two mean-covariance models respectively. As to the cell probabilities, we also discuss how to estimate them under several special structures of　two-way and three-way tables.　More easily evaluated Score test statistics are derived for the DM and EDM models to test the existence of　the intra-cluster correlation.　And the test of homogeneity of intra-cluster correlation among several populations is also derived.

微陣列晶片是發展及應用較為成熟的生物晶片技術。由於微陣列實驗程序複雜，故資料常包含多種不同來源的實驗誤差，為了適當的區分實驗中來自處理、晶片及基因的效應，我們提出混合效應變異數分析模型來調整系統誤差。針對各基因在不同實驗環境的差異性假設檢定問題，利用最小平方法推導出點估計以及對應的檢定統計量。本研究介紹多重檢定問題中的族型一誤差，並證明在此模型下，Sidak調整法為適當的多重檢定方法。在給定族型一誤差率的顯著水準，利用檢定力的公式，運算出在預設檢定力的最低水準下所需最小樣本（晶片）數。最後我們透過電腦模擬，以蒙地卡羅法來估計檢定力與族型一誤差率，由模擬結果發現，採用此最小樣本數結果，其檢定力可達到預期的水準以上，並且其族型一誤差率皆適當地控制在顯著水準以內。

可轉換公司債是一種在持有期間內，投資人可以在規定的時間內將債券轉換為股票，或是到期時得到債券報酬的一種複合式證券。因此，可轉債除了具有債券性質之外，還包含另一部份可視為一美式選擇權的股票選擇權。 本篇論文將可轉換債券評價結合數值分析中的最小蒙地卡羅法(Least square monte carlo)，使得在評價可轉債時，能夠具有更多的彈性處理發行公司自行設計的贖回條款與其他各種不同的契約情況。 此外，本篇論文針對股價考慮跳躍的性質，使用Compound Poisson 過程模擬發生跳躍的次數，導入Merton的跳躍模型(Jump-diffusion Model)，在Merton的假設下，模擬未來股價的動態變化。 信用風險方面，本文採用Duffie提出的風險CIR模型評價。考慮存活函數(Survival Function)和違約強度(Hazard Rate Function)，使用CIR模型描述信用違約強度在可轉債持有期間的動態變化，最後模擬出違約的時點，結合LSMC下的可轉債評價評價法。 最後利率部份，雖然Brennan and Schwartz(1980)認為隨機利率對於可轉換債券的評價，並沒有明顯的效果，反而會降低評價時的效率，但是為了符合評價過程的合理性，本文使用CIR短期利率模型。

This thesis concerns the mathematical modelling of phase transformations with a special emphasis on martensitic phase transformations and their application to the modelling of steels. In Chapter 1, we develop a framework that determines the optimal transformation strain between any two Bravais lattices and use it to give a rigorous proof of a conjecture by E.C. Bain in 1924 on the optimality of the so-called Bain strain. In Chapter 2, we review the Ball-James model and related concepts. We present some simplification of existing results. In Chapter 3, we pose a conjecture for the explicit form of the quasiconvex hull of the three tetragonal wells, known as the three-well problem. We present a new approach to finding inner and outer bounds. In Chapter 4, we focus on highly compatible, so called self-accommodating, martensitic structures and present new results on their fine properties such as estimates on their minimum complexity and bounds on the relative proportion of each martensitic variant in them. In Chapter 5, we investigate the contrary situation when self-accommodating microstructures do not exist. We determine, whether in this situation, it is still energetically favourable to nucleate martensite within austenite. By constructing different types of inclusions, we find that the optimal shape of an inclusion is flat and thin which is in agreement with experimental observation. In Chapter 6, we introduce a mechanism that identifies transformation strains with orientation relationships. This mechanism allows us to develop a simpler, strain-based approach to phase transformation models in steels. One novelty of this approach is the derivation of an explicit dependence of the orientation relationships on the ratio of tetragonality of the product phase. In Chapter 7, we establish a correspondence between common phenomenological models for steels and the Ball-James model. This correspondence is then used to develop a new theory for the (5 5 7) lath transformation in low-carbon steels. Compared to existing theories, this new approach requires a significantly smaller number of input parameters. Furthermore, it predicts a microstructure morphology which differs from what is conventionally believed.