數(shù)字的集群MDR中4FSK地matlab仿真

實用標準文案1.基本概念物理層的定義h=0.27；T是樣期1/4800秒。D的含義是Maximumdeviation最大頻偏4FSK恒包絡調(diào)制，所以在終端率上，應該比較省電。4制的隨機信號產(chǎn)生，并把信強度調(diào)整到規(guī)范規(guī)定的3,1,-1,-3：M=4;Ts=1/4800;%的特率；x=randint(1000,1,M);x=x.*2-3;調(diào)到正確的symbol；plot(x);部分信號繪制如下：3210-1-2-3210220230250260270280290由于是隨機信號，所以是寬帶信。頻譜如下。只繪制到信號的采樣頻率，繪制頻譜：lx=length(x);Freqs=1/Ts;freq=[-Freqs/2:Freqs/lx:Freqs/2-Freqs/lx];S=fftshift(abs(fft(x)));Sdb=10*log10(S);plot(freq,Sdb);精彩文檔

實用標準文案500500進行基帶濾波，濾波后的信號和譜為下圖。注意信號帶寬被限制為根升余弦濾波的3dB帶寬以內(nèi)。這個帶，f=640+1920=2560Hz。一半帶寬：%波器的基本參數(shù)sqrt，0.2。Fs至是，delay默認的。y=rcosflt(x,1,2,'fir/sqrt',0.2,3,0);%繪制頻譜：ly=length(y);Freqs=2/Ts;freq=[-Freqs/2:Freqs/ly:Freqs/2-Freqs/ly];Sy=fftshift(abs(fft(y)));Sdby=10*log10(Sy);figure;plot(freq,Sdby);2520151050-5-10-15-20-5000-2000-1000

0100030005000

3210-1-2-3620640680760接下去，對4FSK制進行仿真精彩文檔

實用標準文案2.4FSK調(diào)制仿真4FSK3種式前余弦濾波后的基帶數(shù)據(jù)進行頻率振蕩器的控制方式實現(xiàn)容易，但是頻率的穩(wěn)定性差方式2，鍵控法。鍵控分頻器的分頻比，以控制2種分頻脈沖的輸出經(jīng)過整形可以得到信在DSP實現(xiàn)中常用到法3用函不管用何種方式DMR規(guī)中，4FSK的調(diào)制的結(jié)果給出了要求，要求在給定序列下，輸出的調(diào)制信號該的號過大偏差的FM調(diào)輸出的信號一致。方1fmmod。是調(diào)制函數(shù)以Fs采樣來看待輸入和輸出的數(shù)據(jù)系帶的帶寬小于可以設計的低一些，以降低采樣頻率，減少計算量樣可以設計Fc=15kHz足夠高。這樣，為了讓輸出的圖像細致些，我們設計采樣頻率Fs=150KHz這個精度應該購了定好了這個調(diào)后，回過頭看函數(shù)。這個函數(shù)中我們感興趣的就是它的Fs設。因為是數(shù)字信號以可以1個其實表1/Ts=4800Hz的目的是讓rcosflt的的數(shù)據(jù)能夠fmmod的的速率匹配上，于是不難算出，每個輸入symbol當于被采樣個據(jù)，就個點。取整后，是32個點。反過來，調(diào)整采樣頻率成為。既然采樣點數(shù)已經(jīng)確定是32，rcosflt的也確定了，調(diào)整為32所以matlab代碼如下（減少些數(shù)為20便于計算M=4;Ts=1/4800;%的特率；x=randint(200,1,M);x=x.*2-3;調(diào)到正確的symbol；%波器的基本參數(shù)sqrt，0.2。Fs至是，取32，delay用認的。y=rcosflt(x,1,32,'fir/sqrt',0.2,3,0);過渡補濾波器增益會下降。y=(3/max(y)).*y;%調(diào)波增益，恢復信號強度。%升余弦濾波后，輸出的就是擬的基帶信號了。%擬FM。Fc=15KHz,Fs=153.6KHz,（為后最值就在1左右）Fc=15000;Fs=153600;freqdev=;ym=fmmod(y,Fc,Fs,freqdev);%制頻譜lx=length(ym);Freqs=Fs;freq=[-Freqs/2:Freqs/lx:Freqs/2-Freqs/lx];S=fftshift(abs(fft(ym)));Sdb=10*log10(S);plot(freq,Sdb);精彩文檔

實用標準文案用上面的方法，根據(jù)標準輸入，以看到輸出的信號頻譜為：M=4;Ts=1/4800;%的特率；fori=1:200if(mod(i,4)<2)3;elsex(i)=-3;endend%波器的基本參數(shù)sqrt，0.2。Fs至是，取32，delay用認的。y=rcosflt(x,1,32,'fir/sqrt',0.2,3,0);過渡補濾波器增益會下降。y=(3/max(y)).*y;%調(diào)波增益，恢復信號強度。%擬FM。Fc=15KHz,Fs=153.6KHz,（為后最值就在1左右）Fc=15000;Fs=153600;freqdev=ym=fmmod(y,Fc,Fs,freqdev);%制頻譜lx=length(ym);Freqs=Fs;freq=[-Freqs/2:Freqs/lx:Freqs/2-Freqs/lx];S=fftshift(abs(fft(ym)));Sdb=10*log10(S);plot(freq,Sdb);方2直寫數(shù)其實就是模擬FM方，對輸入信號進行積分。略。方3fskmod(x,M,freq_sep,nsamp)注意，這個函數(shù)是有問題的。我它進行了修改。對fskmod.m的中修改一行下：phaseIncr(0:nSamp-1)'*(-(M-1):2:(M-1))*2*pifreq_sep/2*samptime/nSamp;對fskdemod的修改一行如：t=[0:1/(nSamp*Fs):-1/(nSamp*Fs)]';另注，fskmod這函中已經(jīng)x的據(jù)整-(M-1):2:(M-1)進處了更新后的函數(shù)我重新命令為已輸入的調(diào)整-(M-1):2:(M-1)了符合的需求。思路還是：通過函數(shù)得到輸出的信號對信進行根升余弦濾波；進行載波調(diào)制。所以過程如下：精彩文檔

實用標準文案M=4;Ts=1/4800;%的特率；%x=randint(100,1,M);x=[0123210123210];Num=length(x);freqsep=0.27;nsamp=10;Fs=1;%Fc=nsamp*1/Ts;載率Fc=10000;%波頻率y=gyzfskmod(x,M,freqsep,nsamp,Fs);注意這里調(diào)制出來的信號復數(shù)表示的模擬信號對擬的IQ信不能夠直接用數(shù)字濾波器正的如下：%計濾波器，用的Nyquist類。%出信號的Fs是nsamp/Ts以Band=nsamp/2;%Frequency設置為歸化為1;%Transition是filterbuilder%輸出Hnyq。個濾波器的時延有直接的控制參數(shù)。通過解調(diào)，可以看到大概時延2個symbol。%出濾波yf=filter(Hnyq,y);%入濾波yfin=filter(Hnyq,y);%調(diào)Z=gyzfskdemod(yfin,M,freqsep,nsamp);%出波形繪圖tlen=length(x)*Ts;%%間分片t=[0:tlen/length(yf):tlen-tlen/length(yf)];%算載波ccos=cos(2*pi*Fc.*t);csin=sin(2*pi*Fc.*t);%算輸出波形S=reshape(real(yf),1,length(yf)).*ccos-reshape(imag(yf),1,length(yf)).*csin;plot(S);%制頻譜精彩文檔

實用標準文案figure;p=rea