matlab编制程序，包括产生输入信号、单位脉冲响应序列的子程序，用filter函数或conv函数求解系统输出响应的主程序。程序中要有绘制信号波形的功能。 列出源程序清单：

Matlab是一种强大的数值计算工具，常用于信号处理和控制系统设计。编写这样的程序通常会涉及以下几个步骤：

首先，你需要创建一个产生输入信号的子程序。例如，你可以使用`sin()`函数生成正弦波作为输入。这是一个简单的例子：

function input_signal = generate_sinusoidal_input(freq, duration)
    sample_rate = 1000; % Hz (假设采样频率)
    time_vec = 0:1/sample_rate:duration; % 时间向量
    input_signal = sin(2*pi*freq*time_vec); % 生成特定频率的正弦信号
end

接下来，定义单位脉冲响应序列的子程序。如果是一个线性滤波器，可以使用`impz()`函数从传递函数得到 impulse response：

function h = get_impulse_response(b, a) % b 和 a 分别是滤波器的系数
    h = impz(b, a);
end

主程序部分会结合这两个子程序，并利用`filter()`或`conv()`函数计算系统的输出响应：

function system_output = main_program(input_freq, filter_coeffs, duration)
    % 生成输入信号
    input_signal = generate_sinusoidal_input(input_freq, duration);
    
    % 获取单位脉冲响应
    h = get_impulse_response(filter_coeffs.b, filter_coeffs.a);
    
    % 使用filter()函数（线性相位滤波）
    if isvector(filter_coeffs.b) % 确保滤波器是一阶系统
        system_output = filter(h, 1, input_signal);
    elseif ismatrix(filter_coeffs.b) % 高阶系统用conv()
        system_output = conv(input_signal, h);
    else
        error('Invalid filter coefficients.');
    end
    
    % 绘制输入信号和输出信号波形
    plot(time_vec, [input_signal, system_output]);
    legend({'Input Signal', 'System Output'});
end

在这个示例中，`filter_coeffs`应该包含滤波器的系数，`input_freq`是输入信号的频率。在主程序中，你可以通过调用`main_program()`并传入相应参数来运行整个过程。