题目:用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');