当前位置：首页 > news >正文

余弦信号的均匀量化和μ律非均匀量化

news 2026/6/10 16:19:33

基于MATLAB实现的余弦信号的均匀量化和μ律非均匀量化的代码。代码包括信号生成、量化过程以及量化误差的计算和可视化。

1. 参数设置

% 信号参数
fs = 1000; % 采样频率 (Hz)
T = 1; % 信号时长 (秒)
t = 0:1/fs:T-1/fs; % 时间向量
f = 5; % 余弦信号频率 (Hz)
A = 1; % 信号振幅% 量化参数
numBits = 8; % 量化位数
mu = 255; % μ律参数

2. 生成余弦信号

% 生成余弦信号
x = A * cos(2 * pi * f * t);

3. 均匀量化

% 均匀量化
L = 2^numBits; % 量化级数
xMax = max(x);
xMin = min(x);
delta = (xMax - xMin) / L; % 量化步长
xQuantizedUniform = floor((x - xMin) / delta) * delta + xMin + delta / 2;

4. μ律非均匀量化

% μ律非均匀量化
% μ律压缩
xCompressed = sign(x) * log(1 + mu * abs(x)) / log(1 + mu);% 均匀量化
xCompressedQuantized = floor((xCompressed - min(xCompressed)) / delta) * delta + min(xCompressed) + delta / 2;% μ律扩张
xQuantizedMuLaw = sign(xCompressedQuantized) * (exp(abs(xCompressedQuantized) * log(1 + mu)) - 1) / mu;

5. 量化误差计算

% 计算量化误差
errorUniform = x - xQuantizedUniform;
errorMuLaw = x - xQuantizedMuLaw;% 计算均方误差 (MSE)
mseUniform = mean(errorUniform.^2);
mseMuLaw = mean(errorMuLaw.^2);

6. 可视化结果

% 绘制原始信号和量化信号
figure;
subplot(3, 1, 1);
plot(t, x);
title('原始余弦信号');
xlabel('时间 (s)');
ylabel('幅度');subplot(3, 1, 2);
plot(t, xQuantizedUniform);
title('均匀量化信号');
xlabel('时间 (s)');
ylabel('幅度');subplot(3, 1, 3);
plot(t, xQuantizedMuLaw);
title('μ律非均匀量化信号');
xlabel('时间 (s)');
ylabel('幅度');% 绘制量化误差
figure;
subplot(2, 1, 1);
plot(t, errorUniform);
title('均匀量化误差');
xlabel('时间 (s)');
ylabel('误差');subplot(2, 1, 2);
plot(t, errorMuLaw);
title('μ律非均匀量化误差');
xlabel('时间 (s)');
ylabel('误差');% 显示均方误差
disp(['均匀量化均方误差 (MSE): ', num2str(mseUniform)]);
disp(['μ律非均匀量化均方误差 (MSE): ', num2str(mseMuLaw)]);

参考代码实现对一个余弦信号的均匀量化和μ律非均匀量化 www.youwenfan.com/contentcng/81290.html