y的运行结果：
程序：n=0:30; %自变量n的取值范围为0到50的整数值
h=(0.9.^n).*[n>=0]; %函数h用逻辑表达式和函数值表示 x=[n<=10]; %x为单位阶跃函数 y=conv(double(x),h);
结果
y =
Columns 1 through 13
1.0000 1.9000 2.7100 3.4390 4.0951 6.1258 6.5132 6.8619 6.1757 5.5581
Columns 14 through 26
5.0023 4.5021 4.0519 3.6467 3.2820 2.1533 1.9380 1.7442 1.5698 1.4128
Columns 27 through 39
1.2715 1.1444 1.0299 0.9269 0.8342 0.4161 0.3364 0.2646 0.2000 0.1418
Columns 40 through 52
0.0895 0.0424 0 0 0 0 0 0 0 0
Columns 53 through 61
0 0 0 0 0 0
stem(n,x,'.') %将x、h、y绘图显示 stem(n,h,'.') n=0:60;
stem(n,y,'.')
运行结果：
x(n)
- 16 -
4.6856 5.2170 2.9538 2.6584 0.7127 0.6033 0 0 0 0 5.6953 2.3926 0.5048 0 0
h(n)：
y(n)：
- 17 -
2.求因果线性移不变系统
y(n)?0.81y(n?2)?x(n)?x(n?2)
j?的单位抽样响应h(n)，并绘出H(e
)的幅频及相频特性曲线。
源程序：
a=[1 0 -0.81]; %a为y的系数序列 b=[1 0 -1]; %b为x的系数序列 n=0:50; %n为自变量取值范围
x=[n<=0]; %输入函数x为单位抽样函数 y=filter(b,a,x) %求解差分方程的零状态响应
运行结果： y的求解结果：
y =
Columns 1 through 9
- 18 -
1.0000 0 -0.1900 0 -0.1539 0 -0.1247 0 -0.1010
Columns 10 through 18
0 -0.0818 0 -0.0662 0 -0.0537 0 -0.0435 0
Columns 19 through 27
-0.0352 0 -0.0285 0 -0.0231 0 -0.0152
Columns 28 through 36
0 -0.0123 0 -0.0099 0 -0.0081 0
Columns 37 through 45
-0.0053 0 -0.0043 0 -0.0035 0 -0.0023
Columns 46 through 51
0 -0.0018 0 -0.0015 0 -0.0012
源程序：
n=[0:50]; %将y绘图显示 stem(n,y,'.')
运行结果：
- 19 -
-0.0187 0 -0.0028 0 -0.0065 0
源程序：
c=zeros(1,50); %产生一个n=0时取1，其余取0的序列c，为下一步准备 c(1)=1
w=0:2*pi/30:2*pi %w在0到2*pi内取等间隔离散值 X=freqz(y,c,w) %X为所求频域函数H
运行结果：
将X的求解结果序列显示：
X =
Columns 1 through 5
0.0052 0.8900 + 0.4403i 1.0458 + 0.2476i 1.0820 + 0.1556i 1.0959 + 0.1034i
Columns 6 through 10
1.1009 + 0.0673i 1.1032 + 0.0375i 1.1052 + 0.0118i 1.1052 - 0.0118i 1.1032 - 0.0375i
- 20 -