通信系統(tǒng)課群綜合訓(xùn)練與設(shè)計(jì)(報(bào)告參考模板)-資料下載頁
2025-04-14 01:50本頁面
【正文】 else code(i,1)=0。 endif abs(A(i))=0amp。amp。abs(A(i))16 code(i,2)=0。code(i,3)=0。code(i,4)=0。step=1。start=0。elseif abs(A(i))=16amp。amp。abs(A(i))32 code(i,2)=0。code(i,3)=0。code(i,4)=1。step=1。start=16。 elseif abs(A(i))=32amp。amp。abs(A(i))64 code(i,2)=0。code(i,3)=1。code(i,4)=0。step=2。start=32。 elseif abs(A(i))=64amp。amp。abs(A(i))128 code(i,2)=0。code(i,3)=1。code(i,4)=1。step=3。start=64。 elseif abs(A(i))=128amp。amp。abs(A(i))256 code(i,2)=1。code(i,3)=0。code(i,4)=0。step=4。start=128。 elseif abs(A(i))=256amp。amp。abs(A(i))512 code(i,2)=1。code(i,3)=0。code(i,4)=1。step=5。start=256。 elseif abs(A(i))=512amp。amp。abs(A(i))1024 code(i,2)=1。code(i,3)=1。code(i,4)=0。step=6。start=512。 elseif abs(A(i))=1024amp。amp。abs(A(i))2048 code(i,2)=1。code(i,3)=1。code(i,4)=1。step=7。start=1024。 endB=floor((abs(A(i))start)/step)。%循環(huán)碼編碼m=3。n=7。k=4。p=cyclpoly(n,k)。code=encode(code,n,k,39。cyclic39。,p)%HBB3編碼ami=code。num=0。for k=1:length(ami) if ami(k)==1 num=num+1。 if num/2==fix(num/2) ami(k)=1。 else ami(k)=1。 end endendnum=0。hd=ami。sign=0。V=zeros(1,length(ami))。B=zeros(1,length(ami))。for k=1:length(ami) if ami(k)==0 num=num+1。 if num==4。 num=0。 hd(k)=1*hd(k4)。 V(k)=hd(k)。 if hd(k)==sign hd(k)=1*hd(k)。 hd(k3)=hd(k)。 B(k3)=hd(k)。 V(k)=hd(k)。 hd(k+1:length(ami))=1*hd(k+1:length(ami))。 end sign=hd(k)。 end else num=0。 endend%psk調(diào)制cp=[]。mod1=[]。f=2*2*pi。t=0:2*pi/99:2*pi。for n=1:length(hd) if hd(n)==0 A=zeros(1,100)。 elseif hd(n)==1 A=ones(1,100)。 elseif hd(n)==1 A=1.*ones(1,100)。 end cp=[cp A]。 c=cos(f*t)。 mod1=[mod1 c]。endfigure subplot(4,2,1)plot(cp)。grid onaxis([0 100*length(hd) 2 2])。title(39。二進(jìn)制信號(hào)序列39。)cm=[]。mod=[]。for n=1:length(hd) if hd(n)==0 B=ones(1,100)。 c=cos(f*t+pi/2)。 elseif hd(n)==1 B=ones(1,100)。 c=cos(f*t+pi)。 elseif hd(n)==1 B=ones(1,100)。 c=cos(f*t)。 end cm=[cm B]。 mod=[mod c]。endtiaozhi=cm.*mod。subplot(4,2,2)plot(tiaozhi)grid onaxis([0 2000 2 2])title(39。psk調(diào)制序列39。)subplot(4,2,3)plot(abs(fft(cp)))axis([0 100*length(hd) 0 400])。title(39。原信號(hào)頻譜39。)subplot(4,2,4)plot(abs(fft(tiaozhi)))axis([0 100*length(hd) 0 400])。title(39。調(diào)制信號(hào)頻譜39。)%瑞利衰落fd=0。%多普勒頻移為0 td=1/fs。chan=rayleighchan(td,fd)%瑞麗衰落信道psk1=filter(chan,tiaozhi)。%信號(hào)通過瑞麗衰落信道subplot(4,2,5)plot(abs(psk1))axis([0 1000 2 2])。title(39。調(diào)制信號(hào)通過衰落信道后幅度特性圖39。) %psk解調(diào)jiet=2*mod1.*psk1。subplot(4,2,6)plot(angle(jiet))。grid ontitle(39。相乘后信號(hào)波形39。)axis([0 1000 2 2])。%低通濾波器fp=500。fs=700。rp=3。rs=20。fn=11025。ws=fs/(fn/2)。wp=fp/(fn/2)。[n,wn]=buttord(wp,ws,rp,rs)。[b,a]=butter(n,wn)。jt=filter(b,a,jiet)。subplot(4,2,7)plot(abs(jt))grid onaxis([0 100*length(hd) 2 2])。title(39。經(jīng)過低通濾波器信號(hào)波形39。)for m=1:100*length(hd) if jt(m) jt(m)=1。 elseif =jt(m)=。 jt(m)=0。 elseif jt(m) jt(m)=1。 endendsubplot(4,2,8)plot(abs(jt))axis([0 100*length(hd) 2 2])。title(39。經(jīng)過抽樣判決后信號(hào)波形39。)jt1=zeros(1,length(hd))。for i=1:length(hd) if (jt((i1)*100+50)) jt1(i)=1。 elseif (jt((i1)*100+50)) jt1(i)=0。(jt((i1)*100+50))= jt1(i)=1。 elseif (jt((i1)*100+50))= jt1(i)=1。 endendjt1%PCM譯碼encode=(reshape(encode39。,8,41))39。l=size(encode,1)。a=[0,16,32,64,128,256,512,1024]。b=[1 1 2 4 8 16 32 64]。c=[0 :]。 for i=1:l x=encode(i,1)。 T=bin2dec(num2str(encode(i,(2:4))))+1。 Y=bin2dec(num2str(encode(i,(5:8))))。 if Y==0。 k(i)=a(T)/2048*。 else k(i)=(a(T)+b(T)*c(Y))/2048*。 end if x==0 s(i)=k(i)。elses(i)=k(i)。end end32
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1