如何在MATLAB中实现微分算子与指定多项式的对应运算

核心思路 #

我们可以把运算拆解为「多项式存储-算子解析-逐阶微分计算-结果合并」四个步骤实现：

• 用MATLAB原生的多项式系数向量（按降幂排列）存储y(x)及各阶微分结果
• 解析用户输入的a(x)表达式，提取每一项的微分阶数和对应系数
• 调用polyder函数计算指定阶数的多项式微分
• 所有项加权求和后得到最终结果，可转符号表达式直观展示

完整实现代码 #

% 1. 预定义固定多项式 y(x) = x^3 + x^2 + x + 1，系数按降幂排列
y_coeff = [1 1 1 1];

% 2. 解析用户输入的a(x)，示例输入为 a(x) = D2 + D + 6，这里可替换为用户输入
a_expr = "D2 + D + 6";
term_list = split(a_expr, "+"); % 拆分所有项
result_coeff = zeros(1, length(y_coeff)); % 初始化结果系数向量

for i = 1:length(term_list)
    term = strip(term_list(i)); % 去除项两端的空白
    if startsWith(term, "D") % 判断是不是微分项
        % 提取微分阶数，D后面没有数字就是1阶，有数字就是对应阶数
        order_str = extractAfter(term, "D");
        if isempty(order_str)
            order = 1;
        else
            order = str2double(order_str);
        end
        % 微分项的系数默认是1，可扩展支持带系数的项比如3D2这类
        coeff = 1;
        % 计算对应阶数的微分
        der_coeff = y_coeff;
        for k = 1:order
            der_coeff = polyder(der_coeff);
            % 微分阶数超过多项式最高阶时结果为0，提前退出
            if isempty(der_coeff)
                der_coeff = 0;
                break;
            end
        end
    else % 常数项，直接乘以y(x)
        coeff = str2double(term);
        der_coeff = y_coeff;
    end
    % 对齐系数长度后加权累加
    len_diff = length(result_coeff) - length(der_coeff);
    if len_diff > 0
        der_coeff = [zeros(1, len_diff), der_coeff];
    end
    result_coeff = result_coeff + coeff * der_coeff;
end

% 3. 转符号表达式输出方便查看
syms x
result_poly = poly2sym(result_coeff, x);
disp("运算结果为：");
disp(result_poly);

示例运行结果 #

输入a(x) = D2 + D + 6时，运行代码输出结果为：
6*x^3 + 9*x^2 + 9*x + 9
和手动计算的结果一致：6*(x³+x²+x+1) + (3x²+2x+1) + (6x+2) = 6x³+6x²+6x+6 +3x²+2x+1 +6x+2 = 6x³+9x²+9x+9。

扩展说明 #

如果需要支持带系数的微分项（比如3D2 - 2D + 5这类包含负号、前置系数的输入），只需要调整项解析部分的逻辑，拆分系数和算子部分即可。

内容的提问来源于stack exchange，提问作者C-137