网站大量收购闲置独家精品文档,联系QQ:2885784924

数字频带系统matlab程序.doc

  1. 1、本文档共10页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数字频带系统matlab程序

实验八 clear all; close all; echo on; t0=3;%信号时长 dt1=0.001 %系统时域采样间隔 fs=1/dt1;%系统采样频率 df=0.2; %所需的频率分辨率 t=0:dt1:t0; % 产生模拟调制信号 m=[ones(1,t0/(3*dt1)),-2*ones(1,t0/(3*dt1)),zeros(1,t0/(3*dt1)+1)]; %画出调制信号波形及频谱 clf figure(1); subplot(121); plot(t,m(1:length(t)));% 画出调制信号波形 R=4.2; axis([0 t0 -R/2 R/2]); xlabel(t); ylabel(调制信号); subplot(122); [M,m,df1,f]=T2F4(m,dt1,df,fs);%求出调制信号频谱 plot(f,fftshift(abs(M)))% 画出调制信号频谱 xlabel(f); ylabel(调制信号频谱); axis([-5,5,0,2.5]); %A/D %低通抽样定理 t=0:dt1:t0; xt=m; [f,xf]=T2F(t,xt); dt=0.01; %抽样信号,抽样频率为fs fs=10; sdt=1/fs; st=[]; for t1=sdt:sdt:t0+sdt st=[st m(ceil(t1/dt1))]; end t1=0:sdt:t0; [f1,sf]=T2F(t1,st); figure(2); subplot(211); t=0:dt1:t0; plot(t,xt(1:length(t))); axis([0 t0 -2.5 2.5]); title(原始信号); subplot(212); stem(t1,st); title(抽样信号); axis([0 t0 -2.5 2.5]); %基带信号形成器 figure(3); subplot(211); plot(t,xt(1:length(t))); title(原始模拟信号); axis([0 t0 -2.5 2.5]); out1=pcm_encode(st); out=2*out1-1;%将编码结果由单极性不归0转换为双极性不归0码 dt=0.1; y=[]; for k=1:length(out) y=[y ones(1,1/dt)*out(k)]; end subplot(212); t=0:t0/(length(y)-1):t0; plot(t,y); xlabel(t); title(PCM编码结果); axis([0 t0 -1.5 1.5]); A=1; %载波振幅 fc=2; %载波频率(Hz) snr=15; %信噪比dB N_sample=8;% 基带信号中每个码元的的采样点数 N=length(out1); % 码元数为编码码元数 Ts=1; % 码元宽度 df=0.01;%频率分辨率 B=1/Ts; f_start=fc-B; f_cutoff=fc+B; fs=fc*N_sample%系统采样频率,即考虑载波后,一个码元内的采样点数 ts=Ts/fs; % 系统采样间隔 t=0:ts:N*Ts-ts; Lt=length(t); d=out; dd=sigexpand(out1,fc*N_sample); gt=ones(1,fc*N_sample); % NRZ波形 d_NRZ=conv(dd,gt); d_NRZ1=d_NRZ(1:Lt); d_NRZ1=mabianhuan(d_NRZ1,fs);%进行2DPSK调制的码变换,获得相对码 figure(4); out=2*d_NRZ1-1;%将相对码由单极性不归0转换为双极性不归0码 y=[]; dt=0.1; for k=1:length(out) y=[y ones(1,1/dt)*out(k)]; end subplot(211); t=0:t0/(length(y)-1):t0; plot(t,y); xlabel(t); title(相对码); axis([0 t0 -1.5 1.5]); d_sjx=2*d_NRZ1-1;%生成的双极性NRZ信号 d_sjx1=d_sjx(1:Lt); subplot(212); [d_sjx1f,d_sjx1,df1,f]=T2F4(d_sjx1,ts,df,fs);%求出双极性NRZ信号功率谱 plot(f,10*log10(abs(fftshift(d_sjx1f).^2/length(f))));% 画出双极性NRZ信号功率谱 axis([-3*B 3*B -100

文档评论(0)

dajuhyy + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档