24小时热门版块排行榜    

查看: 664  |  回复: 3

woshishishi

铁虫 (初入文坛)

[求助] matlab程序运行不出来,求大神指点一二 已有2人参与

function u=b()
dt=0.1;      %时间间隔
n=201;       %网格数
h=2/(n-1);   %空间步长
m=2/dt;      %迭代次数
u=ones(m,n);
u(:,1)=0;
u(:,n)=0;

    for j=2:n-1
       u(1,j)=sin(pi*j*h);
    end  
    for i=1:m
        for j=2:n-1
            u(i+1,j)=u(i,j)-u(i,j)*dt/2*h*(u(i,j+1)-u(i,j-1))+0.01*(dt/h^2)*(u(i,j+1)-2*u(i,j)+u(i,j-1));
        end
    end
end


运行出来的结果:(部分)为什么第二行算不出来?
Columns 1 through 7

         0    0.0628    0.0941    0.1253    0.1564    0.1874    0.2181
         0         0         0         0         0         0         0
         0         0         0         0         0         0         0
         0         0         0         0         0         0         0
         0         0         0         0         0         0         0
         0         0         0         0         0         0         0
         0         0         0         0         0         0         0
         0         0         0         0         0         0         0
         0         0         0         0         0         0         0
         0         0         0         0         0         0         0
         0         0         0         0         0         0         0
         0         0         0         0         0         0         0
         0         0         0         0         0         0         0
         0         0         0         0         0         0         0
         0         0         0         0         0         0         0
         0         0         0         0         0         0         0
         0         0         0         0         0         0         0
         0         0         0         0         0         0         0
         0         0         0         0         0         0         0
         0         0         0         0         0         0         0
         0         0         0         0         0         0         0
回复此楼

» 收录本帖的淘帖专辑推荐

matlab编程绘图

» 猜你喜欢

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

好好科研 好好生活
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +1
数组运算 用点乘“.*”和点除“./”。。。。运算符号前没有“.”代表矩阵运算。。。
数值计算
2楼2015-07-12 09:14:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

woshishishi

铁虫 (初入文坛)

引用回帖:
2楼: Originally posted by 独孤神宇 at 2015-07-12 09:14:20
数组运算 用点乘“.*”和点除“./”。。。。运算符号前没有“.”代表矩阵运算。。。

我的运算不是数组,是元素,也需要点乘点除吗,我加上了,还是没有结果
好好科研 好好生活
3楼2015-07-12 21:12:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

littleMaybe

铁虫 (初入文坛)

【答案】应助回帖

我觉得是不是你for语句中的顺序写错了?最好step一步步运行一下 ==看看是这么个变化过程。笔者也编写过类似的差分方程,算的不对往往就是因为for语句中的条件给弄错了。
嘟嘟噜~
4楼2015-08-07 08:41:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 woshishishi 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 一志愿吉林大学材料学硕321求调剂 +4 Ymlll 2026-03-18 6/300 2026-03-18 22:15 by li123456789.
[考研] 321求调剂 +3 何润采123 2026-03-18 3/150 2026-03-18 21:27 by li123456789.
[考研] 267一志愿南京工业大学0817化工求调剂 +8 SUICHILD 2026-03-12 8/400 2026-03-18 20:55 by winsuccess
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +7 步川酷紫123 2026-03-13 7/350 2026-03-18 17:12 by 尽舜尧1
[考研] 286求调剂 +6 lemonzzn 2026-03-16 9/450 2026-03-18 15:31 by vgtyfty
[考研] 311求调剂 +6 26研0 2026-03-15 6/300 2026-03-18 14:43 by haxia
[考研] 304求调剂 +12 小熊joy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
信息提示
请填处理意见