24小时热门版块排行榜    

查看: 800  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 长年满屏的广告,版主太不责任了。基金也等的急 +3 gltch 2026-08-12 3/150 2026-08-12 13:25 by dcqxinyang
[基金申请] filecode +11 documentary 2026-08-10 12/600 2026-08-12 12:42 by kudofaye
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 19/950 2026-08-12 10:55 by 要乐观耀哥
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +5 Tide man 2026-08-11 5/250 2026-08-12 09:24 by kangban
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +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 布布和一二
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 基金中了 +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
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 这样的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 仁砚薪传
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
信息提示
请填处理意见