q705.m

来自「matlab在电子信息中的应用和实例中的源代码」· M 代码 · 共 38 行

%《MATLAB及在电子信息课程中的应用》第七章例7.5程序q705
% 时域离散系统线性性质验证
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
%N: 输入信号长度 
%B: 差分方程系数
%A: 差分方程系数
clear;close all;
N=64;n=0:N-1;m=20;
B=1;
A=[1,-0.8];
x1=0.8.^n;				                % 产生输入信号x1(n)
x=[(n>=0)&(n<32)];
x1=x1.*x;
y1=filter(B,A,x1); 	                    % 对x1(n)的响应y1(n)
x2=[(n-m)==0];
y2=filter(B,A,x2); 	                    % 对x2(n)的响应y2(n)
x3=5*x1+3*x2;
y3=filter(B,A,x3); 	                    % 对5x1(n)+3x2(n)的响应y3(n)
y=5*y1+3*y2;       	                    % y(n)=5y1(n)+3y2(n)
subplot(3,2,1);stem(n,y1,'.');
line([0,N],[0,0])
set(gcf,'color','w')                    % 设置图形背景色为白色
axis([0,N,min(y1),max(y1)]);
ylabel('y1(n)]')
subplot(3,2,2);stem(n,y2,'.');
line([0,N],[0,0])
axis([0,N,min(y2),max(y2)]);
ylabel('y2(n)]')
subplot(3,2,3);stem(n,y3,'.');
line([0,N],[0,0])
axis([0,N,min(y3),max(y3)]);
xlabel('n');ylabel('y3(n)]')
subplot(3,2,4);stem(n,y,'.');
line([0,N],[0,0])
axis([0,N,min(y),max(y)]);
xlabel('n');ylabel('y(n))')