1回答
卡尔曼滤波:一些疑问
、、、、
我有几个问题：
在给出的openCV文档中：
/*生成测量*/ cvMatMulAdd(卡尔曼->测量矩阵、状态、测量、测量)；
这是正确的吗？在本教程中:韦尔奇和毕晓普在方程1.2中的表示测量=H*状态+测量噪声
似乎两者不一样。
我试着为一个球实现。我尝试了以下几点：(请指出我是否做错了。)
对于测量，我测量两件事：( a) x，b) y，球的质心。
我只是提到与opencv文档中给出的示例不同的行。
CvKalman* kalman = cvCreateKalman( 5, 2, 0 );
const float A[] = { 1, 0, 1, 0, 0,
浏览 0提问于2011-05-06得票数 4
假设我有一个矩阵A，B和C，当我将矩阵初始化为
A = np.array([[]])
B = np.array([[1,2,3]])
C = np.array([[1,2,3],[4,5,6]])
然后
B.shape[0]
C.shape[0]
分别给予1和2 (如预期的)，但
A.shape[0]
给出1，就像B.shape。
获得给定矩阵的行数的最简单方法是什么，但仍然可以确保像A这样的空矩阵给出的值为零。
搜索堆栈溢出一段时间后，我找不到答案，所以我在下面发布我自己的答案，但是如果你能想出一个更干净、更一般的答案，我会接受你的答案。谢谢!
浏览 6提问于2016-04-04得票数 0
回答已采纳
1回答
我有一个100000x8矩阵，我需要将每一行转换成一些符号表示。我有这样的工作没有问题，但做我的矩阵的大小，我想要矢量化的解。
以下是我所做的工作：
cutoff = [-inf -1.15 -0.67 -0.32 0 0.32 0.67 1.15];
for k = 1 : length(matrix)
tmp = matrix(k,:);
for i = 1 : length(tmp)
string(i) = sum( (cutoff <= tmp(i)), 2 );
end
end
例如，如果我矩阵的第1行是0.1053，0.1053，0.1053，
浏览 3提问于2014-01-02得票数 0
回答已采纳
1回答
我想存储一个矩阵列表。为此，我使用matrix构造矩阵，并使用list将其连接起来。本质上，假设我有：
A = matrix(1,2,2)
B = list(A,A)
那么我如何仅从B获得A的尺寸？我可以使用B1获取矩阵元素，但无法提取维度。请帮帮忙。
浏览 2提问于2016-08-31得票数 0
1回答
我得到了一个一万乘六十四个整数矩阵1，2，…，16，通过我之前写的一些程序来识别。这个矩阵表示10，08乘8幅矢量化的图像。
我的问题不在于如何识别它们(我已经讨论了这个部分)，而是更多关于我正在查看的数据的问题。一个带有矩阵的简短的.txt文件给了我一点解释：
32×32位图被划分为4乘4的不重叠块，每个块中的像素数被计数。这将生成一个8乘8的输入矩阵，其中每个元素都是范围为0的整数。。。、16。
不过我还是有点困惑。32乘32位图与8乘8的图像是如何对应的？每个整数代表什么？如何将给定的图像转换为矢量图像，然后将其创建为数组？在这种情况下，位图和图片之间的关系是什么？
感谢您提供的
浏览 2提问于2015-05-03得票数 0
四元数(复数的四维扩展)可以用来表示三维向量的旋转和缩放，而四元数在三维向量上的应用需要两个四元数乘法，因此所需的运算比通过相应的变换矩阵进行乘法所需的运算少。然而，通常使用线性和仿射变换矩阵，特别是在着色码中。
(由于速度、稳定性等原因)什么时候合适和更好？用四元数来表示三维的尺度和旋转，而不是相应的变换矩阵？
浏览 0提问于2015-08-09得票数 16
1回答
我试图将两个复矩阵(b+c*i)相乘，但没有得到结果。
??? Error using ==> mtimes
Inner matrix dimensions must agree.
Error in ==> impedancaZ at 14
l=mtimes(R1,h)
我不理解这个错误，因为矩阵维数是相同的(2 coluoms和9行)
你能帮帮我吗?
浏览 1提问于2014-03-14得票数 0
回答已采纳
1回答
矩阵与向量的乘积
、、
我试图在Matlab中乘2个矩阵，但是它们没有相同的维数。事实上，我的乘法不是标准的矩阵乘法。我有一个31乘1矩阵(或向量)和31乘512乘512矩阵.我想取我的向量的第一个元素，然后用它乘以第一个512乘以512，得到一个新的31乘512乘以512的数组。但我不想使用for循环来解决性能问题。
浏览 1提问于2015-10-06得票数 0
2回答
解一个约束线性方程组
、、、、
我有一个形式为y=Ax+b的方程组，其中y、x和b为n×1向量，A<code>E 29</code>是n×n (对称)矩阵。
这是皱纹。并非所有的x都是未知的。指定x的某些行，并未知y的相应行。下面是一个例子
| 10 | | 5 -2 1 | | * | | -1 |
| * | = | -2 2 0 | | 1 | + | 1 |
| 1 | | 1 0 1 | | * | | 2 |
其中*指定未知的数量。
我已经在Fortran中为上述问题建立了一个解决程序，但是我想知道是否是一个很好的健壮的解决方案--作为La
浏览 3提问于2017-04-06得票数 5
我对tensorflow有点陌生，我只想知道它是否已经在它的操作符中执行了向量化，比如乘法矩阵等等。如果不可行的话。
举个例子，假设我们有两个矩阵，我们想要把它们相乘。我正在寻找一种矢量化的方法来完成这个操作，而不是for循环。
提前谢谢。
浏览 0提问于2017-01-20得票数 0
回答已采纳
1回答
我正在尝试用C + +编写一个算法，我使用了特征库。
问题是我需要用C++解决下面的方程：
Find the coefficients: s1,..., sn ∈ R such that a(-1)ψi = s1ψ1 +...+snψn
ψi:矩阵是已知的。
浏览 0提问于2013-05-13得票数 2
回答已采纳
1回答
我有一个矩阵，它被分配给变量a。
矩阵的大小小于112* 92 uint8>。
因此，我运行了以下代码
b = a'; c = b .*a;
但是，据说维数不匹配，我无法解决原因。
浏览 2提问于2014-03-05得票数 0
回答已采纳
2回答
我有一个数据集，包括到约1800个融合雷达触点的距离和方位，以及到这些触点的实际距离和方位，我需要开发一个校正方程，以使感知的值尽可能接近实际值。 在可视化时，错误似乎有一种趋势，所以在我看来，应该有一个稍微简单的等式来纠正它。 这是~1800方程式的形式： actual_distance = perceived_distance + X(percieved_bearing) + Y(speed_over_ground) + Z(course_over_ground) +A(标题) 求解X、Y、Z和A的最佳方法是什么？ 此外，我不相信所有这些因素都是必要的，所以我完全愿意省略其中的一两个因素
浏览 22提问于2019-04-16得票数 1
我有一个立方体m乘以k，这意味着我有k个矩阵m乘以n，我想把它重塑成一个大的矩阵，例如p行和q列(假设每个m乘以n个矩阵作为一个元素)。
我怎么能这么做？我可以使用整形功能吗？例如，这些矩阵组合在一起时形成一个立方体：
1 1 , 2 2 , 3 3 , ... , 16 16
1 1 2 2 3 3 16 16
在上面的示例中，k=16，m=n=2。
我想像这样重塑他们：
1 1 5 5 9 9 13 13
1 1 5 5 9 9 13 13
2 2 6 6 10 10 14 14
2 2 6 6 10 10 14 14
3 3 7 7 11 11 15 15
浏览 0提问于2014-04-13得票数 2
回答已采纳
我的目标是求解满足A*B=C的矩阵A，其中C是已知的，B是随机生成的。以下是一个例子：
C=[1/3 1/3 1/3]'*[1/3 1/6 1/6 1/6 1/6];
B=rand(5,5);
A=C*pinv(B);
A*B=C_test;
norm(C-C_test);
ans =
4.6671e-16
在这里，C_test的元素与原始C在1e-15之间，但是当B的行数少于列时，错误急剧增加(不确定norm()是否是显示错误的最佳方式，但我认为它说明了问题)。例如：
B=rand(4,5);
A=C*pinv(B);
A*B=C_test;
norm(C-C_test);
ans
浏览 3提问于2014-10-02得票数 1
回答已采纳
假设我们有一个4乘4的矩阵A，它表示某种变换.
我们可以用这个矩阵对向量进行两种方式的变换。
Ax假设x是4×1列向量.
xA假设x是1乘4行向量.
直到昨天，我还没有看到这个xA变体，我认为它是不必要的。然后我遇到了此DirectX代码来自微软。
该函数只接受投影矩阵，计算逆投影矩阵，然后执行这个xA变量(而不是Ax)来转换(或取消项目)一组坐标。
一开始，我没有看，做了Ax，这是荒谬的，现在我试图理解为什么。
我注意到了Ax = xA^T，但是当我看到矩阵乘法的执行方式时。组成结果的成分是非常不同的，这取决于我是使用Ax还是xA (应该如此，矩阵乘法不是可交换的)。
Ax (3×3矩阵和行
浏览 0提问于2016-11-25得票数 1
回答已采纳
1回答
我想使用SciLab工具在大图像中隐藏一个图像，下面是我使用的代码片段
S1_diag = diag(s1);
S2_diag = diag(s2);
S1_diag(1:length(s1), :) = S2_diag(1:length(s1), :);
其中s1和s2是图像1和2的奇异对角线矩阵。
相同的代码在Matlab中工作，但在SciLab中生成一个“无效索引”错误(21)。我错过了什么？
我是SciLab语法的新手，所以无法理解如何在SciLab中解决这个问题。
任何帮助都是非常感谢的。
浏览 3提问于2014-09-04得票数 0
我试图计算出在一个C函数的disas中的什么地方实现了% 4。这是一个未经优化编译的程序。下面是C语言中的函数：
uint function (uint len, int* intArray) {
uint i; uint sum = 0;
for (i = 0; i < len-2; i++)
switch (i % 4)
{
case 0 :
sum += countPairsUpTo(len-2,intArray,8,6);
break;
case 1 :
sum += coun
浏览 3提问于2014-04-03得票数 0
2回答
图形行与列的主要转换
、、、、
这使我困惑了很长时间。为什么大多数图形API(带有HLSL和OpenGL的DirectX)通常表示位置和法线向量都是行向量？这与数学课上的典型形式正好相反。两者在数学上都是正确的，我只是好奇为什么图形学的人选择这种表示法。让我解释一下..。
向量和矩阵乘法的效率：
在我看来，行向量x和变换矩阵'A‘被相乘为x*A，从而导致`y=x*a的存在。
y[0] = x[0]*M[0][0] + x[1]*M[1][0] + x[2]*M[2][0] + x[3]*M[3][0];
y[1] = x[0]*M[0][1] + x[1]*M[1][1] + x[2]*M[2][1] + x[3]*
浏览 9提问于2014-12-09得票数 1
回答已采纳
我试着用FFT算法来绘制一个简单的音频语音文件的图形，它的大小与频率成正比，但我试图通过窗口来消除主频率旁瓣上的频率泄漏。只是给了我一个错误
下面是不加窗口的代码(这是工作的)：
[x,fs]=audioread('E:\voice_hum.wav');
x1=x;
plot_points = 100000;
h=fft(x1,plot_points);
freq = h(1:50001); %dump data after half FS
mag = abs(freq);
freq = 0.5 * fs *(0:50000)/50000; %frequency
浏览 1提问于2013-12-07得票数 0
我正在以下列方式实现具有投影失真的图像的公制校正：
从原始图像中，我找到了两组平行线，并找到了它们的交点(无穷远处的消失点)。
我选择一个圆上的五个非共线点来拟合一个圆锥曲线，然后用上面提到的点来检查这个圆锥曲线在无穷远处的交点位置。
利用这些点求出扭曲的对偶退化二次曲线。
理论上，由于畸变二次曲线是由C*'=HC*H‘(C*是对偶退化二次曲线，’是转置的，H是我的同形)确定的，那么我应该能够运行SVD来确定H.不失真，C*是一个3x3恒等矩阵，在对角线上有最后一个元素。但是，如果我运行SVD，就不会在对角线矩阵中得到它们。对于某些矩阵，我可以通过使用Cholesk
浏览 1提问于2012-09-16得票数 1
回答已采纳
我有一个GLSL着色器，工作很好，直到我添加以下的线条到它。
"vec4 coord = texture2D(Coordinates,gl_TexCoord[0].st);"
"mat4x3 theMatrix;"
"theMatrix[1][1] = 1.0;"
"vec3 CampProj = theMatrix * coord;"
当我检查错误日志时，我被告知：
"ERROR: 0:2: '=' : cannot convert from '4-component vector of flo
浏览 7提问于2014-11-26得票数 5
回答已采纳
1回答
使用CUDA C计算下列表达式的最有效方法是什么？
(A (D^-1)B^T )^-1
其中，D是一个很大的对称矩阵，A是一个小对称矩阵，使得B和B^T中等大小的矩形非对称矩阵。当然，(^-1)和(^T)分别是逆运算和转置运算。
浏览 4提问于2013-08-17得票数 0
我已经知道了为了执行旋转而必须使用的矩阵。如果我必须在z轴上旋转，然后在x轴上旋转，我会用2步完成。我的问题是，是否有可能将两个旋转组合成一个矩阵？我会感谢你的反馈。
浏览 0提问于2015-10-09得票数 9
3回答
是的，它已经做过无数次了，但该死的，我还想再做一次。我正在为C++编写一个简单的矩阵库，目的是实现right。我遇到了一些在数学上很明显的东西，但是对于强类型的系统来说就不那么明显了-- 1x1矩阵只是一个数字。为了避免这种情况，我开始沿着矩阵的多毛路径作为向量的组合，但也偶然发现两个向量相乘可以是一个数，也可以是一个dyad，这取决于两个向量的方向。
我的问题是，在强类型语言(如或)中处理这种情况的方法是什么？
浏览 6提问于2010-06-05得票数 2
从跟踪系统，我可以访问与相机(NeckWRTCameraMatrix)相关的颈4x4变换矩阵。
我需要在three.js动画的颈部骨，有作为父母的其他骨头，但以一种方式，颈部旋转与相机有关。
var m = new THREE.Matrix4();
m.fromArray(NeckWRTCameraMatrix);
var quaternion = new THREE.Quaternion();
quaternion.setFromRotationMatrix(m);
head.boneByName("Joint_Neck").quaternion.copy(quatern
浏览 1提问于2016-06-12得票数 0
1回答
问题1-R中的矩阵幂没有获取矩阵幂的内置命令.编写一个函数matrixpower()，其中包含两个参数mat和k，它们将取矩阵mat的整数幂k。
我的尝试解决方案.
matrixMul <- function(mat1, mat2)
{
rows <- nrow(mat1)
cols <- ncol(mat2)
if(rows == cols)
{
matOut <- matrix(nrow = rows, ncol = cols) # empty matrix
for (i in 1:rows)
浏览 0提问于2019-04-24得票数 5
3回答
将矩阵应用于函数
、、
试图将矩阵应用于函数，但使用mapply却没有成功。
我试图为不同的参数求解一组方程。在一种更简单的函数集合形式中，我试图将一个函数传递给一个矩阵常量-
a b c
[1,] 1 4 7
[2,] 2 5 8
[3,] 3 6 9
并试图求解方程3*a + 2*b + 3*c，并返回矩阵中每一行的答案。我已经将原来的函数改成了一个线性的、更简单的函数--这就是为什么我更喜欢使用#mapply，而以前的解释并没有对我有所帮助。
建立矩阵
my_vector <- 1:9
constants <- matrix(my_vector, 3, 3)
colnames(co
浏览 1提问于2019-09-15得票数 7
4回答
我试图循环一个矩阵，但找不到一个简单和优雅的方式，而不是写许多(>10)方程.有人能帮我吗？
我的母体是这样的：
我想计算以下几点：
(0 *0* 4/24) + (0 *1* 6/24) + (0 *2* 3/24) + (1 *0* 3/24) + (1 *1* 4/24) + (1 *2* 4/24)
而不是使用
__顺便说:我的矩阵代码
vals<- c(4/24, 6/24, 3/24, 3/24, 4/24, 4/24)
x <- c(0,1)
y <- c(0,1,2)
df <- matrix(vals, byrow = TRUE, nr
浏览 6提问于2022-07-10得票数 1
回答已采纳
1回答
如何在三维数组中运行迭代互相关？
、、、、
我正在Matlab中编写代码，但我想听听关于如何在其他语言中实现这一点的建议。
例如，我有一个3D ( M，n，t)矩阵M。我想运行一个迭代for循环，它将遍历每一个i (1 : m)和每一个j (1 : n)，然后计算循环中当前(i，j)对的沿t的互相关，与矩阵中所有可能的(i，j)组合。因此，每一个(i，j)都有一个n*m关联的关联。
下面我已经把它写成了一种伪代码，但是我想不出如何正确地实现它。不幸的是，在Matlab中，“corrcoef”函数不允许在函数本身中迭代1:m，1:n。
任何帮助都会很好！
M = randn(m,n,t);
Xcor_vals = zeros(m,n,m*
浏览 1提问于2018-04-30得票数 0
回答已采纳
2回答
查找最优旋转
、、
我有一个应用程序，我必须找到从一组15个排序器和索引的3D点(X1，X2，...，X15)到另一组具有相同索引的15个点(1个初始点对应1个终结点)的旋转。
我读过很多关于使用Euler角(对某些人来说是邪恶的)、四元数或将向量投影到基轴上的旋转。但是我有一个额外的约束:我的最终集合中的一些点可能是错误的(即有错误的坐标)，所以我想区分要求旋转的点离中值旋转非常远。
我的问题是:对于每组3个点(不是对齐的)和它们的图像，我可以计算四元数(根据转换矩阵不是纯旋转的事实，我有一些额外的计算，但这是可以做到的)。所以我得到了一组四元数(最多455个)，我想去掉那些错误的。
有没有办法找出哪些点的旋转
浏览 0提问于2012-07-03得票数 4
回答已采纳
1回答
基本上，我需要查询用户的矩阵。我需要找到A^2，A^3，A^4.以此类推(其中A是矩阵)。我需要找出A+ A^2 + A^3 .等至6。
到目前为止，这就是我所做的
public class Driver {
public static void main(String[] args) {
int i, j, l, k, sum = 0 ;
int matrixAColumnSize ;
int matrixARowSize ;
double numberOfNodes
浏览 3提问于2014-04-03得票数 0
回答已采纳
1回答
我无法理解这些问题的正确答案。有人能帮我说：“我们的矩阵向量乘法公式假定矩阵M是平方的，把算法推广到M是一些行r和列c的r-by-c矩阵的情况下。”
浏览 5提问于2017-10-15得票数 0
4回答
位姿绕原点的三维旋转
、、、
我想改变一个关于起源的姿势。我的姿势包含一个位置x，y，z和一个以四元数表示的方位x，y，z，w。我有一个4x4变换矩阵，我用它来平移和旋转姿态在原点周围的位置。
np.dot(transforma_mat, np.transpose(np.append(pose.position, 1)))
因此，当原点移动和旋转时，姿态的起点也会移动和旋转。然而，我不知道我如何能够旋转的姿态的方向。
我可以使用将四元数转换为欧拉角和旋转矩阵(反之亦然)。
故事:有一个物体可以移动和旋转。机器人末端执行器在物体移动和旋转之前有一个姿势.我想更新末端执行器的姿态，使其与物体一起移动和旋转。
浏览 4提问于2015-05-13得票数 1
回答已采纳
此语句检查数字是否为32位。
0 < res <= (1 << 31) -1
我似乎不明白为什么，有人能帮助理解这个位移位语法吗？
浏览 0提问于2017-11-26得票数 0
回答已采纳
1回答
左矩阵向量除法
、、、、
如何显式地在c1中求解最右边的方程，特别是((x-1)\y)
我很清楚当您使用一个矩阵时会发生什么，例如A\b (其中A是矩阵，b是列向量)，但是如果在两个相等行的向量上使用反斜杠会发生什么呢？
问题是：
x = (1:3)';
y = ones(3,1);
c1 = ((x-1)\y) % why does this =0.6?
浏览 0提问于2019-01-10得票数 3
回答已采纳
我得到了一个正定矩阵。然后，我修改了矩阵，使该矩阵的所有值都在0.8到1或-0.8到-1之间。但在此之后，矩阵不再是正定矩阵。我用过
[R,p]=chol(NewMat)
函数来验证这一点。如何使修正矩阵正定？
clc
clf
clear all
close all
Matrix=csvread('new1.csv'); %reading the csv file.
fid1 = fopen('new2.csv');
X = textscan(fid1, '%s%s%s%s%f%f%f', 'Delimiter', ',&
浏览 7提问于2014-05-01得票数 0
回答已采纳
1回答
我在3D游戏编程中更深入地研究了四元数的使用(是的，我知道矩阵，它们很好，但我们总是需要学习一些新的东西)，所以，我们可以旋转像这个Pout = q*Pin*conjug(q)这样的对象，其中q是四元数，Pin是一个对象(假设我们使用的是某个框架，这个类是为我们定义的) Vector3类，conjug(q)是它共轭后的Q四元数，最后，Pout是新的向量3‘s对象，在旋转初始向量3’s对象之后，我们得到了一个角阿尔法(或θ，)。你喜欢什么都行)。另外，我知道有一种组合旋转的方法，像这样：q_final = q2*q1 (这代表由alfa1旋转，然后由alfa2角度旋转)。最后，点积表示两个四元数
浏览 6提问于2012-10-30得票数 5
回答已采纳
地狱，我是C++ OpenCV编程的新手。我有一个彩色图像，我希望这个图像乘以几个矩阵(掩蔽)。这是代码的大纲：
#include <opencv2/core/core.hpp>
#include <opencv2/highgui/highgui.hpp>
#include <opencv2/imgproc/imgproc.hpp>
#include <iostream>
using namespace std;
using namespace cv;
char* hai;
int main()
{
hai = "c:/Dalton/g
浏览 2提问于2014-12-22得票数 0
回答已采纳
JAMA: Matrix.times()和Matrix.arrayTimes()在JAMA (一个用于矩阵计算的java库)之间有什么区别？
如果我有一个d维向量x和一个k维向量z，并且我想得到xz^T (x到z转置)，我应该使用Matrix.times还是Matrix.arrayTimes？
如何用JAMA计算这个乘法？
浏览 1提问于2016-05-07得票数 1
回答已采纳
1回答
图上的困惑
、、、、
给定一个无向图G=(V，E)，每个节点i与对象的'Ci‘数相关联。在每一步中，对于每个节点i，Ci对象在I的邻居之间被平均分配。经过K步，输出对象最多的前五个节点的对象数。
下面是一个步骤中发生的事情的一个例子：
A的对象被B和C平分。
B的对象被A和C平分。
C的对象被A和B平分。
一些约束：|V|<10^5，|E|<2*10^5，K<10^7，Ci<1000
我目前的想法是:用矩阵表示每个步骤中的转换。将该问题转化为矩阵幂的计算。但是这个解太慢了，考虑到x-V可以是10^5。
有什么更快的方法吗？
浏览 3提问于2012-10-14得票数 1
请给我。参见模式识别第4.1.3节-主教：“分类的最小二乘”：
在一个2类线性判别系统中，当y( \mathbf{x} )>0时，我们将向量\bf{x}分类为\mathcal{C}_1，而\mathcal{C}_2则将其分类为D3。在4.1.3节的推广中，我们定义了\mathcal{K}线性判别方程--每个类一个：
y_{k}(\bf{x}) = \bf{w_k^Tx} + \mathit{w_{k0}} \tag {4.13}
将领先的1添加到向量\bf{x}中会产生\tilde{\mathbf{x}}。
给出了\mathcal{K}类的线性判别函数：\bf y(x) = \widet
浏览 0提问于2019-08-05得票数 1
我试图用特征值中的JacobiSVD来求解一个线性方程。JacobiSVD有两个分解选项: ComputeThinU和ComputeFullU。它们的含义解释为
computed.ComputeThinU：in ComputeFullU:在JacobiSVD中，表示方阵U为 JacobiSVD，表示要计算薄矩阵U。
“薄矩阵”是什么意思？以及如何确定使用哪种选项？
浏览 1提问于2020-02-09得票数 1
2回答
运算A\B的结果是什么，其中A(1，m)和B (1，m)？
在手册中写道：
A\B returns a least-squares solution to the system of equations A*x= B.
所以它意味着x= inv (A'*A)*A'*B？然而，矩阵A'*A是奇异的.
让我们假设：
A=[1 2 3]
B=[6 7 6]
A\B
0 0 0
0 0 0
2.0000 2.3333 2.0000
如果我们使用MLS：
C = inv (A'*A) sin
浏览 0提问于2012-11-14得票数 9
回答已采纳
1回答
马尔可夫链自由度
、、、
我有一组长度为4的5000字符串，其中字符串中的每个字符都可以是A、B、C或D。
0阶马尔可夫链(不依赖)，使A，B，C，D列的4×1数组。
1阶马尔可夫链(pos j依赖于以前的pos i)，构造了一个4×4列的A、Bi、Ci、Di，以及Aj、Bj、Cj、Dj的列。
2阶马尔可夫链(pos k依赖于pos j和pos i)，构造了一个4×4×4维矩阵: Ai，Bi，Ci，Di；Aj，Bj，Cj，Dj；Ak，Bk，Ck，Dk，或由此形成一个16×4维矩阵Aij，Bij，Cij，Dij；Ak，Bk，Ck，Dk。
3阶马尔可夫链(pos l依赖于pos k，pos j和pos i)
浏览 1提问于2013-04-20得票数 1
3回答
我们可以用几种方法进行矩阵乘法。A和B是两个矩阵，大小为1000 * 1000。
使用matlab内建矩阵算子。A*B
显式地写两个巢。
显式地写三个巢。
动态链接C代码来完成这三个嵌套。
结果是
*算子是非常快速的。大约需要1秒。
巢非常慢，特别是三个巢，500秒。
C循环比matlab嵌套快得多，10秒。
有谁能解释一下。为什么*速度这么快，为什么C循环比在matlab中编写循环更有效。
我用另一个软件，而不是matlab，来做计算。但我认为类似的理由应该适用。
浏览 6提问于2014-05-25得票数 1
回答已采纳
我有一个大小为Nx2 (包含(x，y)坐标)的矩阵x和一个大小为Px1 (P<=N)的矩阵c，它包含我感兴趣的某些x坐标。例如：
x = [10 3; 21 9; 98 54; 4 30; 37 12];
c = [4 98];
我想要获取c中元素的坐标(在上面的例子中是[4 3])。我该怎么做呢？我只找到了当c是1x1矩阵(即标量)时的方法。
浏览 0提问于2012-05-13得票数 1
回答已采纳
1回答
利用变换矩阵绕任意点旋转
、、、、
我想用变换矩阵围绕任意点旋转点。
所以方程是：
p'x = cos(theta) * (px-ox) - sin(theta) * (py-oy) + ox
p'y = sin(theta) * (px-ox) + cos(theta) * (py-oy) + oy
但是我想以矩阵变换的形式重写它，它看起来应该是(c=cos(θ)，s=sin(θ))：
默认情况下，旋转矩阵使用原点作为旋转的中心。若要围绕任意点旋转，必须使用平移矩阵减去到原点的距离，进行旋转，然后再平移。
如果我正确地理解了它，我们首先从原来的点减去t (0,0)，然后围绕点(0,0)旋转，然后加到结
浏览 2提问于2013-12-23得票数 4
回答已采纳
2回答
我一直在试图计算一个36乘36矩阵的约旦范式，该矩阵只有三个不同的条目：1、1/2和0。矩阵是一个概率转移矩阵，因此，给定这些项，矩阵是明显稀疏的。
我一直面临的问题是:每当我试图计算
[V, J] = jordan(A),
或
[V, J] = jordan(sym(A)),
我收到以下错误消息：
使用mupadmex时出错
MuPAD命令中的错误:相似性矩阵太大。
sym/mupadmexnout中的错误(第1546行)
(B).= mupadmex(fcn，args{:})；
sym/jordan中的错误(第32行)
(B).=.=.
我在MATLAB中读到过，Jordan形式的
浏览 1提问于2014-02-14得票数 7
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
实时音视频活动推荐
腾讯云开发者
