这个用matlab怎么编程

% A律压缩函数
function y = A_compress(x, Jhta, rs, A)
% x: 输入信号向量
% Jhta: 压缩特性参数
% rs: 压缩信号与原始信号幅度比值
% A: 电信号振幅最大值
x = x / A; % 标准化幅度值
y = sign(x) .* (log(1+Jhta*abs(x))/log(1+Jhta)) * rs * A; % A律压缩公式
end

% 绘制A律13折线压缩特性曲线
function plot_A_law(Jhta, A)
% Jhta: 压缩特性参数
% A: 电信号振幅最大值
x = -A:0.01:A; % 选取输入信号范围
y = A_compress(x, Jhta, 1, A); % 压缩后的信号值
subplot(2,1,1);
plot(x, y, 'b'); grid on;
title('A律13折线压缩特性曲线');
xlabel('输入信号幅度'); ylabel('压缩后信号幅度');

% 仿真压缩后的信号
rs = 0.5; % 压缩后信号与原始信号幅度比值
x = linspace(-A, A, 500); % 生成输入信号
y = A_compress(x, Jhta, rs, A); % 仿真压缩后的信号
subplot(2,1,2);
plot(x, y, 'r'); grid on;
title(sprintf(['仿真A律压缩特性曲线\nJhta=%.2f, rs=%.2f'], Jhta, rs));
xlabel('输入信号幅度'); ylabel('压缩后信号幅度');
end

% 示例：绘制A律13折线压缩特性曲线并仿真
Jhta = 0.1; % 压缩参数
A = 87.6; % 电信号振幅最大值
plot_A_law(Jhta, A);