sin.m

正弦波信号的均匀量化和线性编码 A/D 变换

close all
clear
b=input('请输入线性编码位数')
if b==[],b=8;end
m=2^(b-1)-1;
pause
disp('最大量化电平') %显示量化电平
disp(m)
pause
de=input('取样精度')
if de==[],de=0.1;end
pause
figure(1)
set(1,'Position',[60,60,290,350]) %设置图形的大小
figure(2)
set(2,'Position',[360,60,400,350])
x=(-pi:de:pi); %sinx 的自变量范围
y=m*sin(x); %输入电压
ym=max(y); %将输入最大值ym 转换为2^(b-1)-1
y11=y./ym %归一化
y12=y11.*(2^(b-1)-1); %反归一化
y1=round(y12); %求量化输出电压y
a=length(y);
aa=zeros(a,b);
for i=1:a
c=adc(b,y1(i));
aa(i,:)=c; %纪录量化编码
y2=dac(b,c);
Y(i)=y2./(2^(b-1)-1).*ym
end
figure(1)
plot(x,y,'m')
ylabel('电压');
hold on
stairs(x,Y) %画出阶梯波形
通信原理（软件）实验报告 2005.4 page 32 of 37
title('均匀量化过程')
grid on
hold off
figure(2)
t=y-Y;
plot(y,t)
xlabel('输入电压');ylabel('误差');
title('均匀量化误差分析')
grid on
aa %输出量化编码