数字频带系统matlab程序.doc

实验八 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