基于matlab的qpsk調制解調仿真設計與研究設計說明書-資料下載頁
2025-06-27 18:01本頁面
【正文】 波data_I = 2*bit_I+1。 % 將bit_I中的1變成1，0變成1。 注意data_I是500點data_Q = 2*bit_Q+1。 %將bit_Q中的1變成1，0變成1data_I1=repmat(data_I39。,20,1)。 %將500行的列向量data_I的共軛轉置data_I’復制為20*500的矩陣，20行數據是相同的。data_Q1=repmat(data_Q39。,20,1)。for i=1:1e4 %data_I2是將data_I1這個20*500的矩陣拉長為1*10000的行向量data_I2(i)=data_I1(i)。 data_Q2(i)=data_Q1(i)。end。f=0::1。xrc=+*cos(pi*f)。 data_I2_rc=conv(data_I2,xrc)/。 % data_I2_rc就是Idatadata_Q2_rc=conv(data_Q2,xrc)/。 % data_Q2_rc就是Qdataf1=1。t1=0::1e3+。 % 10010個數據，長度和data_I2_rc以及data_Q2_rc相同n0=rand(size(t1))。 %n0是1*10010的均勻分布的噪聲I_rc=data_I2_rc.*cos(2*pi*f1*t1)。Q_rc=data_Q2_rc.*sin(2*pi*f1*t1)。QPSK_rc=(sqrt(1/2).*I_rc+sqrt(1/2).*Q_rc)。QPSK_rc_n0=QPSK_rc+n0。I_demo=QPSK_rc_n0.*cos(2*pi*f1*t1)。 % 解調(相干解調，與載波相乘)Q_demo=QPSK_rc_n0.*sin(2*pi*f1*t1)。I_recover=conv(I_demo,xrc)。 % 低通濾波 Q_recover=conv(Q_demo,xrc)。I=I_recover(11:10010)。 Q=Q_recover(11:10010)。t2=0::。 %t2有20000個點t3=0::。 %t3有10000個點data_recover=[]。 % 抽樣判決for i=1:20:10000 %data_recover是待判決的20000點data_recover=[data_recover I(i:1:i+19) Q(i:1:i+19)]。 end。bit_recover=[]。for i=1:20:20000 if sum(data_recover(i:i+19))0 %20點為同一個值，20點數據疊加后與閾值0比較 data_recover_a(i:i+19)=1。 %data_recover_a是并/串轉換后的20000點 bit_recover=[bit_recover 1]。 %bit_recover是1000點數據 else data_recover_a(i:i+19)=1。 bit_recover=[bit_recover 1]。 endenderror=0。dd = 2*bit_in+1。 % 將bit_in中的1變成1，0變成1ddd=[dd39。]。 %ddd是1表示0，1表示1的原始序列，1000個點ddd1=repmat(ddd,20,1)。 %ddd1是20*1000的矩陣for i=1:2e4 ddd2(i)=ddd1(i)。 %將ddd1拉直為1*20000的行向量ddd2endfor i=1:1e3 if bit_recover(i)~=ddd(i) error=error+1。 endendp=error/1000。figure(1)subplot(4,1,1)。plot(t2,ddd2)。axis([0 100 2 2])。title(39。原序列39。)。subplot(4,1,2)。plot(t1,I_demo)。axis([0 100 2 2])。title(39。I 支路解調39。)。subplot(4,1,3)。plot(t1,Q_demo)。axis([0 100 2 2])。title(39。Q 支路解調39。)。subplot(4,1,4)。plot(t2,data_recover_a)。axis([0 100 2 2])。title(39。解調后序列39。)。24
環(huán)評公示相關推薦
鄂ICP備17016276號-1