利用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   大湿Mastwet  阅读(1659)  评论(0编辑  收藏  举报

编辑推荐:
· AI与.NET技术实操系列(二):开始使用ML.NET
· 记一次.NET内存居高不下排查解决与启示
· 探究高空视频全景AR技术的实现原理
· 理解Rust引用及其生命周期标识(上)
· 浏览器原生「磁吸」效果!Anchor Positioning 锚点定位神器解析
阅读排行:
· 全程不用写代码,我用AI程序员写了一个飞机大战
· DeepSeek 开源周回顾「GitHub 热点速览」
· 记一次.NET内存居高不下排查解决与启示
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· .NET10 - 预览版1新功能体验(一)

导航

< 2025年3月 >
23 24 25 26 27 28 1
2 3 4 5 6 7 8
9 10 11 12 13 14 15
16 17 18 19 20 21 22
23 24 25 26 27 28 29
30 31 1 2 3 4 5
点击右上角即可分享
微信分享提示