matlab教程源程序代碼仿真fm系統(tǒng)

和FM系統(tǒng)調(diào)制解調(diào)的基本過程，利用仿真FM信號的產(chǎn)生與接收，包括加深對FM系統(tǒng)的理解。、系統(tǒng)，如電視信號的呃傳輸以及中。頻率調(diào)制是一種模擬角度調(diào)制方、設信號為m(t)，則FMfi(t)kFMFM(t)2kFM

可知，F(xiàn)M

fmax

FM

fmax 則有調(diào)頻指數(shù)為 =????????下面，考慮調(diào)頻信號的接收。FM

檢波以及隔直濾波后得到s(t)=k*m(t).由此可以解調(diào)得到模擬信號，由于解調(diào)設備簡單是基帶頻譜的搬移，根據(jù)計算，有卡森規(guī)則可以得到FM信號的頻譜為BT2fmax2B2kFMmt

頻偏，大于兩倍調(diào)制信號帶寬。雖然FM一般具有較大帶寬，但是較好的抗噪性??????根據(jù)FM信號系統(tǒng)原理，仿真FM信號的產(chǎn)生與接收。調(diào)制信號：單音信號m(t=cos(2πft)f=20hz。載頻信號：x(t)=cos(2π????t)，其中????=500hz。調(diào)頻指數(shù)：??????=10，由此可以計算出??????=200。 。值，實驗中取小信噪比為15dB，大信噪比為35dB，觀察解況。一、FM %t0秒的信號，t01 一秒內(nèi)20hz單音信號1幅度0幅度 時間一個周期20hz單音信號1幅度0幅度 時間 s=A*cos(2*pi*fc*t+2*pi*k*thelta);%（35dBr0=awgn(s,35)35db 150hz500700hzfsamp=fs;fcuts=[100250750mags=[01devs=[0.010.05hh=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale'); 乎和調(diào)制信號一樣。二者較大區(qū)別在于邊界處，是由于包絡求解使用的方法存在的加35dB噪聲的FM信號21幅度0幅度 時間不加噪聲FM信號21幅度0幅度 時間

微分加噪聲FM信號幅度0幅度 時間/s解調(diào)信號幅度幅度0 時間/s調(diào)制信號1幅度0幅度 時間（15dB % %帶通濾波,10，fm20B=440hz,500hz，[250750]為150hz500700hzfsamp=fcuts=[100250750mags=[01devs=[0.010.05hh=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');15dB噪聲接收FM信號21幅度0幅度 時間不加噪聲FM信號21幅度0幅度 時間

微分加噪聲接收FM信號幅度0幅度 時間/s解調(diào)信號幅度幅度0 時間/s調(diào)制信號1幅度0幅度 時間 % FM信號的方式，其解調(diào)方法簡單，不需要載波同步，受到了親睞。

微分接收FM信號幅度0幅度 時間/s解調(diào)信號幅度幅度0 時間/s調(diào)制信號1幅度0幅度 時間150hz500700hzfsamp=fs;fcuts=[100250750mags=[01devs=[0.010.05hh=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');為通帶頻率，200hzfsamp=fs;fcuts=[50200];mags=[1devs=[0.05hh=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale'); 10FM系統(tǒng)，相干解調(diào)不能得到信號，但是對于窄帶調(diào)頻信號（1，可以使用相干解調(diào)方法。0 （NFFT=ym= %fft fy=kym/NFFT*fs;圖2為已調(diào)信號的頻譜，根據(jù)卡森可以得到，帶寬為440hz，并且以500hz對稱，351比較，可以看到基本一致，有零頻信號，隔直后則可。-通過學習FM系統(tǒng)，我們分析仿真了FM信號的產(chǎn)生與接收的過程。接收方FM系程序%FM f=20t0=0.1t0秒的信號，t01title('20hz單音信號');title('20hz單音信號');%FM % %35dbfigure(3),subplot(2,1,1),plot(t,r0);title('35dBFM信號');xlabel('時間/s');ylabel('');axis([0t0-2- fsamp=fs; mags=[01devs=[0.010.05

500hz，[250750]為通帶頻率，150hz500，頻帶700hzhh=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');mm=[mzeros(1,length(t1)-length(t))];%15db');axis([0t0-2- fsamp=fs;fcuts=[100250750900]; %10，fm20B=440hz,，500hz，[250750]為通帶頻率，150hz500，頻帶700hzmags=[01devs=[0.010.05hh=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');微分加噪聲接收mm=[mzeros(1,length(t1)-length(t))]; % %fsamp=fs;fcuts=[100250750900]; 500hz，[250750]150hz500，頻帶寬度700hzmags=[01devs=[0.010.05hh=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');fsamp=fs;fcuts=[50200]; 500hz，[250750]150hz500700hzmags=[10];devs=[0.05hh=fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');NFFT=2^nextpow2(length(m));ym=fftshift(fft(m,NFFT));fy=NFFT=2^nextpow2(length(s));ym=fy=subplot(5,1,2),plot(fy,abs(ym));title('FM信號頻譜');xlabel('頻率/Hz');ylabel('|ym|');NFFT=2^nextpow2(length(r0));ym=fftshift(fft(r0,NFFT));fy=subplot(5,1,3),plot(fy,abs(ym)