資源描述:
《直接序列擴(kuò)頻.docx》由會員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、一、實(shí)驗(yàn)?zāi)康模?、熟悉MATLAB語言的基本用法;2、掌握MATLAB語言中數(shù)據(jù)信號的產(chǎn)生;3、掌握直接序列擴(kuò)頻信號的產(chǎn)生;4、掌握直接序列擴(kuò)頻信號的解擴(kuò)方法;5、掌握MATLAB語言中信號頻譜的繪制方法。二、實(shí)驗(yàn)原理:(1)利用MATLAB隨機(jī)產(chǎn)生數(shù)據(jù)比特;利用MATLAB隨機(jī)生成PN序列;將數(shù)據(jù)比特與PN序列相乘完成信號擴(kuò)頻;將擴(kuò)頻信號與PN序列再次相乘完成解擴(kuò)。圖1直接序列擴(kuò)頻系統(tǒng)原理圖(2)最大相關(guān)峰值解擴(kuò),就是設(shè)置一個(gè)標(biāo)準(zhǔn)門限,然后通過條件判斷將信號值變?yōu)?,-1,從而將原先加噪的擴(kuò)頻
2、信號解擴(kuò)。三、實(shí)驗(yàn)步驟:程序代碼:Wave_test2.mclearall;clc;%***********************%基本要求內(nèi)容len=5;%要產(chǎn)生的數(shù)據(jù)長度;N=10;%一個(gè)符號用10個(gè)采樣點(diǎn)表示signal=randi([01],1,len);%randi產(chǎn)生的數(shù)據(jù)范圍[0,1],產(chǎn)生的大小1*lensignal2=2*signal-1;%將數(shù)據(jù)范圍變?yōu)閇-1,1]signal3=repmat(signal2,N,1);%將數(shù)據(jù)進(jìn)行復(fù)制,將signal2復(fù)制為N行1列sig
3、nal4=reshape(signal3,1,N*len);%將signal3變形為1行N*len列figure(1);subplot(411);stem(signal4);title('原始采樣信號');%***********%產(chǎn)生PN碼,長度為為N*LenPN1=randi([01],1,N*len);PN2=2*PN1-1;subplot(412);stem(PN2);title('PN序列');%************ds1=signal4.*PN2;%進(jìn)行擴(kuò)頻subplot(413
4、);stem(ds1);title('擴(kuò)頻信號');%************ds2=ds1.*PN2;subplot(414);stem(ds2);title('解擴(kuò)信號');%********************************%畫頻譜圖NS=512;fft_signal4=fft(signal4,NS);fft_PN2=fft(PN2,NS);fft_ds1=fft(ds1,NS);fft_ds2=fft(ds2,NS);figure(2);subplot(411);plot(
5、(-NS/2+1):NS/2,abs(fftshift(fft_signal4)));title('原始采樣信號頻譜');subplot(412);plot((-NS/2+1):NS/2,abs(fftshift(fft_PN2)));title('PN碼頻譜');subplot(413);plot((-NS/2+1):NS/2,abs(fftshift(fft_ds1)));title('擴(kuò)頻信號頻譜');subplot(414);plot((-NS/2+1):NS/2,abs(fftshif
6、t(fft_ds2)));title('解擴(kuò)信號頻譜');%********************************************%提高要求內(nèi)容。coef=[011];%反饋系數(shù)reg=[001];%移位寄存器初始狀態(tài)len=50;%需要生產(chǎn)的m序列的長度m_code=zeros(1,len);m=3;%移位寄存器級數(shù)fori=1:lenm_code(i)=reg(end);reg_next=mod(sum(reg.*coef),2);reg=[reg_nextreg(1:m-
7、1)];endfigure(3);subplot(311);stem(signal4);title('原始信號');%******subplot(312);stem(m_code);title('m序列');%*********PN3=2*m_code-1;%將m序列變?yōu)閇-11]形成PN3碼subplot(313);stem(PN3);title('用m序列產(chǎn)生的PN3');%*****************************%畫頻譜NS=512;fft_signal4=fft(sig
8、nal4,NS);fft_m_code=fft(m_code,NS);fft_PN3=fft(PN3,NS);figure(4);subplot(311);plot((-NS/2+1):NS/2,abs(fftshift(fft_signal4)));title('原始采樣信號頻譜');subplot(312);plot((-NS/2+1):NS/2,abs(fftshift(fft_m_code)));title('m序列頻譜');subplot(313);plot((-NS/2+1):NS/