How to coding and plotting

2014 2 26
1 个回答
回答已采纳
更新时间：2014 2 26
5 次查看（30 天）

0 个投票

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 ?

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

 采纳的回答

Mischa Kim
Mischa Kim 2014-2-26
编辑：Mischa Kim 2014-2-26

0 个投票

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更早的评论

dab483
dab483 2014-2-26
eg: now the result for pp=
pp= 0 0 0 0 0.3 0.4 0.5 0.5 0.5 0.5
i want to make where pp>0 it will become
new pp = 0 0 0 0 1 1 1 1 1 1 %%(and plot this)
using
if pp>0
new pp=1;
else
new pp=0;
end
but all i got is all zeros. any advice? Thanks
Mischa Kim
Mischa Kim 2014-2-26
编辑：Mischa Kim 2014-2-26
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
dab483 2014-2-26
ok. i think i got it! thanks for your help.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Annotations 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by