第二章—微分方程答案.doc
文本预览下载声明
1. a)
b)
c)
% a)
a=.5; b=[0 2];
n = 0:30;
y0 = 0;
x0 = 0;
x = zeros(size(n));
x(1) = 1;
y=recur(a,b,n,x,x0,y0);
subplot(221),stem(n,y)
title(’1a)’)
xlabel(’n’)
ylabel(’y[n]’)
y(1:5) % displays y for n=0 to 5, used to double-check hand answers
% b)
a=2; b=[0 2];
n = 0:30;
y0 = 0;
x0 = 0;
x = zeros(size(n));
x(1) = 1;
y=recur(a,b,n,x,x0,y0);
subplot(222),stem(n,y)
title(’1b)’)
xlabel(’n’)
ylabel(’y[n]’)
y(1:5)
% c)
a=[1.2 .32]; b = [1 -1];
y0 = [1 2];
x0 = 0;
x = ones(size(n));
y=recur(a,b,n,x,x0,y0);
subplot(223),stem(n,y)
title(’1c)’)
xlabel(’n’)
ylabel(’y[n]’)
y(1:5)
显示全部