一个循环程序，求指点

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

function Y=fic(E,N)
%find initial condition 寻求微粒的初值条件
%总能量为E，动能为T，势能为V，并且e=t+v
Y=zeros(N,4);
for i=1:N
T=unifrnd(0,E)
Y(i,2)=abs((2*T*cos(unifrnd(-pi,pi)))^0.5)
Y(i,4)=abs((2*T-Y(i,2)^2))^0.5
Y(i,1)=unifrnd(-pi,pi)
test=acos(E-T-2.5-1.5*cos(Y(i,1))/(1.5+0.5*cos(Y(i,1))))
   if isreal(test)==1
      Y(i,3)=test;
   else （这一块怎么写）
   end
end


我想要的结果是，如果test不是实数，那就重新计算第i步，否则才进行下一步，
求指点，谢谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：求指点 condition function Initial dition function 下一步 程序

function Y=fic(E,N)
%find initial condition 寻求微粒的初值条件
%总能量为E，动能为T，势能为V，并且e=t+v
Y=zeros(N,4);
for i=1:N
I_Flag = 0
Do while I_Flag = 0
T=unifrnd(0,E)
Y(i,2)=abs((2*T*cos(unifrnd(-pi,pi)))^0.5)
Y(i,4)=abs((2*T-Y(i,2)^2))^0.5
Y(i,1)=unifrnd(-pi,pi)
test=acos(E-T-2.5-1.5*cos(Y(i,1))/(1.5+0.5*cos(Y(i,1))))
   if isreal(test)==1
      Y(i,3)=test;
      I_Flag = 1
  end if
End Do While
end

postcardh 发表于 2013-8-20 04:11
function Y=fic(E,N)
%find initial condition 寻求微粒的初值条件
%总能量为E，动能为T，势能为V，并且e ...

你编写的是matlab程序吗，怎么提示错误

postcardh 发表于 2013-8-20 04:11
function Y=fic(E,N)
%find initial condition 寻求微粒的初值条件
%总能量为E，动能为T，势能为V，并且e ...

对了  非常感谢您的回答！思路很好！！！！

function Y=fic(E,N)
%find initial condition 寻求微粒的初值条件
%总能量为E，动能为T，势能为V，并且e=t+v
Y=zeros(N,4);
i=0;
while (i<N)
    i=i+1;
    T=unifrnd(0,E);
   Y(i,2)=abs((2*T*cos(unifrnd(-pi,pi)))^0.5);
   Y(i,4)=abs((2*T-Y(i,2)^2))^0.5;
   Y(i,1)=unifrnd(-pi,pi);
    test=acos(E-T-2.5-1.5*cos(Y(i,1))/(1.5+0.5*cos(Y(i,1))));
        if isreal(test)==1
            Y(i,3)=test;
        else
          i=i-1;
   end
end

bluethky 发表于 2013-8-20 20:22
function Y=fic(E,N)
%find initial condition 寻求微粒的初值条件
%总能量为E，动能为T，势能为V，并且e ...

运行成功，正在测试数据是否正确！