文章目录
一、卷积 与 " 线性常系数差分方程 "
" 线性常系数差分方程 " 不能使用 卷积函数 conv 函数进行求解 , 因为卷积的右侧没有 y ( n ) y(n) y(n) , 卷积公式如下 :
y ( n ) = ∑ m = − ∞ + ∞ x ( m ) h ( n − m ) = x ( n ) ∗ h ( n ) y(n) = \sum^{+\infty}_{m = -\infty} x(m) h(n-m) = x(n) * h(n) y(n)=m=−∞∑+∞x(m)h(n−m)=x(n)∗h(n)
而 " 线性常系数差分方程 " 如下 :
y ( n ) = ∑ i = 0 M b i x ( n − i ) − ∑ i = 1 N a i y ( n − i ) n ≥ M y(n) = \sum_{i = 0}^M b_i x(n - i) - \sum_{i = 1}^N a_i y(n - i) \ \ \ \ \ \ \ n \geq M y(n)=i=0∑Mbix(n−i)−i=1∑Naiy(n−i) n≥M
在 " 线性常系数差分方程 " 公式的右侧比 卷积 公式中 , 多了一个 ∑ i = 1 N a i y ( n − i ) \sum_{i = 1}^N a_i y(n - i) ∑i=1Naiy(n−i) 项 , 其中有 y ( n ) y(n) y(n) 序列 , 这样就无法使用 conv 卷积函数求解 " 线性常系数差分方程 " ;
二、使用 matlab 求解 " 线性常系数差分方程 "
matlab 中 , 使用 filter 函数, 求解 " 线性常系数差分方程 " ;
参考文档 :
- filter 函数 : https://ww2.mathworks.cn/help/matlab/ref/filter.html
- filtic 函数 : https://ww2.mathworks.cn/help/signal/ref/filtic.html
filter 函数语法如下 :
yn = filter(B, A, xn, xi)
xi = filtic(B, A, ys, xs)
" 线性常系数差分方程 " 公式如下 :
y ( n ) = ∑ i = 0 M b i x ( n − i ) − ∑ i = 1 N a i y ( n − i ) n ≥ M y(n) = \sum_{i = 0}^M b_i x(n - i) - \sum_{i = 1}^N a_i y(n - i) \ \ \ \ \ \ \ n \geq M y(n)=i=0∑Mbix(n−i)−i=1∑Naiy(n−i) n≥M
matlab 中的 filter 函数中的参数 与 " 线性常系数差分方程 " 公式项的对应关系 :
① B 参数 : filter 函数中的 B 向量
B = [ b 0 , b 1 , ⋯ , b M ] B = [b_0 , b_1, \cdots,b_M] B=[b0,b1,⋯,bM]
就是公式中的 b i b_i bi , 注意 i i i 范围是 [ 0 , M ] [0,M] [0,M] ;
② A 参数 : filter 函数中的 A 向量
A = [ a 1 , a 2 , ⋯ , a N ] A = [a_1 , a_2, \cdots,a_N] A=[a1,a2,⋯,aN]
就是公式中的 a i a_i ai , 注意 i i i 范围是 [ 1 , N ] [1,N] [1,N] ;
③ xn 参数 : 输入序列 对应的 向量 ;
④ xi 参数 : 该参数 与 ys 和 xs 条件有关 , ys 和 xs 是初始条件向量 , 分别是 :
y s = [ y ( − 1 ) , y ( − 2 ) , ⋯ , y ( − N ) ] y_s = [y(-1), y(-2), \cdots , y(-N)] ys=[y(−1),y(−2),⋯,y(−N)]
x s = [ x ( − 1 ) , x ( − 2 ) , ⋯ , x ( − N ) ] x_s = [x(-1), x(-2), \cdots , x(-N)] xs=[x(−1),x(−2),⋯,x(−N)]
xi 是通过 filtic 函数 计算出来的 , 需要传入 A , B A , B A,B 向量 , 和 ys 和 xs 条件 ;
标签：方程,系数,卷积,差分,filter,线性 来源： https://blog.csdn.net/han1202012/article/details/123165905
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。