用matlab求解一个积分方程组,求解这样一个方程组,其中M0=2400,u=1/6的重力加速度.fx fy z 是变量.fy z 的值. 谢谢.>> syms fx fy z t f m u;>> F1=fx/(m-(f*t)/2940);>> F2=u-(fy/(m-(f*t)/2940));>> d1=int(F1,t,0,z);>> d2=int(

来源:学生作业帮助网 编辑:作业帮 时间:2024/05/06 18:01:00

用matlab求解一个积分方程组,求解这样一个方程组,其中M0=2400,u=1/6的重力加速度.fx fy z 是变量.fy z 的值. 谢谢.>> syms fx fy z t f m u;>> F1=fx/(m-(f*t)/2940);>> F2=u-(fy/(m-(f*t)/2940));>> d1=int(F1,t,0,z);>> d2=int(
用matlab求解一个积分方程组,

求解这样一个方程组,其中M0=2400,u=1/6的重力加速度.fx fy  z  是变量.fy  z 的值. 谢谢.
>> syms fx fy z t f m u;
>> F1=fx/(m-(f*t)/2940);
>> F2=u-(fy/(m-(f*t)/2940));
>> d1=int(F1,t,0,z);
>> d2=int(F2,t,0,z);
>> d1
d1 =
-2940*fx*(log(m-1/2940*f*z)-log(m))/f
>> d2
d2 =
(u*z*f+2940*fy*log(m-1/2940*f*z)-2940*fy*log(m))/f
>> e1=sym('1692=d1');
>> e2=sym('57=d2');
>> e3=sym('fx^2+fy^2=7500^2');
>> e4=sym('f^2=7500^2');
接下来就不知道怎么解了

用matlab求解一个积分方程组,求解这样一个方程组,其中M0=2400,u=1/6的重力加速度.fx fy z 是变量.fy z 的值. 谢谢.>> syms fx fy z t f m u;>> F1=fx/(m-(f*t)/2940);>> F2=u-(fy/(m-(f*t)/2940));>> d1=int(F1,t,0,z);>> d2=int(
刚好回答了另外一个类似的问题(1830658703514901060,替换网页地址question后面的那串数字即可).
 
这个问题应该很难求出解析解,应该使用fsolve求数值解.
参考代码(题中所给条件u=1/6疑为u=g/6之误):
m = 2400;
Q = 7500/2490;
g = 9.8;
a = g/6;
% Tx = x1,Ty = x2,t = x3
f = @(x) [quadl(@(t)x(1)./(m-Q*t),0,x(3))-1692;
    quadl(@(t)x(2)./(m-Q*t)-a,0,x(3))-57;
    sqrt(x(1)^2+x(2)^2)-7500];
x = fsolve(f,[5304,5304,700]);
fprintf('Tx = %.2f\nTy = %.2f\nt = %.2f\n',x)
 
求出的结果:
Tx = 6873.97
Ty = 2999.75
t = 417.17

同学 你也是啊。。。