|
|
|
|
|
|
| n=1(Total_x/deta_x)+1); %% |
|
|
| omega=2/(1+sqrt(1-
| ((cos(pi/nn)+cos(pi/nn))/2 ).^2) ); %% |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| for
| iteration=1:iteration_max |
|
|
| u_next(i,j)=omega.*sigma.*(u(i+1,j)+u(i,j+1)+u(i-1,j)+u(i,j-1))... |
|
|
|
|
| if( max(max(u-u_last))< convergence) |
|
|
|
|
|
|
|
| clabel(C,h,V,\'FontSize\',12,\'Color\',\'red\',\'LabelSpacing\',1200); |
|
|
|
|
| 5.结果分析
| 容许误差为1e-8,最大迭代次数10000次,实际迭代次数335次。雅各比迭代次数或高斯赛德尔迭代相比较收敛较快。将结果绘制成云图如图3,较符合物理规律。
| 6.注1:椭圆型方程也常用于有限差分法的网格变换。 |
注2:初值条件打错了,是x=1时,u=0....... |
|
|