求关于q2不定积分的Matlab运算代码

Question

已知m12和m11想要对q2积分的完整matlab代码以及求积分运算时那个q2的倒数是怎么处理的，感谢

Answer 1

你可以试一试用微分方程ode来做

a1 = 1; a2 = 2; a3 = 3; a4 = 4; a5 = 5; a6 = 6;
M11 = @(q2, q3) a1 + a2 + a4 + 2*a3*cos(q2) + 2*a6*cos(q3) + 2*a5*cos(q2+q3);
M12 = @(q2, q3) a2 + a4 + a3*cos(q2) + 2*a6*cos(q3) + a5*cos(q2+q3);
ode = @(t, y) -M12(y(1), y(2))/M11(y(1), y(2));
q2_initial = 0; 
q3_initial = 0; 
y_initial = [q2_initial, q3_initial];
[t, y] = ode45(ode, [t_start, t_end], y_initial);
q2_values = y(:, 1);
q3_values = y(:, 2);

integral_value = trapz(t, -M12(q2_values, q3_values)./M11(q2_values, q3_values));