基于MATLAB的数字频带通信系统的仿真 - 图文

北京邮电大学世纪学院毕业设计（论文）

s2_delay=[-ones(1,N_sample*fc/2) s2(1:end-N_sample*fc/2)]; st_oqpsk=s1.*cos(2*pi*fc*t)-s2_delay.*sin(2*pi*fc*t);

%经过带通后，在经过非线性电路 [f y1f]=T2F(t,st_qpsk); [f y2f]=T2F(t,st_oqpsk);

[t y1]=bpf(f,y1f,fc-1/Ts,fc+1/Ts); [t y2]=bpf(f,y2f,fc-1/Ts,fc+1/Ts); subplot(221);

plot(t,y1);xlabel('t');ylabel('QPSK波形');

axis([5 15 -1.6 1.6]);title('经过带通后的波形');

subplot(222);

plot(t,y2);xlabel('t');ylabel('OQPSK波形'); axis([5 15 -1.6 1.6]);title('经过带通后的波形');

%经过非线性电路 y1=1.5*tanh(2*y1); y2=1.5*tanh(2*y2); [f y1f]=T2F(t,y1); [f y2f]=T2F(t,y2);

py1f=py1f+abs(y1f).^2/T; %QPSK不同段信号功率谱密度相加 py2f=py2f+abs(y2f).^2/T; %OQPSK不同段信号功率谱密度相加 end

py1f=py1f/100 %QPSK100段功率谱密度平均 py2f=py2f/100 %OQPSK100段功率谱密度平均

subplot(223);

plot(f,10*log10(py1f));xlabel('f');ylabel('QPSK功率谱密度(dB/Hz)'); title('经过非线性电路后的功率谱密度');axis([-15 15 -30 10]);

subplot(224);

plot(f,10*log10(py2f));xlabel('f');ylabel('OQPSK功率谱密度(dB/Hz)'); title('经过非线性电路后的功率谱密度');axis([-15 15 -30 10]);

figure(2) x=-2:0.1:2;

y=1.5*tanh(2*x);

plot(x,y);title('非线性电路的输入输入函数')

38