新手小张 2021-07-31 15:54 采纳率: 50%
浏览 78
已结题

关于用Python的Sympy库求解微分方程解析解时返回值是空列表的问题

1、对于一个微分方程,使用MATLAB可以很轻松得到其解析解,但是使用Python的sympy库却返回一个空列表。
2、
MATLAB代码如下:

clear;   
clc;
syms T(x) x m n q A;
x1 = 339.5;
x2 = 410.5;
T1 = 255;
T2 = 25;
m=2.4772
n = 0.0071
q = (m+n*(T2+T1)/2)*((T2-T1)/(x2-x1))
e = diff(T,x) == q/(m+n*T(x));  %待求解方程
cond = [T(x1)==T1];     %边界条件
T = dsolve(e,cond)

结果

T =

(2^(1/2)*(7259410169/2 - 7983760*x)^(1/2))/71 - 24772/71

Python代码如下

import sympy as sy
def solveequ(x1,x2,T1,T2,m=2.4772,n=0.0071):
    T_e = sy.symbols('T_e', cls=sy.Function)
    x = sy.symbols("x")
    Q = ((T2-T1)/(x2-x1))*(m+n*(T2+T1)/2)
    eq = sy.Eq((m+n*T_e(x))*T_e(x).diff(x),Q)
    edge = {}
    edge[T_e(x1)] = T1
    getT = sy.dsolve(eq,T_e(x), ics=edge)
    return getT
T0 = 25
T1 = 175
T2 = 195
T3 = 235
T4 = 255
T_e = solveequ(339.5,410.5,T4,T0)
[]
  • 写回答

1条回答 默认 最新

  • 新手小张 2021-07-31 20:41
    关注

    Python3.7
    Sympy1.8

    评论

报告相同问题?

问题事件

  • 系统已结题 8月8日
  • 创建了问题 7月31日

悬赏问题

  • ¥15 有两个非常“自以为是”烦人的问题急期待大家解决!
  • ¥30 STM32 INMP441无法读取数据
  • ¥100 求汇川机器人IRCB300控制器和示教器同版本升级固件文件升级包
  • ¥15 用visualstudio2022创建vue项目后无法启动
  • ¥15 x趋于0时tanx-sinx极限可以拆开算吗
  • ¥15 pyqt信号槽连接写法
  • ¥500 把面具戴到人脸上,请大家贡献智慧,别用大模型回答,大模型的答案没啥用
  • ¥15 任意一个散点图自己下载其js脚本文件并做成独立的案例页面,不要作在线的,要离线状态。
  • ¥15 各位 帮我看看如何写代码,打出来的图形要和如下图呈现的一样,急
  • ¥30 c#打开word开启修订并实时显示批注