Statistical Interpolation Method of Turbulent Phase Screen

A relative displacement between the grid points of optical fields and those of phase screens may occur in the simulation of light propagation through the turbulent atmosphere. A statistical interpolator is proposed to solve this problem in this paper. It is evaluated by the phase structure function and numerical experiments of light propagation through atmospheric turbulence with/without adaptive optics (AO) and it is also compared with the well-known linear interpolator under the same condition. Results of the phase structure function show that the statistical interpolator is more accurate in comparison with the linear one, especially in the high frequency region. More importantly, the long-exposure results of light propagation through the turbulent atmosphere with/without AO also show that the statistical interpolator is more accurate and reliable than the linear one. (C) 2009 Optical Society of America.

改进叠前深度偏移的计算效率＆叠前高分辨率成像

An high-resolution prestack imaging technique of seismic data is developed in this thesis. By using this technique, the reflected coefficients of sheet sands can be gained in order to understand and identify thin oil reservoirs. One-way wave equation based migration methods can more accurately model seismic wave propagation effect such as multi-arrivals and obtain almost correct reflected energy in the presence of complex inhomogeneous media, and therefore, achieve more superiorities in imaging complex structure. So it is a good choice to apply the proposed high-resolution imaging to the presatck depth migration gathers. But one of the main shorting of one-way wave equation based migration methods is the low computational efficiency, thus the improvement on computational efficiency is first carried out. The method to improve the computational efficiency of prestack depth migration is first presented in this thesis, that is frequency-dependent varying-step depth exploration scheme plus a table-driven, one-point wavefield interpolation technology for wave equation based migration methods; The frequency-dependent varying-step depth exploration scheme reduces the computational cost of wavefield depth extrapolation, and the a table-driven, one-point wavefield interpolation technology reconstructs the extrapolated wavefield with an equal, desired vertical step with high computational efficiency. The proposed varying-step depth extrapolation plus one-point interpolation scheme results in 2/3 reduction in computational cost when compared to the equal-step depth extrapolation of wavefield, but gives the almost same imaging. The frequency-dependent varying-step depth exploration scheme is presented in theory by using the optimum split-step Fourier. But the proposed scheme can also be used by other wave equation based migration methods of the frequency domain. The proposed method is demonstrated by using impulse response, 2-D Marmousi dataset, 3-D salt dataset and the 3-D field dataset. A method of high-resolution prestack imaging is presented in the 2nd part of this thesis. The seismic interference method to solve the relative reflected coefficients is presented. The high-resolution imaging is obtained by introducing a sparseness- constrained least-square inversion into the reflected coefficient imaging. Gaussian regularization is first imposed and a smoothed solution is obtained by solving equation derived from the least-square inversion. Then the Cauchy regularization is introducing to the least-square inversion , the sparse solution of relative reflected coefficients can be obtained, that is high-resolution solution. The proposed scheme can be used together with other prestack imaging if the higher resolution is needed in a target zone. The seismic interference method in theory and the solution to sparseness-constrained least-square inversion are presented. The proposed method is demonstrated by synthetic examples and filed data.

基于显式分形插值的高保真地震数据重建方法研究

Based on the fractal theories, contractive mapping principles as well as the fixed point theory, by means of affine transform, this dissertation develops a novel Explicit Fractal Interpolation Function（EFIF）which can be used to reconstruct the seismic data with high fidelity and precision. Spatial trace interpolation is one of the important issues in seismic data processing. Under the ideal circumstances, seismic data should be sampled with a uniform spatial coverage. However, practical constraints such as the complex surface conditions indicate that the sampling density may be sparse or for other reasons some traces may be lost. The wide spacing between receivers can result in sparse sampling along traverse lines, thus result in a spatial aliasing of short-wavelength features. Hence, the method of interpolation is of very importance. It not only needs to make the amplitude information obvious but the phase information, especially that of the point that the phase changes acutely. Many people put forward several interpolation methods, yet this dissertation focuses attention on a special class of fractal interpolation function, referred to as explicit fractal interpolation function to improve the accuracy of the interpolation reconstruction and to make the local information obvious. The traditional fractal interpolation method mainly based on the randomly Fractional Brown Motion (FBM) model, furthermore, the vertical scaling factor which plays a critical role in the implementation of fractal interpolation is assigned the same value during the whole interpolating process, so it can not make the local information obvious. In addition, the maximal defect of the traditional fractal interpolation method is that it cannot obtain the function values on each interpolating nodes, thereby it cannot analyze the node error quantitatively and cannot evaluate the feasibility of this method. Detailed discussions about the applications of fractal interpolation in seismology have not been given by the pioneers, let alone the interpolating processing of the single trace seismogram. On the basis of the previous work and fractal theory this dissertation discusses the fractal interpolation thoroughly and the stability of this special kind of interpolating function is discussed, at the same time the explicit presentation of the vertical scaling factor which controls the precision of the interpolation has been proposed. This novel method develops the traditional fractal interpolation method and converts the fractal interpolation with random algorithms into the interpolation with determined algorithms. The data structure of binary tree method has been applied during the process of interpolation, and it avoids the process of iteration that is inevitable in traditional fractal interpolation and improves the computation efficiency. To illustrate the validity of the novel method, this dissertation develops several theoretical models and synthesizes the common shot gathers and seismograms and reconstructs the traces that were erased from the initial section using the explicit fractal interpolation method. In order to compare the differences between the theoretical traces that were erased in the initial section and the resulting traces after reconstruction on waveform and amplitudes quantitatively, each missing traces are reconstructed and the residuals are analyzed. The numerical experiments demonstrate that the novel fractal interpolation method is not only applicable to reconstruct the seismograms with small offset but to the seismograms with large offset. The seismograms reconstructed by explicit fractal interpolation method resemble the original ones well. The waveform of the missing traces could be estimated very well and also the amplitudes of the interpolated traces are a good approximation of the original ones. The high precision and computational efficiency of the explicit fractal interpolation make it a useful tool to reconstruct the seismic data; it can not only make the local information obvious but preserve the overall characteristics of the object investigated. To illustrate the influence of the explicit fractal interpolation method to the accuracy of the imaging of the structure in the earth’s interior, this dissertation applies the method mentioned above to the reverse-time migration. The imaging sections obtained by using the fractal interpolated reflected data resemble the original ones very well. The numerical experiments demonstrate that even with the sparse sampling we can still obtain the high accurate imaging of the earth’s interior’s structure by means of the explicit fractal interpolation method. So we can obtain the imaging results of the earth’s interior with fine quality by using relatively small number of seismic stations. With the fractal interpolation method we will improve the efficiency and the accuracy of the reverse-time migration under economic conditions. To verify the application effect to real data of the method presented in this paper, we tested the method by using the real data provided by the Broadband Seismic Array Laboratory, IGGCAS. The results demonstrate that the accuracy of explicit fractal interpolation is still very high even with the real data with large epicenter and large offset. The amplitudes and the phase of the reconstructed station data resemble the original ones that were erased in the initial section very well. Altogether, the novel fractal interpolation function provides a new and useful tool to reconstruct the seismic data with high precision and efficiency, and presents an alternative to image the deep structure of the earth accurately.

电离层参量M(3000)F2和hmF2的经验正交函数方法建模研究

The ionospheric parameter M(3000)F2 (the so-called transmission factor or the propagation factor) is important not only in practical applications such as frequency planning for radio-communication but also in ionospheric modeling. This parameter is strongly anti-correlated with the ionospheric F2-layer peak height hmF2，a parameter often used as a key anchor point in some widely used empirical models of the ionospheric electron density profile (e.g., in IRI and NeQuick models). Since hmF2 is not easy to obtain from measurements and M(3000)F2 can be routinely scaled from ionograms recorded by ionosonde/digisonde stations distributed globally and its data has been accumulated for a long history, usually the value of hmF2 is calculated from M(3000)F2 using the empirical formula connecting them. In practice, CCIR M(3000)F2 model is widely used to obtain M(3000)F2 value. However, recently some authors found that the CCIR M(3000)F2 model has remarkable discrepancies with the measured M(3000)F2, especially in low-latitude and equatorial regions. For this reason, the International Reference Ionosphere (IRI) research community proposes to improve or update the currently used CCIR M(3000)F2 model. Any efforts toward the improvement and updating of the current M(3000)F2 model or newly development of a global hmF2 model are encouraged. In this dissertation, an effort is made to construct the empirical models of M(3000)F2 and hmF2 based on the empirical orthogonal function (EOF) analysis combined with regression analysis method. The main results are as follows: 1. A single station model is constructed using monthly median hourly values of M(3000)F2 data observed at Wuhan Ionospheric Observatory during the years of 1957–1991 and compared with the IRI model. The result shows that EOF method is possible to use only a few orders of EOF components to represent most of the variance of the original data set. It is a powerful method for ionospheric modeling. 2. Using the values of M(3000)F2 observed by ionosondes distributed globally, data at grids uniformly distributed globally were obtained by using the Kriging interpolation method. Then the gridded data were decomposed into EOF components using two different coordinates: (1) geographical longitude and latitude; (2) modified dip (Modip) and local time. Based on the EOF decompositions of the gridded data under these two coordinates systems, two types of the global M(3000)F2 model are constructed. Statistical analysis showed that the two types of the constructed M(3000)F2 model have better agreement with the observational M(3000)F2 than the M(3000)F2 model currently used by IRI. The constructed models can represent the global variations of M(3000)F2 better. 3. The hmF2 data used to construct the hmF2 model were converted from the observed M(3000)F2 based on the empirical formula connecting them. We also constructed two types of the global hmF2 model using the similar method of modeling M(3000)F2. Statistical analysis showed that the prediction of our models is more accurate than the model of IRI. This demonstrated that using EOF analysis method to construct global model of hmF2 directly is feasible. The results in this thesis indicate that the modeling technique based on EOF expansion combined with regression analysis is very promising when used to construct the global models of M(3000)F2 and hmF2. It is worthwhile to investigate further and has the potential to be used to the global modeling of other ionospheric parameters.

一种测量光刻机工件台方镜不平度的新方法

提出一种新的步进扫描投影光刻机工件台方镜不平度测量方法。以方镜平移补偿量与旋转补偿量为测量目标,使用两个双频激光干涉仪分别测量工件台在x和y方向的位置和旋转量；将方镜不平度的测量按照一定的偏移量分成若干个序列,每一个序列包括对方镜有效区域的若干次往返测量；根据所有序列的测量结果计算出方镜的旋转补偿量；为每一个序列建立临时边界条件,并据此计算出每一序列所测得的方镜粗略平移补偿量；采用三次样条插值与最小二乘法建立每一个序列间的关系,以平滑连接所有测量序列得到精确的方镜平移补偿量。结果表明,该方法用于测量方镜平

激光辐照引起Ge2Sb2Te5非晶态薄膜的电/光性质变化

研究了激光辐照引起Ge2Sb2Te5非晶态薄膜的电/光性质变化,当激光功率为580mW时薄膜的方块电阻有四个数量级（10^7～10^3Ω/□）的突变;对电阻发生突变前、中、后的三个样品进行了XRD测试,结果表明,随着激光功率的增大,薄膜由非晶态向晶态转变,用椭偏仪测试了结构转变前、中、后三个样品的光学常数,在可见光范围内薄膜的光学常数在波长相同情况下有：n（非晶态）〉n（中间态）〉n（晶态）,k（晶态）〉k（中间态）〉k（非晶态）,α（晶态）〉α（中间态）〉α（非晶态）,结合电阻变化曲线和XRD图谱讨论了激光辐照Ge2Sb2Te5非晶态薄膜的电/光性质变化同激光功率和结构转变之间的关系.

Structural change of laser-irradiated Ge2Sb2Te5 films studied by electrical property measurement

Sheet resistance of laser-irradiated Ge2Sb2Te5 thin films prepared by magnetron sputtering was measured by the four-point probe method. With increasing laser power the sheet resistance undergoes an abrupt drop from 10(7) to 10(3) Omega/square at about 580 mW. The abrupt drop in resistance is due to the structural change from amorphous to crystalline state as revealed by X-ray diffraction (XRD) study of the samples around the abrupt change point. Crystallized dots were also formed in the amorphous Ge2Sb2Te5 films by focused short pulse laser-irradiated, the resistivities at the crystallized dots and the non-crystallized area are 3.375 x 10(-3) and 2.725 Omega m, sheet resistance is 3.37 x 10(4) and 2.725 x 10(7) Omega/square respectively, deduced from the I-V Curves that is obtained by conductive atomic force microscope (C-AFM). (C) 2008 Elsevier B.V. All rights reserved.

样点法在森林鸟类调查中的运用

下载PDF阅读器样点法由于具有易于实施、更易做到随机化或系统化以及适合于复杂及斑块化生境等优点,已成为目前使用最为广泛的森林鸟类调查方法.本文在介绍样点法的假设条件、类型和影响鸟类调查准确性因素的基础上,重点介绍了样点法在调查森林鸟类需要考虑的样本大小、样点布设、调查时间、每点的停留时间和调查周期等.同时也分析了样点法的不足,并提出样点法与网捕法结合是提高调查结果准确性的有效途径.

Analysis of interdigital electro-optic Bragg diffraction grating

Some integrated optics devices can be made based on the interdigital electro-optic Bragg diffraction grating. The point-matching method is extended to the analysis of interdigital electro-optic Bragg diffraction gratings. This method provides a simple and fast analytic expression of the electric field in the structure. The field distributions are used to calculate the optical and electrical characteristic parameters of the gratings. The effects of finite conductor thickness have been taken into account in the analysis. It is shown that the simulation results agree well with the measured data.

Information geometry theory of high-dimension space and application for speaker independent continuous digit speech recognition

In the light of descriptive geometry and notions in set theory, this paper re-defines the basic elements in space such as curve and surface and so on, presents some fundamental notions with respect to the point cover based on the High-dimension space (HDS) point covering theory, finally takes points from mapping part of speech signals to HDS, so as to analyze distribution information of these speech points in HDS, and various geometric covering objects for speech points and their relationship. Besides, this paper also proposes a new algorithm for speaker independent continuous digit speech recognition based on the HDS point dynamic searching theory without end-points detection and segmentation. First from the different digit syllables in real continuous digit speech, we establish the covering area in feature space for continuous speech. During recognition, we make use of the point covering dynamic searching theory in HDS to do recognition, and then get the satisfying recognized results. At last, compared to HMM (Hidden Markov models)-based method, from the development trend of the comparing results, as sample amount increasing, the difference of recognition rate between two methods will decrease slowly, while sample amount approaching to be very large, two recognition rates all close to 100% little by little. As seen from the results, the recognition rate of HDS point covering method is higher than that of in HMM (Hidden Markov models) based method, because, the point covering describes the morphological distribution for speech in HDS, whereas HMM-based method is only a probability distribution, whose accuracy is certainly inferior to point covering.

Electrical and Optical Phase Transition Properties of Nano Vanadium Dioxide Thin Films

Nano-vanadium dioxide thin films were prepared through thermal annealing vanadium oxide thin films deposited by dual ion beam sputtering. The nano-vanadium dioxide thin films changed its state from semiconductor phase to metal phase through heating by homemade system. Four point probe method and Fourier transform infrared spectrum technology were employed to measure and anaylze the electrical and optical semiconductor-to-metal phase transition properties of nano-vanadium dioxide thin films, respectively. The results show that there is an obvious discrepancy between the semiconductor-to-metal phase transition properties of electrical and optical phase transition. The nano-vanadium dioxide thin films' phase transiton temperature defined by electrical phase transiton property is 63 degrees C, higher than that defined by optical phase transiton property at 5 mu m, 60 degrees C; and the temperature width of electrical phase transition duration is also wider than that of optical phase transiton duration. The semiconductor-to-metal phase transiton temperature defined by optical properties increases with increasing wavelength in the region of infrared wave band, and the occuring temperature of phase transiton from semiconductor to metal also increases with wavelength increasing, but the duration temperature width of transition decreases with wavelength increasing. The phase transition properties of nano-vanadium dioxide thin film has obvious relationship with wavelength in infrared wave band. The phase transition properties can be tuned through wavelength in infrared wave band, and the semiconductor-to-metal phase transition properties of nano vanadiium dioxide thin films can be better characterized by electrical property.

Dislocation movement in nitrogen-doped Czochralski silicon

Dislocation movement in N-doped Czochralski silicon (Cz-Si) was surveyed by four point bend method. Dislocation movement velocities in Cz-Si doped with nitrogen, with both nitrogen and antimony, and with only antimony were investigated. The order of measured dislocation movement velocities, at 700 degrees C less than or equal to T less than or equal to 800 degrees C and under resolved stress sigma=4.1 kg/mm(2), was V-Sb.O > V-n.Sb.O>V-N.O. The experiments showed that nigtrogen doping could retard the movement of dislocations.

地形可视性分析中观察点设置问题研究

观察点设置问题是地形可视性分析中的一类重要问题，对该问题的研究可以在空间信息辅助决策、通信、旅游、野生动物保护等领域发挥重大作用。本文在对地形可视性分析中观察点设置问题现有研究成果总结和分析基础上对该问题展开深入研究。 首先，针对现有解决方法只从智能算法或地形数据表示方法单一角度进行分析和研究的局限性，提出了一种问题相关的智能算法和数据表示方法相结合的解决问题新框架。该框架考虑了解决观察点设置问题时智能算法的优点和数据表示方式的特点相互配合问题，目的是充分发挥二者各自的优势以提高观察点设置问题解决的准确度与效率。 其次，在深入分析观察点设置问题本身特点的基础上，结合隶属云理论的基本理论和方法，对经典模拟退火算法从退温函数设计、温度产生过程、状态生成过程三方面进行了问题相关的改进，提出了一种适于观察点设置问题的改进模拟退火算法(Improved Simulated Annealing algorithm, ISA)。该算法一方面保持了经典模拟退火算法的稳定倾向特性，保证了算法满足伴随退火温度的不断下降，对恶化的新状态越来越难于接受这一模拟退火算法的最基本特征；另一方面其退火温度的连续性随机变化特性和隐含的“回火升温”过程，则有利于算法有效拒绝恶化解，加速算法收敛，能够更好地满足观察点设置问题对于算法收敛速度的要求。 再次，在分析地形数据的精度、误差等因素对观察点设置问题的解决准确性和解决效率影响程度的基础上，提出了一种基于离散余弦变换的地形数据内插方法(Discrete Cosine Transformation Interpolation method, DCTI)。新方法将传统空域上的地形内插转换到变换域上进行，同时充分利用了离散余弦变换的熵保持特性和能量压缩特性，简化了变换域上的内插过程，提高了地形数据内插的效率和精度。DCTI方法与其他现有典型地形数据内插方法相比，对地形可视性信息获取的准确性和效率影响最小，为平衡观察点设置问题解决过程中时间效率和准确度之间的关系，最终有效地解决观察点设置问题提供了数据基础。 最后，从智能算法和地形数据相结合的角度出发，提出了一种基于ISA和DCTI相结合的观察点设置问题多分辨率处理方法(Multi-Resolution Processing method, MRP)。新方法将模拟退火算法的逐次退火特点和地形数据的多分辨率表示充分结合，达到了发挥算法数据相结合的综合优势的目的。与现有单纯基于模拟退火算法的解决方法相比，在问题解决准确度保持不变的前提下，基于MRP方法的观察点设置问题解决的平均耗时减少85%~95%，为实际工程应用问题的解决提供了一条重要途径。

区域水土流失研究中降雨量空间插值方法选取

在区域水土流失模型研究中,空间插值可提供每个计算栅格的气象要素资料。考虑到研究区域降雨与高程相关性很弱,不宜采用梯度距离反比法(GIDS),故采用距离反比法(IDW)和普通克里格法(Kriging),对延安示范区及其周围共50个站点2000—2003年的5—10月逐月降雨量进行插值。交叉验证结果表明:对2种插值方法,二者经过对数变换后平均相对误差(MRE)为8.30%和7.67%,分别比原始数据插值后的MRE下降了23.17%和23.50%,说明插值精度得到了提升,对研究区域某一年逐月降水的插值Kriging方法比IDW方法更加精确。

Insulin nanoparticle preparation and encapsulation into poly(lactic-co-glycolic acid) microspheres by using an anhydrous system

Insulin has been encapsulated in poly(lactic-co-glycolic acid) (PLGA) microspheres by solid-in-oil-in-oil (S/O/O) emulsion technique using DMF/corn oil as new solvent pairs. To get better encapsulation efficiency, insulin nanoparticles were prepared by the modified isoelectric point precipitation method so that it had good dispersion in the inner oil phase. The resulting microspheres had drug loading of 10% (w/w), while the encapsulation efficiency could be up to 90-100%. And the insulin release from the microspheres could last for 60 days. Microspheres encapsulated original insulin with the same method had lower encapsulation efficiency, and shorter release period. Laser scanning confocal microscopy indicated the insulin nanoparticle and original insulin had different distribution in microspheres. The results suggested that using insulin nanoparticle was better than original insulin for microsphere preparation by S/O/O method.