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

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

目录 目录 ……………………………………………………………………………………………… 正文………………………………………………………………………………………………… 1 基本波形……………………………………………………………………………… 1.1 正弦波………………………………………………………………………… 1.2 三角波和锯齿波……………………………………………………………… 1.3 方波和不对称方波…………………………………………………… 1.4伪码…………………………………………………………………………………… 2 调幅……………………………………………………………………………… 2.1常规双边带调幅(AM)…………………………………………………………… 2.2抑制载波双边带调幅(DSB)…………………………………………………………… 2.3 单边带调幅(SSB)………………………………………………………………… 3.4振幅键控(ASK)…………………………………………………………………… 3 调频……………………………………………………………………………… 3.1单音频调频…………………………………………………………………………… 3.2锯齿波调频…………………………………………………………………………… 3.3 三角波调频………………………………………………………………………… 3.4 移频键控(FSK)………… …………………………………………………………… 4 调相…………………………………………………………………………………………… 4.1 方波调相………………………………………………………………………………… 4.2 伪码调相………………………………………………………………………………… 一、 基本波形 (1)正弦波 为正弦波频率，为初始相位 理论频谱： 为正弦波角频率，且 时域波形： 图1-1-1 正弦波的时域波形(横坐标为时间t，纵坐标为幅度) 图1-1-1为正弦波的时域波形。从图中可以看出，正弦波的频率为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); grid on xlabel( t ); ylabel( 正弦波 ); N=1024; %Plot amplitude spectrum sf=fft(y,N); f=-fs/2:fs/N:(fs/2-fs/N); figure(2) plot(f,fftshift(abs(sf))/max(abs(sf))); axis([-100 100 0 1]); grid on; title( Amplitude Spectrum of正弦波 ) xlabel( Frequency (Hz) ) ylabel( |F(f)| ) (2)三角波和锯齿波(sawtooth) 三角波和锯齿波时域表达式：f(t)=2*fm*twidth-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为三角波锯齿波的时域波形。从图中可以看出，频率为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(2*pi*fm*t,0); %锯齿波 x=sawtooth(2*pi*fm*t,1); %锯齿波 figure(1); subplot(3,1,1); plot(t,y); grid on; xlabel( t ); ylabel( 三角波 ); subplot(3,1,2); plot(t,z); grid on; xlabel( t ); ylabel( 锯齿波1 ); subplot(3,1,3); plot(t,x); grid on; xlabel( t ); ylabel( 锯齿波2 ); N=2^20; %Plot amplitude spectrum Y=fft(y,N); Z=fft(z,N); X=fft(x,N); f=-fs/2:fs/N:(fs/2-fs/N); figure(2) subplot(3,1,1); plot(f,fftshift(abs(Y))/max(abs(Y))); axis([-300 300 0 1]); grid on; title( Amplitude Spectrum of 三角波 ) xlabel( Frequency (Hz) ) ylabel( |F(f)| ) subplot(3,1,2); plot(f,fftshift(abs(Z))/max(abs(Z))); axis([-300 300 0 1]); grid on; title( Amplitude Spectrum of 锯齿波1 ) xlabel( Frequency (Hz) ) ylabel( |F(f)| ) subpl