qq_32109565的博客

原创 DQM法求解微分方程

求解方程 𝑦′′(𝜓,𝑡)+𝑎1𝑦′(𝜓,𝑡)+𝑎2𝑦(𝜓,𝑡)+𝑎3cos⁡(𝑎4𝑡)𝑦(𝜓,𝑡)=0y′′(ψ,t)+a1​y′(ψ,t)+a2​y(ψ,t)+a3​cos(a4​t)y(ψ,t)=0 在 𝜓=0ψ=0 处的动力系统响应。最后，我们绘制了中间节点（即 𝜓=0ψ=0 处）的系统响应随时间的变化。我们首先定义了问题的常数和初始条件，然后创建了切比雪夫节点。接下来，我们应用了边界条件，并初始化了解矩阵。以下是运行的代码，特别注意矩阵尺寸的匹配。我们手动构建了一阶导数近似矩阵。