线性方程组迭代解的随机模型测试研究

本文讨论大型线性方程组迭代解的随机模型测试评价问题。给出了常用迭代解法ＣＧ、ＬＳＱＲ、ＳＩＲＴ、ＳＡＲＴ、ＳＡＳＩＲＴ等的测试结果。结果表明：（１）方程组系数矩阵的特性（条件数）及解结构都对解精度有重要影响。解模型越粗糙,解的精度越低。（２）各种求解算法都有一定的平滑效应,同时各种算法也都会产生误差大于２００％的奇异解,奇异解元素数一般约占１０％。（３）数据的拟合残差一般不能真实反映解的精度。（４）对含误差数据的求解问题,较好的求解算法是ＤＬＳＱＲ与ＳＡＳＩＲＴ。     

The structure of the mid- and high-latitude ionosphere during the November 2004 storm event obtained from GPS observations

GPS data from the International GNSS Service (IGS) network were used to study the development of the severe geomagnetic storm of November 7–12, 2004, in the total electron content (TEC) on a global scale. The TEC maps were produced for analyzing the storm. For producing the maps over European and North American sectors, GPS measurements from more than 100 stations were used. The dense network of GPS stations provided TEC measurements with a high temporal and spatial resolution. To present the temporal and spatial variation of TEC during the storm, differential TEC maps relative to a quiet day (November 6, 2004) were created. The features of geomagnetic storm attributed to the complex development of ionospheric storm depend on latitude, longitude and local time. The positive, as well as negative effects were detected in TEC variations as a consequence of the evolution of the geomagnetic storm. The maximal effect was registered in the subauroral/auroral ionosphere during substorm activity in the evening and night period. The latitudinal profiles obtained from TEC maps for Europe gave rise to the storm-time dynamic of the ionospheric trough, which was detected on November 7 and 9 at latitudes below 50°N. In the report, features of the response of TEC to the storm for European and North American sectors are analyzed.     

基于类柔度差曲率和频率摄动的结构损伤识别

为提高梁式结构损伤诊断的效率,提出一种基于类柔度差曲率和频率摄动的结构损伤识别方法。首先根据结构振动理论,研究广义柔度矩阵计算公式;再利用模态柔度对结构损伤灵敏性高的优点,改进基于柔度差曲率的损伤定位指标,定义类柔度差曲率LCFC损伤指标,并初步识别损伤;最后基于矩阵摄动进行结构损伤识别结果确认。考虑多种损伤工况,对一简支梁结构进行损伤识别数值模拟验证。结果表明:仅使用一阶模态,建立的类柔度差曲率LCFC指标对梁式结构损伤定位具有良好的诊断效果,且计算工作量小;对于含边界损伤单元的多损伤工况,当损伤程度大于10%时,LCFC指标识别有效;当损伤程度不大于25%时,各工况二阶摄动识别结果精度较高,相对误差较一阶摄动结果明显降低,证明了该方法的实用性、有效性和精确性。     

利用二维DLT及光束法平差进行数字摄像机标定

利用共线方程和二维DLT之间的对应关系导出了由二维DLT的8个参数表达的主纵线方程式,讨论了主点初值的求解方法,给出了单张像片摄像机参数分解不惟一性及临界序列的证明。详细推导了利用二维DLT参数分解摄像机外方位元素初值的实用算法,论述了利用光束法平差进行摄像机标定的数学模型。实际图像数据实验取得了很好的结果,验证了本文所提出的摄像机标定算法的可行性。     

The effect of tidal inlets on open coast storm surge hydrographs

Bridge scour modeling requires storm surge hydrographs as open ocean boundary conditions for coastal waters surrounding tidal inlets. These open coast storm surge hydrographs are used to accurately determine both horizontal and vertical circulation patterns, and thus scour, within the inlet and bay for an extreme event. At present, very little information is available on the effect that tidal inlets have on these open coast storm surge hydrographs. Furthermore, current modeling practice enforces a single design hydrograph along the open coast boundary for bridge scour models. This study expands on these concepts and provides a more fundamental understanding on both of these modeling areas.     

A Compact Explicit Difference Scheme of High Accuracy for Extended Boussinesq Equations

Presented here is a compact explicit difference scheme of high accuracy for solving the extended Boussinesq equations.For time discretization,a three-stage explicit Runge-Kutta method with TVD property is used at predicting stage,a cubic spline function is adopted at correcting stage,which made the time discretization accuracy up to fourth order;For spatial discretization,a three-point explicit compact difference scheme with arbitrary order accuracy is employed.The extended Boussinesq equations derived by Beji and Nadaoka are solved by the proposed scheme.The numerical results agree well with the experimental data.At the same time,the comparisons of the two numerical results between the present scheme and low accuracy difference method are made,which further show the necessity of using high accuracy scheme to solve the extended Boussinesq equations.As a valid sample,the wave propagation on the rectangular step is formulated by the present scheme,the modelled results are in better agreement with the experimental data than those of Kittitanasuan.     

Numerical simulation of wave-induced currents combined with parabolic mild-slope equation in curvilinear coordinates

Researches on breaking-induced currents by waves are summarized firstly in this paper. Then, a combined numerical model in orthogonal curvilinear coordinates is presented to simulate wave-induced current in areas with curved boundary or irregular coastline. The proposed wave-induced current model includes a nearshore current module established through orthogonal curvilinear transformation form of shallow water equations and a wave module based on the curvilinear parabolic approximation wave equation. The wave module actually serves as the driving force to provide the current module with required radiation stresses. The Crank-Nicolson finite difference scheme and the alternating directions implicit method are used to solve the wave and current module, respectively. The established surf zone currents model is validated by two numerical experiments about longshore currents and rip currents in basins with rip channel and breakwater. The numerical results are compared with the measured data and published numerical results.     

关于二阶伴随模型的理论研究

Hesse矩阵-目标函数关于控制变量的二阶偏导数形成的矩阵,在变分数据同化过程中以及敏感性分析等方面起着重要的作用;它可以通过建立数学模型的一阶和二阶伴随模型求得.以浅水方程模型为例,利用泛函的Gâteaux微分和Hilbert空间上伴随算子的概念,导出了一阶和二阶伴随模型并由此得到Hesse矩阵.改进了Zhi Wang等(1992)建立的二阶伴随模型理论.     

斜压海洋动力学的一种三维数值模式──Ⅱ.温度、盐度和垂直涡动粘性系数的计算

根据海水温度和盐度平流扩散方程给出一种数值计算方案并采用混合长度理论给出垂直涡动粘性系数的计算方法。对于温度和盐度方程,其平流过程采用了Lax－Wendroff格式,水平扩散采用显格式,垂直扩散采用隐格式。时间步长主要受平流过程的Courant-Friedrichs-Lewy条件限制。垂直涡动粘性系数计算依据Prandtl混合长度理论,并考虑了海水层化的抑制作用,因而其数值与流场及密度场结构有关。温度、盐度及垂直涡动粘性系数的计算与动力方程中内模态的计算同步进行。应用本模式模拟渤、黄、东海由潮流、密度流和风海流迭加而成的综合海流,得出了良好的结果。     

New finite-gap solutions for the coupled Burgers equations engendered by the Neumann systems

On the tangent bundle TS^N-1 of the unit sphere S^N-1, this paper reduces the coupled Burgers equations to two Neumann systems by using the nonlinearization of the Lax pair, whose Liouville integrability is displayed in the scheme of the r-matrix technique. Based on the Lax matrix of the Neumann systems, the Abel--Jacobi coordinates are appropriately chosen to straighten out the restricted Neumann flows on the complex torus, from which the new finite-gap solutions expressed by Riemann theta functions for the coupled Burgers equations are given in view of the Jacobi inversion.