24小时热门版块排行榜    

查看: 795  |  回复: 4

小果果199128

新虫 (小有名气)

[求助] MATLAB初学者告急!!!程序运行不了?是哪里错了?卡在这动不了了,忘大神指点!!! 已有2人参与

clear all;
clc;
bx=[0 0 0 0 0 0 ;0 0 0  0 0 0; 1 0 0 0 0 0; 0 1 0 0 0 0;0 0 sqrt(2) 0 0 0;0 0 0 sqrt(2) 0 0]; %Creation
b=bx';%Anniation
n=bx*b;
imu=0;
dat=0;
for mu=-0.4:0.01:-1;
    mu
    imu=imu+1;
z=4;
y=1;%k
v=0;
im=eye(6);
   t=0.1;
    g=1;
    id=0;
    w=0;
    ppa=1;%psi
    roa=1;
    while true
        id=id+1;
        p=ppa;
       ro=roa;
         h=[0 0 -z*t*p 0 0 0;0 y-mu  g -z*t*p 0 0;-z*t*p g w-mu 0 -sqrt(2)*z*t*p 0;0 -z*t*p 0 w+y-2*mu sqrt(2)*g -sqrt(2)*z*t*p;0 0 -sqrt(2)*z*t*p sqrt(2)*g 2*w-2*mu 0;0 0 0 -sqrt(2)*z*t*p 0 2*w+y-3*mu]+z*t*p*p*im;
    [va,da]=eig(h);
        p=va(:,1)'*bx*va(:,1);%
        ro=va(:,1)'*n*va(:,1);
    if(id>200);
        break
    end
    end
    dat(imu,1)=z*t;
    dat(imu,2)=ppa;
    dat(imu,3)=roa;
    dat(imu,4)=mu;
end
plot(dat(:,4),dat(:,2),'-rs',dat(:,4),dat(:,3),'-gs')
错误提示如下:
Index exceeds matrix dimensions.

Error in jcm3 (line 57)
plot(dat(:,4),dat(:,2),'-rs')%,dat(:,4),dat(:,3),'-gs')
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsy12

银虫 (初入文坛)

【答案】应助回帖


感谢参与,应助指数 +1
jjdg: 金币+1, 感谢参与 2015-05-04 13:38:39
你的dat变量定义的有问题,应该定义成矩阵,而不应该赋值为0;不知道你的dat是几行几列的矩阵,我也不知道该改成什么样
2楼2015-05-03 22:48:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ziying007

新虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
jjdg: 金币+1, 感谢参与 2015-05-04 13:38:52
for循环的循环条件有问题。如果你在command window里检查一下的话,你会发现 mu=-0.4:0.01:-1会输出:
mu =
   Empty matrix: 1-by-0
所以,你的for没有起作用。而且你要把for的那个分号;去掉。
把for那里改成如下所示:
for mu=-0.4:-0.01:-1 %nBegin:nStep:nEnd  nStep要改成-1(负一)
    mu;
    imu=imu+1;
改后程序就能出图了。出了一个不知道是什么含义的图.....但是,我又不知道你要画什么
3楼2015-05-04 04:49:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小果果199128

新虫 (小有名气)

引用回帖:
2楼: Originally posted by lsy12 at 2015-05-03 22:48:07
你的dat变量定义的有问题,应该定义成矩阵,而不应该赋值为0;不知道你的dat是几行几列的矩阵,我也不知道该改成什么样

谢谢你啊,dat=0是定义初始值是0,后面我找见问题了,嘿嘿,就是呢个mu循环应该是-0.01。总之谢谢啦!
4楼2015-05-05 23:45:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小果果199128

新虫 (小有名气)

引用回帖:
3楼: Originally posted by ziying007 at 2015-05-04 04:49:47
for循环的循环条件有问题。如果你在command window里检查一下的话,你会发现 mu=-0.4:0.01:-1会输出:
mu =
   Empty matrix: 1-by-0
所以,你的for没有起作用。而且你要把for的那个分号;去掉。
把for那里改成 ...

嗯嗯,谢谢啦,老师说小木虫很好,果然啊,很温暖。
5楼2015-05-05 23:49:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小果果199128 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 12/600 2026-08-12 01:27 by 南京中西医
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +4 Tide man 2026-08-11 4/200 2026-08-12 00:39 by liuqidi
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 9/450 2026-08-11 22:12 by 医学老男孩
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] filecode +9 documentary 2026-08-10 9/450 2026-08-11 10:08 by wk7465
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
信息提示
请填处理意见