y=[1;1]; %y(1)=y(2)=1,根据问题，自行设置。n=20; %% 自己设置y=[y;zeros(n-2,1)]; % 这一行是开辟空间，可以不要。 for i=3:ny(i)=5.775+1.04*y(i-1)+-0.48*y(i-2);endplot((1:n)',y);
