📄 chx3.m
字号:
%量化对开根号升余弦信号频谱的影响
L=64;
Rb=1000;
a=0.5;
t=-10/Rb:0.0003:10/Rb;
t1=cos(a*Rb*pi*t);
t2=sin(pi*Rb*t)./(pi*Rb*t);
t3=1-4*a*a*Rb*Rb*t.*t;
P=t2.*t1./(t3);
for N=2:14 %量化比特数
partition=[-1:1/2.^(N-1):1-1/2.^(N-1)]; %设定等距量化间隔
codebook=[-1:1/2.^(N-1):1]; %设定等差量化码数
[index,quants]=quantiz(P,partition,codebook); %进行量化
Y=fft(quants,L);
e(N-1)=sum(abs(fft(P,L)-Y).^2); %误差能量计算
end
Q=sqrt(Y);
M=ifft(Q,L);
subplot(2,2,1);
plot(M);
grid;
xlabel('Time');
ylabel('Amplitude');
subplot(2,2,2);
plot(quants,'.');
grid;
xlabel('Time');
ylabel('Amplitude');
subplot(2,2,3);
plot(abs(Q(1:L/2)));
grid;
xlabel('k');
ylabel('Amplitude');
subplot(2,2,4);
plot(e,'*');
grid;
xlabel('量化比特数');
ylabel('能量误差e');
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -