matlab带参数的方程组求解
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/05 14:49:18
matlab带参数的方程组求解
matlab带参数的方程组求解
matlab带参数的方程组求解
>> syms x y v t1 t2;
>> L1=(((600-x)^2+y^2)^(1/2))-(x^2+y^2)^(1/2)-v*t1;
>> L2=(((600-x)^2+y^2)^(1/2))-((x^2+(450-y)^2)^(1/2))-v*t2;
>> [x,y]=solve(L1,L2)
x =
(t1^2*v^2)/1200 + (t1*v*(12*((t1*v - 600)*(t1*v + 600)*(t2*v - 750)*(t2*v + 750)*(t1*v - t2*v + 450)*(t2*v - t1*v + 450))^(1/2) - 16*t2^3*v^3 - 25*t1^3*v^3 + 3240000*t2*v + 16*t1*t2^2*v^3 + 16*t1^2*t2*v^3))/(1200*(25*t1^2*v^2 - 32*t1*t2*v^2 + 16*t2^2*v^2 - 3240000)) + 300
(t1^2*v^2)/1200 - (t1*v*(25*t1^3*v^3 + 16*t2^3*v^3 + 12*((t1*v - 600)*(t1*v + 600)*(t2*v - 750)*(t2*v + 750)*(t1*v - t2*v + 450)*(t2*v - t1*v + 450))^(1/2) - 3240000*t2*v - 16*t1*t2^2*v^3 - 16*t1^2*t2*v^3))/(1200*(25*t1^2*v^2 - 32*t1*t2*v^2 + 16*t2^2*v^2 - 3240000)) + 300
y =
(t1^2*v^2)/900 - (t2^2*v^2)/900 + (t1*v*(12*((t1*v - 600)*(t1*v + 600)*(t2*v - 750)*(t2*v + 750)*(t1*v - t2*v + 450)*(t2*v - t1*v + 450))^(1/2) - 16*t2^3*v^3 - 25*t1^3*v^3 + 3240000*t2*v + 16*t1*t2^2*v^3 + 16*t1^2*t2*v^3))/(900*(25*t1^2*v^2 - 32*t1*t2*v^2 + 16*t2^2*v^2 - 3240000)) - (t2*v*(12*((t1*v - 600)*(t1*v + 600)*(t2*v - 750)*(t2*v + 750)*(t1*v - t2*v + 450)*(t2*v - t1*v + 450))^(1/2) - 16*t2^3*v^3 - 25*t1^3*v^3 + 3240000*t2*v + 16*t1*t2^2*v^3 + 16*t1^2*t2*v^3))/(900*(25*t1^2*v^2 - 32*t1*t2*v^2 + 16*t2^2*v^2 - 3240000)) + 225
(t1^2*v^2)/900 - (t2^2*v^2)/900 - (t1*v*(25*t1^3*v^3 + 16*t2^3*v^3 + 12*((t1*v - 600)*(t1*v + 600)*(t2*v - 750)*(t2*v + 750)*(t1*v - t2*v + 450)*(t2*v - t1*v + 450))^(1/2) - 3240000*t2*v - 16*t1*t2^2*v^3 - 16*t1^2*t2*v^3))/(900*(25*t1^2*v^2 - 32*t1*t2*v^2 + 16*t2^2*v^2 - 3240000)) + (t2*v*(25*t1^3*v^3 + 16*t2^3*v^3 + 12*((t1*v - 600)*(t1*v + 600)*(t2*v - 750)*(t2*v + 750)*(t1*v - t2*v + 450)*(t2*v - t1*v + 450))^(1/2) - 3240000*t2*v - 16*t1*t2^2*v^3 - 16*t1^2*t2*v^3))/(900*(25*t1^2*v^2 - 32*t1*t2*v^2 + 16*t2^2*v^2 - 3240000)) + 225
运行结果好复杂.但是确实是这个.
这个使用solve函数就行了
MATLAB中的SOLVE函数不但能解直接求解,还能求解带参数的方程: clear all 方程组可化为矩阵形式,(没有办法加公式,下面为一个4x4矩阵乘以4x1矩阵=