24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

chenyg07

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1370.9
散金: 6
帖子: 28
在线: 47.8小时
虫号: 1162865
注册: 2010-12-05
性别: GG
专业: 等离子体物理

[求助] 求助急！！！谢谢

题干：
我的程序为什么不行

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

初值及边界设置

算法

问题方程

回复此楼

» 猜你喜欢

西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向已经有3人回复
救命帖已经有6人回复
存款400万可以在学校里躺平吗已经有37人回复
英文综述是否需要润色及查重已经有5人回复
限项规定已经有5人回复
为什么nbs上溴没有产物点出现呢已经有9人回复
招博士已经有3人回复
最失望的一年已经有18人回复

» 本主题相关价值贴推荐，对您同样有帮助:

很高兴认识你

1楼 2011-11-02 19:40:41

已阅关注TA 给TA发消息送TA红花 TA的回帖

chenyg07

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1370.9
散金: 6
帖子: 28
在线: 47.8小时
虫号: 1162865
注册: 2010-12-05
性别: GG
专业: 等离子体物理

★
xiegangmai(金币+1): 鼓励讨论交流！ 2011-11-02 20:43:24

谢谢，别人帮我解决了，把
m(1)=m(N+1);
m(N)=m(2);
%

   for i=2:N+2
      u(i)=m(i);
   end
换成
m(1)=m(N+1);
m(N+2)=m(2);
%

   for i=2:N+2
      u(i)=m(i);
   end
  m(1)=u(1);