How to coding and plotting

Question

dab483 2014-2-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/119182-how-to-coding-and-plotting

评论： dab483 2014-2-26

Hi i have this Page hinkley test

sumPH=0
  for t=1:Nt
       mew2=(inv(Nt)*sum((y2-xhat2)));
       vari2=(inv(Nt-1)*sum((xhat2-mew2).^(2)));
       sumPH= sumPH + vv/vari2*(rg12 - mean(mew2) -vv/2);
       pp=max(0,sumPH);
    end

Now i want to plot (t,pp) using either 1 or zero. How do i coding if i want to make maximum of pp become 1 instead its own value of sumPH ?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Mischa Kim 2014-2-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/119182-how-to-coding-and-plotting#answer_126202

编辑：Mischa Kim 2014-2-26

在 MATLAB Online 中打开

You mean normalize the plot?

 sumPH=0
 for t=1:Nt
    mew2  = (inv(Nt)*sum((y2-xhat2)));
    vari2 = (inv(Nt-1)*sum((xhat2-mew2).^(2)));
    sumPH = sumPH + vv/vari2*(rg12 - mean(mew2) -vv/2);
    pp(t) = max(0,sumPH);         % save each data point
 end
 plot(t,pp/max(pp))               % plot normalized data

If you want to plot only 1s and 0s, then you would use

    pp(t) = sumPH>0;              % save each data point
    ...
 plot(t,pp)                       % plot data

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

Mischa Kim 2014-2-26

编辑：Mischa Kim 2014-2-26

在 MATLAB Online 中打开

Yep, as outlined above, use

 sumPH=0
 for t=1:Nt
    mew2  = (inv(Nt)*sum((y2-xhat2)));
    vari2 = (inv(Nt-1)*sum((xhat2-mew2).^(2)));
    sumPH = sumPH + vv/vari2*(rg12 - mean(mew2) -vv/2);
    pp(t) = sumPH>0;     % save each data point
 end
 plot(t,pp,'*')          % plot data

...and make sure that t and pp are both row or both column vectors. Note, with '*' parameter the data points are plotted individually using a star as the marker symbol.

dab483 2014-2-26