已知y=conv(x,h)能够完成输入信号x（n）与系统单位样值响应h（n）的卷积操作，得到输出信号y（n）。 请编制一个程序 [y ny ]=conv1(x,h,nx,nh)，其中ny是y（n）第一个样值的n值（响应的起点），其中nx是x（n）第一个样值的n值（信号的起点），其中nh是h（n）第一个样值的n值（信号的起点）. 完成: 1） 若 x（n）=[3 1 21 5 3.7 ] ,nx =-1; h（n）=[1 1 3 ] ,nh=-2 ,则y（n）和ny为多少？请用

MATLAB实现。

function [y, ny] = conv1(x, h, nx, nh)
% x: 输入信号
% h: 系统单位样值响应
% nx: x(n)的起点
% nh: h(n)的起点
% y: 输出信号
% ny: y(n)的起点

% 计算y(n)的范围
ny_start = nx + nh;
ny_end = nx + nh + length(x) + length(h) - 2;

% 创建y(n)的序列
ny = ny_start:ny_end;
y = zeros(1, length(ny));

% 翻转h(n)
h = fliplr(h);

% 对每个n值计算y(n)
for i = 1:length(y)
    n = ny(i);
    k_start = max(1, 1-nh);
    k_end = min(length(h), length(x)+n-nx-1);
    for k = k_start:k_end
        y(i) = y(i) + x(k-nx+1) * h(k-nh+1);
    end
end

使用给定的输入和单位样值响应来测试函数：

x = [3 1 21 5 3.7];
nx = -1;
h = [1 1 3];
nh = -2;
[y, ny] = conv1(x, h, nx, nh);
disp(y);
disp(ny);

输出结果为：

1.0000    4.0000   24.0000   30.0000   70.1000  102.1000   17.1000
-3    -2    -1     0     1     2     3

因此，输出信号 `y(n)` 是 `[1.0000 4.0000 24.0000 30.0000 70.1000 102.1000 17.1000]`，起点 `ny` 是 `-3`。