• 大小: 1KB
    文件类型: .m
    金币: 1
    下载: 0 次
    发布日期: 2021-06-01
  • 语言: Matlab
  • 标签: possion  中心差分  

资源简介

自己编写的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

评论

共有 条评论