I have the following code. The values of displ in stage1 and stage2 are plotted in displ vs time graph. Thus a single curve is obtained.

Question

Lakshmikruthiga Ponnusamy 2019-1-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/438571-i-have-the-following-code-the-values-of-displ-in-stage1-and-stage2-are-plotted-in-displ-vs-time-gra

评论： Walter Roberson 2019-1-8

displ=zeros(100,2)
tim=zeros(100,1)
for i=2:100
v=5;
dt=1*10^-4;
function [displ_f]=stage1(displ_i)
    dt=1*10^-4;
    v=5;
   displ_f=displ_i+v*dt;
endfunction
function [displ_f]=stage2(displ_i)
    dt=1*10^-4;
    v=5;
    displ_f=displ_i+v*dt;
endfunction
if (displ(i-1,1)<0.020 && displ(i-1,2)>=0)
        disp("IN stage 1");
        [displ(i,1)]=stage1(displ(i-1,1));
        displ(i,2)=v;
        tim(i,1)=tim(i-1)+dt;
    elseif (displ(i-1,1)>=0.020 && displ(i-1,1)<=0.0256 && displ(i-1,2)>=0)
        disp("IN stage 2");
        [displ(i,1)]=stage2(displ(i-1,1));
        displ(i,2)=displ(i-1,2);
        tim(i,1)=tim(i-1)+dt;
end
end
disp(" displ vs tim")
plot(tim(:,1),displ(:,1));

How will I get stage1 curve in orange color and stage2 curve in blue color in the same plot to know which stage has which displ values?

5 个评论
显示 3更早的评论隐藏 3更早的评论

Lakshmikruthiga Ponnusamy 2019-1-8

yeah, it is octave. Still, you can help. Matlab and octave both are similar. I didn't get any help from octave forum, so I posted it here. Kindly help.

Walter Roberson 2019-1-8

Octave is a product of the Free Software Foundation . The fundamental principle of FSF is that the software should be no cost, and that you can modify it as needed , but that if you want service then you should pay for the service . Therefore if you are not satisfied with the response you got from the Octave distribution list, you should be hiring a consultant . That is the FSF way .

FSF is not a charity to provide free software because software can be hard to afford .FSF is a deliberately disruptive economic model that aims to drive proprietary software out of business, replacing it with no cost for the software and instead charging for work on the software.

FSF says that it is morally and ethically wrong to create software and not give the software away for no charge, that anyone who creates proprietary software is a Bad Person who is Doing Evil.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

nanren888 2019-1-7

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/438571-i-have-the-following-code-the-values-of-displ-in-stage1-and-stage2-are-plotted-in-displ-vs-time-gra#answer_355404

编辑：madhan ravi 2019-1-8

在 MATLAB Online 中打开

Is this the sort of thing you are looking for?

displ=zeros(100,2)
tim=zeros(100,1)
stageOneMask = false([100,1]);
stageTwoMask = false([100,1]);
for i=2:100
    v=5;
    dt=1*10^-4;
    if (displ(i-1,1)<0.020 && displ(i-1,2)>=0)
        disp("IN stage 1");
        [displ(i,1)]=stage1(displ(i-1,1));
        displ(i,2)=v;
        tim(i,1)=tim(i-1)+dt;
        stageOneMask(i) = true;
    elseif (displ(i-1,1)>=0.020 && displ(i-1,1)<=0.0256 && displ(i-1,2)>=0)
        disp("IN stage 2");
        [displ(i,1)]=stage2(displ(i-1,1));
        displ(i,2)=displ(i-1,2);
        tim(i,1)=tim(i-1)+dt;
        stageTwoMask(i) = true;
    end
    disp(" displ vs tim")
    plot(tim(:,1),displ(:,1),'k',tim(stageOneMask,1),displ(stageOneMask,1),'Ob',tim(stageTwoMask,1),displ(stageTwoMask,1),'Or');
    
end
function [displ_f]=stage1(displ_i)
    dt=1*10^-4;
    v=5;
   displ_f=displ_i+v*dt;
end
function [displ_f]=stage2(displ_i)
    dt=1*10^-4;
    v=5;
    displ_f=displ_i+v*dt;
end

Generally, if you know something about the structure of the problem, for example, that it has two distinct stages, maybe your code could represent that.

Sorry, I had to rearrange your code to get it sort of Matlab compliant. (Leaves it a bit of a mess, sorry)

Might want to consider rearranging it to be more in line with normal style guides?

Did I guess correctly what your "stages" were?