文章目錄

• 一、使用 matlab 求解 “ 線性常系數(shù)差分方程 “ 示例
• • 1、B 向量元素 : x(n) 參數(shù)
  • 2、A 向量元素 : y(n) 參數(shù)
  • 3、輸入序列
  • 4、matlab 代碼

一、使用 matlab 求解 “ 線性常系數(shù)差分方程 “ 示例

---

描述 某個(gè) " 線性時(shí)不變系統(tǒng) " 的 " 線性常系數(shù)差分方程 " 如下 :

$$y(n)=1.5x(n)+0.7y(n?1)$$

輸入序列 :

$$x(n)=\delta (n)$$

邊界條件 / 初始條件 :

$$y(?1)=1$$

求該 LTI 系統(tǒng)的 輸出序列 ;

線性常系數(shù)差分方程 公式 :

$$y(n)=\sum _{i=0}^M{b}_{i}x(n?i)?\sum _{i=1}^N{a}_{i}y(n?i)n\ge M$$

1、B 向量元素 : x(n) 參數(shù)

討論 $B$ 向量 , $B$ 向量是 $x(n)$ 的參數(shù) , 有幾個(gè) $x(n)$ 項(xiàng) , $B$ 向量 就有幾個(gè)元素 ;

上式中 $M=0$ , $x(n)$ 的項(xiàng)只有 $1$ 項(xiàng) , ${\sum }_{i=0}^M{b}_{i}x(n?i)$ 只有一項(xiàng) , 加和式只有一項(xiàng) , 因此對應(yīng)的 $B$ 向量 , 只有 $1$ 個(gè)元素 ;

B = [1.5];

2、A 向量元素 : y(n) 參數(shù)

下面討論 $A$ 向量 , $A$ 向量是 $y(n)$ 的參數(shù) , 有幾個(gè) $y(n)$ 項(xiàng) , $A$ 向量 就有幾個(gè)元素 ;

線性常系數(shù)差分方程 :

$$y(n)=1.5x(n)+0.7y(n?1)$$

將 $0.7y(n?1)$ 移到左邊 , 得到 :

$$y(n)?0.7y(n?1)=1.5x(n)$$

這里有 $2$ 個(gè) $y(n)$ 項(xiàng) , $A$ 向量的元素有兩個(gè) , $1,?0.7$ ;

A = [1, -0.7];

3、輸入序列

輸入序列 :

$$x(n)=\delta (n)$$

輸入序列 的元素個(gè)數(shù) , 等于 輸出序列 的元素個(gè)數(shù) ;

$n=0$ 時(shí) , $x(n)=1$ , 然后再次生成 $30$ 個(gè) $0$ 元素 , 放到 輸入序列 中 ;

輸入序列為 $\{1,\underset{30個(gè)0}{0,0,?,0}\}$ , 共 $31$ 個(gè)元素 ;

對應(yīng)的 matlab 代碼為

xn=[1,zeros(1,30)];

4、matlab 代碼

matlab 代碼 :

% 邊界條件 y(-1) = 1 , 這里設(shè)置 ys = 1 ys = 1;% 輸入序列 為 單位脈沖序列 xn=[1,zeros(1,30)]; % 線性常系數(shù)差分方程 中的 x(n) 項(xiàng)系數(shù) B=1.5;% 線性常系數(shù)差分方程 中的 y(n) 項(xiàng)系數(shù) A=[1, -0.7];% 等效 初始條件 的 輸入序列 xi xi=filtic(B,A,ys);% 輸出序列 yn=filter(B,A,xn,xi); %建立幕布 figure; %繪制 "輸出序列" 圖像 , 點(diǎn)用上三角表示 plot(yn, '^');% 打開網(wǎng)格 grid on;

繪圖效果 :

總結(jié)

以上是生活随笔為你收集整理的【数字信号处理】线性常系数差分方程 ( 使用 matlab 求解 “ 线性常系数差分方程 “ 示例 | A 向量分析 | B 向量分析 | 输入序列分析 | matlab 代码 )的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。