差分方程求解

题目：用Matlab仿真。设一个离散时间系统的差分方程为：

                  y(n)+0. 5y(n-1)-0.3y(n-2)+0.2y(n-3)=0.75x(n)+x(n-1)-0.6x(n-2)-0.4x(n-3) ,

        分别输入三个不同的序列x1(n) , x2(n)和x(n)=4x1(n)-3x2(n) ,并计算对于输入x1(n)，x2(n)和x(n)的输出y1(n)，y2(n)和y(n)，判断该系统是否为线性系统。

% 定义系统参数  
a = [1, 0.5, -0.3, 0.2];  % 反馈系数  
b = [0.75, 1, -0.6, -0.4];  % 输入系数  
  
% 定义输入信号  
x1 =ones(50,1);  % 输入信号1 
x2 =[zeros(1,3) ones(1,50)];   % 输入信号2  
x = 4*x1 - 3*x2;  % 通过x1和x2生成的信号作为输入信号  
    
y1=filter(b,a,x1);    %filter函数基本用法 y=filter (b,a,x); 
y2=filter(b,a,x2);
y=filter(b,a,x);
 
% 绘制结果  
figure;  
subplot(3, 1, 1);  
stem(x1(1:20));  
title('输入信号x1');  
subplot(3, 1, 2);  
stem(x2(1:20));  
title('输入信号x2');  
subplot(3, 1, 3);  
stem(x(1:20));  
title('输入信号x');  
  
figure;  
subplot(3, 1, 1);  
stem(y1(1:20));  
title('输出信号y1');  
subplot(3, 1, 2);  
stem(y2(1:20));  
title('输出信号y2');  
subplot(3, 1, 3);  
stem(y(1:20));  
title('输出信号y');