利用python、scipy求解差分方程的解
文档:https://docs.scipy.org/doc/scipy/reference/generated/scipy.signal.lfilter.html#scipy.signal.lfilter
import matplotlib.pyplot as plt import numpy as np from scipy import signal nmin = 0 nmax = 8 n = np.arange(nmin,nmax+1,1) nl = len(n) den = np.array([1,0.6]) num = np.array([1,2]) xn = (0.1*7)**n y0 = np.array([4,2]) #0状态 x01 = np.array([0]); zi1 = signal.lfilter_zi(num,den) y1,_ = signal.lfilter(num,den,xn,zi=zi1*x01) #0输入 x02 = np.zeros(nl) zi2 = signal.lfiltic(num,den,y0) y2,_ = signal.lfilter(num,den,x02,zi = zi2) #全响应 y3,_ = signal.lfilter(num,den,xn,zi=zi1) #计算单位冲激响应 t4,y4 = signal.dimpulse((num,den,1),n=nl) plt.subplot(411) plt.stem(n,y1) plt.ylim(-5,5) plt.subplot(412) plt.stem(n,y2) plt.ylim(-5,5) plt.subplot(413) plt.stem(n,y3) plt.ylim(-3,3) plt.subplot(414) plt.stem(t4,np.squeeze(y4)) plt.ylim(-5,5) plt.show()
posted on 2021-04-03 17:18 大湿Mastwet 阅读(1660) 评论(0) 编辑 收藏 举报
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· AI与.NET技术实操系列(二):开始使用ML.NET
· 记一次.NET内存居高不下排查解决与启示
· 探究高空视频全景AR技术的实现原理
· 理解Rust引用及其生命周期标识(上)
· 浏览器原生「磁吸」效果!Anchor Positioning 锚点定位神器解析
· 全程不用写代码,我用AI程序员写了一个飞机大战
· DeepSeek 开源周回顾「GitHub 热点速览」
· 记一次.NET内存居高不下排查解决与启示
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· .NET10 - 预览版1新功能体验(一)