filtic函数详解（附3个实例）

科技前沿 • 2025-02-28 08:55 • 阅读 36

大家好，我是讯享网，很高兴认识大家。

filtic函数

目的是为了求得初始条件

Z = filtic( B, A, Y, X )

讯享网

讯享网Z = filtic( B, A, Y )

表示在过去X=0

举3个例子

例1

求y(n) + 1.12y(n-1) = 0.1x(n) + 0.2x(n-1)的零输入响应，初始条件y(-1)=1

讯享网

 clc; clear all; close all; B = [0.1, 0.2]; % Numerator Coefficients A = [1, 1.12]; % Denominator Coefficients Y = [1]; % Initial conditions for output xic = filtic(B,A,Y) % Finding initial conditions for the system yzi = filter(B,A,zeros(1,20),xic) % Zero Input response stem(yzi)

在这里插入图片描述

例2

讯享网 clc; clear all; close all; A=[1,-0.6,0.08]; B=[1]; Y=[0,0]; xic=filtic(B,A,Y); yni=filter(B,A,[1,zeros(1,19)],xic); n1=0:length(yni)-1; subplot(2,1,1); stem(n1,yni,'.'); title('unit impulse response'); xlabel('n1');ylabel('yni'); yns=filter(B,A,ones(1,20),xic); n2=0:length(yns)-1; subplot(2,1,2); stem(n2,yns,'.'); title('unit step response'); xlabel('n2');ylabel('yns');

在这里插入图片描述

例3

 clc; clear all; close all; A=[1,-0.7,0.1]; B=[2,0,-1]; Y=[0,0]; xic=filtic(B,A,Y); yni=filter(B,A,[1,zeros(1,19)],xic); n1=0:length(yni)-1; subplot(2,1,1); stem(n1,yni,'.'); title('unit impulse response'); xlabel('n1');ylabel('yni'); yns=filter(B,A,ones(1,20),xic); n2=0:length(yns)-1; subplot(2,1,2); stem(n2,yns,'.'); title('unit step response'); xlabel('n2');ylabel('yns');

在这里插入图片描述

filtic函数详解（附3个实例）

filtic函数

例1

例2

例3

相关推荐