matlab解拉普拉斯方程ppt课件

资源ID：57076945 资源大小：276.50KB 全文页数：25页
资源格式： PPT 下载积分：20金贝

快捷下载

账号登录下载

微信登录下载

三方登录下载：

微信扫一扫登录

下载资源需要20金贝

邮箱/手机：
温馨提示：	快捷下载时，用户名和密码都是您填写的邮箱或者手机号，方便查询和重复下载（系统自动生成）。如填写123，账号就是123，密码也是123。
支付方式：
验证码：	换一换

账号：
密码：
验证码：	换一换
当日自动登录忘记密码？

1、金锄头文库是“C2C”交易模式，即卖家上传的文档直接由买家下载，本站只是中间服务平台，本站所有文档下载所得的收益全部归上传人（卖家）所有，作为网络服务商，若您的权利被侵害请及时联系右侧客服；
2、如你看到网页展示的文档有jinchutou.com水印，是因预览和防盗链等技术需要对部份页面进行转换压缩成图而已，我们并不对上传的文档进行任何编辑或修改，文档下载后都不会有jinchutou.com水印标识，下载后原文更清晰；
3、所有的PPT和DOC文档都被视为“模板”，允许上传人保留章节、目录结构的情况下删减部份的内容；下载前须认真查看，确认无误后再购买；
4、文档大部份都是可以预览的，金锄头文库作为内容存储提供商，无法对各卖家所售文档的真实性、完整性、准确性以及专业性等问题提供审核和保证，请慎重购买；
5、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确)，网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据；
6、如果您还有什么不清楚的或需要我们协助，可以点击右侧栏的客服。

下载须知 | 常见问题汇总

1、会员注册 2、如何支付与充值 3、个人资料信息修改 4、我的收藏和“收藏文件夹” 5、我的读者群和加入读者群 6、我的书签 7、金锄头文库批量上传工具（绿色版）V1.0 8、下载文档（资源）相关问题整理 9、解决下载文档时，自动弹出迅雷的问题 10、下载时为什么支付不成功？

matlab解拉普拉斯方程ppt课件

拉普拉斯方程有限差分解,均匀介质情况下,均匀介质情况下,clear all clc m=10 for k=1:mfor j=1:mU(j,k)=0;end end U(5,5)=100.0,for i=1:200for k=2:m-1for j=2:m-1Unew(j,k)=1/4*(U(j+1,k)+U(j-1,k)+U(j,k+1)+U(j,k-1);U(j,k)=Unew(j,k);U(5,5)=100.0;endend end contour(U) axis equal,clear all %清除变量 Clc %清除command window m=10 %网格数目为10 for k=1:mfor j=1:mU(j,k)=0;%赋初值end end,U =0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0,U(5,5)=100.0,U =0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 0 0 100 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0,for i=1:200for k=2:m-1for j=2:m-1Unew(j,k)=1/4*(U(j+1,k)+U(j- 1,k)+U(j,k+1)+U(j,k-1);%拉普拉斯方程有限差分解U(j,k)=Unew(j,k);U(m/2,m/2)=100.0;endend end,for i=1for k=2:m-1for j=2:m-1Unew(j,k)=1/4*(U(j+1,k)+U(j-1,k)+U(j,k+1)+U(j,k-1);U(j,k)=Unew(j,k);U(m/2,m/2)=100.0;endend end,i=1,i=2,i=3,(for i=1:200) i=200,(for i=1:201) i=201,均匀介质情况下,clear all clc m=10 for k=1:mfor j=1:mU(j,k)=0;end end U(5,5)=100.0,for i=1:200for k=2:m-1for j=2:m-1Unew(j,k)=1/4*(U(j+1,k)+U(j-1,k)+U(j,k+1)+U(j,k-1);U(j,k)=Unew(j,k);U(5,5)=100.0;endend end contour(U) axis equal,均匀介质情况下,clear all clc m=10 for k=1:mfor j=1:mU(j,k)=0;end end U(5,5)=100.0 d=1,While d<0.01 for k=2:m-1for j=2:m-1Unew(j,k)=1/4*(U(j+1,k)+U(j-1,k)+U(j,k+1)+U(j,k-1);U(j,k)=Unew(j,k);U(5,5)=100.0;d=abs(U(j,k)- Unew(j,k);)endend end contour(U) axis equal,均匀介质情况下,clear all clc m=10 for k=1:mfor j=1:mU(j,k)=0;end end U(5,5)=100.0,for i=1:200for k=2:m-1for j=2:m-1Unew(j,k)=1/4*(U(j+1,k)+U(j-1,k)+U(j,k+1)+U(j,k-1);U(j,k)=Unew(j,k);U(5,5)=100.0;endend end contour(U) axis equal,2018/10/18,非均匀介质情况下,ep： 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 10 10 10 1 1 1 1 1 10 10 10 1 1 1 1 1 10 10 10 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1,2018/10/18,非均匀介质情况下,2018/10/18,非均匀介质情况下,ep： for k=1:mfor j=1:mep(j,k)=1;end endfor k=4:6for j=4:6ep(j,k)=10;end end,2018/10/18,非均匀介质情况下,Unew(j,k)=1/4*(U(j+1,k)+U(j-1,k)+U(j,k+1)+U(j,k-1);修改为：Unew(j,k)=1/4*(ep(j+1,k)*U(j+1,k)+ep(j-1,k)*U(j-1,k)+ep(j,k+1)*U(j,k+1)+ep(j,k-1)*U(j,k-1);,2018/10/18,非均匀介质情况下,U(5,5)=100.0for i=1:200for k=2:m-1for j=2:m-1Unew(j,k)=1/4*(ep(j+1,k)*U(j+1,k)+ep(j-1,k)*U(j-1,k)+ep(j,k+1)*U(j,k+1)+ep(j,k-1)*U(j,k-1);U(j,k)=Unew(j,k);U(5,5)=100.0;endend end contour(U) axis equal,clear all clc m=10 for k=1:mfor j=1:mU(j,k)=0;end endfor k=1:mfor j=1:mep(j,k)=1;end endfor k=4:6for j=4:6ep(j,k)=10;end end,2018/10/18,非均匀介质情况下,2018/10/18,非均匀介质情况下,ep：1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100,2018/10/18,第一类边界条件,clear all clc m=10 for k=1:mfor j=1:mU(j,k)=0;end endfor j=1for k=1:mU(j,k)=10;end end,2018/10/18,第二类边界条件,for i=1:200for k=2:m-1for j=2:m-1Unew(j,k)=1/4*(ep(j+1,k)*U(j+1,k)+ep(j-1,k)*U(j-1,k)+ep(j,k+1)*U(j,k+1)+ep(j,k-1)*U(j,k-1);U(j,k)=Unew(j,k);U(5,5)=100.0;endend j=100for k=2:m-1Unew(j,k)=1/4*(ep(j+1,k)*U(j+1,k)+2*ep(j-1,k)*U(j-1,k)+ep(j,k+1)*U(j,k+1)+ep(j,k-1)*U(j,k-1);U(j,k)=Unew(j,k);endend,2018/10/18,本门课讲解的内容,梯度、散度和旋度介绍拉普拉斯方程和泊松方程均匀介质情况下非均匀介质情况下边界条件第一类边界条件第二类边界条件解拉普拉斯方程数值方法（有限差分法、有限元法）解析法（电像法、分离变量法、保角变换、格林函数等）,

注意事项

本文（matlab解拉普拉斯方程ppt课件）为本站会员（ZJ****3）主动上传，金锄头文库仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对上载内容本身不做任何修改或编辑。若此文所含内容侵犯了您的版权或隐私，请立即阅读金锄头文库的“版权提示”【网址:https://www.jinchutou.com/h-59.html】，按提示上传提交保证函及证明材料，经审查核实后我们立即给予删除！

温馨提示：如果因为网速或其他原因下载失败请重新下载，重复下载不扣分。