MATLAB编译独立应用时vpasolve等符号函数受限问题求解

问题背景 #

以下MATLAB代码在IDE中运行正常，但编译为独立应用时触发报错：

"syms, subs, vpasolve" are excluded from packaging for the MATLAB Runtime environment according to the MATLAB Compiler license. Either remove the file or function from your code, or use the MATLAB function "isdeployed" to ensure the function is not invoked in the deployed component.

原代码依赖符号工具箱函数，无法打包到MATLAB Runtime：

function [Y1] = roots_2021a(q)
syms r
Y=airy(1,-r)-q*airy(-r);
x=0:0.1:20;
y=double(subs(Y,r,x));
d=diff(y>0);
k=num2cell(find(d~=0)+1);
X=double(cellfun(@(z)vpasolve(Y==0,r,x(z)),k));
Y1 = unique(X*exp(1j*pi/3));
end

用户尝试用matlabFunction转换表达式，但仍需替换vpasolve的逻辑以完成编译。

解决方案：纯数值方法替代符号计算 #

核心思路是完全避开符号工具箱，直接使用MATLAB基础数值函数实现根求解：

修改后的代码 #

function [Y1] = roots_2021a(q)
    % 构造数值计算函数，直接调用airy的数值版本
    f = @(r) airy(1, -r) - q * airy(-r);
    
    % 生成采样点，检测函数符号变化的区间
    x = 0:0.1:20;
    y = f(x);
    
    % 定位符号变化的位置（对应根所在的区间）
    sign_change_indices = find(diff(sign(y)) ~= 0);
    
    % 遍历每个区间，用fzero求精确根
    X = zeros(size(sign_change_indices));
    for i = 1:length(sign_change_indices)
        interval_left = x(sign_change_indices(i));
        interval_right = x(sign_change_indices(i) + 1);
        % 在区间内查找根，fzero支持MATLAB Runtime打包
        X(i) = fzero(f, [interval_left, interval_right]);
    end
    
    % 完成转换并去重
    Y1 = unique(X * exp(1j*pi/3));
end

关键说明 #

1. 替换符号函数：直接用匿名函数f调用airy的数值实现，无需syms和subs，完全脱离符号工具箱
2. 根区间定位：通过diff(sign(y))检测函数符号变化的位置，确保每个区间内存在且仅有一个根（符合原方程根的分布特性）
3. 数值求根：用MATLAB基础函数fzero替代vpasolve，fzero支持打包到MATLAB Runtime，能在指定区间内高效求解精确根
4. 无需matlabFunction：原尝试的matlabFunction转换步骤多余，因为airy本身就是可直接调用的数值函数

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