% 生成时间序列

fs = 5000;      % [Hz] 信号采样频率

T = 1;          % [s] 信号长度

x = 0:1/fs:T;   % [s] 时间序列

% 生成信号序列

f = 440;        % [Hz] 信号频率

y = 1*sin(2*pi*f*x); 

% 输出音频文件

fname = 'Asin1.wav';        % 设定文件名称 注意格式

audiowrite(fname,y,fs);     % 输出文件

% 音频文件测试

clear y Fs                  % 初始化工作区

[y,Fs] = audioread(fname);  % 读取音频文件

sound(y,Fs);                % 收听音频

% 生成时间序列

fs = 12000;      % [Hz]	信号采样频率

T = 1;           % [s]	信号长度

x = 0:1/fs:T;    % [s]  时间序列

% 生成信号序列

f = 440;         % [Hz]  信号频率

A = [0.5 0.7 1 2 5 10];             % 设了6个信号幅度

y1 = zeros(6,length(x));            % 未截断信号初始化

for n = 1:6

    y1(n,:) = A(n)*sin(2*pi*f*x);   % 未截断信号赋值

end

y2 = y1;                            % 截断信号初始化

for n = 1:6

    for m = 1:length(x)             % 截断信号赋值

        if y1(n,m)>1

            y2(n,m) = 1;

        else if y1(n,m)<-1

                y2(n,m) = -1;

            end

        end

    end

end

% 测试控制

k = 5;                      % 选择第k组进行听音比较，k取1至6

y = [y1(k,:),y2(k,:)];      % 生成目标信号

% 输出音频文件

fname = 'Asin2.wav';        % 设定文件名称 注意格式

audiowrite(fname,y,fs);     % 输出文件

% 音频文件测试

clear y Fs                  % 初始化工作区

[y,Fs] = audioread(fname);  % 读取音频文件

sound(y,Fs);                % 收听音频