sy2.txt

几个Matlab编程示例

%离散系统的差分方程、冲激响应和卷积分析

%2-1
N=21;
a=[1 -1];
b=[1 0.75 0.125];
x1=[1 zeros(1,N-1)];
x2=ones(1,N);
k=0:1:N-1;
y1=filter(a,b,x1);
y2=filter(a,b,x2);
subplot(2,1,1);stem(k,y1);
xlabel('n');ylabel('幅度');title('单位冲激响应');
subplot(2,1,2);stem(k,y2);
xlabel('n');ylabel('幅度');title('阶跃响应');

%2-1
N=21;
a=0.25*[1 1 1 1];
b=[1];
x1=[1 zeros(1,N-1)];
x2=ones(1,N);
k=0:1:N-1;
y1=filter(a,b,x1);
y2=filter(a,b,x2);
subplot(2,1,1);stem(k,y1);
xlabel('n');ylabel('幅度');title('单位冲激响应');
subplot(2,1,2);stem(k,y2);
xlabel('n');ylabel('幅度');title('阶跃响应');