【正文】 it(length(bit))]。w =300。 %角頻率for i=1:length(bit)/4 decode(i,[1:4])=bit((i1)*4 + 1 : i*4)。endbit2=xianxi(decode)。 %線性分組碼編碼 Sign = [bit2,sheng]。L=length (Sign )。 %信號長度st=zeros (1,m*L)。sign1=zeros (1,m*L)。t = 0:1/fs:Ts*L1/fs。for i=1:L if Sign(i)==1 sign1((i1)*m +1:i*m) = High。 else sign1((i1)*m +1:i*m) = Low。 endendfor i=1:L if Sign(i) == 1 st((i1)*m +1:i*m) =cos(2* pi*w *t((i1)*m + 1:i*m)+(pi/ 2))。 else st((i1)*m +1:i*m) =cos(2* pi*w *t((i1)*m + 1:i*m))。 endendst1=st。st1=st1+cos(3000*t)。figureplot(t, sign1)。axis( [0,Ts*(L+1),(1/2),3/2] )。title ( 39。編碼后的信號39。 )。gridfiguresubplot (2,1,1)。plot (t,st )。axis([0,Ts*(L+1), 3/ 2,3/2] )。title ( 39。數(shù)字調(diào)制后的信號39。 )。gridsubplot (2,1,2)。plot ( t, st1 )。axis( [0,Ts*(L+1)+1, 5/2,5/2])。title ( 39。加噪后的信號39。 )。griddt = st1.*cos(2*pi*w*t )。figureplot (t,dt)。axis( [0,Ts*(L+1),2,5/2] )。title ( 39。與相干波相乘后的波形39。 )。grid[N,Wn] = buttord( 2*pi*50, 2*pi*100,3,25,39。s39。)。 %臨界頻率采用角頻率表示[b,a]=butter(N,Wn,39。s39。)。[bz,az]=impinvar(b,a,fs)。 %映射為數(shù)字的dt = filter(bz,az,dt)。figureplot (t,dt)。axis( [0,Ts*(L+1), 3/2,3/2] )。title ( 39。低通濾波后的波形39。 )。 ss=dt。 %存放抽樣判決后的值bit1=[]。tt=length(ss)/L。for i=tt/2:tt:length(ss) if ss(i)(1/4) ss((itt/2+1):(i+tt/2))=zeros ( 1,3000)。 bit1=[bit1,0]。 else ss((itt/2+1):(i+tt/2))=ones ( 1,3000)。 bit1=[bit1,1]。 endendfigureplot (t,ss,39。r39。)。axis( [0,Ts*(L+1), 1/2,3/2] )。title ( 39。信道譯碼信號39。 )。gridbb=[]。 U1=rem(length(bit1),7)。 l=length(bit1)。sheng1=[bit1(length(bit1))U1+1:bit1(length(bit1))]。bit1([lU1+1:bit1(l)])=[]。bb=xianxiyima(bit1)。 %線性分組碼譯碼。bb=[bb,sheng1]。此函數(shù)的功能是將信源編碼后的0、1序列再經(jīng)過信道編碼（7，4）線性分組碼編碼，將信道編碼后的信號進行BPSK調(diào)制，由于考慮到在信道中傳輸是會有噪聲的影響，所以在仿真時我人為的加入一定的噪聲。在接收端進行相干解調(diào)，用一個與調(diào)制信號同頻的信號與接收到的信號（加入噪聲后的已調(diào)信號）相乘，再用一個數(shù)字濾波器進行濾波，去掉了信號中的高頻成份和噪聲，取出直流分量，得到方波信號，對這個信號進行抽樣判決將其轉(zhuǎn)換成0、1比特流后，再經(jīng)過信道譯碼。子函數(shù)5：function a3=yima(y,n)m=2/n。for i=1:n+1 q(i)=(1)+(i1)*m。end tt=[]。 tt(1)=1。 for x=2:length(q) tt(x)=(q(x1)+q(x))/2。 %存放第幾個電平段的電平值 end tt=[tt,1]。nu=ceil(log2(n))。decode=[]。for i=1:length(y)/(nu+1) decode(i,[1:(nu+1)])=y((i1)*(nu+1) + 1 : i*(nu+1))。endab=[] 。 % 存放的十進制數(shù)，到tt中去找對應(yīng)值 [w,z]=size(decode)。 for j=1:w L=nu。 value=0。 for i=0:L value=decode(j,i+1)*(2^(Li))+ value。 end ab(j)=value。 end decode。 w。 a3=[] 。 % 存放信源譯碼值 for i=1:w a3(i)=tt(ab(i))。 end此函數(shù)的功能是將接收到的0、1比特流分為八位二進制一組一組的，再將每八位二進制轉(zhuǎn)換為十進制數(shù)，再到每段中去尋到與發(fā)關(guān)端有相同約定段號值。即可恢復(fù)原信號的量化值。子函數(shù)6：function y5=Ayi(y,A)for i=1:length(y)if((abs(y(i))=0)amp。(abs(y(i))=(1/(1+log(A))))) y1(i)=y(i)*(1+log(A))/A。else if((abs(y(i))=1)amp。(abs(y(i))(1/(1+log(A))))) y1(i)=(1/A)*sign(y(i))*exp(abs(y(i))*(1+log(A))1)。 endendendy5=y1。end此函數(shù)的功能是將恢復(fù)的量化值通過A律反轉(zhuǎn)換公式得到對應(yīng)的值，是用13折線轉(zhuǎn)換的反過程。因為在當A＝，可以用13折線用13條折線近似A律的平滑曲線，在實際工程中A律的平滑曲線不易實現(xiàn)，而用13折線法可以降低工程實現(xiàn)中的難度。主函數(shù)：clcclear allt=0::6*pi。a1=sin(t)。 %模擬信源figure(1)plot(t,a1)。title ( 39。原始模擬信號39。 )。t1=0:2*pi/32:6*pi。 %抽樣a=sin(t1)。figureplot(t1,a)。stem(t1,a)。title ( 39。抽樣信號39。 )。y4=zhexian(a) 。 %13折線 n=256。 %量化電平數(shù)。bitstream=bianma(y4,n)。bit1=BPSK(bitstream)。y2=yima(bit1,n)。y5=Ayi(y2,)。figureplot(y5)。stem(y5,39。r39。)axis([0,100,1,])figureplot(y5,39。r39。)。 hold onplot(y5,39。.b39。)。axis([0,100,1,])主函數(shù)的功能是是產(chǎn)生模擬信號并且進行抽樣得到抽樣值，然后將各個子函數(shù)銜接到一塊，形成一個比較完整的通信系統(tǒng)，進行仿真畫出各個階段重要的圖形，以便觀察和分析。