Octave实战：求解函数极值、定积分误差分析与微分方程可视化

# Octave 解决方案 ## 1. 求函数 f(x)=x^(2/3).*(x^2-8) 的极值 ```octave % 定义函数 f = @(x) x.^(2/3).*(x.^2 - 8); % 绘制函数图像观察极值点 x = linspace(-3, 3, 1000); plot(x, f(x)); grid on; title('f(x) = x^{2/3}(x^2 - 8)'); xlabel('x'); ylabel('f(x)'); % 计算导数并寻找临界点 syms x_sym; f_sym = x_sym^(2/3)*(x_sym^2 - 8); df = diff(f_sym); critical_points = double(solve(df == 0)); % 计算二阶导数判断极值性质 d2f = diff(df); for i = 1:length(critical_points) x0 = critical_points(i); if isreal(x0) second_deriv = double(subs(d2f, x_sym, x0)); if second_deriv > 0 fprintf('在 x = %.4f 处有极小值: %.4f\n', x0, f(x0)); elseif second_deriv < 0 fprintf('在 x = %.4f 处有极大值: %.4f\n', x0, f(x0)); else fprintf('在 x = %.4f 处可能是拐点\n', x0); end end end ``` ## 2. 计算定积分 ∫₀^{π/4} 1/(1-sinx) dx 并与精确值 √2 比较 ```octave % 精确值 exact_value = sqrt(2); % 定义被积函数 f = @(x) 1./(1 - sin(x)); % 积分区间 a = 0; b = pi/4; % (1) 矩形法 n = 1000; h = (b - a)/n; x = a:h:b; rectangular = h * sum(f(x(1:end-1))); error_rect = abs(rectangular - exact_value); % (2) 梯形法 trapezoidal = trapz(x, f(x)); error_trap = abs(trapezoidal - exact_value); % (3) Simpson法 simpson = quad(f, a, b); error_simp = abs(simpson - exact_value); % (4) Gauss积分 gauss = quadgk(f, a, b); error_gauss = abs(gauss - exact_value); % 显示结果 fprintf('精确值: %.8f\n', exact_value); fprintf('矩形法结果: %.8f, 误差: %.8f\n', rectangular, error_rect); fprintf('梯形法结果: %.8f, 误差: %.8f\n', trapezoidal, error_trap); fprintf('Simpson法结果: %.8f, 误差: %.8f\n', simpson, error_simp); fprintf('Gauss积分结果: %.8f, 误差: %.8f\n', gauss, error_gauss); ``` ## 3. 求解微分方程 y'' - y' + y = 3cos(t) ```octave % 定义微分方程 ode = @(t, y) [y(2); y(2) - y(1) + 3*cos(t)]; % 初始条件和时间范围 tspan = [0, 2*pi]; y0 = [0; 0]; % 使用ode45求解 [t, y] = ode45(ode, tspan, y0); % 绘制结果 figure; plot(t, y(:,1), 'b-', 'LineWidth', 2); % y(t) 蓝色实线 hold on; plot(t, y(:,2), 'r--', 'LineWidth', 2); % y'(t) 红色虚线 grid on; legend('y(t)', "y'(t)"); title('微分方程 y'''' - y'' + y = 3cos(t) 的解'); xlabel('t'); ylabel('y 和 y'''); ``` 说明： 1. 对于第一个问题，我们通过符号计算找到导数并求解临界点，然后使用二阶导数判断极值性质。 2. 对于积分问题，我们比较了四种数值积分方法（矩形法、梯形法、Simpson法和Gauss积分）的结果与精确值√2的误差。 3. 对于微分方程，我们将其转化为一阶方程组并使用ode45求解，然后绘制y和y'的曲线，使用不同颜色和线型区分。 运行这些代码需要Octave环境，并可能需要安装符号计算包（对于第一个问题）。 ######[AI问答 | 1057点数解答 | 2025-04-01 19:42:28]