通信原理課程設(shè)計(jì)1-展示頁
2024-10-25 03:24本頁面
【正文】 lot(f,10*log10(abs(st3f).^2/T))。)。xlabel(39。雙極性sinc波形39。grid axis([0 202 2])。)。xlabel(39。plot(f,10*log10(abs(st2f).^2/T))。)。xlabel(39。subplot(323)plot(t,[st2(1:length(t))])。單極性NRZ功率譜密度(dB/Hz)39。grid axis([5 540 10])。grid subplot(322)。單極性NRZ波形39。axis([0 ])。[f,st5f]=T2F(t,[st5(1:length(t))])。[f,st3f]=T2F(t,[st3(1:length(t))])。[f,st1f]=T2F(t,[st1(1:length(t))])。st4=conv(data,gt4)。% 變成雙極性序列data=sigexpand(d,N_sample)。%調(diào)用Matlab的卷積函數(shù)conv st2=conv(data,gt2)。gt5=[gt5 zeros(1,N_sample/2)]。%插入N_sample1個0 gt4=ones(1,N_sample)。%截段取10個碼元d=(sign(randn(1,N))+1)/2。%1x8 mt3=sinc((t5)/Ts)。%1x8 gt2=ones(1,N_sample/2)。%1x8000 T=t(end)。%抽樣時(shí)間間隔 N=1000。%每個碼元的長度N_sample=8。close all。out=reshape(out,1,M*N)。out=zeros(M,N)。sf=T/N*fftshift(sf)。f=N/2*df:df:N/2*dfdf。df=1/T。第一篇：通信原理課程設(shè)計(jì)1附錄Ⅰ：數(shù)字基帶信號的仿真代碼與結(jié)果 function[f,sf]=T2F(t,st)dt=t(2)t(1)。T=t(end)。N=length(st)。sf=fft(st)。end %在把序列d插入到序列M中,得到一個新序列 function [out]=sigexpand(d,M)N=length(d)。out(1,:)=d。%數(shù)字基帶信號的功率譜密度 clear all。Ts=1。%每個碼元的抽樣點(diǎn)數(shù) dt=Ts/N_sample。%碼元數(shù)t=0:dt:(N*N_sample1)*dt。gt1=ones(1,N_sample)。% 1x4?1x8? gt2=[gt2 zeros(1,N_sample/2)]。% 1x8000 sin(pi*t/Ts)/(pi*t/Ts)波形gt3=mt3(1:10*N_sample)。%1x1000 d=?? data=sigexpand(d,N_sample)。gt5=ones(1,N_sample/2)。st1=conv(data,gt1)。d=2*d1。st3=conv(data,gt3)。st5=conv(data,gt5)。[f,st2f]=T2F(t,[st2(1:length(t))])。[f,st4f]=T2F(t,[st4(1:length(t))])。figure(1)subplot(321)plot(t,[st1(1:length(t))])。xlabel(39。)。plot(f,10*log10(abs(st1f).^2/T))。xlabel(39。)。grid axis([0 ])。單極性RZ波形39。subplot(324)。grid axis([5 540 10])。單極性RZ功率譜密度(dB/Hz)39。subplot(325)plot(t5,[st3(1:length(t))])。ylabel(39。)。t/Ts39。subplot(326)。gridaxis([5 540 10])。sinc波形功率譜密度(dB/Hz)39。xlabel(39。)。axis([0 ])。雙極性NRZ波形39。grid subplot(222)。grid axis([5 540 10])。雙極性NRZ功率譜密度(dB/Hz)39。subplot(223)subplot(224)。grid plot(f,10*log10(abs(st5f).^2/T))。axis([5 540 10])。雙極性RZ波形39。xlabel(39。)。N=100。N_sample=8。Ts=1。dt=Ts/N_sample。ht1=gt。dd=rt1(N_sample:N_sample:end)。ddd=sigexpand(dd,N_sample)。stem(t,ddd(1:length(t))/8)。figure(1)xlabel(39。)。plot(t,st(1:length(t)))。axis([0 ])。xlabel(39。)。subplot(322)。t/Ts39。stem(t,a)。理想低通濾波后輸出39。axis([0 ])。xlabel(39。)。ddd=sigexpand(dd,N_sample)。stem(t5,ddd(1:length(t))/8)。axis([0 ])。xlabel(39。)。方波濾波后輸出39。ylabel(39。)。附錄Ⅲ： 升余弦滾降系統(tǒng)仿真代碼與結(jié)果%升余弦滾降系統(tǒng)示意圖 Ts=1。dt=Ts/N_sample。t=10*Ts:dt:10*Ts。alpha=[0,1]。elseif abs(f(k))Xf(n,k)=Ts。end end xt(n,:)=sinc(t/Ts).*(cos(alpha(n)*pi*t/Ts))./(14*alpha(n)^2*t.^2/Ts^2+eps)。axis([1 1 0 ])。f/Ts39。ylabel(39。)。axis([10 ])。t39。ylabel(39。)。N=length(S)。tmp2=N*ifft(S)/T。s(N/2+1:N)=tmp2(1:N/2)。endfunction S=t2f(s,fs)%s為輸入信號，S為s的頻譜，close all N=256。L=N/M。fs=N/T。fs為采樣率； N=length(s)。f=[N/2:(N/21)]/T。tmp2=N*ifft(s)/fs。S(N/2+1:N)=tmp1(1:N/2)。end %升余弦滾降的PAM信號 clear allt=[(T/2):1/fs:(T/21/fs)]。G=cos(pi*f/2).^2。a=12*(rand(1,M))。delta(L/2:L:N)=a*fs。s=real(f2t(S,fs))。T=N/fst=[(T/2):1/fs:(T/21/fs)]。
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1