请教各位大神,我想解这个方程,最后将α表示出来。
mathematica吧
全部回复
仅看楼主
level 3
雨默默地_ 楼主
Solve[A*Cos[2 \[Alpha] - \[CurlyPhi] - a]*Sin[\[Alpha]]*
Cos[\[CurlyPhi] - \[Alpha] + \[Omega] + \[Xi]] -
A*Sin[\[Alpha] - \[CurlyPhi]]*Cos[\[Alpha] - \[Omega]]*
Cos[2 \[Alpha] - \[CurlyPhi] - \[Omega] - \[Xi]] -
B*Sin[\[Alpha]]*Sin[\[Alpha]] + C*Sin[\[Alpha]]*Cos[\[Alpha]] -
F*Sin[\[Alpha] - \[CurlyPhi]]*
Cos[\[CurlyPhi] - \[Alpha] + \[Omega] + \[Xi]] == 0, \[Alpha]]
这是我的代码,可是运行后一直出不来结果,想问问大腿们这个是什么问题呢?
2017年12月08日 11点12分 1
level 10
语法没错,但是你的方程太复杂了。
Reduce[A*Cos[2 \[Alpha] - \[CurlyPhi] - a]*Sin[\[Alpha]]*
Cos[\[CurlyPhi] - \[Alpha] + \[Omega] + \[Xi]] -
A*Sin[\[Alpha] - \[CurlyPhi]]*Cos[\[Alpha] - \[Omega]]*
Cos[2 \[Alpha] - \[CurlyPhi] - \[Omega] - \[Xi]] -
B*Sin[\[Alpha]]*Sin[\[Alpha]] + C*Sin[\[Alpha]]*Cos[\[Alpha]] -
F*Sin[\[Alpha] - \[CurlyPhi]]*
Cos[\[CurlyPhi] - \[Alpha] + \[Omega] + \[Xi]] == 0]
部分解是可以得到你要的表达式的!
2017年12月08日 13点12分 5
我试了下,得到的是很长的一串,看起来好复杂。还有这里为什么用reduce呢?谢谢啦~老哥
2017年12月10日 12点12分
level 9
没解析解吧,试NSolve
2017年12月08日 22点12分 6
level 8
不想要那么长的结果再加一个 Quartics -> False
2017年12月09日 08点12分 7
这个数学上的处理很漂亮!
2017年12月09日 11点12分
大腿,我想问问最开始我的那个为什么会发生一直在运行呢?
2017年12月10日 12点12分
@雨默默地_ 不用管为什么,这是软件的局限性,其他软件也会有(可能在其他地方),只要知道怎么能算出来就行了
2017年12月11日 07点12分
@孺子剑牛不群 嗯,也对哈,感谢大腿[真棒]
2017年12月12日 03点12分
吧务
level 15
这类问题比较稳健的做法其实是用万能公式换元:
eq = A*Cos[2 \[Alpha] - \[CurlyPhi] - a]*Sin[\[Alpha]]*
Cos[\[CurlyPhi] - \[Alpha] + \[Omega] + \[Xi]] -
A*Sin[\[Alpha] - \[CurlyPhi]]*Cos[\[Alpha] - \[Omega]]*
Cos[2 \[Alpha] - \[CurlyPhi] - \[Omega] - \[Xi]] - B*Sin[\[Alpha]]*Sin[\[Alpha]] +
C*Sin[\[Alpha]]*Cos[\[Alpha]] -
F*Sin[\[Alpha] - \[CurlyPhi]]*Cos[\[CurlyPhi] - \[Alpha] + \[Omega] + \[Xi]] == 0;
neweq = TrigExpand@eq /.
Tan[\[Alpha]] :> Sin@\[Alpha]/Cos@\[Alpha] /. {Sin@\[Alpha] :> (2 tan)/(1 + tan^2),
Cos@\[Alpha] :> (1 - tan^2)/(1 + tan^2)} // Simplify
sol = Solve[neweq, tan]; // AbsoluteTiming
Solve[#, \[Alpha]] & /@ (Tan[\[Alpha]/2] == tan /. sol) // Simplify
不过这个方法解出来的解里面含了Root,并且好像很难继续化简。
2018年01月06日 12点01分 11
这个方法虽然稳,但是这样换元把方程升次了(4次变8次),难以写成根式的形式
2018年01月08日 04点01分
1