| 查看: 633 | 回复: 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 请教一下规定数组这么写合适吗？ |
回复此楼» 猜你喜欢
filecode
已经有14人回复
-
关于Filecode分析方法
已经有4人回复
-
FileCode能看出啥？
已经有20人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有35人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有3人回复
-
不应该看fileCode
已经有4人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有5人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有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
