用MATLAB程序开发设计PCM编码调制系统教学总结.docx
《用MATLAB程序开发设计PCM编码调制系统教学总结.docx》由会员分享,可在线阅读,更多相关《用MATLAB程序开发设计PCM编码调制系统教学总结.docx(19页珍藏版)》请在咨信网上搜索。
1、用MATLAB程序开发设计PCM编码调制系统目录一、作业要求1二、PCM编解码原理简述1三、MATLAB程序及仿真结果6一、作业要求用MATLAB程序开发设计PCM编码调制系统。系统参数设置如下:1 输入模拟信源的最高频率限制在4KHz以内(可以是正弦信号,也可随机产生信号);2 信号平均功率的动态范围30dB(自行设定幅度区间)。要求完成以下工作:1 采用均匀量化,若要求信噪比不应低于25dB,不考虑线路衰减和损耗,则设计该均匀量化器,并显示原始输入信号,量化输出信号,以及编码结果。2 若采用A律13折线PCM编码,显示原始输入信号图及解码后的信号图;打印输出编码结果,并比较译码后的量化误差
2、。二、PCM编解码原理简述均匀量化的基本原理在脉冲编码调制中,模拟信号首先以高于奈奎斯特的速率采样,然后将所的样本量化。假设模拟信号是以-max,max表示的区间内分布的,而量化电平数很大。量化电平可以是相等的或是不相等的;前者就属于均匀PCM,而后者就是非均匀PCM。关于量化的几个基本概念,量化间隔;量化误差;量化信噪比。 (1)相邻量化电平间距离称量化间隔, 用“”表示。(2)设抽样值为,量化后的值为, xq(kTs)与x(kTs)的误差称为量化误差,又称为量化噪声;量化误差不超过/2,而量化级数目越多,值越小,量化误差也越小。(3)衡量量化的性能好坏最常用指标是量化信噪比(Sq/Nq),
3、其中Sq表示量化信号值xq(kTs)产生的功率,Nq表示量化误差功率,量化信噪比越大,则量化性能越好。在均匀PCM中,长度为2Xmax的区间-max,max被划分为N个相等的子区间,每一子区间长度为=2Xmax/N。如果N足够大,那么在每一子区间内输入的密度函数就能认为是均匀的,产生的失真为D=2/12。如果N是2的幂次方即,那么就要求用比特来表示每个量化电平。这就意味着,如果模拟信号的带宽是,采样又是在奈奎斯特率下完成的,那么传输PCM信号所要求的带宽至少是(实际上1.5比较接近于实际)。这时失真由下式给出,如果模拟信号的功率用表示,则信号/量化噪声的比(SQNR)由下式给出式中表示归一化输
4、入,定义为以分贝(dB)计的SQNR为量化以后,这些已量化的电平用比特对每个已量化电平进行编码.编码通常使用自然二进制码(NBC),即最低电平映射为全0序列,最高电平映射为全1序列,全部其余的电平按已量化值的递增次序映射。非均匀量化的基本原理正变换:,其中是归一化输入(),是一个参数,在标准律的非线性中它等于255。反变换:在非均匀量化PCM中,输入信号首先通过一非线性环节以减小输入的动态范围,再将输出加到某一均匀PCM系统上。在接收端,输出再通过另一非线形环节,该环节是在发送端所用的非线性环节的逆特征。这样,总的效果就等效于一个在量化电平之间具有非均匀间隔的PCM系统。非均匀量化时,量化器随
5、着输入信号的大小采用不同的量化间隔,大信号时采用大的量化间隔,小信号时采用小的量化间隔,可以以较少的量化电平数达到输入动态范围的要求一般对语音信号传输来说,所使用的非线性可以是律的非线性,或是A律的非线性,中国和欧洲采用A率压缩特性(A=87.56),北美和日本采用律压缩特性(=255),压缩特性分别如下:律:A律:式中,x为归一化输入,y为归一化输出,A、为压缩系数,上式也可以表示为,由于A律在工程上不好实现,所以我们经常用近似的13折线压缩法去代替A率压缩下表列出了计算值与13折线时的值的比较。计算值与A律13折线时值的比较0101按折线分段时的01段落12345678斜率16168421
6、表中第二行的值是根据时计算得到的,第三行的值是13折线分段时的值。可见,13折线各段落的分界点与曲线十分逼近,同时按2的幂次分割有利于数字化。A律压扩特性是连续曲线,A律不同压扩特性也不同,在电路上实现这样的函数规律是相当复杂的。实际中,往往采用近似于A律函数规律的13折线()的压扩特性。这样,它基本上保持了连续压扩特性曲线的优点,。本设计中所用到的PCM编码正是采用这种压扩特性来编码的。在13折线法中,无论输入信号是正是负,均按8段折线(8个段落)进行编码。若用8位折叠二进制码来表示输入信号的抽样量化值,其中用第一位表示量化值的极性,其余七位(第二位至第八位)则表示抽样量化值的绝对大小。具体
7、的做法是:用第二至第四位表示段落码,它的8种可能状态来分别代表8个段落的起点电平。其它四位表示段内码,它的16种可能状态来分别代表每一段落的16个均匀划分的量化级。这样处理的结果,8个段落被划分成128个量化级。段落码和8个段落之间的关系如表1所示;段内码与16个量化级之间的关系见表2所示。表1 段落码 表2 段内码段落序号段落码段落范围量化间隔段内码量化间隔段内码81111024-2048151111701117110512-1024141110601106101256-512131101501015100128-25612110040100401164-128111011300113010
8、32-6410101020010200116-32910011000110000-168100000000在13折线中,用8位的折叠二进制码表示信号量化值的具体步骤为:用第2到4位表示段落码,8个段落的起点电平由它的8种可能状态来分别表示。其他四位表示段内码,每一个段落它的16个均匀的划分量化级由它的16种可能状态来分别表示。这样就使得8个段落被划分为128个量化级。再加上负的,相当于一共有256种量化电平数。解码时,码字电平=段落起始电平+8M5+4M6+2M7+M8*k 解码电平=码字电平+k2三、MATLAB程序及仿真结果注:1、 为便于阅读,本节仅仅给出核心代码,详细代码请见本文附录。
9、2、 文中的图片若不够清晰,子文件夹中还有高清输出的图片(600ppi)1、 生成题中要求的动态信号题中要求输入模拟信源的最高频率限制在4KHz以内(可以是正弦信号,也可随机产生信号);信号平均功率的动态范围30dB(自行设定幅度区间)。为了能够适应不同条件,本文生成的是频率范围在04k之间变化,动态范围在3035dB之间的随机信号。fana=104;%analog freqT=1/fana;%analog intervalSignalLength=64;tana=(0:SignalLength-1)*T;RealDynamicRange=0;while(RealDynamicRange=35
10、) while(RealDynamicRange=35) signal=randi(-300,300,1,SignalLength)/10;%raise accuracy idx=find(signal=0); % find all 0 signal(idx)=1; % set 0.1 to these indexes RealDynamicRange=20*log10(max(abs(signal)/ min(abs(signal); end%check freq range-NFFT=2nextpow2(SignalLength);ini_fft=fft(signal,NFFT)/NFFT
11、*2;%show reality amplitudef=fana/NFFT*(0:1:NFFT-1);abs_ini_fft=abs(ini_fft);%limit band below 4kk_stop_freq=floor(4000*NFFT*T);ini_fft(:,k_stop_freq:NFFT-k_stop_freq)=0;%lpfinisignal=real(ifft(ini_fft)*NFFT/2);%get the input signalRealDynamicRange=20*log10(max(abs(inisignal)/ min(abs(inisignal);end%
12、check freq range againinisignal_fft=fft(inisignal,NFFT)/NFFT*2;%show reality amplitudeabs_inisignal_fft=abs(inisignal_fft);原始产生的随机信号,按照题目要求处理过之后的信号以及频谱如下图所示:分别为原始随机信号、原始随机信号的带宽;处理之后带限在04kHz之内的信号,处理之后的信号的带宽。生成信号的动态范围RealDynamicRange本次随机产生结果为34.8dB。2、 信号采样本节对上述产生的信号进行采样,采样的频率为10000Hz。%sampling%fsam=10
13、5;%sampling freqTsam=1/fsam;tsam=0:Tsam:max(tana);sigsampling=interp1(tana,inisignal,tsam,linear);subplot(2,1,1);plot(tana,inisignal);title(analog signal);xlabel(t/s);ylabel(f(t);subplot(2,1,2);stem(tsam,sigsampling,Marker,none);title(sampling signal);xlabel(t/s);ylabel(fs(t);print -r600 -djpeg Anal
14、ogAndSamplingSignal;采样前后信号为3、信号的均匀量化%uniform quantization%N_uniquan=round(25-20*log10(min(abs(inisignal)/max(abs(inisignal)-4.77)/6.02);M_uniquan=2N_uniquan;sig_uniquan=quant(inisignal,M_uniquan,max(abs(inisignal);plot(tana,inisignal,b);hold on;plot(tana,sig_uniquan,r);print -r600 -djpeg UniQuan;del
15、ta=2*(max(abs(sig_uniquan)/(M_uniquan);sig_uniquancode=round(sig_uniquan-min(sig_uniquan)/delta);sig_uniquancode_s=dec2bin(sig_uniquancode,N_uniquan);fid=fopen(sig_uniquancode_s.txt,w);for i=1:length(inisignal) fprintf(fid,%srn,sig_uniquancode_s(i,:);endfclose(all);量化之前的信号与量化之后恢复的信号分别为,本次信号均匀量化所需的量化
16、位数为9。量化信噪比小于25dB。均匀编码结果如附2所示。4、PCM编解码%A_PCM%UniTo1=inisignal/max(abs(inisignal);%jxm dlm dnmjxm=zeros(1,length(UniTo1);dlm=jxm;dnm=jxm;quannumdlm=jxm;for i=1:length(UniTo1) %jxm if(UniTo1(i)0 jxm(i)=0; else jxm(i)=1; end %dlm tmpsig=abs(UniTo1(i); quannumdlm(i)=round(tmpsig/(1/2048); startquanV=0,16
17、,32,64,128,256,512,1024; quaninterval=1,1,2,4,8,16,32,64; if(quannumdlm(i)=16&quannumdlm(i)=32&quannumdlm(i)=64&quannumdlm(i)=128&quannumdlm(i)=256&quannumdlm(i)=512&quannumdlm(i)1024) dlm(i)=6; else dlm(i)=7; end %dnm dnm(i)=floor(quannumdlm(i)-(startquanV(dlm(i)+1)/quaninterval(dlm(i)+1); if(dnm(i
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB 程序 开发 设计 PCM 编码 调制 系统 教学 总结
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【人****来】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【人****来】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。