PCM源码程序 - matlab数学建模 - 谷速源码
下载频道> 资源分类> matlab源码> 数学建模> PCM源码程序

标题：PCM源码程序
分享到：

所属分类: 数学建模 资源类型： 文件大小: 1.67 KB 上传时间: 2016-01-25 20:15:55 下载次数: 6 资源积分：1分 提 供 者: xiaopeng2 PCM源码程序
内容：
PCM源码程序，程序员在编程的过程中可以参考学习使用，希望对IT程序员有用，此源码程序简单易懂、方便阅读，有很好的学习价值！
close all
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%UNIVERSIDAD NACIONAL PEDRO RUIZ GALLO
%CURSO: TELECOMUNICAICONES II  SEMESTRE 2013-I
%ING. FRANCISCO SEGURA ALTAMIRANO
%LAMBAYEQUE 2013
%Simulacion de modulacion PCM con cuatizacion Uniforme de una se馻l con una
%frecuencias Fmax.
Fmax=100;   %Frecuencia de la Se馻l Analogica
Fs=2*Fmax;  %Frecuencia de muestreo
Ts=1/Fs;    %Periodo de muestreo
N=100;      %Numero de muestras
nb=5;       %Tama駉 de la palabra de la codificacion para cada muestra
L=2^nb-1      %Numero de niveles permitidos
x=randi(L,N,1); %Se genera las muestras aleatoria de la se馻l
xpcm=dec2bin(x)-48; %se codifica en binario natural
xpcm=xpcm(:); %se coloca la se馻l PCM en un vector
%Aca es necesario calcular la frecuencia de muestreo de la se馻l PCM
%resultante que esta relacionada con la Fs de la se馻l y el numero de bit y
%el numero de puntos que se usara para construir la representacion digital
%de la cadena PCM.
Np=20   %Cada pulso se conformara por 20 puntos.
Ts1=(Ts/nb)*(1/Np);
Fs1=1/Ts1;
%Con la informaci髇 de Fs1 y Ts1 y sabiendo el numero de muestras y la
%cantiadad de bits por muestra se realiza el barrido de tiempo
t=0:Ts1:N*nb*Np*Ts1-Ts1;
%Tambien para mostrar correctamente el espectro es necesario realizar el
%barrido de frecuencias necesarios asi
f=linspace(-Fs1/2,Fs1/2,length(t));
%Se usara codificacion NRZ para esta parte
pnrz=ones(1,Np);    %Pulso basico para transmision
%Realizamos la se馻lizacion 
senpcm=pnrz'*xpcm';
%conertimos a un vector
senpcm=senpcm(:);
%graficamos
subplot(321)
plot(t,senpcm)
title('Se馻l PCM en el tiempo')
%Podemos obtener mas informacion si  observamos esta se馻l en el dominio de
%la frecuecia
SENPCM=fft(senpcm)/length(senpcm);
%graficamos
subplot(322)
plot(f,fftshift(abs(SENPCM))/max(abs(SENPCM)))
title('Se馻l PCM en la frecuencia y filtro ideal')
%Segun nuestros calculos vtx=Fs*5=1000bps, entonces BWtx=vtx/2 (NRZ)
%BWtx=500Hz
hold on
fc_ideal=500;  %Frecuencia de corte ideal, Puede cambiarse para probar
H=(abs(f)<=fc_ideal);  %se construye un filtro ideal
plot(f,H,'r','linewidth',3);
legend(' PCM','Filtro Ideal')
SENPCM_FIL=fftshift(H').*(SENPCM);  %se realiaza el fitrado de la se馻l
subplot(323)
plot(f,fftshift(abs(SENPCM_FIL))/max(abs(SENPCM_FIL)));
title(' PCM filtrada en la frecuencia')
%Calculamos la transformada inversa y observamos la se馻l PCM
senpcm_fil=ifft(SENPCM_FIL)*length(SENPCM_FIL);
subplot(324)
plot(t,real(senpcm_fil))
title(' PCM filtrada en el tiempo')
hold on
plot(t,senpcm,'r','linewidth',2)
legend(' filtrada',' Original')
%Ahora consideraremos el caso mas real al usar un filtro RC con frecuencia
%de corte igual a la del ancho de banda de transmission calculado
fcorte=500;
Hrc=1./(1+j*f/fcorte);   %Funcion de transferencia del filtro
SENPCM_FILrc=fftshift(Hrc').*(SENPCM);
subplot(325)
plot(f,fftshift(abs(SENPCM_FILrc)));
title(' PCM filtrada con filtro RC')
senpcm_filrc=ifft(SENPCM_FILrc)*length(SENPCM_FILrc);
subplot(326)
plot(t,real(senpcm_filrc))
title(' PCM filtrada en el tiempo')
hold on
plot(t,senpcm,'r','linewidth',2)
legend(' filtrada',' Original')
 

文件列表（点击上边下载按钮，如果是垃圾文件请在下面评价差评或者投诉）:

PCM源码程序/
PCM源码程序/modpcm.m

关键词： 源码 程序

相关推荐

matlab源码下载排行

Top_arrow
回到顶部

友情链接

matlab源码 网站源码 编程语言 网络编程 手机通讯编程 KUNAG 官网 KUNAG 电商 更多
联系方式| 版权声明| 招聘信息| 广告服务| 银行汇款| 法律顾问| 兼职技术| 付款方式| 关于我们|
网站客服网站客服 程序员兼职招聘 程序员兼职招聘
沪ICP备19040327号-3
公安备案号：沪公网安备 31011802003874号
库纳格流体控制系统（上海）有限公司 版权所有
Copyright © 1999-2014, GUSUCODE.COM, All Rights Reserved