【最新资料】五点差分法(matlab)解椭圆型偏微分方程

《【最新资料】五点差分法(matlab)解椭圆型偏微分方程》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、用差分法解椭圆型偏微分方程-（Uxx+Uyy)=(pi*pi-1)e^xsin(pi*y)0

2、误差，p为精确解symstemp;u=zeros(n+1,m+1);x=0+(0:m)*h;y=0+(0:n)*h;for(i=1:n+1)u(i,1)=sin(pi*y(i));u(i,m+1)=exp(1)*exp(1)*sin(pi*y(i));endfor(i=1:n)for(j=1:m)f(i,j)=(pi*pi-1)*exp(x(j))*sin(pi*y(i));endendt=zeros(n-1,m-1);for(k=1:kmax)for(i=2:n)for(j=2:m)temp=h*h*f(i,j)

3、/4+(u(i,j+1)+u(i,j-1)+u(i+1,j)+u(i-1,j))/4;t(i,j)=(temp-u(i,j))*(temp-u(i,j));u(i,j)=temp;endendt(i,j)=sqrt(t(i,j));if(k>kmax)break;endif(max(max(t))

4、;endEnd在命令窗口中输入：[peuxyk]=wudianchafenfa(0.1,20,10,10000,1e-6)k=147surf(x,y,u);xlabel(‘x’);ylabel(‘y’);zlabel(‘u’);Title(‘五点差分法解椭圆型偏微分方程例1’)就可以得到下图ｓｕｒｆ（ｘ，ｙ，ｐ）ｓｕｒｆ（ｘ，ｙ，ｅ）[peuxyk]=wudianchafenfa(0.05,40,20,10000,1e-6)[peuxyk]=wudianchafenfa(0.025,80,40,10000,1e-6

5、)为什么分得越小，误差会变大呢？我们试试运行：[peuxyk]=wudianchafenfa(0.025,80，40,10000,1e-8)K=2164surf(x,y,e)误差变小了吧还可以试试[peuxyk]=wudianchafenfa(0.025,80，40,10000,1e-10)K=3355误差又大了一点再试试[peuxyk]=wudianchafenfa(0.025,80,40,10000,1e-11)k=3952误差趋于稳定总结：最终的误差曲面与网格数有关，也与设定的迭代前后两次差值（ep,看程序）

6、有关；固定网格数，随着设定的迭代前后两次差值变小，误差由大比变小，中间有一个最小值，随着又增大一点，最后趋于稳定。也许可以去研究一下那个误差最小的地方或者研究趋于稳定时的临界值。以下是附加文档，不需要的朋友下载后删除，谢谢顶岗实习总结专题13篇第一篇:顶岗实习总结为了进一步巩固理论知识，将理论与实践有机地结合起来，按照学校的计划要求，本人进行了为期个月的顶岗实习。这个月里的时间里，经过我个人的实践和努力学习，在同事们的指导和帮助下，对村的概况和村委会有了一定的了解，对村村委会的日常工作及内部制度有了初步的认识，同时

7、，在与其他工作人员交谈过程中学到了许多难能可贵经验和知识。通过这次实践，使我对村委会实务有所了解，也为我今后的顺利工作打下了良好的基础。一、实习工作情况村是一个（此处可添加一些你实习的那个村和村委会的介绍）我到村村委会后，先了解了村的发展史以及村委会各个机构的设置情况，村委会的规模、人员数量等，做一些力所能及的工作，帮忙清理卫生，做一些后勤工作；再了解村的文化历史，认识了一些同事，村委会给我安排了一个特定的指导人；然后在村委会学习了解其他人员工作情况，实习期间我努力将自己在学校所学的理论知识向实践方面转化，尽量做到

8、理论与实践相结合。在实习期间我遵守了工作纪律，不迟到、不早退，认真完成领导交办的工作。我在村委会主要是负责管理日常信件的工作，这个工作看似轻松，却是责任重大，来不得办点马虎。一封信件没有及时收发，很有可能造成工作的失误、严重的甚至会造成巨大的经济损失。很感谢村委会对我这个实习生的信任，委派了如此重要的工作给我。在实习过程中，在信件收发管理上，我一直亲力亲为，