function[sqnr,a_quan,code]=u_pcm(a,n)
%u_pcm 一个序列的均匀PCM编码
%[sqnr,a_quan,code]=u_pcm(a,n)
%a=输入序列
%n=量化级数的数目(偶数)
%sqnr=输入信号量化噪声比(以dB为单位)
%a_quan=编码前的量化输出
%code=编码后的输出
amax=max(abs(a));
a_quan=a/amax;%归一化
b_quan=a_quan;%
d=2/n;%量化间隔
q=d.*[0:n-1];
q=q-((n-1)/2)*d;%量化电平
for i=1:n%
a_quan(find((q(i)-d/2<=a_quan)&(a_quan<=q(i)+d/2)))=...
q(i).*ones(1,length(find((q(i)-d/2<=a_quan)&(a_quan<=q(i)+d/2))));
b_quan(find(a_quan==q(i))) =(i-1).*ones(1,length(find(a_quan==q(i))));
end
a_quan=a_quan*amax;%恢复原值
nu=ceil(log2(n));%
code=zeros(length(a),nu);
for i=1:length(a)
for j=nu:-1:0
if (fix(b_quan(i)/(2^j))==1)
code(i,(nu-j))=1;%
b_quan(i)=b_quan(i)-2^j;%
end
end
end
sqnr=20*log10(norm(a)/norm(a-a_quan));%求量化信噪比
模拟信号的数字化传输_实现均匀和非均匀量化器_matlab源码
版权申诉
5星 · 超过95%的资源 83 浏览量
2022-03-22
22:10:50
上传
评论 1
收藏 3KB RAR 举报
阿里matlab建模师
- 粉丝: 3223
- 资源: 2782
最新资源
- 基于STM32F103C8T6单片机蓄电池在线监测系统主板硬件(原理图+PCB)工程文件.zip
- mysql大纲资料.txt
- c++大纲资料.txt
- 效率工具bat脚本实现日志提取
- MyBatis 中动态 SQL 的示例
- STM8L101F3P6单片机+CC1100模块433M遥控器设计硬件(原理图+PCB)工程文件.zip
- 上传下载铁人下载系统 Liuxing 1.0-liuxing1.0.rar
- 南京邮电大学数学实验实力雄厚,凭借其优秀的师资力量、丰富的实践教学资源和卓越的科研成果,成为国内一流的数学实验教学和科研基地
- 【火爆朋友圈的今天吃什么源码 v1.0】随机的为用户带来每一天的用餐选择和推荐.rar
- MPU6050中文版数据手册
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
- 1
- 2
- 3
- 4
前往页