[已解决]python可以数学建模吗？

Simoncolosi

我想做的数学模型是把一组复杂周期函数的时域数据写出它的傅里叶级数的三角展开式，并绘制光谱。
比如：t=t1,t2,t3,t4.....
         a=a1,a2,a3,a4.....
求大佬指点

最佳答案

月排行榜 / 总排行榜

facevoid

2019-9-19 18:56:17

我找了下好像没有那种功能完善的包，只有别人写好的代码块，只要求级数的话应该够用了

https://blog.csdn.net/Night_MFC/article/details/84669194

跳转到最佳答案楼层

facevoid

在计算这一块，你可以用python的math数学包来自己写出傅里叶级数表达式，或者你可以用别人写好的傅里叶级数包。

画图的话可以用matplotlib包实现

Simoncolosi

facevoid 发表于 2019-9-19 16:17
在计算这一块，你可以用python的math数学包来自己写出傅里叶级数表达式，或者你可以用别人写好的傅里叶级数 ...

傅里叶级数包？有现成的呀？

facevoid

我找了下好像没有那种功能完善的包，只有别人写好的代码块，只要求级数的话应该够用了

https://blog.csdn.net/Night_MFC/article/details/84669194

Simoncolosi

facevoid 发表于 2019-9-19 18:56
我找了下好像没有那种功能完善的包，只有别人写好的代码块，只要求级数的话应该够用了

https://blog.csd ...

好像可以用给~谢谢啦！

Simoncolosi

facevoid 发表于 2019-9-19 18:56
我找了下好像没有那种功能完善的包，只有别人写好的代码块，只要求级数的话应该够用了

https://blog.csd ...

from symfit import parameters, variables, sin, cos, Fit
import numpy as np
import matplotlib.pyplot as plt

def fourier_series(x, f, n=0):
    a0, *cos_a = parameters(','.join(['a{}'.format(i) for i in range(0, n + 1)]))
    sin_b = parameters(','.join(['b{}'.format(i) for i in range(1, n + 1)]))
    # Construct the series
    series = a0 + sum(ai * cos(i * f * x) + bi * sin(i * f * x)
                     for i, (ai, bi) in enumerate(zip(cos_a, sin_b), start=1))
    return series

x, y = variables('x, y')
w, = parameters('w')
model_dict = {y: fourier_series(x, f=w, n=3)}
print(model_dict)
xdata = np.linspace(-np.pi, np.pi)
ydata = np.zeros_like(xdata)
ydata[xdata > 0] = 1
fit = Fit(model_dict, x=xdata, y=ydata)
fit_result = fit.execute()
print(fit_result)
plt.plot(xdata, ydata)
plt.plot(xdata, fit.model(x=xdata, **fit_result.params).y, color='green', ls=':')

大佬~可不可以给我讲解一下每行大概是什么意思啊。虽然我这样有点过分。。。

facevoid

Simoncolosi 发表于 2019-9-19 19:31
from symfit import parameters, variables, sin, cos, Fit
import numpy as np
import matplotlib.pyp ...

我似乎也不太理解。。。或者你可以先搜一下 parameters和variables的用途