weixin_57972756 2021-11-02 11:55 采纳率: 50%
浏览 1472
已结题

如何用python求解三角函数方程呢?

想求解这样一个三角函数方程(其他数据都是已知的,只用求角度就行),使用sympy模块得不出解

img

img

img

想问问大家,python有什么方法可以求解三角函数方程吗?
代码如下:

from sympy import*
import math
x=Symbol('x')
T=float(input('T为多少:'))
K=float(input('K为多少:'))
a=0.196
b=0.36
A=T*(math.sqrt(math.pi*a))/K
print(solve((-4+8*b)*sin(x)+2*sin(2*x)+A*a*((-3+8*b)*sin(x/2)-5*sin(5*x/2))),x)
  • 写回答

2条回答 默认 最新

  • 「已注销」 2021-11-02 16:00
    关注

    math中的sin()函数有错误,用sympy的sin()函数;
    另外,solve(fx,x)你的x的位置写到solve外面去了;
    另外,我改好了以后发现程序也是执行不出来,拆解成两部分运行后,发现后半部分的式子解不出来,查了下资料可以是要nsolve解决
    但是nsolve解出来的是近似值,精度应该也是很高的,但是,函数的x要换成0,代码在下面
    另外,你的k和t被我写死了,自己改下哈~

    img

    img

    import sympy,math
    x=sympy.Symbol('x')
    T=3
    K=5
    a=0.196
    b=0.36
    A=T*(math.sqrt(math.pi*a))/K
    fx=(-4+8*b)*sympy.sin(x)+2*sympy.sin(2*x)+A*a*((-3+8*b)*sympy.sin(x/2)-5*sympy.sin(5*x/2))
    print(sympy.nsolve(fx,0))
    
    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论
查看更多回答(1条)

报告相同问题?

问题事件

  • 系统已结题 11月10日
  • 已采纳回答 11月2日
  • 修改了问题 11月2日
  • 创建了问题 11月2日

悬赏问题

  • ¥15 如何在scanpy上做差异基因和通路富集?
  • ¥20 关于#硬件工程#的问题,请各位专家解答!
  • ¥15 关于#matlab#的问题:期望的系统闭环传递函数为G(s)=wn^2/s^2+2¢wn+wn^2阻尼系数¢=0.707,使系统具有较小的超调量
  • ¥15 FLUENT如何实现在堆积颗粒的上表面加载高斯热源
  • ¥30 截图中的mathematics程序转换成matlab
  • ¥15 动力学代码报错,维度不匹配
  • ¥15 Power query添加列问题
  • ¥50 Kubernetes&Fission&Eleasticsearch
  • ¥15 報錯:Person is not mapped,如何解決?
  • ¥15 c++头文件不能识别CDialog