我正在尝试部分跟踪一个以分数作为MATLAB中元素的符号矩阵,但是Qetlab中的部分跟踪函数似乎存在一些问题。
我尝试查找源代码,但似乎没有找到解决方案。 Qetlab的(完整)跟踪功能(trace
)可以正常工作。以下是我要制作的简化示例:
N_sym=sym(1/32*eye(8)); % I define my symbolic matrix
sys=[1,2]; % Defines the subsystems that I want to trace out
dims=2*[1 1 1]; % Defines the dimensions of each subsystem (here I have defined 3 in total).
PartialTrace(N_sym,sys,dims)
我得到了错误:
Error using symengine Arguments must be 2-dimensional. Error in sym/privUnaryOp (line 1019) Csym = mupadmex(op,args{1}.s,varargin{:}); Error in sym/sum (line 55) s = privUnaryOp(A, 'symobj::prodsumdim', num2str(dim), '_plus'); Error in PartialTrace (line 93) Xpt = sum(Xpt(:,:,1:sub_sys_vec(1)+1:sub_sys_vec(1)^2),3);
我该如何做?