数字频带传输系统

1、数字调制系统的基本结构调制原理调制原理二进制符号序列：其中，2ASK信号：2ASK信号产生方法信号产生方法F=200; %数字信号的带宽fc=800; %正弦载波信号的频率Tc=1/fc; %正弦载波信号的周期Q=fc/F; %频率比M=500; %一个正弦周期内的采样点数t=0:Tc/M:Q/F; %单码元周期时间 carry_cyc=cos(2*pi*fc*t); %单码元周期内的载波信号m=randint(1,10); %基带序列dsig=kron(m,ones(size(carry_cyc); %数字信号carry=kron(ones(size(m),carry_cyc); %载波信号

2、ask=kron(m,carry_cyc); %已调信号00.020.040.060.080.10.120.140.160.180.200.51数 字 信 号00.020.040.060.080.10.120.140.160.180.2-101载 波 波 形00.020.040.060.080.10.120.140.160.180.2-101经 过 2ASK调 制 后 的 波 形Figure & Code解调原理解调原理ask_mul_coh=ask.*carry; %已调信号与相干载波相乘time=0:Tc/M:(length(ask)/M)*Tc; %时间time=time(1:end-1

3、);%获得 Butterworth 模拟低通原型滤波器的阶数及3dB 截止频率 Ord,omega_c=buttord(4*pi*fc*0.6,4*pi*fc*0.8,2,30,s); %原型滤波器向实际滤波器转换，获得滤波器的分子，分母系数num,den=butter(Ord,omega_c,s); h=tf(num,den); %获得滤波器传递函数滤x=lsim(h,ask_mul_coh,time); %运用模拟滤波器对信号进行滤波th=0.25; %抽样判决的阈值设置 %抽样判决点的选取t_judge=(0:length(m)-1)*length(carry_cyc)+round(le

4、ngth(carry_cyc)/2); y=(x(t_judge); %抽样判决时刻时的信号值y_judge=1*(y=th)+0*(y=st2(i) st(i)=1; elsest(i)=0;endend012345678910-2-1012经 过 信 道 后 的 2fsk波 形幅度的大小t012345678910-2-1012经 过 带 通 滤 波 器 f1后 的 波 形幅度012345678910-2-1012经 过 带 通 滤 波 器 f2后 的 波 形幅度tFigure 01234567891000.20.40.60.811.21.4经 过 相 乘 器 h1后 的 波 形幅度0123

5、4567891000.20.40.60.811.21.4经 过 相 乘 器 h2后 的 波 形幅度t012345678910-0.100.10.20.30.40.50.6经 过 低 通 滤 波 器 sw1后 的 波 形幅度012345678910-0.100.10.20.30.40.50.6经 过 低 通 滤 波 器 sw2后 的 波 形幅度Figure 01234567891000.20.40.60.811.2经 过 抽 样 判 决 器 后 解 调 出 的 波 形幅度01234567891000.20.40.60.81原 始 数 字 序 列 的 波 形幅度t % 2ASK/2FSK解调误码率 snr=1:0.1:100; %信噪比 DBsnr=10*log10(snr); %归一化信噪比E/N %ASK Peno_Ask=0.5*exp(-snr/4); %非相干解调 Pehe_Ask=0.5*erfc(sqrt(snr)/2); %相干解调 %FSK Peno_Fsk=0.5*exp(-snr/2); %非相干解调 Pehe_Fsk=0.5*erfc(sqrt(snr)/2); %相干解调 M Code二进制数字调制系统误码率公式表M Figure0246810121416182000.0