资源简介
自己编写的possion离散方法,给大家提供一些思路
代码片段和文件信息
function [p e v x y k]=possion(hMxMykmaxep) % g-s迭代法解五点差分法问题
%kmax为最大迭代次数
%mn为xy方向的网格数,例如(2-0)/0.01=200;
%e为误差,p为精确解
syms temp;
v=zeros(Mx+1My+1);
u1=zeros(Mx+1My+1);
u2=zeros(Mx+1My+1);
u3=zeros(Mx+1My+1);
f=zeros(Mx+1My+1);
x=0+(0:Mx)*h;
y=0+(0:My)*h;
w=4/pi*atan(0.2144*Mx+0.4822);
w=1;
for i=Mx/4:Mx*3/4
v(i1)=1500;
end
t=zeros(My-1Mx-1);
for k=1:kmax
for i=2:Mx
for j=2:My
%temp=(1-w)*u(ij)+w*(h*h*f(ij)/4+(u(ij+1)+u(ij-1)+u(i+1j)+u(i-1j)))/4;
temp=-h*h*f(ij)/4+(u(ij+1)+u(i
评论
共有 条评论