matlab 拉普拉斯求解,matlab解拉普拉斯方程.ppt

均匀介质情况下 非均匀介质情况下 ep： 11111111 11111111 11111111 11110101011 11110101011 11110101011 11111111 11111111 11111111 11111111 非均匀介质情况下 非均匀介质情况下 ep： for k=1:mfor j=1:mep(j,k)=1;end end for k=4:6for j=4:6ep(j,k)=10;end end 非均匀介质情况下 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)); 非均匀介质情况下 非均匀介质情况下 非均匀介质情况下 ep： 1111111111 1111111111 1111111111 10101010101010101010 10101010101010101010 10101010101010101010 10101010101010101010 100100100100100100100100100100 100100100100100100100100100100 100100100100100100100100100100 第一类边界条件 clear all clc m=10 for k=1:mfor j=1:mU(j,k)=0;end end for j=1for k=1:mU(j,k)=10;end end 第二类边界条件 for i=1:200 for 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;end endj=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);end end 本门课讲解的内容 梯度、散度和旋度 介绍拉普拉斯方程和泊松方程 均匀介质情况下 非均匀介质情况下 边界条件 第一类边界条件 第二类边界条件 解拉普拉斯方程 数值方法(有限差分法、有限元法) 解析法(电像法、分离变量法、保角变换、格林函数等) 拉普拉斯方程有限差分解 均匀介质情况下 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:200 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;end end 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 =0000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000 U(5,5)=100.0 U =00000000000000000000000000000000000000000000 100 0000000000000000000000000000000000000000000000000000000 for i=1:200 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)