编写一个[y,y1,y2]=mwave(f1,m1,f2,m2)函数，实现以下功能，并绘出y1、y2、y在t∈[0，2π]区间500个样点的图形。（其中调用参数2 ≤ f1、f2 ≤ 20 HZ；0.5 ≤ m1、m2 ≤ 2）matlab

`mwave` 函数是一个用于生成两个简谐波合成信号的MATLAB函数，通常用于模拟声波或其他周期性物理现象。这个函数需要四个输入参数：两个频率 `f1` 和 `f2`，以及它们对应的振幅 `m1` 和 `m2`。在这个函数中，我们假设频率是离散的，因此会创建一系列时间点 `t`，然后计算每个时间点上 y1 (第一个波) 和 y2 (第二个波) 的值。最后，我们将这三者相加得到总信号 `y`。

以下是该函数的基本实现：

function [y, y1, y2] = mwave(f1, m1, f2, m2)
    % 确保频率和振幅在给定范围内
    assert(2 <= f1 && f1 <= 20, 'Frequency f1 must be between 2 and 20 Hz');
    assert(0.5 <= m1 && m1 <= 2, 'Amplitude m1 must be between 0.5 and 2');
    assert(2 <= f2 && f2 <= 20, 'Frequency f2 must be between 2 and 20 Hz');
    assert(0.5 <= m2 && m2 <= 2, 'Amplitude m2 must be between 0.5 and 2');

    % 创建时间轴，范围从0到2*pi，步长为1/500
    t = linspace(0, 2*pi, 500);
    
    % 计算每个简谐波的值
    y1 = m1 * sin(2 * pi * f1 * t); % 第一个正弦波
    y2 = m2 * cos(2 * pi * f2 * t); % 第二个余弦波

    % 合成信号
    y = y1 + y2;

    % 绘制图形
    plot(t, y1, 'b', 'LineWidth', 1.5, 'DisplayName', 'Waveform 1');
    hold on;
    plot(t, y2, 'r', 'LineWidth', 1.5, 'DisplayName', 'Waveform 2');
    plot(t, y, 'g', 'LineWidth', 1.5, 'DisplayName', 'Total Signal');
    legend('show');
    xlabel('Time (s)');
    ylabel('Amplitude');
    title(sprintf('Sine and Cosine Waves with Frequencies %d Hz and %d Hz', f1, f2));
end

运行此函数后，它将返回三个向量 `y`, `y1`, 和 `y2`，并在当前工作空间显示相应的图形。