Question

matlab求解含积分的方程组

这个可以用solve函数求解吗？？？新手，求大佬帮忙解答，最好能提供代码，感谢

Answer 1

用matlab求解含积分的方程组，可以这样实现：

1、先用integral函数，求解定积分，即

integral(@(t)sqrt(a^2.*sin(t).^2+b^2.*cos(t).^2),0,pi/2)

2、再用vpasolve函数，求解β、θ值，即

[beta1,theta1]=vpasolve(eq1,eq2, [beta1,theta1])

这里：eq1是指方程一，eq2是指方程二

3、详细的运行代码为

a=15;b=14;r=14.5;

syms theta1 beta1

eq1=integral(@(t)sqrt(a^2.*sin(t).^2+b^2.*cos(t).^2),0,pi/2)+beta1*pi*r/180;

eq2=(a*sin(theta1))/(b*cos(theta1));

[beta1,theta1]=vpasolve(eq1==100,eq2==tan(beta1), [beta1,theta1])

这里:beta1是指β，theta1是指θ

运行结果

追问

感谢，可是方程1下限是0，所以不能等于100的吧，大佬指导一下

感谢，可是方程1下限是0，所以不能等于100的吧，大佬指导一下

追答

没有错吧，按你写的公式写的。

eq1=integral(@(t)sqrt(a^2.*sin(t).^2+b^2.*cos(t).^2),0,pi/2)+beta1*pi*r/180;

这里的积分下限是0，上限是pi/2