| 查看: 2970 | 回复: 16 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
haixing2008荣誉版主 (文坛精英)
|
[求助]
高手进，用matlab求解一维抛物方程的初边值问题
|
|
本人初学者，请各位帮帮忙，若附上具体的matlab程序，本人将不胜感激。 [ 来自科研家族 流体数值模拟 ] |
回复此楼» 猜你喜欢
职称评审没过，求安慰
已经有54人回复
-
垃圾破二本职称评审标准
已经有19人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有5人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教高手...醌类物质显色问题
已经有12人回复
- charge的意义 高分求解
已经有5人回复
- 求高手出山...这三个极性大小大概是个怎样的关系
已经有5人回复
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- 有关磁流体方程计算的询问
已经有8人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】用MATLAB求解偏微分方程的数值解 悬赏20
已经有8人回复
- 【求助】使用Matlab预估动力学方程问题
已经有13人回复
haixing2008
荣誉版主 (文坛精英)
- 应助: 22 (小学生)
- 贵宾: 2.252
- 金币: 32603.8
- 散金: 2020
- 红花: 44
- 沙发: 1
- 帖子: 25068
- 在线: 429.1小时
- 虫号: 540264
- 注册: 2008-04-06
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 信息科学
4楼2012-05-07 22:25:24
【答案】应助回帖
★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
lovibond: 金币+1, 鼓励交流 2012-05-07 20:29:34
haixing2008: 金币+5, ★有帮助, 边界处理格式没有加进去，这是我的主要问题，呵呵 2012-05-07 22:26:54
感谢参与，应助指数 +1
lovibond: 金币+1, 鼓励交流 2012-05-07 20:29:34
haixing2008: 金币+5, ★有帮助, 边界处理格式没有加进去，这是我的主要问题，呵呵 2012-05-07 22:26:54
|
下面这个程序我运行过，可以了 %Finite different method for 1D parabolic equation %Backward Euler for Time, Center Different for Space function u=parabolicFD(m,n) h=1/m;%space size tao=1/n;%time step r=tao/h^2; u0=zeros(m+1,1); u1=zeros(m+1,1); A=zeros(m+1,m+1); F=zeros(m+1,1); for j=1:m+1 u0(j)=cos(pi*(j-1)*h); end %assembling the Matrix A for i=1:m+1 for j=1:m+1 if (abs(i-j)==1) A(i,j)=-1*r; end if (i==j) A(i,j)=1+2*r; end end end A(1,2)=-2*r; A(m+1,m)=-2*r; % for i=1:n for j=1:m+1 F(j)=u0(j)+tao*sin(tao*i); end u1=A\F; for j=1:m+1 u0(j)=u1(j); end end for j=1:m+1 u(j)=u1(j); end end 最后在命令窗口分别运行parabolicFD(40,1600)，parabolicFD(80,3200) |
2楼2012-05-07 19:13:26
cool_smile
木虫 (著名写手)
- 应助: 23 (小学生)
- 贵宾: 0.01
- 金币: 4916.4
- 散金: 3067
- 红花: 28
- 帖子: 1033
- 在线: 966.9小时
- 虫号: 503053
- 注册: 2008-02-15
- 性别: GG
- 专业: 计算数学与科学工程计算
3楼2012-05-07 19:44:01
5楼2012-05-07 23:02:21