24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

柱子234

新虫 (小有名气)

应助: 6 (幼儿园)
金币: 94.8
散金: 20
红花: 4
帖子: 154
在线: 61.4小时
虫号: 3284156
注册: 2014-06-20
性别: GG
专业: 应用数学方法

[求助] 5点差分格式解边值问题。。。求大神帮忙已有1人参与

有没有类似的程序，哪位大神给我发一下，或者帮我写以下程序。

IMG_20151229_193059.jpg

回复此楼

» 猜你喜欢

1楼 2015-12-29 19:34:32

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

米米米BO

金虫 (初入文坛)

应助: 6 (幼儿园)
金币: 2065.3
红花: 4
帖子: 48
在线: 12.5小时
虫号: 2814649
注册: 2013-11-20
性别: GG
专业: 土壤学

引用回帖:

6楼: Originally posted by 米米米BO at 2015-12-30 10:13:51
clear,clc;
n=9;          % grid number
h=0.25;       % step length
u=zeros(n,n);
u(:,1)=1;    % top boundary
u(:,n)=0;    % bottom boundary
for k=1:1000 % maximum iteration
...

上下边界那个标注反了

赞一下

回复此楼

7楼2015-12-30 10:16:32

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 8 个回答

米米米BO

金虫 (初入文坛)

应助: 6 (幼儿园)
金币: 2065.3
红花: 4
帖子: 48
在线: 12.5小时
虫号: 2814649
注册: 2013-11-20
性别: GG
专业: 土壤学

把差分方程写出来，直接可以在Excel里面算

发自小木虫Android客户端

赞一下

回复此楼

2楼2015-12-29 19:55:20

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

柱子234

新虫 (小有名气)

应助: 6 (幼儿园)
金币: 94.8
散金: 20
红花: 4
帖子: 154
在线: 61.4小时
虫号: 3284156
注册: 2014-06-20
性别: GG
专业: 应用数学方法

引用回帖:

2楼: Originally posted by 米米米BO at 2015-12-29 19:55:20
把差分方程写出来，直接可以在Excel里面算

下面别人的程序，我该如何利用我的问题里面的边值条件呢？

function [ u x y k]=wudianchafenfa1(h,m,n,kmax,ep)
% g-s迭代法解五点差分法问题
%kmax为最大迭代次数
%m,n为x,y方向的网格数，例如（2-0）/0.01=200;aa
syms temp;
u=zeros(n+1,m+1);
x=-1+(0:m)*h;
y=-1+(0:n)*h;
for(i=1:n+1)
u(i,1)=sin(3*pi/4+pi*y(i)/4);
u(i,m+1)=sin(3*pi/4+pi*y(i)/4);
end
for(i=1:n)
for(j=1:m)
      f(i,j)=32*u(i,j);
end
end
t=zeros(n-1,m-1);
for(k=1:kmax)
for(i=2:n)
      for(j=2:m)
         temp=h*h*f(i,j)/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;
      end
end
t(i,j)=sqrt(t(i,j));
if(k>kmax)
      break;
end
if(max(max(t))<ep)
      break;
end
end

赞一下

回复此楼

3楼2015-12-29 21:42:43

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖