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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 求助 急！！！谢谢
31/1返回列表
查看: 290  |  回复: 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返回列表
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 12:24 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 12:24 by BCdvBAEpt3zm
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 15/750 2026-11-01 12:11 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 11:14 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 10:53 by BCdvBAEpt3zm
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 6/300 2026-11-01 10:53 by BCdvBAEpt3zm
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 10:27 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 07:01 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[基金申请] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 05:41 by YXM1FDCSE8Bn
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 05:41 by YXM1FDCSE8Bn
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 04:04 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 3/150 2026-11-01 04:04 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +8 zuesli 2026-10-27 14/700 2026-10-30 22:44 by 好运连连公主
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
关闭
请填处理意见
关闭