摘要 |
本发明提供一种特殊鞍点问题的高效预处理方法;对鞍点问题(3)的结构特征,设计结构化的免增广和免Schur余的块三角预处理子(11),(12);对设计的预处理子进行谱分析,分析当选择合适参数时预处理矩阵将有更好的特征值聚集性,同时弥补原块对角预处理子只能解决k<sup>2</sup><1的不足,扩大其适用范围;对具有奇异(1,1)块的非对称鞍点问题(6)结构特征,设计更加广义的结构化块三角预处理子;对设计的预处理子(13)进行谱分析,研究设计的预处理矩阵特征值分布、相应的特征向量和最小多项式,并分析选择最优参数时特征值的聚集性。 |
主权项 |
特殊鞍点问题的高效预处理方法,其特征在于,(1)对由离散化混合型时谐Maxwell方程离散产生的鞍点问题<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mover><mi>A</mi><mo>‾</mo></mover><mn>1</mn></msub><mi>x</mi><mo>=</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><mrow><mi>A</mi><mo>-</mo><msup><mi>k</mi><mn>2</mn></msup><mi>M</mi></mrow></mtd><mtd><msup><mi>B</mi><mi>T</mi></msup></mtd></mtr><mtr><mtd><mi>B</mi></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mfenced open = "(" close = ")"><mtable><mtr><mtd><mi>u</mi></mtd></mtr><mtr><mtd><mi>v</mi></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><mi>f</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>=</mo><mi>b</mi><mo>,</mo></mrow>]]></math><img file="FDA0000869015700000011.GIF" wi="758" he="158" /></maths>根据此鞍点问题的结构特征,设计两种结构化免增广和免Schur余的块三角预处理子:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>M</mi><mrow><mi>ξ</mi><mo>,</mo><mi>η</mi></mrow></msub><mo>=</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><mrow><mi>A</mi><mo>+</mo><mrow><mo>(</mo><mi>ξ</mi><mo>-</mo><msup><mi>k</mi><mn>2</mn></msup><mo>)</mo></mrow><mi>M</mi></mrow></mtd><mtd><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>ξ</mi><mi>η</mi><mo>)</mo><msup><mi>B</mi><mi>T</mi></msup></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mi>η</mi><mi>L</mi></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000869015700000012.GIF" wi="1318" he="158" /></maths>和<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>H</mi><mrow><mi>ξ</mi><mo>,</mo><mi>η</mi></mrow></msub><mo>=</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><mrow><mi>A</mi><mo>-</mo><msup><mi>k</mi><mn>2</mn></msup><mi>M</mi><mo>+</mo><msup><mi>ξB</mi><mi>T</mi></msup><msup><mi>L</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>B</mi></mrow></mtd><mtd><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>ξ</mi><mi>η</mi><mo>)</mo><msup><mi>B</mi><mi>T</mi></msup></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><mi>η</mi><mi>L</mi></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000869015700000013.GIF" wi="1331" he="159" /></maths>(2)对设计的预处理子(11)和(12),分析其构造及应用代价和已有的免增广和免Schur余的块对角预处理子相当,分析参数最优选取时特征值的聚集性,同时弥补原块对角预处理子只能解决k<sup>2</sup><1的不足,扩大其适用范围;(3)对具有奇异(1,1)块的非对称鞍点问题<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mover><mi>A</mi><mo>‾</mo></mover><mn>2</mn></msub><mi>x</mi><mo>=</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><mi>A</mi></mtd><mtd><msup><mi>B</mi><mi>T</mi></msup></mtd></mtr><mtr><mtd><mi>C</mi></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mfenced open = "(" close = ")"><mtable><mtr><mtd><mi>u</mi></mtd></mtr><mtr><mtd><mi>v</mi></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><mi>f</mi></mtd></mtr><mtr><mtd><mi>g</mi></mtd></mtr></mtable></mfenced><mo>=</mo><mi>b</mi><mo>,</mo></mrow>]]></math><img file="FDA0000869015700000014.GIF" wi="630" he="159" /></maths>根据此鞍点问题的结构特征,研究设计出一种广义的结构化块三角预处理子:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>T</mi><mrow><mi>ξ</mi><mo>,</mo><mi>η</mi></mrow></msub><mo>=</mo><mfenced open = "(" close = ")"><mtable><mtr><mtd><mrow><mi>A</mi><mo>+</mo><msup><mi>ξB</mi><mi>T</mi></msup><msup><mi>W</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mi>C</mi></mrow></mtd><mtd><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>ξ</mi><mi>η</mi><mo>)</mo><msup><mi>B</mi><mi>T</mi></msup></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>η</mi><mi>W</mi></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000869015700000015.GIF" wi="1287" he="158" /></maths>(4)对设计的结构化的预处理子(13),进行理论谱分析,研究预处理矩阵<img file="FDA0000869015700000016.GIF" wi="125" he="79" />等的特征值分布、相应的特征向量和最小多项式,最后分析,参数满足一定条件时,预处理子<img file="FDA0000869015700000017.GIF" wi="126" he="78" />等的特征值更加聚集。 |