| 查看: 660 | 回复: 1 | ||
百里道木虫 (初入文坛)
|
[求助]
有限差分的实现
|
|
n=200; L=1; t=0.1; Kn=1; %输入值 dx=L/n; w=dx/Kn; dt=w; m=t/dt; %规定数组 x=linspace(0,1,n); e1=[n,m]; e2=[n,m]; e=[n,m]; %计算过程 %初始值 for j=2:n; e1(j,1)=0; e2(j,1)=0; end %边界值 for k=1:m; e1(1,k)=0.5; e2(1,k)=0.5; e1(n,k)=0; e2(n,k)=0; end %内部过程 for j=2:n-1; for k=2:m; e1(j+1,k+1)=(1-w/2)*e1(j,k)+w/2*e2(j,k); e2(j-1,k+1)=(1-w/2)*e2(j,k)+w/2*e1(j,k); end end for j=1:n; for k=1:m; e(j,k)=e1(j,k)+e2(j,k); end end plot(x,e,'k-'); grid on 请教一下规定数组这么写合适吗？ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有7人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Matlab求解二阶偏微分方程组，希望能给出相应的求解代码
已经有10人回复
- 请问这个4阶的非线性偏微分方程组在PDE中怎么解？等待高手解答。
已经有17人回复
- 有限差分离散二阶偏微分方程后如何编译matlab程序解方程组？
已经有3人回复
- 求助，含有常数项的差分方程怎么转换成状态空间方程？
已经有5人回复
- matlab有限差分法求解复杂的偏微分方程组
已经有19人回复
- matlab 积分求助
已经有8人回复
- 积分方程和微分方程
已经有4人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】有人用有限差分方法计算过G点声子谱吗？
已经有10人回复
- 【求助】求微弱差分电流转差分电压芯片型号
已经有4人回复
- 100金币下载论文求助
已经有3人回复
- 【求助】请教虫友关于matlab差分方程的求解和画图程序
已经有4人回复
ltcc2000
木虫 (著名写手)
- 应助: 16 (小学生)
- 金币: 4074.5
- 散金: 365
- 红花: 3
- 帖子: 2330
- 在线: 154.5小时
- 虫号: 1647979
- 注册: 2012-02-27
- 性别: GG
- 专业: 控制理论与方法
2楼2012-04-09 08:37:43
50