有限差分法实验报告材料(参考)

《有限差分法实验报告材料(参考)》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实用标准文案工程电磁场实验报告——有限差分法文档大全实用标准文案用超松弛迭代法求解接地金属槽内电位的分布=V100j一、实验要求按对称场差分格式求解电位的分布已知：给定边值：如图1-7示图1-7接地金属槽内半场域的网格给定初值误范围差:计算：迭代次数，，将计算结果保存到文件中二、实验思想有限差分法有限差分法（FiniteDifferentialMethod）是基于差分原理的一种数值计算法。其基本思想：将场域离散为许多小网格，应用差分原理，将求解连续函数的泊松方程的问题转换为求解网格节点上文档大全实用标准文案的差分方程组的问题。泊松方程的五点差分格式当场域中得到拉普拉

2、斯方程的五点差分格式图1-4高斯——赛德尔迭代法差分方程组的求解方法(1)高斯——赛德尔迭代法(1-14)式中：·迭代顺序可按先行后列，或先列后行进行。·迭代过程遇到边界节点时，代入边界值或边界差分格式，直到所有节点电位满足为止。（2）超松弛迭代法(1-15)式中：——加速收敛因子可见：迭代收敛的速度与有明显关系二、程序源代码#include#include#includedoubleA[5][5];voidmain(void){doubleBJ[5][5];//数组B用于比较电势ints[100];//

3、用于储存迭代次数文档大全实用标准文案doubled[100];//用于记录所有的加速因子d[0]=1.0;inti,j,N=0,M=0,x;for(i=0;i<100;i++)d[i]=0.01*i+d[0];//加速因子从1.0到2.0之间的20个数！doublew[100][10];intP,Q;for(P=0;P<4;P++)for(Q=0;Q<5;Q++)A[P][Q]=0;for(P=0;P<5;P++)A[4][P]=100;cout<<"数组A的所有元素是:"<

4、[i][j]<

5、[i][j]-BJ[i][j])<1e-5)pp++;}N++;}while(pp<=9);pp=0;for(i=0;i<3;i++)文档大全实用标准文案w[M][i+1]=A[1][i+1];for(i=3;i<6;i++)w[M][i+1]=A[2][i-2];for(i=6;i<9;i++)w[M][i+1]=A[3][i-5];s[M]=N;M++;N=0;intP,Q;for(P=0;P<4;P++)for(Q=0;Q<5;Q++)A[P][Q]=0;for(P=0;P<5;P++)A[4][P]=100;}intmin=s[0];intp,q;cout<

6、<"输出所有的加速因子的迭代次数:"<<'';for(q=1;q<100;q++){//cout<s[q]){min=s[q];p=q;}}cout<

7、==0)cout<<'';}cout<<'';}文档大全实用标准文案二、程序框图启动赋边界节点已知电位值赋予场域内各节点电位初始值累计迭代次数N=0N=N+1按超松弛法进行一次迭代，求所有内点相邻二次迭代值的最大误差是否小于停机打印NY文档大全实用标准文案迭代解程序框图二、结果分析迭代收敛的速度与的关系收敛因子（）1.01.71.81.831.851.871.902.0迭代次数（）>1000269174143122133171发散最佳收敛因子的经验公式：（正方形场域、正方形网格）（矩形场域、正方形网格）程序执行结果如下文档大全实用标准文案文档大全