20060401

％ 奇异值分解 (sigular value decomposition,SVD) 是另一种正交矩阵分解法；SVD是最可靠的分解法， ％ 但是它比QR 分解法要花上近十倍的计算时间。[U,S,V]=svd(A)，其中U和V代表二个相互正交矩阵， ％ 而S代表一对角矩阵。 和QR分解法相同者， 原矩阵A不必为正方矩阵。 ％ 使用SVD分解法的用途是解最小平方误差法和数据压缩。用svd分解法解线性方程组，在Quke2中就用这个来计算图形信息，性能相当的好。在计算线性方程组时，一些不能分解的矩阵或者严重病态矩阵的线性方程都能很好的得到解( Singular value decomposition (sigular value decomposition, SVD) is another orthogonal matrix decomposition method SVD decomposition is the most reliable method, but it takes more than QR decomposition near ten times the computing time. [U, S, V] = svd (A), in which U and V on behalf of two mutually orthogonal matrix, and the S on behalf of a diagonal matrix. And QR decomposition are the same, the original matrix A is no need for the square matrix. The use of SVD decomposition method are used as a solution of least squares error method and data compression. Using SVD decomposition solution of linear equations, in Quke2 on to use this information to calculate the graphics performance quite good. In the calculation of linear equations, some indecomposable matrix or serious pathological matrix of linear equations can be a very good solution)

CT迭代重建算法中3维投影矩阵的生成算法，该算法对学习计算机断层层析成像很有帮助。(CT image reconstruction projection ,3D)

Orthogonal Polynomials Approximation 数值分析，计算正交基多项式的系数 (Given a function f and a set of m >0 distinct points . You are supposed to write a function to approximate f by an orthogonal polynomial using the exact function values at the given m points with a weight assigned to each point . The total error must be no larger than a given tolerance. Format of function int OPA( double (*f)(double t), int m, double x[], double w[], double c[], double*eps ) where the function pointer double (*f)(double t) defines the function f int m is the number of points double x[] contains points double w[] contains the values of a weight function at the given points x[] double c[] contains the coefficients of the approximation polynomial double*eps is passed into the function as the tolerance for the error, and is supposed to be returned as the value of error. The function OPA is supposed to return the degree of the approximation polynomial. Note: a constant Max_n is defined so that if the total error is still not small enough when n = Ma)

一圆形游泳池如图所示，现在需在其周围建一圆形过道，并在其四周围上栅栏。 栅栏价格为35元/米，过道造价为20元/平方米。 过道宽度为3米， 游泳池半径由键盘输入。 要求编程计算并输出过道和栅栏的造价。(A circular pool as shown in figure, now in its built around to a circular corridor, and in its 4 around on the fence. The fence for 35 yuan/m price, corridor cost for 20 yuan/square metre. Aisle width is 3 meters, The radius of the swimming pool by keyboard input. Demand and output programming calculation of corridor and the cost of the fence. )

MFC实现快速傅里叶变换（FFT）MFC实现快速傅里叶变换（FFT）(MFC MFC Fast Fourier Transform (FFT) fast Fourier transform (FFT))

交错网格PML弹性波模拟地震波形的c语言程序(Staggered grid PML elastic wave seismic waveform simulation)

FFT变换，有硬件说明架构和框图，说明傅里叶变换的原理和意义，同时附有128点的傅里叶变换(FFT )

fft code in visual c++

Concert c++ example of facility location problem. C++ concert technology is a registered trademark of ILOG CPLEX

实现重力式热管换热器设计计算功能，简单按照设计任务书，给出所需热管的结构尺寸，涉及的热管蒸发段和冷凝段长度分别为1m，绝热段长度30cm，内径20cm，外径25cm，热管工质为水。(C language for thermosyphon s design.Realize the gravity heat pipe heat exchanger design and calculation functions , simple design in accordance with the mission statement , given the size of the required structure of the heat pipe , heat pipe evaporator and condenser section lengths were involved 1m, adiabatic length 30cm, inner diameter 20cm, 25cm outside diameter heat pipe working fluid is water.)