错误使用 symengine Unable to convert expression containing symbolic variables into double array. Apply 'subs' function first to substitute values for variables. 出错 sym/double (第 868 行) Xstr = mupadmex('symobj::double', S.s, 0); 出错 MAINONE (第 26 行) y0 = double([u_steadystate, v_steadystate]);

这个错误提示表明你试图将一个包含符号变量的表达式转换为双精度数组，但是符号变量的值还没有被替换成数值。在Matlab中，符号计算和数值计算是分开的，符号变量和数值变量是不同的类型。如果你想将符号变量转换为数值变量，需要使用subs函数将符号变量替换为数值，然后再进行转换。因此，可以通过下面的代码来解决错误：

% 将符号变量替换为数值变量
u_ss = double(u_steadystate);
v_ss = double(v_steadystate);

% 定义时间t和初值条件
tspan = [0 1440];
y0 = [u_ss, v_ss];

% 调用ode45函数求解ODE方程
[t, y] = ode45(odefun, tspan, y0);

% 绘制u和v随时间变化的图像
figure;
plot(t, y(:, 1), 'r-', t, y(:, 2), 'b--', 'LineWidth', 2);
xlabel('Time (s)');
ylabel('Value');
legend('u', 'v');
title('u and v vs. Time');

这样就可以顺利地将符号变量替换为数值变量，然后使用ode45函数求解ODE方程，并绘制出u和v随时间变化的图像。