光学Cherenkov差频产生太赫兹技术研究.ppt

有限差分迭代及数值拟合 2014科学计算选讲课程报告 报告人：刘宏翔1013202027 组 员：宋寰宇1013202028 王茂榕1013202029 电磁场问题求解 理论/实验值对比 数值 计算 引言 有限差分 简单迭代与超松弛迭代 最小二乘法 曲线拟合 差分 简单、直观。无论是常微分方程还是偏微分方程，初值问题或边界问题，椭圆形、双曲形、抛物形二阶线性方程，以及高阶或非线性方程，通常都可以利用此法将其转化为代数方程组，再借助计算机做数值计算。 前向差分 后向差分 中心差分 差分 做Taylor级数展开 中心差分截断误差最小 二阶导数 离散系统求解电磁场连续域问题 差分 将求解域内连续场分布用网格节点上的离散数值解代替。当网格分的充分细时，才能达到足够的精度。 迭代 简单迭代 用围绕它的四个点的电位的平均值作为它的新值，当所有点计算完后，用它们的新值代替旧值，完成一次迭代计算，再进行下一次迭代，直到每一点计算的新值与旧值之差小于指定的范围为止。 迭代 超松弛法 引入松弛因子，加快收敛速度 迭代 分为10*10格 迭代次数k=290 简单迭代 迭代 分为10*10格 迭代次数k=154 松弛法 迭代 分为10*10格 迭代次数k=154 未引入松弛因子 迭代 分为10*10格 迭代次数k=37 松弛因子w=1.5279 引入松弛因子w 迭代 分为50*50格 迭代次数k=3037 未引入松弛因子 迭代 分为50*50格 迭代次数k=178 松弛因子w=1.8818 引入松弛因子w 0 0 0 0 0 0 0 0 0 0 0 0 1.10749900509585 2.09934359632906 2.87750174144730 3.37156898764052 3.54066662379080 3.37156910794312 2.87750191719825 2.09934388381103 1.10749901908117 0 0 2.33065193757706 4.41237422378446 6.03909443092522 7.06810779228716 7.41952850753098 7.06810794655764 6.03909473497658 4.41237461794051 2.33065218623347 0 0 3.80273482474707 7.18040727713608 9.79839417314094 11.4422392993502 12.0012317260702 11.4422394594873 9.79839447648837 7.18040767357934 3.80273508958600 0 0 5.69988028158639 10.7081261716114 14.5318357879330 16.9012235933887 17.7009196978387 16.9012237367560 14.5318360543472 10.7081264985037 5.69988048464073 0 0 8.28866016584190 15.4203815154595 20.7195992732075 23.9298996256100 24.9999997869452 23.9298997501224 20.7195994942518 15.4203817677753 8.28866033925234 0 0 12.0343789704918 21.9651405492133 28.9962802278297 33.0987758786947 34.4392800873920 33.0987759660732 28.9962803846733 21.9651407313991 12.0343790933102 0 0 17.8837152124614 31.4095215592073 40.2016052483881 45.0296435807337 46.5595687108899 45.0296436319817 40.2016053448796 31.4095216706439 17.8837152902679 0 0 28.0909603495355 45.5876252515968 55.3709756326577 60.2586244708383 61.7397075377044 60.2586244994186 55.3709756824080 45.5876253031910 28.0909603901734 0 0 48.8925009553839 67.479043442685