24小时热门版块排行榜    

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

金虫 (初入文坛)


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的回帖
查看全部 3 个回答

chenyg07

金虫 (初入文坛)

感觉边界条件不对 一段时间就出现锯齿了  还请各位大神相救
很高兴认识你
2楼2011-11-02 19:42:04
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 澳大利亚 Murdoch University 全奖博士招生(3个名额)地质化工冶金领域 +9 AI8RaGaPaSCR 2026-08-07 9/450 2026-08-13 15:45 by GGlceF2k3sb4
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +5 医学老男孩 2026-08-13 11/550 2026-08-13 15:41 by sdfapple719
[基金申请] 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
[基金申请] FileCode能看出啥? +8 要乐观耀哥 2026-08-10 26/1300 2026-08-13 12:01 by Tide man
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 11/550 2026-08-13 10:40 by ZH远山
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +9 家与远方 2026-08-10 14/700 2026-08-13 09:41 by 臭臭不臭01
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 奇怪,两个人的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 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
信息提示
请填处理意见