信号波形的产生与频谱的matlab程序

《信号波形的产生与频谱的matlab程序》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、目录目录………………………………………………………………………………………………正文…………………………………………………………………………………………………1基本波形………………………………………………………………………………1.1正弦波…………………………………………………………………………1.2三角波和锯齿波………………………………………………………………1.3方波和不对称方波……………………………………………………1.4伪码……………………………………………………………………………………2调幅………………………………………………

2、………………………………2.1常规双边带调幅（AM）……………………………………………………………2.2抑制载波双边带调幅(DSB)……………………………………………………………2.3单边带调幅（SSB）…………………………………………………………………3.4振幅键控（ASK）……………………………………………………………………3调频………………………………………………………………………………3.1单音频调频……………………………………………………………………………3.2锯齿波调频…………………………………………………………………………

3、…3.3三角波调频…………………………………………………………………………3.4移频键控（FSK）………………………………………………………………………4调相……………………………………………………………………………………………4.1方波调相…………………………………………………………………………………4.2伪码调相…………………………………………………………………………………一、基本波形（1）正弦波为正弦波频率，为初始相位理论频谱：为正弦波角频率，且时域波形：图1-1-1正弦波的时域波形（横坐标为时间t，纵坐标为幅度）图1-1-1为正

4、弦波的时域波形。从图中可以看出，正弦波的频率为50Hz，峰值为1，为1，所以初始相位为。设采样率=1000Hz，一个周期采样20个点。频域波形：图1-1-2正弦波的频域波形（横坐标为频率HZ，纵坐标为归一化幅度）从图中可以看出，正弦波在正负50HZ的地方有二根谱线，与理论相符。MATLAB程序：%正弦波波形fm=50;%正弦波频率fs=1000;%采样频率dt=1/fs;phase=0;%初相t=0:dt:0.1;figure(1)y=sin(2*pi*fm*t+phase);%正弦波plot(t,y);gridonxlabel('t'

5、);ylabel('正弦波');N=1024;%Plotamplitudespectrumsf=fft(y,N);f=-fs/2:fs/N:(fs/2-fs/N);figure(2)plot(f,fftshift(abs(sf))/max(abs(sf)));axis([-10010001]);gridon;title('AmplitudeSpectrumof正弦波')xlabel('Frequency(Hz)')ylabel('

6、F(f)

7、')（2）三角波和锯齿波（sawtooth）三角波和锯齿波时域表达式：f(t)=2*fm*twi

8、dth-1,0≤t≤T2-2*fm*twidth+1+width1-width,T2≤t≤T为三角波频率。周期三角波的傅立叶级数ft=4Eπ21∞1n2sin2nπ2cosnω1t周期锯齿波的傅立叶级数ft=Eπ1∞(-1)n+11nsinnω1t产生一个周期为T、幅度在-1到+1之间的周期性三角波信号。其中表示最大幅度出现的位置：即在一个周期内，信号从t=0到×T时函数值从-1到+1线性增加，而从×T到T又是从+1到-1线性下降。取值在0~1之间。 图1-2-1三角波和锯齿波的时域波形（横坐标为时间t，纵坐标为幅度）图1-2-1为三角

9、波锯齿波的时域波形。从图中可以看出，频率为50Hz，设采样率=10000Hz，一个周期采样200个点。频域波形：图1-2-2三角波和锯齿波的时域波形（横坐标为频率HZ，纵坐标为归一化幅度）。由上图可以看出三角波的频谱含有奇次谐波（50HZ的奇数倍）的余弦分量，幅度以1n2的规律收敛。锯齿波的频谱只包含正弦分量，幅度以1n的规律收敛。MATLAB程序：%锯齿波和三角波波形fm=50;%锯齿波频率fs=10000;%采样频率dt=1/fs;t=0:dt:0.1;y=sawtooth(2*pi*fm*t,0.5);%三角波z=sawtooth

10、(2*pi*fm*t,0);%锯齿波x=sawtooth(2*pi*fm*t,1);%锯齿波figure(1);subplot(3,1,1);plot(t,y);gridon;xlabel('t');yla