Matlab CVX环境下subs()符号替换失效触发类型报错排查

• subs()本身是Matlab符号工具箱函数，仅负责完成符号表达式内的变量替换操作，不会自动修改返回值类型：哪怕你把表达式里所有符号变量都替换成数值，返回结果依然是sym类型。打印时看到数值只是sym类型对纯数值常量的显示效果，本质不是双精度浮点数，用whos/class()查询就能看到类型仍为sym，这个行为和CVX无关，是符号工具箱本身的默认设定。
• 报错本质是两类工具箱的类型系统不兼容：进入cvx_begin到cvx_end的建模块后，CVX会重载所有矩阵运算、赋值操作的规则，仅支持double类型普通数值、CVX自身定义的优化变量/表达式两类数据。当你把sym类型的rhot和CVX声明的优化变量X_ele做乘法、求迹运算时，Matlab会尝试在sym类和cvx类之间做隐式类型转换，但两个工具箱没有做适配接口，最终抛出类型转换错误。
• 你把double()转换写在CVX块内可以临时解决问题，但属于不规范写法：一旦后续代码调整时，在转换语句前就调用了sym变量和CVX变量做运算，会再次触发相同报错。

所有和符号计算相关的操作（求导、符号替换、常数矩阵构造）必须全部放在cvx_begin语句之前完成，在进入CVX建模块前就把所有参与优化建模的常数矩阵显式转为double类型，不要把任何sym类型变量带入CVX块。
针对你给出的代码，只需要调整符号替换部分的代码即可：

% 符号计算、替换、类型转换全放在CVX块外完成
dim = size(rhot,1);
for i = 1:p
    drhot(:,:,i) = diff(rhot,t(i));
end
% subs替换后立刻转double，从根源上消除sym类型
drhot = double(subs(drhot,t,value));
rhot = double(subs(rhot,t,value));
rhott = kron(eye(p),rhot);

% 再进入CVX建模流程
cvx_begin sdp
    variable L_ele(dim,dim,(1+p)*p/2) hermitian;
    variable X_ele(dim,dim,p) hermitian;
    cnt = 1;
    for i = 1:p
        for j = i:p
            if i == j
                L((i-1)*dim+1:i*dim,(j-1)*dim+1:j*dim) = L_ele(:,:,cnt);
                cnt = cnt + 1;
            else
                L((i-1)*dim+1:i*dim,(j-1)*dim+1:j*dim) = L_ele(:,:,cnt);
                L((j-1)*dim+1:j*dim,(i-1)*dim+1:i*dim) = L_ele(:,:,cnt);
                cnt = cnt + 1;
            end
        end
    end
    for i = 1:p
        X((i-1)*dim+1:i*dim,1:dim) = X_ele(:,:,i);
    end
    for i = 1:p
        X2(1:dim,(i-1)*dim+1:i*dim) = X_ele(:,:,i);
    end

    % 块内不需要再做任何double转换，直接使用预处理好的数值矩阵即可
    minimize(trace(rhott*L));
    subject to
        for i =1:p
            trace(rhot*X_ele(:,:,i)) == value(i);
        end
        for i = 1:p
            for j = 1:p
                if i == j
                    trace(drhot(:,:,i)*X_ele(:,:,j)) == 1;
                else
                    trace(drhot(:,:,i)*X_ele(:,:,j)) == 0;
                end
            end
        end
        [L X; X2 eye(dim)] >= 0;
cvx_end
res = cvx_optval;

• 不要靠命令行打印的输出判断变量类型：存纯数值的sym类型打印效果和double类型完全一致，判断类型必须用class(变量名)或者whos查询。
• CVX建模块内禁止执行任何符号工具箱相关的操作（包括sym变量运算、subs调用等），两类工具箱完全不兼容。
• 所有优化问题里的常数项、系数矩阵，都应该在进入CVX块前完成计算和类型校验，避免建模过程中出现非预期的类型冲突。

内容的提问来源于stack exchange，提问作者narip