MATLAB FM仿真報告

1、實驗二FM仿真一實驗題目假設(shè)基帶信號，載波頻率為20kHz， FM的頻率偏移常數(shù)為5kHz/V。仿真產(chǎn)生FM信號，觀察已調(diào)信號的波形和頻譜。二 基本原理單音頻信號經(jīng)FM調(diào)制后的表達式為其中調(diào)制指數(shù)。同實驗一中相仿，定義必要的仿真參數(shù)，在此基礎(chǔ)上可得到載波信號和調(diào)制信號。根據(jù)可得到頻偏，由此可寫出最終的FM信號的表達式進行仿真計算。對FM信號進行傅里葉變換可得頻譜特性，變換依舊使用實驗一中給出的t2f.m函數(shù)。三仿真方案 開始設(shè)置參數(shù)調(diào)制信號時域范圍 t1采樣頻率ts載波頻率 fc調(diào)制信號頻率 fm頻偏常數(shù) kf生成載波s=cos（2*pi*fc*t）生成調(diào)制信號x=sin(2*pi*4*fm*

2、t)+2*cos(2*pi*2*fm*t)+4*sin(2*pi*fm*t+pi/3)調(diào)制結(jié)束四 仿真源代碼cleart1=0.1; %調(diào)制信號的時域范圍 fs=600000; %抽樣頻率 ts=1/fs; %采樣率 t=-t1:ts:t1; fc=20000; %設(shè)定載波頻率 fc=20kHzs=cos(2*pi*fc*t); %生成載波 figureplot(t,s);xlabel(時間);ylabel(幅度);title(載波波形);grid onaxis(-0.0001 0.0001 -2 2);fm=250; %設(shè)定調(diào)制信號頻率 kf=5000;x=sin(2*pi*4*fm*t)+

3、2*cos(2*pi*2*fm*t)+4*sin(2*pi*fm*t+pi/3);%生成調(diào)制信號 figure;plot(t,x);grid onxlabel(t);ylabel(x);axis(0 0.01 -10 10)title(調(diào)制信號波形)figurea=fftshift(fft(x);f=linspace(-fs/2,fs/2,length(t);plot(f,abs(a);grid onxlabel(ffrequence(Hz);ylabel(powerspectrum(x);axis(-1500 1500 0 100);title(調(diào)制信號頻譜圖);figurephi=2*pi

4、*kf*cumsum(x)*ts;y=cos(2*pi*fc*t+phi);plot(t,y);grid onxlabel(t(s);ylabel(y);axis(0 0.01 -5 5);b=fft(y,1024);f=(0:length(b)-1)*fs/length(b)-fs/2;title(已調(diào)信號波形);figureplot(f,abs(b);grid onxlabel(frequence(Hz);ylabel(powerspectrum(x);axis(-300000 300000 0 200);title(已調(diào)信號頻譜圖);grid on五實驗結(jié)果及分析1.調(diào)制信號波形極其頻譜圖根據(jù)調(diào)制信號表達式可知信號x由三個正弦波相加而成，其中三個正弦波頻率符合頻率頻譜圖所示。2.載波波形載波波形符合程序設(shè)計的20kHz的正弦波。3.仿真FM信號波形，其形狀為疏密波，最大頻偏5kHz/V。由圖可讀出帶寬為70kHz，符合公式計算B2（+1）fm=2（6.15+4）=69kHz4.已調(diào)信號頻譜圖六實驗心得通過這次通信原理實驗，我發(fā)現(xiàn)自己的基礎(chǔ)知識和應(yīng)用能力都很差。剛拿到題目時，在如何實現(xiàn)FM調(diào)制的地方就卡住了。最后通過各種網(wǎng)上的資料解決了問題。通過了這次實驗，