24小时热门版块排行榜    

查看: 260  |  回复: 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

初值及边界设置



算法



问题方程
回复此楼
很高兴认识你
已阅   关注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);
很高兴认识你
3楼2011-11-02 20:19:14
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chenyg07 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +11 医学老男孩 2026-08-13 29/1450 2026-08-14 00:04 by 大快人心啦
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +9 Tide man 2026-08-10 10/500 2026-08-13 20:33 by 温淼淼淼
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 28/1400 2026-08-13 17:51 by 医学老男孩
[基金申请] filecode +14 documentary 2026-08-10 16/800 2026-08-13 16:39 by 润生71
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +9 家与远方 2026-08-10 14/700 2026-08-13 09:41 by 臭臭不臭01
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
信息提示
请填处理意见