长期观测资料的分析.ppt

式（29）又可以写为其中，各项系数可以写为以及方程组（30）称为矛盾方程组（26）的法方程或正规方程。法方程中方程的个数与未知数的个数相等，可按普通方程组求解，而且法方程的系数行列式是对称的，即法方程系数cij的求解，只要把矛盾方程组（26）的第i列和第j列系数取出来，一一对应求积，然后再求和；如计算fi，则把矛盾方程的第第i列系数和等号右边的一列数做相同的运算即可（32）等时间间隔连续观测序列法方程的系数式（31）和（32）可以求解任意时间间隔连续观测序列法方程的各元素对于等时间间隔连续观测序列（与无关）的法方程的系数，假定观测记录数为，取中间时刻（），计算可以进一步简化。令则观测时刻依次为，…，，0，，…，，对应水位依次为，…，，，，…，，则法方程可写为方程组(25）的法方程为根据法方程系数运算规则可以算出CBAK数学归纳法以及以及（37a）（38）（37b）由于所有的系数C和D均为零，法方程（33）实际可以分为两个方程组和（39a）（39b）1上面两个方程组的系数行列式依然是对称的2方程组（39a）包含K+1个方程，可以解出S0,x1,x2,﹍,xK；方程组（39b）包含K个方程，可以解出y1,y2,﹍,yK。然后再根据式（24a，b）和式（22）可以求出调和常数3对于潮流的分析可采取相同的方法，只要先将流速和流向的观测值和分解为北、东的分流和，然后对和像对一样处理，可以求出北分流和东分流的调和常数，进而还可以求出潮流椭圆要素一些结论傅立叶系数的计算：公式（37b）和（38）中的和分别包含NK个余弦和正弦值，直接计算会影响程序的效率格尔策（Goertzer）迭代法由于削减了大部分三角函数的计算过程可以显著提高计算效率，首先令则有（40）（41）同理，对数列f0,f1,f2,…,fK按如下迭代方式计算则有（42b）（42a）比较式（37b）和(42b）以及（38）和（42a）如果，且令另外假定则有取样：在离散的时刻观测或从曲线读取观测值01样本：读取下来用于分析的观测序列02潮汐分析中任何一个样本都是离散的，且时间长度是有限的03第4节潮汐分潮的选取与观测记录04的时间间隔以及观测时段长05度的关系06如果观测时间段的长度，两个观测记录之间的时间间隔与选取的分潮配合不好，就可能得到不准确的结果，甚至可能计算不出结果来。设两个分潮角速率，则理论潮位和实际潮位分别为以和为算例，得到两个分潮的法方程（45）（46）（47）法方程系数和右端项可以由（34）和（37）给出*第三章长期观测资料的分析第1节实际潮汐分潮的调和常数实际潮位的可预测部分：引力潮（天文潮）、气象潮、天文气象复合潮、浅水潮引力潮：直接由引潮力作用产生的分潮气象潮：由气象因素(如风、气压、降水和蒸发等)所引起的海面振动现象迎岸风可以引起水位水位上升，离岸风引起水位下降，高气压能使水位降低而低气压则使水位升高反映水位周年变化的分潮是Sa和Ssa，另外气象条件的周日变化也可引起相应水位的变化，引入一个平太阳日分潮S1，上述分潮就是气象分潮气象天文复合潮：气象条件改变不仅可以直接引起气象分潮，还可以通过改变天文潮波在海洋中的运动状况，从而使得海洋对引潮力的响应发生改变，主要表现为实际分潮的振幅和迟角产生季节变化，即形成天文气象复合潮，如MB2，SA2浅水潮：天文潮进入浅水区由于海底摩擦或天文分潮的非线性耗散所产生的分潮，如M4，MS4噪声：气象扰动引起水位的不规则变化，会造成分析结果的误差从属分潮：由天文或气象源分潮复合得到的天文气象复合潮或浅水分潮统称为从属分潮字母下标代表周期，如a、sa、m、f分别代表一年、半年、一月、半月01大部分从属分潮的源分潮用单独一个字母表示，如M2、Sa、O1，例外的是2Q1、2N2、OO1三个小分潮02从属分潮中倍潮只用一个字母表示，如M4、O3、M6，复合分潮则由构成成它的所有源分潮表示，如MS4、MSN6，03天文气象复合潮由源分潮加上字母A或B构成04从属分潮的命名天文潮可以展开为许多余弦振动之和，对应某一频率做周期变化的引潮力分潮，海洋也要产生这一频率的震荡，因此某一定地