python求解一维弦振动方程

1、一维弦振动方程

2、差分公式及求解

python 复制代码
import numpy as np
import matplotlib.pyplot as plt
l=1    #长度
nx = 101  #x方向网格节点数量
dx = l/(nx-1)  #空间网格步长
t=6    #总时间
dt = 0.01  # 时间步长
t_set=np.array([0.00,0.10,0.20,0.30,0.40,0.50,0.60]) #需要观察的时间
t_get=t_set/dt #需要观察的时间的步数

# 指定初始条件
xp=np.linspace(0,l,num=nx)
u = np.sin(np.pi*xp)
utp=xp*(1-xp)*dt
u1=np.zeros(nx)
u1=u+utp
un1=u.copy()
un2=u1.copy()
# 计算
u_get=[] #需要观察的时间的波形
for n in range(int(t/dt)+1):
    if n in t_get.astype(int):
        if n==0:
            u_get.append(u)
            continue
        if n==1:
            u_get.append(u1)
            continue
        u_get.append(un2)
    un = np.zeros(nx)
    for i in range(1,nx-1):
       un[i] = dt**2/dx**2*(un2[i+1]+un2[i-1])+2*(1-dt**2/dx**2)*un2[i] - un1[i]
    un1=un2.copy()
    un2 = un.copy()


fig = plt.figure()
plt.rcParams['xtick.direction'] = 'in'  # 将x周的刻度线方向设置向内
plt.rcParams['ytick.direction'] = 'in'  # 将y轴的刻度方向设置向内
clist = ['blue', 'red', 'green', 'black', 'darkgreen', 'lime', 'gold', 'purple', 'green', 'cyan', 'salmon', 'grey',
             'mediumvioletred', 'darkkhaki', 'gray', 'darkcyan', 'violet', 'powderblue']
#绘制选的时间步的结果
for i in range(len(t_get)):
    plt.plot(xp, u_get[i], c=clist[i],lw =3,label= f't={t_set[i]}s')
plt.legend()
plt.show()

**注意:**时间不长和空间步长之比dt/dx<=1时,次解法是稳定且收敛的

相关推荐
用户8356290780513 分钟前
用Python高效处理Excel数据:Excel数据读取指南
后端·python
我星期八休息26 分钟前
深入理解跳表(Skip List):原理、实现与应用
开发语言·数据结构·人工智能·python·算法·list
蒋星熠35 分钟前
如何在Anaconda中配置你的CUDA & Pytorch & cuNN环境(2025最新教程)
开发语言·人工智能·pytorch·python·深度学习·机器学习·ai
合作小小程序员小小店1 小时前
机器学习介绍
人工智能·python·机器学习·scikit-learn·安全威胁分析
JavaEdge在掘金1 小时前
掌握Spring IoC容器和Bean作用,轻松实现依赖注入!
python
flysh052 小时前
pyAutoGUI 模块主要功能介绍-(2)键盘功能
python·pyautogui
强盛小灵通专卖员2 小时前
闪电科创 SCI专业辅导
python·深度强化学习·研究生·ei会议·导师·sci期刊
跟橙姐学代码3 小时前
自动化邮件发送的终极秘籍:Python库smtplib与email的完整玩法
前端·python·ipython
扯淡的闲人3 小时前
多语言编码Agent解决方案(2)-后端服务实现
开发语言·python·深度学习
蒋星熠3 小时前
深度学习实战指南:从神经网络基础到模型优化的完整攻略
人工智能·python·深度学习·神经网络·机器学习·卷积神经网络·transformer