一道matlab解方程

来源:百度知道 编辑:UC知道 时间:2024/07/13 07:23:10
syms Ay Bx By Cx Cy Dx Dy Ex Ey Ax
eq1=sym('Ax=0')
eq2=sym('Ay=0')
eq3=sym('Bx=Ax+lab*cos(wt)')
eq4=sym('By=Ay+lab*sin(wt)')
eq5=sym('Cy=0')
eq6=sym('Cx=(lbc^2-(Cy-By)^2)^(1/2)+Bx')
eq7=sym('Ex=Dx-led*(Cy-By)/(Cy-By+Cx-Bx)')
eq8=sym('Ey=(By*Cx-Cy*Bx+Ex*(Cy-By))/(Cx-Bx)')
eq9=sym('Dx=m')
eq10=sym('Dy=Ey-(led^2-(Dx-Ex)^2)^(1/2)')
[Dy]=solve(eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9,eq10)

想把Dy用 lab lbc led m wt表示
这样能行吗?

我解出来是无解
而实际上我用vb编程的话,按照上面的关系,运动图都可以画出来。
上面大打错一个应该是这个
syms Ay Bx By Cx Cy Dx Dy Ex Ey Ax
eq1=sym('Ax=0')
eq2=sym('Ay=0')
eq3=sym('Bx=Ax+lab*cos(wt)')
eq4=sym('By=Ay+lab*sin(wt)')
eq5=sym('Cy=0')
eq6=sym('Cx=(lbc^2-(Cy-By)^2)^(1/2)+Bx')
eq7=sym('Ex=Dx-led*(Cy-By)/(Cy-By+Cx-Bx)')
eq8=sym('Ey=(By*Cx

syms Ay Bx By Cx Cy Dx Dy Ex Ey Ax
eq1=sym('Ax=0')
eq2=sym('Ay=0')
eq3=sym('Bx=Ax+lab*cos(wt)')
eq4=sym('By=Ay+lab*sin(wt)')
eq5=sym('Cy=0')
eq6=sym('Cx=(lbc^2-(Cy-By)^2)^(1/2)+Bx')
eq7=sym('Ex=Dx-led*(Cy-By)/(Cy-By+Cx-Bx)')
eq8=sym('Ey=(By*Cx-Cy*Bx+Ex*(Cy-By))/(Cx-Bx)')
eq9=sym('Dx=m')
eq10=sym('Dy=Ey-(led^2-(Dx-Ex)^2)^(1/2)')
[Dy]=solve(eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9,eq10)

想把Dy用 lab lbc led m wt表示
这样能行吗?

我解出来是无解
问题补充:而实际上我用vb编程的话,按照上面的关系,运动图都可以画出来。
上面大打错一个应该是这个
syms Ay Bx By Cx Cy Dx Dy Ex Ey Ax
eq1=sym('Ax=0')
eq2=sym('Ay=0')
eq3=sym('Bx=Ax+lab*cos(wt)')
eq4=sym('By=Ay+lab*sin(wt)')
eq5=sym('Cy=0')
eq6=sym('Cx=(lbc^2-(Cy-By)^2)^(1/2)+Bx')
eq7=sym('Ex=Dx-led*(Cy-By)/(Cy-By+Cx-Bx)')
eq