| 查看: 797 | 回复: 4 | ||
[求助]
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') |
回复此楼» 猜你喜欢
FileCode能看出啥？
已经有19人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有32人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有5人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有13人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
-
帮忙看看fileCode
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
2楼2015-05-03 22:48:07
【答案】应助回帖
★
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2015-05-04 13:38:52
感谢参与，应助指数 +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
4楼2015-05-05 23:45:21
5楼2015-05-05 23:49:06
