通信原理課程設(shè)計(jì)(編輯修改稿)
2025-02-09 19:19 本頁面
【文章內(nèi)容簡介】 N=1/dt。F=X(1:N/2+1)。 % F(k)=X(k)(k=1:N/2+1)f=f*(0:N/2)/N。 % 使頻率軸f從零開始figure(4)subplot(3,1,1),plot(f,abs(F)),title(39。dsb調(diào)制信號頻譜圖39。)。xlim([0,fc*2])。xlabel(39。Frequency39。)。ylabel(39。|F(k)|39。) 調(diào)制與解調(diào)的波形，頻譜（2）SSB方式 clear all。close all。clc。t0=。dt=。t=[0:dt:1]。fc=phi0=0。ct=cos(2*pi*fc*t+phi0)。 %載波信號fm=10。fun = @(t)1.*(t=0amp。t=t0/3)+(2).*(t=t0/3amp。t=2*t0/3)+0.*(t=2*t0/3)。mt=fun(t)。mh=imag(hilbert(mt))。b=sin(2*pi*fc.*t)ussb=mt.*ctimag(hilbert(mt)).*b。lssb=mt.*ct+imag(hilbert(mt)).*b。figure(1)subplot(3,1,1),plot(t,mt),title(39。mt調(diào)制信號39。)subplot(3,1,2),plot(t,ussb),title(39。ussb載波信號39。)subplot(3,1,3),plot(t,lssb),title(39。lssb已調(diào)信號39。) dt=t(2)t(1)。 % 采樣周期f=1/dt。 % 采樣頻率(Hz)X1=fft(ussb)。 % 計(jì)算x的快速傅立葉變換XX2=fft(lssb)。N=1/dt。F1=X1(1:N/2+1)。 % F(k)=X(k)(k=1:N/2+1)F2=X2(1:N/2+1)。f=f*(0:N/2)/N。 % 使頻率軸f從零開始figure(4)subplot(3,1,1),plot(f,abs(F1)),hold on,title(39。ussb調(diào)制信號頻譜圖39。)。plot(f,abs(F1))。xlim([fc*2,fc*2])。subplot(3,1,2),plot(f,abs(F2)),hold on,title(39。lssb39。)。plot(f,abs(F2))。xlim([fc*2,fc*2])。xlabel(39。Frequency39。)。ylabel(39。|F(k)|39。)（3）高斯白噪聲n0=30。fc=t=[0::]。u=cos(2*pi*fc.*t)。M=length(u)。x=sqrt(4*n0)*randn(1,M)。plot(t,x(1:length(t)))。 數(shù)字基帶通信系統(tǒng)的仿真(1)一個(gè)升余弦頻譜的濾波器， t=1/100+eps:1/10000:1/100。alfa=。ts=1/1000。h_t=sin(pi*t/ts)./(pi*t/ts).*(cos(alfa*pi*t/ts)./(14*alfa^2*t.^2/ts^2))。plot(t,h_t)。 (2) 利用matlab的SIMULINK功能建立一個(gè)基帶傳輸模型，采用單極性, 或雙極性碼作為基帶信號，發(fā)送濾波器為上述升余弦濾波器，發(fā)送數(shù)據(jù)率為1000bps，分別觀察輸出信號在無噪聲干擾及有噪聲干擾下波形及眼圖。無噪聲時(shí)眼圖有噪聲時(shí)的眼圖注意：必須首先運(yùn)行實(shí)驗(yàn)步驟1 中的程序得出h(t)后，才能運(yùn)行該模型。(3)建立基帶接收機(jī)模型，觀察判決輸出，與發(fā)送數(shù)碼進(jìn)行比較。（選做）（1）試做單極性歸零碼、雙極性非歸零碼、單極性非歸零碼、雙極性歸零碼，占空比50%（選擇其中2種）；（2）雙相碼、AMI碼、HDB3碼等的仿真（選其中2種）。單極性非歸零碼數(shù)字調(diào)制與解調(diào)的仿真(選擇兩種調(diào)function y=snrz(x)%本函數(shù)實(shí)現(xiàn)將輸入的一段二進(jìn)制代碼編為相應(yīng)的單極性非歸零碼輸出%輸入x為二進(jìn)制碼，輸出y為編好的碼%給出計(jì)算每一個(gè)碼元的點(diǎn)數(shù)，因?yàn)槲覀冎挥杏秒x散的點(diǎn)來得出連續(xù)的函數(shù)表示。grid=300。t=0:1/grid:length(x)。for i=1:length(x), %計(jì)算碼元的值 if(x(i)==1), %如果信息為1 for j=1:grid, %該碼元對應(yīng)的點(diǎn)值取1 y((i1)*grid+j)=1。 end else for j=1:grid。 %反之，信息為0，該碼元對應(yīng)的點(diǎn)值取0 y((i1)*grid+j)=0。 end。 end。 endy=[y,x(i)]。M=max(y)。m=min(y)。subplot(2,1,1)。plot(t,y)。axis([0,i,M+])。%采用title命令來實(shí)現(xiàn)標(biāo)記出各碼元對應(yīng)的二進(jìn)制信息title(39。單極性非歸零碼39。)。
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1