畢業(yè)論文msk調(diào)制解調(diào)系統(tǒng)的設(shè)計-資料下載頁
2025-06-20 13:05本頁面
【正文】 subplot(4,1,3)。 plot(tI,I,39。k39。)。axis([Ts,2,2])。 title(39。同相調(diào)制信號I39。)。 grid on subplot(4,1,4)。 plot(tQ,Q,39。k39。)。axis([Ts,2,2])。 title(39。正交調(diào)制信號Q39。)。 grid on figure subplot(3,1,1)。 plot(tI,U,39。k39。)。axis([Ts,2,2])。 title(39。Icos(wct)cos(wt)39。)。grid on subplot(3,1,2)。 plot(tQ,V,39。k39。)。axis([Ts,2,2])。 title(39。Qsin(wct)sin(wt)39。)。grid on subplot(3,1,3)。 plot(tI,W,39。k39。)。axis([Ts,2,2])。 title(39。MSK信號39。)。grid on noise=*randn(1,n)。 S=W+noise。 [S1,s,df1]=fftseq(S,ts,df)。 S1=S1/fs。 f=[0:df1:df1*(length(s)1)]fs/2。 n_cutoff=100。 H=zeros(size(f))。 for i=floor((fcn_cutoff)/df1):floor((fc+n_cutoff)/df1) H(i)=1。 end。 for i=(65536floor((fcn_cutoff)/df1)):1:(65536floor((fc+n_cutoff)/df1)), H(i)=1。 end。 DEM=S1.*H。 dem=real(ifft(DEM))*fs。 dm=dem(1:length(tQ))。 ds1=dm.*cos(2*pi*fc*tQ+pi*tQ/(2*Ts))。 ds2=dm.*cos(2*pi*fc*tQpi*tQ/(2*Ts))。 [DW1,ww,df1]=fftseq(ds1,ts,df)。 [DW2,ww,df1]=fftseq(ds2,ts,df)。 DW1=DW1/fs。 DW2=DW2/fs。 n_cutoff1=floor(200/df1)。 f1=[0:df1:df1*(length(ww)1)]fs/2。 H1=zeros(size(f1))。 H1(1:n_cutoff1)=2*ones(1,n_cutoff1)。 H1(length(f1)n_cutoff1+1:length(f1))=2*ones(1,n_cutoff1)。 DEM1=H1.*DW1。 DEM2=H1.*DW2。 s1=real(ifft(DEM1))*fs。 s11=s1(1:length(tQ))。 s2=real(ifft(DEM2))*fs。 s22=s2(1:length(tQ))。 for i=1:n if (s11(i)=)||(s11(i)=),sa(i)=1。 else sa(i)=1。 end。 end。 for i=1:n if (s22(i)=)||(s22(i)=),sb(i)=1。 else sb(i)=1。 end。 end。 for i=0:N1, u1(i+1)=sa(i*count+count/2)。 u2(i+1)=sb(i*count+count/2)。 end。 for j=1:N, if (mod(j,2)==0) out(j)=u2(j)。 else out(j)=u1(j)。 end。 end。 for i=1:n, OUT(i)=out(((i1)mod((i1),count))/count+1)。 end。 figure subplot(3,1,1)。 plot(tQ,S,39。k39。)。axis([0,2,2])。 title(39。 加噪聲后的信號S39。)。axis([0,2,2])。 grid on subplot(3,1,2)。 plot(tQ,dm,39。k39。)。 title(39。 進(jìn)過帶通后的信號dm39。)。axis([0,2,2])。 grid on subplot(3,1,3)。 plot(tQ,ds1,39。k39。)。 title(39。 乘以載波后的信號ds139。)。axis([0,2,2])。 grid on figure subplot(3,1,1)。 plot(tQ,ds2,39。k39。)。 title(39。 乘以載波后的信號ds239。)。axis([0,2,2])。 grid on subplot(3,1,2)。 plot(tQ,s11,39。k39。)。 title(39。 經(jīng)過低通后的信號s1139。)。axis([0,2,2])。 grid on subplot(3,1,3)。 plot(tQ,s22,39。k39。)。 title(39。 經(jīng)過低通后的信號s2239。)。axis([0,2,2])。 grid on figure subplot(3,1,1)。 plot(f,fftshift(abs(H)),39。k39。)。 title(39。帶通濾波器H39。)。grid on subplot(3,1,2)。 plot(f1,fftshift(abs(H1)),39。k39。)。 title(39。低通濾波器H139。)。grid on subplot(3,1,3)。 plot(tQ,sa,39。k39。)。 title(39。 經(jīng)過低通后的信號sa39。)。axis([0,2,2])。 grid on figure subplot(3,1,1)。 plot(tQ,sb,39。k39。)。 title(39。 經(jīng)過低通后的信號sb39。)。axis([0,2,2])。 grid on subplot(3,1,2)。 plot(tQ_R,R,39。k39。)。axis([0,1,2,2])。 title(39。調(diào)制信號39。)。grid on subplot(3,1,3)。 plot(tQ_R,OUT(count+1:n),39。k39。)。 title(39。解調(diào)信號39。)。axis([0,1,2,2])。 grid on count1=0。 for i=1:N1, if(r(i)~=out(i+1)) count1=count1+1。 end。 end。 Pe=count1/(N1)*100。 sprintf(39。MSK錯誤概率為：％g％％39。,Pe) %頻率間隔的圖 n=10000。 count=100。 N=n/count。 r=rand(1,N1)。 for i=1:N1, if (r(i)=) r(i)=1。 else r(i)=1。 end。 end。 r([1:10])=[1 1 1 1 1 1 1 1 1 1]。 for i=1:ncount, R(i)=r(((i1)mod((i1),count))/count+1)。 end。 a(1)=1。 for i=2:N if r(i1)==1 a(i)=a(i1)。 else a(i)=a(i1)。 end。 end。 for i=1:n A(i)=a(((i1)mod((i1),count))/count+1)。 end。 for i=1:N if(mod(i,2)) l(i)=a(i)。 l(i+1)=a(i)。 else q(i)=a(i)。 q(i1)=a(i)。 end。 end。 for i=1:n I(i)=l(((i1)mod((i1),count))/count+1)。 end。 for i=1:n Q(i)=q(((i1)mod((i1),count))/count+1)。 end。 ts=。 Ts=ts*count。 fs=1/ts。 fc=3/(2*Ts)。 t0=Ts*Nts。 f=1/(4*Ts)。 df=。 tI=[Ts:ts:t0Ts]。 tQ=[0:ts:t0]。 tQ_R=[0:ts:t0Ts]。 c=2*pi*f。 u=I.*cos(c*tI)。 v=Q.*sin(c*tQ)。 U=u.*cos(2*pi*fc*tI)。 V=v.*sin(2*pi*fc*tQ)。 for i=count+1:n W(i)=U(i)+V(icount)。 end。 figure subplot(3,1,1)。 plot(tQ_R,R,39。k39。)。axis([0,2,2])。 %title(39。數(shù)字信號39。)。 grid on subplot(3,1,2)。 plot(tI,W,39。k39。)。axis([0,2,2])。 %title(39。MSK信號39。)。 gri
規(guī)章制度相關(guān)推薦
鄂ICP備17016276號-1