bianma_jiema.m

来自「PCM编码把模拟信号转换为数字编码」· M 代码 · 共 25 行

function [sqnr,Y]=bianma_jiema(a)
% 对编码并输出码值
% 画出输入——输出图
% 画出输入信号时域与频谱图
% 画出解码后的信号时域与频谱图
%  mu是的值
% n 是量化电平数

[f1,fx1]=sp_1(a);
wmax=max(abs(a));
[y,wmax]=mulaw(a,255) ;     %mulaw  非线性PCM的U律非线性
[sqnr,a_quan,code]=u_pcm(y,256);  %函数实现u率的均匀编码
[x,lmax]=fanmulaw(a_quan,255);
 [f2,fx2]=sp_1(x);
figure(1);
subplot(2,1,1)  ,plot(a,'b'); grid, ylabel('原信号时域图');            %原信号图
subplot(2,1,2)  ,plot(x,'g');   ylabel('解码时域图');           %解码后的复原图
figure(2);
subplot(2,1,1)  , plot(f1,abs(fx1));grid,xlabel('频率f');ylabel('|w(f)|'); title('原信号频谱图'); %原信号的频谱图
subplot(2,1,2)  ,plot(f2,abs(fx2));grid,xlabel('频率f');ylabel('|w(f)|'); title('反u率频谱图');  %解码后的频谱图
figure(3);
shuru_shuchu(a)