| 查看: 250 | 回复: 2 | ||
chenyg07金虫 (初入文坛)
|
[求助]
求助 急！！！谢谢
|
|
题干： 我的程序为什么不行 clc clear sigma=5; a=1; D=1/(sqrt(2*pi)*sigma); L=2*pi; T=100; dett=0.1; N=100; detx=L/N; miu=30; for i=1:N+3 u(i)=100*D*exp(-(i-miu)*(i-miu)/(2*sigma*sigma)); end for t=1:7 for j=2:1:N+2 i=j; m(j)=u(i)-a*dett/(2*detx)*(u(i+1)-u(i-1))+... a*a*dett*dett/(2*detx*detx)*(u(i+1)-2*u(i)+u(i-1)); if (m(j)<0.001) m(j)=0.0; end end m(1)=m(N+1); m(N)=m(2); % for i=2:N+2 u(i)=m(i); end t plot(m) pause(0.2) end  初值及边界设置  算法  问题方程 |
回复此楼» 猜你喜欢
关于filecode
已经有6人回复
-
关于代码变化问题，想知道的进来
已经有8人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有7人回复
-
关于豆爷回答的JTJC与%2F数量
已经有6人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有31人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助/交流】急，求助关于琼脂糖凝胶的问题，谢谢！
已经有8人回复
- 求助，投稿要求中的表格问题，头大，急，谢谢帮助！
已经有6人回复
- 求助翻译 急谢谢
已经有2人回复
- 急。求助翻译。质量一般就行，不要机译。谢谢
已经有4人回复
- 【求助】关于国内外导师介绍签名问题（急，谢谢大家啦！）
已经有16人回复
- 急求助英语翻译，谢谢!
已经有1人回复
- 翻译求助 急用 谢谢
已经有3人回复
- 【求助】请教图像处理问题，急，谢谢各位了
已经有7人回复
- 【求助】求2个化合物的合成，急，谢谢各位大虾了
已经有11人回复
chenyg07
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1370.9
- 散金: 6
- 帖子: 28
- 在线: 47.8小时
- 虫号: 1162865
- 注册: 2010-12-05
- 性别: GG
- 专业: 等离子体物理
2楼2011-11-02 19:42:04
chenyg07
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1370.9
- 散金: 6
- 帖子: 28
- 在线: 47.8小时
- 虫号: 1162865
- 注册: 2010-12-05
- 性别: GG
- 专业: 等离子体物理
3楼2011-11-02 20:19:14
100