• Python动态建模（4）

    上一篇：

    Python动态建模-（3）_Unconquerable&Llxy的博客-CSDN博客modefunc：函数（就是上面定义的write()._write_t、write()._write_t_x、write()._write_t_y、write()._write_t_z四选一）bool布尔值（True、False）------------默认为True，真。...https://blog.csdn.net/html_finder/article/details/126161824

    感谢以下文章和资料作为参考：

    用 Matplotlib 库生成动画图表 - UCloud云社区动画是一种展示现象的有趣方式。相对于静态图表，人类总是容易被动画和交互式图表所吸引。在描述多年来的股票价格、过去十年的气候变化、季节性和趋势等时间序列数据时，动画更有意义，因为我们可以看到特定的参...https://www.ucloud.cn/yun/43671.htmlvideo - FFMPEG No such filter: 'palettegen' - Stack Overflowhttps://stackoverflow.com/questions/34134854/ffmpeg-no-such-filter-palettegenhttps://www.lfd.uci.edu/~gohlke/pythonlibs/#pythonmagickhttps://www.lfd.uci.edu/~gohlke/pythonlibs/#pythonmagick【Python】图片处理之“PythonMagick”库简易安装笔记_蛙鳜鸡鹳狸猿的博客-CSDN博客_pythonmagick“ImageMagick”可以说是开源的代码/命令行版PS，它支持包括创建、编辑、构图再到格式转换的图片处理功能。更重要的是，它几乎完整支持大部分编程语言，其中Python的接口库就是“PythonMagick”了。具体参考http://www.imagemagick.org/script/index.php。以下整理了“PythonMagick”库在各操作环境的简易安装方法。①RPM式https://blog.csdn.net/sweeper_freedoman/article/details/52994690【Python】批量直接修改图片存储大小脚本_蛙鳜鸡鹳狸猿的博客-CSDN博客对图片的处理有的情况下是对存储大小而非纵横的“width&height”有要求，这种对图片文件磁盘存储大小修改的工作也往往是批量的。借助Python通过“PythonMagick”库（参考：http://blog.csdn.net/sweeper_freedoman/article/details/52994690）可以实现需求。脚本简单如下。# !/usr/bin/python#https://blog.csdn.net/sweeper_freedoman/article/details/53000520Python调用“ImageMagick”：图片格式转换、尺寸修改、属性重构及加水印_蛙鳜鸡鹳狸猿的博客-CSDN博客      “ImageMagick”的Python库是“PythonMagick”，如之前的博文（http://blog.csdn.net/sweeper_freedoman/article/details/53000520，http://blog.csdn.net/sweeper_freedoman/article/details/53000145），有借用“PythonMagick”写过...https://blog.csdn.net/sweeper_freedoman/article/details/69789307

    现在明白了：缺少animation.save函数，并且在调用ImageMagick和ffmpeg库（matplotlib依赖）时出错。

    让我们解决一下...

    1）关于上次报的Warning：

    使用animation.save函数保存动态图表。

    所以马上就能解决了

    问题非常多，而且难以解决。

    至此为止，我们的项目结束了破产了

    这是完整的代码，其实基本功能都已经实现，关键在导出为图表的部分。各位大神可以以此为基础改进完善。

    1. # coding=utf-8
    2.  
    3. __author__ = "Unconquerable&Llxy"
    4. __name__ = "3DGraphic_With_Formula"
    5.  
    6. import sympy
    7. from tkinter.messagebox import *
    8. from tkinter import *
    9. from math import *
    10.  
    11. import matplotlib.pyplot as plt
    12. from mpl_toolkits.mplot3d import Axes3D
    13. import matplotlib.animation as animation
    14.  
    15.  
    16. def no_ans():
    17.     Tk().geometry("9x9+990000+999999")
    18.     showerror("警告", "方程无解析解")
    19.  
    20.  
    21. x = []
    22. y = []
    23. z = []
    24. t = []
    25.  
    26.  
    27. def toLeg(fx):
    28.     v = str(fx).split("\n")
    29.     val = []
    30.     for i in v:
    31.         tv = i.split("=")
    32.         length = len(tv)
    33.         last = None
    34.         while length != 0:
    35.             if length == 2:
    36.                 val.append('='.join(tv))
    37.                 break
    38.             if length == 1:
    39.                 val.append('='.join([tv[0], str(last)]))
    40.                 break
    41.             val.append('='.join([tv[-2], tv[-1]]))
    42.             length = length - 1
    43.             last = tv[-1]
    44.             tv.pop()
    45.     return '\n'.join(val)
    46.  
    47.  
    48. def legalSet(fx):
    49.     v = str(fx).split("\n")
    50.     nw = []
    51.     for l in v:
    52.  
    53.         tmpstr = ""
    54.         tmpstr2 = ""
    55.         if l == '': continue
    56.         k = l.split("=")
    57.         tmpstr += "((" + k[0] + ")-(" + k[1] + "))"
    58.         tmpstr2 += "((" + k[1] + ")-(" + k[0] + "))"
    59.         tmpstrs = [tmpstr, tmpstr2]
    60.         nw.append(tmpstrs)
    61.     v = nw
    62.     del nw
    63.  
    64.     index = -1
    65.     for i in v:
    66.         index += 1
    67.         for j in range(2):
    68.             if ":" in i: v[index][j] = i[j].replace(":", "/")
    69.             if "^" in i: v[index][j] = i[j].replace("^", "**")
    70.             if "÷" in i: v[index][j] = i[j].replace("÷", "/")
    71.             if "×" in i: v[index][j] = i[j].replace("×", "*")
    72.             if "[" in i: v[index][j] = i[j].replace("[", "(")
    73.             if "]" in i: v[index][j] = i[j].replace("]", ")")
    74.             if "【" in i: v[index][j] = i[j].replace("【", "(")
    75.             if "】" in i: v[index][j] = i[j].replace("】", ")")
    76.             if "{" in i: v[index][j] = i[j].replace("{", "(")
    77.             if "}" in i: v[index][j] = i[j].replace("}", ")")
    78.             if "（" in i: v[index][j] = i[j].replace("（", "(")
    79.             if "）" in i: v[index][j] = i[j].replace("）", ")")
    80.     newlt = []
    81.     lt = []
    82.     for j in v:
    83.  
    84.         lt__ = []
    85.  
    86.         news = []
    87.         for i_ in j:
    88.             new = ""
    89.             pos = -1
    90.  
    91.             funcs = []
    92.             for i in i_:
    93.  
    94.                 pos += 1
    95.                 tmpos = pos + 0
    96.                 string = ""
    97.                 while i_[tmpos].isalpha():
    98.                     string += i_[tmpos]
    99.                     tmpos = tmpos + 1
    100.  
    101.                 if string in sympy.__all__:
    102.                     funcs += [i for i in range(pos, tmpos + 1)]
    103.  
    104.                 if ((i.isalpha() or i == "(") and (i_[pos - 1].isnumeric() or i_[pos - 1].isalpha())
    105.                         and pos != 0 and pos not in funcs):
    106.                     new += "*" + " " + i
    107.  
    108.                 else:
    109.  
    110.                     if (i.isalpha() and pos not in funcs):
    111.                         new += " " + i
    112.                     else:
    113.                         new += i
    114.                 if i.isalpha() and pos not in funcs:
    115.                     lt__.append(i)
    116.  
    117.             news.append(new)
    118.  
    119.         lt__ = list(set(lt__))
    120.         lt.append(lt__)
    121.         newlt.append(news)
    122.     return newlt, lt
    123.  
    124.  
    125. def sqrt(base, times=2):
    126.     if base < 0:
    127.         ans = (-base) ** (1 / times)
    128.         return eval(f"{ans}j")
    129.     return base ** (1 / times)
    130.  
    131.  
    132. class write:
    133.     def _write_t_x(self, formula, tmin_=-10, tmax_=10, tnums=100,
    134.                    min_=-10, max_=10, nums=100, printout=True):
    135.  
    136.         global x, y, z, t
    137.  
    138.         self.step = (tmax_ - tmin_) / tnums
    139.         self.aa = tmin_ - self.step
    140.         self.xstep = (max_ - min_) / nums
    141.         r"""
    143.         print(un)
    144.         print(lf,toLeg(formula),sep="\n\n",end="\n\n")
    145.         """
    146.         solves = []
    147.         for iVal in range(tnums + 1):
    148.  
    149.             self.aa = self.aa + self.step
    150.             self.aa = float("%.10f" % self.aa)
    151.             xValue = min_ - self.xstep
    152.             for XVAL in range(nums + 1):
    153.                 lf, un = legalSet(toLeg(formula))
    154.                 xValue += self.xstep
    155.                 xValue = float("%.10f" % xValue)
    156.  
    157.                 def key(item):
    158.                     if "t" in item[0] or "x" in item[0]:
    159.                         return 1
    160.                     else:
    161.                         return 999
    162.  
    163.                 bf = dict(enumerate(lf))
    164.                 lf = sorted(lf, key=key)
    165.                 nun = []
    166.                 for i in range(len(lf)):
    167.                     nun.insert(lf.index(bf[i]), un[i])
    168.                 un = nun
    169.                 del nun
    170.                 idx = -1
    171.                 for imp in un:
    172.                     idx += 1
    173.                     for __imp in imp:
    174.                         if __imp == 't':
    175.                             un[idx].remove("t")
    176.                         if __imp == 'x':
    177.                             un[idx].remove("x")
    178.                 del idx
    179.  
    180.                 repl = {"t": [str(self.aa)], "x": [str(xValue)]}
    181.                 is_exit = False
    182.                 times = -1
    183.                 while not is_exit:
    184.                     times += 1
    185.  
    186.                     if times >= 10:
    187.                         # print(f"\nmessage:目前已经循环求解超过10次，目前初始代入值为：x{repl['x']}，t{repl['t']}，\n"
    188.                         #      f"已经求解{repl}，剩{4-len(repl.values())}值不可求解，故不予求解并跳出循环")
    189.                         break
    190.                     indx = -1
    191.                     for lf_ in lf:
    192.                         indx += 1
    193.                         for i in repl.keys():
    194.                             for j in range(2):
    195.                                 lf_[j] = lf_[j].replace(" " + i, str(repl[i][0]))
    196.  
    197.                         idx = -1
    198.                         for imp in un:
    199.                             idx += 1
    200.                             for __imp in imp:
    201.                                 if __imp in repl.keys():
    202.                                     un[idx].remove(__imp)
    203.                         del idx
    204.  
    205.                         if un[indx] == []:
    206.                             del un[indx]
    207.                             del lf[indx]
    208.                             if indx == 0 or (indx - 1) >= len(lf): break
    209.                             indx -= 1
    210.                             continue
    211.  
    212.                         try:
    213.                             ams = sympy.solve(tuple(lf_), tuple(un[indx]))
    214.                             if ams == []:
    215.                                 newlf = []
    216.                                 for i in lf_:
    217.                                     newlf.append(i.replace(str(repl["x"]), str(repl["x"]) + "*sympy.I"))
    218.                                 ams = sympy.solve(tuple(newlf), tuple(un[indx]))
    219.                         except NotImplementedError as err:
    220.                             print(err)
    221.                             no_ans()
    222.                             return
    223.  
    224.                         if isinstance(ams, (dict,)):
    225.                             for km, jj in ams.items():
    226.                                 if not ('x' in str(jj) or
    227.                                         'y' in str(jj) or
    228.                                         'z' in str(jj) or
    229.                                         't' in str(jj)):
    230.                                     if km.name in repl.keys():
    231.                                         repl[km.name].append(str(jj))
    232.                                     else:
    233.                                         repl[km.name] = [str(jj)]
    234.                         else:
    235.                             ix = -1
    236.                             undo = 0
    237.                             for i in ams:
    238.                                 ix += 1
    239.                                 try:
    240.                                     if isinstance(i, (sympy.core.numbers.Float,
    241.                                                       sympy.core.numbers.Integer,
    242.                                                       int, float, complex)):
    243.  
    244.                                         try:
    245.                                             repl[tuple(un[indx])[ix]] = [str(eval(str(i)))]
    246.                                         except IndexError:
    247.  
    248.                                             repl[tuple(un[indx])[undo]].append(str(eval(str(i))))
    249.                                         for i in repl.keys():
    250.                                             if i in lf_:
    251.                                                 lf_ = lf_.replace(i, str(repl[i][0]))
    252.                                     elif isinstance(i, (tuple, list, set)):
    253.                                         inxx = -1
    254.                                         for jj in i:
    255.                                             inxx += 1
    256.                                             if not ('x' in str(jj) or
    257.                                                     'y' in str(jj) or
    258.                                                     'z' in str(jj) or
    259.                                                     't' in str(jj)):
    260.                                                 try:
    261.                                                     repl[tuple(un[indx][inxx])[ix]] = [str(jj)]
    262.                                                 except IndexError:
    263.                                                     repl[tuple(un[indx][inxx])[undo]].append(str(jj))
    264.                                     else:
    265.                                         continue
    266.                                 except TypeError as err:
    267.                                     print("get:err", err)
    268.                                     continue
    269.                                 else:
    270.                                     undo = ix
    271.                         if len(repl.values()) >= 4: is_exit = True;break
    272.                 if is_exit:
    273.                     solves.append(repl)
    274.                 if printout:
    275.                     print(f"\r{(iVal * (nums + 1) + XVAL + 1) / ((nums + 1) * (tnums + 1)) * 100}%", end="",
    276.                           flush=True)
    277.  
    278.         if printout:
    279.             print("\nsucceed in Calculating Values...")
    280.         num = 0
    281.         for i in solves:
    282.             if printout:
    283.                 print(f"\r{(num + 1) / len(solves) * 100}%", end="", flush=True)
    284.             num += 1
    285.             x_ = []
    286.             for j in i["x"]:
    287.                 x_.append(eval(j))
    288.             x.append(x_)
    289.             y_ = []
    290.             for j in i["y"]:
    291.                 y_.append(eval(j))
    292.             y.append(y_)
    293.             z_ = []
    294.             for j in i["z"]:
    295.                 z_.append(eval(j))
    296.             z.append(z_)
    297.             t_ = []
    298.             for j in i["t"]:
    299.                 t_.append(eval(j))
    300.             t.append(t_)
    301.  
    302.         del x_, y_, z_, t_, solves
    303.  
    304.         if printout:
    305.             print("\nsucceed in Saving Values...")
    306.  
    307.     def _write_t_y(self, formula, tmin_=-10, tmax_=10, tnums=100,
    308.                    min_=-10, max_=10, nums=100, printout=True):
    309.  
    310.         global x, y, z, t
    311.  
    312.         self.step = (tmax_ - tmin_) / tnums
    313.         self.aa = tmin_ - self.step
    314.         self.xstep = (max_ - min_) / nums
    315.         r"""
    317.         print(un)
    318.         print(lf,toLeg(formula),sep="\n\n",end="\n\n")
    319.         """
    320.         solves = []
    321.         for iVal in range(tnums + 1):
    322.  
    323.             self.aa = self.aa + self.step
    324.             self.aa = float("%.10f" % self.aa)
    325.             xValue = min_ - self.xstep
    326.             for XVAL in range(nums + 1):
    327.                 lf, un = legalSet(toLeg(formula))
    328.                 xValue += self.xstep
    329.                 xValue = float("%.10f" % xValue)
    330.  
    331.                 def key(item):
    332.                     if "t" in item[0] or "y" in item[0]:
    333.                         return 1
    334.                     else:
    335.                         return 999
    336.  
    337.                 bf = dict(enumerate(lf))
    338.                 lf = sorted(lf, key=key)
    339.                 nun = []
    340.                 for i in range(len(lf)):
    341.                     nun.insert(lf.index(bf[i]), un[i])
    342.                 un = nun
    343.                 del nun
    344.                 idx = -1
    345.                 for imp in un:
    346.                     idx += 1
    347.                     for __imp in imp:
    348.                         if __imp == 't':
    349.                             un[idx].remove("t")
    350.                         if __imp == 'y':
    351.                             un[idx].remove("x")
    352.                 del idx
    353.  
    354.                 repl = {"t": [str(self.aa)], "y": [str(xValue)]}
    355.                 is_exit = False
    356.                 times = -1
    357.                 while not is_exit:
    358.                     times += 1
    359.  
    360.                     if times >= 20:
    361.                         print(f"\nmessage:目前已经循环求解超过20次，目前初始代入值为：y{repl['y']}，t{repl['t']}，\n"
    362.                               f"已经求解{repl}，剩{4 - len(repl.values())}值不可求解，故不予求解并跳出循环")
    363.                         break
    364.                     indx = -1
    365.                     for lf_ in lf:
    366.                         indx += 1
    367.                         for i in repl.keys():
    368.                             for j in range(2):
    369.                                 lf_[j] = lf_[j].replace(" " + i, str(repl[i][0]))
    370.  
    371.                         idx = -1
    372.                         for imp in un:
    373.                             idx += 1
    374.                             for __imp in imp:
    375.                                 if __imp in repl.keys():
    376.                                     un[idx].remove(__imp)
    377.                         del idx
    378.  
    379.                         if un[indx] == []:
    380.                             del un[indx]
    381.                             del lf[indx]
    382.                             if indx == 0 or (indx - 1) >= len(lf): break
    383.                             indx -= 1
    384.                             continue
    385.  
    386.                         try:
    387.                             ams = sympy.solve(tuple(lf_), tuple(un[indx]))
    388.                             if ams == []:
    389.                                 newlf = []
    390.                                 for i in lf_:
    391.                                     newlf.append(i.replace(str(repl["y"]), str(repl["y"]) + "*sympy.I"))
    392.                                 ams = sympy.solve(tuple(newlf), tuple(un[indx]))
    393.                         except NotImplementedError as err:
    394.                             print(err)
    395.                             no_ans()
    396.                             return
    397.  
    398.                         if isinstance(ams, (dict,)):
    399.                             for km, jj in ams.items():
    400.                                 if not ('x' in str(jj) or
    401.                                         'y' in str(jj) or
    402.                                         'z' in str(jj) or
    403.                                         't' in str(jj)):
    404.                                     if km.name in repl.keys():
    405.                                         repl[km.name].append(str(jj))
    406.                                     else:
    407.                                         repl[km.name] = [str(jj)]
    408.                         else:
    409.                             ix = -1
    410.                             undo = 0
    411.                             for i in ams:
    412.                                 ix += 1
    413.                                 try:
    414.                                     if isinstance(i, (sympy.core.numbers.Float,
    415.                                                       sympy.core.numbers.Integer,
    416.                                                       int, float, complex)):
    417.  
    418.                                         try:
    419.                                             repl[tuple(un[indx])[ix]] = [str(eval(str(i)))]
    420.                                         except IndexError:
    421.  
    422.                                             repl[tuple(un[indx])[undo]].append(str(eval(str(i))))
    423.                                         for i in repl.keys():
    424.                                             if i in lf_:
    425.                                                 lf_ = lf_.replace(i, str(repl[i][0]))
    426.                                     elif isinstance(i, (tuple, list, set)):
    427.                                         inxx = -1
    428.                                         for jj in i:
    429.                                             inxx += 1
    430.                                             if not ('x' in str(jj) or
    431.                                                     'y' in str(jj) or
    432.                                                     'z' in str(jj) or
    433.                                                     't' in str(jj)):
    434.                                                 try:
    435.                                                     repl[tuple(un[indx][inxx])[ix]] = [str(jj)]
    436.                                                 except IndexError:
    437.                                                     repl[tuple(un[indx][inxx])[undo]].append(str(jj))
    438.                                     else:
    439.                                         continue
    440.                                 except TypeError as err:
    441.                                     print("get:err", err)
    442.                                     continue
    443.                                 else:
    444.                                     undo = ix
    445.                         if len(repl.values()) >= 4: is_exit = True;break
    446.                 if is_exit:
    447.                     solves.append(repl)
    448.                 if printout:
    449.                     print(f"\r{(iVal * (nums + 1) + XVAL + 1) / ((nums + 1) * (tnums + 1)) * 100}%", end="",
    450.                           flush=True)
    451.  
    452.         if printout:
    453.             print("\nsucceed in Calculating Values...")
    454.         num = 0
    455.         for i in solves:
    456.             if printout:
    457.                 print(f"\r{(num + 1) / len(solves) * 100}%", end="", flush=True)
    458.             num += 1
    459.             x_ = []
    460.             for j in i["x"]:
    461.                 x_.append(eval(j))
    462.             x.append(x_)
    463.             y_ = []
    464.             for j in i["y"]:
    465.                 y_.append(eval(j))
    466.             y.append(y_)
    467.             z_ = []
    468.             for j in i["z"]:
    469.                 z_.append(eval(j))
    470.             z.append(z_)
    471.             t_ = []
    472.             for j in i["t"]:
    473.                 t_.append(eval(j))
    474.             t.append(t_)
    475.  
    476.         del x_, y_, z_, t_, solves
    477.  
    478.         if printout:
    479.             print("\nsucceed in Saving Values...")
    480.  
    481.     def _write_t_z(self, formula, tmin_=-10, tmax_=10, tnums=100,
    482.                    min_=-10, max_=10, nums=100, printout=True):
    483.  
    484.         global x, y, z, t
    485.  
    486.         self.step = (tmax_ - tmin_) / tnums
    487.         self.aa = tmin_ - self.step
    488.         self.xstep = (max_ - min_) / nums
    489.         r"""
    491.         print(un)
    492.         print(lf,toLeg(formula),sep="\n\n",end="\n\n")
    493.         """
    494.         solves = []
    495.         for iVal in range(tnums + 1):
    496.  
    497.             self.aa = self.aa + self.step
    498.             self.aa = float("%.10f" % self.aa)
    499.             xValue = min_ - self.xstep
    500.             for XVAL in range(nums + 1):
    501.                 lf, un = legalSet(toLeg(formula))
    502.                 xValue += self.xstep
    503.                 xValue = float("%.10f" % xValue)
    504.  
    505.                 def key(item):
    506.                     if "t" in item[0] or "z" in item[0]:
    507.                         return 1
    508.                     else:
    509.                         return 999
    510.  
    511.                 bf = dict(enumerate(lf))
    512.                 lf = sorted(lf, key=key)
    513.                 nun = []
    514.                 for i in range(len(lf)):
    515.                     nun.insert(lf.index(bf[i]), un[i])
    516.                 un = nun
    517.                 del nun
    518.                 idx = -1
    519.                 for imp in un:
    520.                     idx += 1
    521.                     for __imp in imp:
    522.                         if __imp == 't':
    523.                             un[idx].remove("t")
    524.                         if __imp == 'z':
    525.                             un[idx].remove("z")
    526.                 del idx
    527.  
    528.                 repl = {"t": [str(self.aa)], "z": [str(xValue)]}
    529.                 is_exit = False
    530.                 times = -1
    531.                 while not is_exit:
    532.                     times += 1
    533.  
    534.                     if times >= 20:
    535.                         print(f"\nmessage:目前已经循环求解超过20次，目前初始代入值为：z{repl['z']}，t{repl['t']}，\n"
    536.                               f"已经求解{repl}，剩{4 - len(repl.values())}值不可求解，故不予求解并跳出循环")
    537.                         break
    538.                     indx = -1
    539.                     for lf_ in lf:
    540.                         indx += 1
    541.                         for i in repl.keys():
    542.                             for j in range(2):
    543.                                 lf_[j] = lf_[j].replace(" " + i, str(repl[i][0]))
    544.  
    545.                         idx = -1
    546.                         for imp in un:
    547.                             idx += 1
    548.                             for __imp in imp:
    549.                                 if __imp in repl.keys():
    550.                                     un[idx].remove(__imp)
    551.                         del idx
    552.  
    553.                         if un[indx] == []:
    554.                             del un[indx]
    555.                             del lf[indx]
    556.                             if indx == 0 or (indx - 1) >= len(lf): break
    557.                             indx -= 1
    558.                             continue
    559.  
    560.                         try:
    561.                             ams = sympy.solve(tuple(lf_), tuple(un[indx]))
    562.                             if ams == []:
    563.                                 newlf = []
    564.                                 for i in lf_:
    565.                                     newlf.append(i.replace(str(repl["z"]), str(repl["z"]) + "*sympy.I"))
    566.                                 ams = sympy.solve(tuple(newlf), tuple(un[indx]))
    567.                         except NotImplementedError as err:
    568.                             print(err)
    569.                             no_ans()
    570.                             return
    571.  
    572.                         if isinstance(ams, (dict,)):
    573.                             for km, jj in ams.items():
    574.                                 if not ('x' in str(jj) or
    575.                                         'y' in str(jj) or
    576.                                         'z' in str(jj) or
    577.                                         't' in str(jj)):
    578.                                     if km.name in repl.keys():
    579.                                         repl[km.name].append(str(jj))
    580.                                     else:
    581.                                         repl[km.name] = [str(jj)]
    582.                         else:
    583.                             ix = -1
    584.                             undo = 0
    585.                             for i in ams:
    586.                                 ix += 1
    587.                                 try:
    588.                                     if isinstance(i, (sympy.core.numbers.Float,
    589.                                                       sympy.core.numbers.Integer,
    590.                                                       int, float, complex)):
    591.  
    592.                                         try:
    593.                                             repl[tuple(un[indx])[ix]] = [str(eval(str(i)))]
    594.                                         except IndexError:
    595.  
    596.                                             repl[tuple(un[indx])[undo]].append(str(eval(str(i))))
    597.                                         for i in repl.keys():
    598.                                             if i in lf_:
    599.                                                 lf_ = lf_.replace(i, str(repl[i][0]))
    600.                                     elif isinstance(i, (tuple, list, set)):
    601.                                         inxx = -1
    602.                                         for jj in i:
    603.                                             inxx += 1
    604.                                             if not ('x' in str(jj) or
    605.                                                     'y' in str(jj) or
    606.                                                     'z' in str(jj) or
    607.                                                     't' in str(jj)):
    608.                                                 try:
    609.                                                     repl[tuple(un[indx][inxx])[ix]] = [str(jj)]
    610.                                                 except IndexError:
    611.                                                     repl[tuple(un[indx][inxx])[undo]].append(str(jj))
    612.                                     else:
    613.                                         continue
    614.                                 except TypeError as err:
    615.                                     print("get:err", err)
    616.                                     continue
    617.                                 else:
    618.                                     undo = ix
    619.                         if len(repl.values()) >= 4: is_exit = True;break
    620.                 if is_exit:
    621.                     solves.append(repl)
    622.                 if printout:
    623.                     print(f"\r{(iVal * (nums + 1) + XVAL + 1) / ((nums + 1) * (tnums + 1)) * 100}%", end="",
    624.                           flush=True)
    625.  
    626.         if printout:
    627.             print("\nsucceed in Calculating Values...")
    628.         num = 0
    629.         for i in solves:
    630.             if printout:
    631.                 print(f"\r{(num + 1) / len(solves) * 100}%", end="", flush=True)
    632.             num += 1
    633.             x_ = []
    634.             for j in i["x"]:
    635.                 x_.append(eval(j))
    636.             x.append(x_)
    637.             y_ = []
    638.             for j in i["y"]:
    639.                 y_.append(eval(j))
    640.             y.append(y_)
    641.             z_ = []
    642.             for j in i["z"]:
    643.                 z_.append(eval(j))
    644.             z.append(z_)
    645.             t_ = []
    646.             for j in i["t"]:
    647.                 t_.append(eval(j))
    648.             t.append(t_)
    649.  
    650.         del x_, y_, z_, t_, solves
    651.  
    652.         if printout:
    653.             print("\nsucceed in Saving Values...")
    654.  
    655.     def _write_t(self, formula, tmin_=-10, tmax_=10, tnums=100, printout=True):
    656.  
    657.         global x, y, z, t
    658.         I = 1j
    659.         oo = inf
    660.  
    661.         self.step = (tmax_ - tmin_) / tnums
    662.         self.aa = tmin_ - self.step
    663.         r"""
    665.         print(un)
    666.         print(lf,toLeg(formula),sep="\n\n",end="\n\n")
    667.         """
    668.         solves = []
    669.         for iVal in range(tnums + 1):
    670.  
    671.             self.aa = self.aa + self.step
    672.             self.aa = float("%.10f" % self.aa)
    673.             lf, un = legalSet(toLeg(formula))
    674.  
    675.             def key(item):
    676.                 if "t" in item[0]:
    677.                     return 1
    678.                 else:
    679.                     return 999
    680.  
    681.             bf = dict(enumerate(lf))
    682.             lf = sorted(lf, key=key)
    683.             nun = []
    684.             for i in range(len(lf)):
    685.                 nun.insert(lf.index(bf[i]), un[i])
    686.             un = nun
    687.             del nun
    688.             idx = -1
    689.             for imp in un:
    690.                 idx += 1
    691.                 for __imp in imp:
    692.                     if __imp == 't':
    693.                         un[idx].remove("t")
    694.             del idx
    695.  
    696.             repl = {"t": [str(self.aa)]}
    697.             is_exit = False
    698.  
    699.             times = -1
    700.             while not is_exit:
    701.                 times += 1
    702.  
    703.                 if times >= 20:
    704.                     print(f"\nmessage:目前已经循环求解超过20次，目前初始代入值为：t{repl['t']}，\n"
    705.                           f"已经求解{repl}，剩{4 - len(repl.values())}值不可求解，故不予求解并跳出循环")
    706.                     break
    707.                 indx = -1
    708.                 for lf_ in lf:
    709.                     indx += 1
    710.                     for i in repl.keys():
    711.                         for j in range(2):
    712.                             if str(i) in str(lf_[j]):
    713.                                 lf_[j] = lf_[j].replace(i, str(repl[i][0]))
    714.  
    715.                     idx = -1
    716.                     for imp in un:
    717.                         idx += 1
    718.                         for __imp in imp:
    719.                             if __imp in repl.keys():
    720.                                 un[idx].remove(__imp)
    721.                     del idx
    722.  
    723.                     if un[indx] == []:
    724.                         del un[indx]
    725.                         del lf[indx]
    726.                         if indx == 0 or (indx - 1) >= len(lf): break
    727.                         indx -= 1
    728.                         continue
    729.  
    730.                     try:
    731.                         ams = sympy.solve(tuple(lf_), tuple(un[indx]))
    732.                     except NotImplementedError as err:
    733.                         print(err)
    734.                         no_ans()
    735.                         return
    736.  
    737.                     if isinstance(ams, (dict,)):
    738.                         for km, jj in ams.items():
    739.                             if not ('x' in str(jj) or
    740.                                     'y' in str(jj) or
    741.                                     'z' in str(jj) or
    742.                                     't' in str(jj)):
    743.                                 if km.name in repl.keys():
    744.                                     repl[km.name].append(str(jj))
    745.                                 else:
    746.                                     repl[km.name] = [str(jj)]
    747.                     else:
    748.                         ix = -1
    749.                         undo = 0
    750.                         for i in ams:
    751.                             ix += 1
    752.                             try:
    753.                                 if isinstance(i, (sympy.core.numbers.Float,
    754.                                                   sympy.core.numbers.Integer,
    755.                                                   int, float, complex)):
    756.  
    757.                                     try:
    758.                                         repl[tuple(un[indx])[ix]] = [str(eval(str(i)))]
    759.                                     except IndexError:
    760.  
    761.                                         repl[tuple(un[indx])[undo]].append(str(eval(str(i))))
    762.                                     for i in repl.keys():
    763.                                         if i in lf_:
    764.                                             lf_ = lf_.replace(i, str(repl[i][0]))
    765.                                 elif isinstance(i, (tuple, list, set)):
    766.                                     inxx = -1
    767.                                     for jj in i:
    768.                                         inxx += 1
    769.                                         if not ('x' in str(jj) or
    770.                                                 'y' in str(jj) or
    771.                                                 'z' in str(jj) or
    772.                                                 't' in str(jj)):
    773.                                             try:
    774.                                                 repl[tuple(un[indx][inxx])[ix]] = [str(jj)]
    775.                                             except IndexError:
    776.                                                 repl[tuple(un[indx][inxx])[undo]].append(str(jj))
    777.                                 else:
    778.                                     continue
    779.                             except TypeError as err:
    780.                                 print("get:err", err)
    781.                                 continue
    782.                             else:
    783.                                 undo = ix
    784.                     if len(repl.values()) >= 4: is_exit = True;break
    785.             if is_exit:
    786.                 solves.append(repl)
    787.             if printout:
    788.                 print(f"\r{(iVal + 1) / (tnums + 1) * 100}%", end="",
    789.                       flush=True)
    790.  
    791.         if printout:
    792.             print("\nsucceed in Calculating Values...")
    793.         num = 0
    794.         for i in solves:
    795.             if printout:
    796.                 print(f"\r{(num + 1) / len(solves) * 100}%", end="", flush=True)
    797.             num += 1
    798.             x_ = []
    799.             for j in i["x"]:
    800.                 x_.append(eval(j))
    801.             x.append(x_)
    802.             y_ = []
    803.             for j in i["y"]:
    804.                 y_.append(eval(j))
    805.             y.append(y_)
    806.             z_ = []
    807.             for j in i["z"]:
    808.                 z_.append(eval(j))
    809.             z.append(z_)
    810.             t_ = []
    811.             for j in i["t"]:
    812.                 t_.append(eval(j))
    813.             t.append(t_)
    814.  
    815.         del x_, y_, z_, t_, solves
    816.  
    817.         if printout:
    818.             print("\nsucceed in Saving Values...")
    819.  
    820. def format():
    821.     global x,y,z,t
    822.     nx={}
    823.     ny={}
    824.     nz={}
    825.     nt=[]
    826.  
    827.     vx=[]
    828.     vy=[]
    829.     vz=[]
    830.  
    831.     for i,j,k,v in zip(x,y,z,t):
    832.         for it,jt,kt,vt in zip(i,j,k,v):
    833.  
    834.             if vt not in nx.keys():
    835.                 nx[vt]=[]
    836.             if vt not in ny.keys():
    837.                 ny[vt]=[]
    838.             if vt not in nz.keys():
    839.                 nz[vt]=[]
    840.  
    841.             nx[vt].append(it)
    842.             ny[vt].append(jt)
    843.             nz[vt].append(kt)
    844.             nt.append(vt)
    845.  
    846.     nt.sort()
    847.  
    848.     for num in set(nt):
    849.         for i,j,k in zip(nx[num],ny[num],nz[num]):
    850.             vx.append(i)
    851.             vy.append(j)
    852.             vz.append(k)
    853.  
    854.     x=vx
    855.     y=vy
    856.     z=vz
    857.     t=nt
    858.  
    859.  
    860.  
    861. def create_img(tnum, fig=None, ax=None):
    862.     global x, y, z, t
    863.  
    864.     if fig is None:
    865.         fig = plt.figure()
    866.     if ax is None:
    867.         ax = Axes3D(fig)
    868.         ax.legend(loc='best')
    869.         ax.axis('off')
    870.     if isinstance(tnum, list):
    871.         tnum = tnum[0]
    872.     ax.cla()
    873.  
    874.     x1 = []
    875.     y1 = []
    876.     z1 = []
    877.     try:
    878.         i = t.index(tnum)
    879.     except (IndexError, ValueError) as err:
    880.         i = 0
    881.     try:
    882.         while t[i] == tnum:
    883.             if isinstance(x[i], (tuple, list, set)):
    884.                 x[i] = x[i][0]
    885.             if isinstance(x[i], complex):
    886.                 i += 1
    887.                 continue
    888.             if isinstance(y[i], (tuple, list, set)):
    889.                 y[i] = y[i][0]
    890.             if isinstance(y[i], complex):
    891.                 i += 1
    892.                 continue
    893.             if isinstance(z[i], (tuple, list, set)):
    894.                 z[i] = z[i][0]
    895.             if isinstance(z[i], complex):
    896.                 i += 1
    897.                 continue
    898.             x1.append(x[i])
    899.             y1.append(y[i])
    900.             z1.append(z[i])
    901.             i += 1
    902.     except IndexError:
    903.         pass
    904.  
    905.     scatter = ax.scatter(x1, y1, z1, c='r')
    906.     return scatter,
    907.  
    908.     # plt.savefig('./1.jpg',)
    909.  
    910.  
    911. def Graphic_main(formula, modefunc, tmin_=-10, tmax_=10, tnums=100,
    912.                  min_=-10, max_=10, nums=100, printout=True):
    913.     global x, y, z, t
    914.     try:
    915.         modefunc(formula=formula, tmin_=tmin_, tmax_=tmax_, tnums=tnums,
    916.                  min_=min_, max_=max_, nums=nums, printout=printout)
    917.     except Exception:
    918.         modefunc(formula=formula, tmin_=tmin_, tmax_=tmax_, tnums=tnums, printout=printout)
    919.     format()
    920.     fig = plt.figure()
    921.  
    922.  
    923.  
    924.     ax = Axes3D(fig, auto_add_to_figure=False)
    925.  
    926.     print(min(x),min(y),min(z),"\n",max(x),max(y),max(z))
    927.  
    928.     ax.set_ylim(min(y),max(y))
    929.     ax.set_xlim(xmin=min(x), xmax=max(x))
    930.     ax.set_zlim(zmin=min(z), zmax=max(z))
    931.  
    932.     fig.add_axes(ax)
    933.  
    934.  
    935.     # ax.legend(loc="best")
    936.  
    937.  
    938.     def demo(num):
    939.         a = tmin_ + ((tmax_ - tmin_) / tnums) * (num % tnums)
    940.         print(f"\rwait:{num-tnums}/{tnums}",end="")
    941.  
    942.         return create_img(a, fig, ax),
    943.  
    944.     ani = animation.FuncAnimation(fig, demo, fargs=[],
    945.                                   interval=int(((tmax_ - tmin_) / tnums) * 1000))
    946.     ani.save('test.mp4')
    947.     # 调试部分
    948.  
    949.  
    950. print("program running")
    951.  
    952. """write()._write_t_x("z=t\n3x+y=z\nx=sin(y)",tnums=100,
    953.       printout=True)
    954. print(t)
    955. create_img([0.0])"""
    956. Graphic_main("sin(x)=y\nz=t\ncos(y)=z", modefunc=write()._write_t_x, tnums=10, nums=10)
    957. print("end...")

    再次感谢各位的观看！

    ---------------------------------【完】--------------------------------------

  • 相关阅读:
    Agile Development
    使用 Selenium 或其他工具模拟浏览器使用及语法代码
    计算机毕业设计Java交通事故档案管理系统(源码+mysql数据库+系统+lw文档）
    PowerCLi 一键批量部署OVA 到esxi 7
    nodejs+vue＋python＋PHP＋微信小程序南七街道志愿者服务平台的设计与实现-计算机毕业设计推荐
    C#中使用python(基于Process-C#自带库)
    视频监控/视频汇聚/视频云存储EasyCVR平台接入华为ivs3800平台提示400报错，该如何解决？
    Linux文件和文件夹命令详解
    java-php-python-ssm网上商城系统计算机毕业设计
    大白话讲清搞好就能“年薪百万”的SpringCloud微服务
  • 原文地址：https://blog.csdn.net/html_finder/article/details/126284437