应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 求助 急！！！谢谢
31/1返回列表
查看: 201  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

初值及边界设置



算法



问题方程
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
很高兴认识你
1楼 2011-11-02 19:40:41
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

chenyg07

金虫 (初入文坛)
感觉边界条件不对 一段时间就出现锯齿了  还请各位大神相救
赞一下 回复此楼
很高兴认识你
2楼2011-11-02 19:42:04
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

chenyg07

金虫 (初入文坛)
★
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);
赞一下(1人) 回复此楼
很高兴认识你
3楼2011-11-02 20:19:14
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chenyg07 的主题更新
31/1返回列表
信息提示
关闭
请填处理意见
关闭