尝试获得具有不同 X-axes 的多图

Question

我需要这样的东西：

但我明白我知道我必须为情节的每个时间段提供一个单独的时间，但我不知道如何。我的情节也有问题我认为它情节不止一次而且它们相互重叠。

更新：

在@Zephyr 的帮助下现在看起来像这样：

但应该是这样的：

代码如下：

import mpmath as mp
import numpy as np
import sympy
beta = 0.25;
L=[0.04980905361573844, 0.0208207352451072, 0.012368753716465475, 0.009117292529338674, 0.007461219338976698, 0.006510364609693688, 0.005899506135250773, 0.005485130537183343, 0.0051898472561455, 0.004961157595209418, 0.004778617403698715, 0.00463084459959999, 0.004510113095117956, 0.004410195593051274, 0.004330450690278247]
Lc=[1.7509008762765992, 0.14986486457338544, 0.03453912303580302, 0.014622269851256788, 0.00831141421008418, 0.005660123321843252, 0.004287823173522503, 0.0034922189865254395, 0.0029879534061896186, 0.0026315863522143363, 0.002367747989524076, 0.0021671535838986545, 0.0020116727106455415, 0.0018885896058416002, 0.0017939246597491803]
for k0 in [0.052,0.12,0.252,0.464,0.792,1.264,1.928,2.824,4,5.600,7.795,10.806,14.928,20.599,28.000]:
    t=np.linspace(k0,30,50)
    for i in range(len(L)):
        j0 = L[i];
        j1 = Lc[i];
        G = []
        def f(s):
            return s**(beta - 1)/(j0*s**beta + j1*sympy.gamma(beta + 1))

        for j in range(len(t)):
            G.append(mp.invertlaplace(f, t[j], method = 'dehoog', dps = 10, degree = 50))

        plt.plot(t,G)

答案应该是这样的，但我得到了这个错误

import numpy as np
import sympy
beta = 0.25;
L=[0.04980905361573844, 0.0208207352451072, 0.012368753716465475, 0.009117292529338674, 0.007461219338976698, 0.006510364609693688, 0.005899506135250773, 0.005485130537183343, 0.0051898472561455, 0.004961157595209418, 0.004778617403698715, 0.00463084459959999, 0.004510113095117956, 0.004410195593051274, 0.004330450690278247]
Lc=[1.7509008762765992, 0.14986486457338544, 0.03453912303580302, 0.014622269851256788, 0.00831141421008418, 0.005660123321843252, 0.004287823173522503, 0.0034922189865254395, 0.0029879534061896186, 0.0026315863522143363, 0.002367747989524076, 0.0021671535838986545, 0.0020116727106455415, 0.0018885896058416002, 0.0017939246597491803]
for k0 in [0.052,0.12,0.252,0.464,0.792,1.264,1.928,2.824,4,5.600,7.795,10.806,14.928,20.599,28.000]:
    t[k0]=np.linspace(k0,30,50)
    for i in range(len(L)):
        j0 = L[i];
        j1 = Lc[i];
        G = []
        def f(s):
            return s**(beta - 1)/(j0*s**beta + j1*sympy.gamma(beta + 1))

        for j in range(len(t[k0])):
            G.append(mp.invertlaplace(f, t[k0][j], method = 'dehoog', dps = 10, degree = 50))

        plt.plot(t[k0],G)

虚线的作用是：

31.9279939766313*exp(-0.18*sympy.sqrt(7)/sympy.sqrt(t))

]

每一步都应该像这样完成。

Answer 1

你循环的次数太多了，这就是你绘制多条重叠线的原因。
您应该以这种方式重新构建您的代码：

import mpmath as mp
import numpy as np
import sympy
import matplotlib.pyplot as plt


beta = 0.25
L = [0.04980905361573844, 0.0208207352451072, 0.012368753716465475, 0.009117292529338674, 0.007461219338976698, 0.006510364609693688, 0.005899506135250773, 0.005485130537183343, 0.0051898472561455, 0.004961157595209418, 0.004778617403698715, 0.00463084459959999, 0.004510113095117956, 0.004410195593051274, 0.004330450690278247]
Lc = [1.7509008762765992, 0.14986486457338544, 0.03453912303580302, 0.014622269851256788, 0.00831141421008418, 0.005660123321843252, 0.004287823173522503, 0.0034922189865254395, 0.0029879534061896186, 0.0026315863522143363, 0.002367747989524076, 0.0021671535838986545, 0.0020116727106455415, 0.0018885896058416002, 0.0017939246597491803]
K = [0.052, 0.12, 0.252, 0.464, 0.792, 1.264, 1.928, 2.824, 4, 5.600, 7.795, 10.806, 14.928, 20.599, 28.000]


def f(s):
    return s**(beta - 1)/(j0*s**beta + j1*sympy.gamma(beta + 1))


for k0, j0, j1 in zip(K, L, Lc):
    t = np.linspace(k0, 30, 50)
    G = []
    for j in range(len(t)):
        G.append(mp.invertlaplace(f, t[j], method = 'dehoog', dps = 10, degree = 50))

    plt.plot(t, G)

plt.show()

感谢 zip，您可以遍历 K、L 和 Lc 列表，同时从每个列表中选取一个元素；不需要 i 和 j 计数器。

尝试获得具有不同 X-axes 的多图

Try to get a multiple plot with different X-axes

python

math

plot

data-visualization

matplotlib