SciPy中傅里叶级数的计算

Question

有绝缘磁头的金属棒传热的一般解决办法是：

​

我试着在给定的时间计算初始温差为100的典型情况下的梯度。

import numpy as np
import matplotlib.pyplot as plt

x = np.linspace(0,3,100)
L = 3
n_max = 20

def bn(n):
    n = int(n)
    if (n%2 != 0):
        return 400/(np.pi**2*n**2)
    else:
        return 0

def wn(n):
    global L
    wn = (np.pi*n)/L
    return wn

def fourierSeries(n_max,x,t):
    a0 = 100/2
    partialSums = a0
    for n in range(1,n_max):
            partialSums = partialSums + bn(n)*np.exp(-.00001*wn(n)**2*t)*np.cos(wn(n)*x)
    return partialSums

u = []
for i in x:
    u.append(fourierSeries(n_max,i,1))

plt.plot(x,u)

但结果不是预期的

​

代码可能有什么问题？

Answer 1

我相信你错过了棒的温度函数：

f(x) =T+ 100/L *x

使用这个，并计算积分将完成这项工作。

import numpy as np
import matplotlib.pyplot as plt

x = np.linspace(0,3,100)
L = 3
n_max = 20

def bn(n):
    b=200/n**2/np.pi**2*(np.cos(n*np.pi)-1)
    return b

def wn(n):
    wn = (np.pi*n)/L
    return wn

def fourierSeries(n_max,x,t):
    a0 = 100/2
    partialSums = a0
    for n in range(1,n_max):
        partialSums = partialSums + bn(n)*np.exp(-.0005*wn(n)**2*t)*np.cos(wn(n)*x)
    return partialSums

u = []
hour = 3600
for i in x:
    u.append(fourierSeries(n_max,i,2*hour))

plt.plot(x,u)

而图是：

​

​