eig(X, 'nobalance') 的八度等效值是多少

发布于 2024-10-03 01:19:03 字数 419 浏览 10 评论 0 原文

我试图找到马尔可夫链的平衡分布，这意味着找到代表它的转移矩阵的特征值，但是，eig函数自动标准化它返回的特征向量，在MatLab中有一个标志你可以传递给函数停止这种行为

eig(X, '无平衡')

其中 X 是矩阵。请参阅 http://www.mathworks.com/help/techdoc/ref/eig。 html。但是，当我在八度音程中尝试此操作时，我收到一个错误：

错误：eig：错误类型参数“sq_string”

是否还有其他我应该调用的函数？

干杯

原文

分享到QQ

分享到微博

如果你对这篇内容有疑问，欢迎到本站社区发帖提问参与讨论，获取更多帮助，或者扫码二维码加入 Web 技术交流群。

发布评论

需要登录才能够评论，你可以免费注册一个本站的账号。

薄情伤 2024-10-10 01:19:03

如果您的目标是计算马尔可夫链的均衡分布，请查看 PMTK3 工具箱中的 rel="nofollow">mcStatDist 函数实现。它显示了计算结果的四种不同方法。示例：

TR = rand(3,3);                          %# random transition matrix
TR = bsxfun(@rdivide, TR, sum(TR,2));    %# normalize so that rows sum to one

[V,D] = eig(TR');                        %'# eigen-decomposition
EQ = V(:,1) ./ sum(V(:,1));              %# state equilibrium distribution

如链接代码的注释中所述，此方法在某些情况下可能在数值上不稳定，因此您可能需要考虑其他选项之一......

If your goal is to compute the equilibrium distribution of a Markov chain, take a look at the mcStatDist function implementation from the PMTK3 toolbox. It shows four different ways to compute the result. Example:

TR = rand(3,3);                          %# random transition matrix
TR = bsxfun(@rdivide, TR, sum(TR,2));    %# normalize so that rows sum to one

[V,D] = eig(TR');                        %'# eigen-decomposition
EQ = V(:,1) ./ sum(V(:,1));              %# state equilibrium distribution

As noted in the comments of the linked code, this method can be numerically unstable for some cases, so you might want to consider one of the other options...

回复收藏 0 原文

~没有更多了~