系统1：y[n]=x[n-3],对系统定义系数向量a1和b1，MATLAB用freqz产生该系统的频率响应h1(e^jw)的幅值和相位

由于系统1是一个时域系统，不能通过系数向量a1和b1直接使用freqz函数产生频率响应。但我们可以手动计算出系统的差分方程，然后使用freqz函数来生成系统的幅度和相位响应。

系统1的差分方程为：

y[n] = x[n - 3]

将其变换为Z域：

Y(z) = z^-3 X(z)

因此，系统1的传递函数为：

H(z) = Y(z) / X(z) = z^-3

现在我们可以使用MATLAB中的freqz函数来生成系统的幅度和相位响应。下面是MATLAB代码：

b1 = [0 0 0 1];  % 系统1的b系数
a1 = 1;          % 系统1的a系数
[h1, w] = freqz(b1, a1);  % 使用freqz函数计算系统1的幅度和相位响应

该代码将生成系统1的幅度响应和相位响应，并将它们存储在h1和w中。h1是系统的复数频率响应，w是频率向量。您可以使用MATLAB的abs和angle函数来分别提取系统的幅度和相位响应：

mag1 = abs(h1);    % 提取系统1的幅度响应
phase1 = angle(h1);  % 提取系统1的相位响应

现在，mag1和phase1就是系统1的幅度响应和相位响应。