MATLAB系统响应求解

MATLAB实验报告

姓名：李金玮

学号：14061114

班级：14184111

实验内容及要求

系统差分方程为：

1．若输入序列为，编程求此系统输出序列y(n) ，并画出其波形。

2．编程得到系统频率响应，并画图。

3．编程得到系统的零极点分布图，分析系统的因果性和稳定性。

clear all

close all

n=1:304

k1=zeros(1,300);

k1(1)=1;

a=[k1,0,0,0,0];

b=[0,k1,0,0,0];

c=[0,0,k1,0,0];

d=[0,0,0,k1,0];

e=[0,0,0,0,k1];

x=a+2*b+3*c+4*d+5*e;

aa=[1,-1.62,1.28]

bb=[0.5,0.1]

y=filter(bb,aa,x);

subplot(3,1,1)

plot(n,y)

subplot(3,1,2)

xk1=fft(y,304)

plot(n,xk1)

subplot(3,1,3)

zplane(bb,aa)

[z,p,k]=tf2zp(bb,aa) sys=zpk(z,p,k)