| 查看: 631 | 回复: 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 请教一下规定数组这么写合适吗？ |
回复此楼» 猜你喜欢
为什么网上很多人说本周 12号出结果
已经有4人回复
-
基金中了
已经有15人回复
-
帮忙看看fileCode
已经有8人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有3人回复
-
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
静等基金结果
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
