Transition probability matrix for markov chain

16 次查看（过去 30 天）
John
John 2011-9-2
编辑： Walter Roberson 2020-11-28
Hi there
I have time, speed and acceleration data for a car in three columns. I'm trying to generate a 2 dimensional transition probability matrix of velocity and acceleration.
The concept is given a particular speed and acceleration I would like to know the next most likely (probable) speed and acceleration.
I have some code below, but cannot fully understand it. Will it generate a 2D transition matrix or a 4D transition matrix?
Thank you
%%First bin data into categories
speedBinN = 5;
aceelBinN = 5;
speed = binit( data(:,2), linspace(min(data(:,2)),max(data(:,2)),speedBinN) ); % bin them into categories
accel = binit( data(:,3), linspace(min(data(:,3)),max(data(:,3)),aceelBinN) );
%%count up transitions
transCountMat = zeros(speedBinN,aceelBinN,speedBinN,aceelBinN);
for ii = 1:size(data,1)-1
transCountMat( speed(ii),accel(ii),speed(ii+1),accel(ii+1) ) = transCountMat( speed(ii),accel(ii),speed(ii+1),accel(ii+1) ) + 1;
end
%%calculate probabilities
sumOverPossibleDestinations = sum( sum(transCountMat, 4), 3);
transMat = bsxfun( @rdivide, transCountMat, sumOverPossibleDestinations );
%%User Interactive stuff
IM = imagesc(squeeze(transMat(1,1,:,:)));
colorbar
set(IM,'ButtonDownFcn',@bdFcn)
set(gca,'ydir','normal')
ylabel speed
xlabel accel
hold on
p = plot(1,1,'w');
updateIndicator(1,1)

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

回答（6 个）

Walter Roberson
Walter Roberson 2011-9-2
Looks to me like it will generate a 2D output for transMat. The count matrix is 4 dimensional, but it is summed twice, which reduces that to 2 dimensions.
Looks to me like binit() is just the second output of histc(). With the linspace nature of the bins, that operation could probably be made more efficient than even histc(). Also the transcount loop could probably be replaced with a single accumarray call.
John
John 2011-9-3
Hi Walter,
Thank you for your reply and for explaining this to me. I'm happy that it generates a 2D output.
The code prints out the values of each element in the matrix. For example
val(:,:,1,3) =
NaN 0.2000 0.7692 0.0043 0
0 0 0.1294 0 NaN
0 0 0 0 NaN
NaN NaN 0 0 NaN
NaN NaN NaN 0 NaN
How can I make it just print out the actual matrix? I understand that this is probably a simple question but I don't know to do this.
I appreciate your help.
Thank you
  2 个评论
Walter Roberson
Walter Roberson 2011-9-3
There is no "val" in the code you show, so I am unsure what you are asking about?
When MATLAB displays a multidimensional matrix (4 dimensions in this case), it displays a "page" at a time, where a "page" is the first 2 dimensions. What format would you like the 4 dimensional matrix printed out in?
Walter Roberson
Walter Roberson 2011-9-3
Hmmm, looking again, it appears that transMat will be 4 dimensional, not the 2 dimensional that I thought.

请先登录，再进行评论。

John
John 2011-9-5
Hi Walter,
Would you know how to change the code to produce a 2D transition probability matrix of velocity and acceleration?
I've uploaded a sample data set to sendspace.com http://www.sendspace.com/file/oqx54r
And here is the full code
Thank you for your help
function transMat = mwExample
load('Sample data set mathworks.mat')
%%First bin data into categories
speedBinN = 5;
aceelBinN = 5;
speed = binit( data(:,2), linspace(min(data(:,2)),max(data(:,2)),speedBinN) ); % bin them into categories
accel = binit( data(:,3), linspace(min(data(:,3)),max(data(:,3)),aceelBinN) );
%%count up transitions
transCountMat = zeros(speedBinN,aceelBinN,speedBinN,aceelBinN);
for ii = 1:size(data,1)-1
transCountMat( speed(ii),accel(ii),speed(ii+1),accel(ii+1) ) = transCountMat( speed(ii),accel(ii),speed(ii+1),accel(ii+1) ) + 1;
end
%%calculate probabilities
sumOverPossibleDestinations = sum( sum(transCountMat, 4), 3);
transMat = bsxfun( @rdivide, transCountMat, sumOverPossibleDestinations );
%%User Interactive stuff
IM = imagesc(squeeze(transMat(1,1,:,:)));
colorbar
set(IM,'ButtonDownFcn',@bdFcn)
set(gca,'ydir','normal')
ylabel speed
xlabel accel
hold on
p = plot(1,1,'w');
updateIndicator(1,1)
Hanane
Hanane 2014-3-14
Hi; what is the role of this function ''binit''?
Harini pushparaj
Harini pushparaj 2017-12-27
@john..HI can i get your sample data file ..I maybe useful for my work ..I don't understand what is binit is ..it will be great if you could help me..Thank you!!
  3 个评论
显示 1更早的评论
Harini pushparaj
Harini pushparaj 2018-1-8
if i use binit ..i am getting error saying undefined binit and if i use histc..there is no error ..how come?
Walter Roberson
Walter Roberson 2018-1-8
You would need to store the above two lines in binit.m in your code directory.

请先登录，再进行评论。

shahab anjum
shahab anjum 2020-3-2
please help me too if i have 1000x286 matrix how can i calculate the transistion and emission probabilites of that matrix plz
help
  12 个评论
显示 10更早的评论
Walter Roberson
Walter Roberson 2020-3-13
编辑：Walter Roberson 2020-11-28
The way you processed your data originally is not compatible with the possibility that your array is an image.
There appear to be a number of techniques available for HMM processing of images for different purposes.
https://www.intechopen.com/books/reviews-refinements-and-new-ideas-in-face-recognition/hidden-markov-models-in-automatic-face-recognition-a-review
Attilio Pittelli
Attilio Pittelli 2020-11-28
you spoke about a list, but in the previous example he got the velocity and acceleration matrix. In my case i've got a double matrix of those inputs: a duration and a sequence of values. I'm trying to code it out but i'm still getting a transition probability matrix of 0 or 1 and the most of the outputs are 0 and it's not working. Do you have any suggestion? I'm really stucked at this point
Thank you

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Matrix Indexing 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by