数字频带传输系统.pptxVIP

数字频带传输系统 ——基于MATLAB仿真 数字调制系统的基本结构2ASK调制解调调制原理?二进制符号序列：其中，2ASK信号：=???2ASK信号产生方法2ASK调制解调Figure CodeF=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);%载波信号ask=kron(m,carry_cyc);%已调信号2ASK调制解调解调原理2ASK调制解调Figure Codeask_mul_coh=ask.*carry;%已调信号与相干载波相乘time=0:Tc/M:(length(ask)/M)*Tc;%时间time=time(1:end-1);%获得 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(length(carry_cyc)/2);y=(x(t_judge));%抽样判决时刻时的信号值y_judge=1*(y=th)+0*(y=th);%抽样判决信号值的0 阶保持 %抽样判决后的数字信号波形y_value=kron(y_judge,ones(size(carry_cyc))); n_tau=time+0.5/F;%抽样判决后的信号对应的时间2FSK调制解调调制原理?2FSK信号：2FSK调制解调Figure Codefs=2000; dt=1/fs;%采样频率、周期f1=50; f2=150;%两个信号的频率base=round(rand(1,10));%产生原始数字随机信号g1=base; g2=~base;%将原始数字信号反转与g1反向g1t=(ones(1,2000))*g1;%进行抽样 ga=g1t(:);%将数字序列变成列向量g2t=(ones(1,2000))*g2;gb=g2t(:);t=0:dt:10-dt;fsk1=ga.*cos(2*pi*f1.*t);%得到频率为f1的fsk1已调信号fsk2=gb.*cos(2*pi*f2.*t);%得到频率为f2的fsk2已调信号fsk=fsk1+fsk2;%已产生2FSK信号2FSK调制解调解调原理2FSK调制解调Figure Codesn=awgn(fsk,30); %产生的随机噪声%设置带通滤波器的参数b1=fir1(101,[45/1000 55/1000]); b2=fir1(101,[145/1000 155/1000]);%经过带通滤波器后的信号H1=filter(b1,1,sn); H2=filter(b2,1,sn); %相干载波相乘的信号sw1=H1.*H1; sw2=H2.*H2;%设置低通滤波器的参数 bn1=fir1(101,[2/1000 55/1000]); bn2=fir1(101,[2/1000 155/1000]);%滤波st1=filter(bn1,1,sw1);st2=filter(bn2,1,sw2);%抽样判决st=zeros(1,t);for i=1:length(t) if(st1(i)=st2(i))st(i)=1; else st(i)=0; endend2FSK调制解调Figure 2FSK调制解调Figure 抗噪声性能的比较M Code%% 2ASK/2FSK解调误码率snr=1:0.1:100;%信噪比DBsnr=10*log10(snr);%归一化信噪比E/N%ASKPeno_Ask=0.5*exp(-snr/4);%非相干解调Pehe_Ask=0.5*erfc(sqrt(snr)/2);%相干解调%FSKPeno_Fsk=0.5*exp(-snr/2);%非相干解调Pehe_Fsk=0.5*erfc(sqrt((snr)/2));%相干