【数字信号调制】基于matlab PM调制与解调【含Matlab源码 2657期】.zip

clc; close all; clear all; %% % 信号的PM调制与解调 % 此程序给出了信号的PM调制的与解调的仿真过程，并分析了在不同信噪比条件下解调的性能 %% 信号参数 fm = 3 ; fc = 10 ; fs = 200 ; Am = 1 ; A = 2 ; N = 1e3 ; K = N -1 ; n = 0 :K; f = n*fs/N; t = 0 : 1 /fs:K/fs; wm = 2 * pi *fm; wc = 2 * pi *fc; thm = 0 ; thc = 0 ; B = 2 *fm; KFM = 5 ; beatFM = KFM*Am; ct = A* cos (wc*t); fft_ct = abs (fft(ct)); figure ( 1 ) subplot(2,1,1); plot (t,ct); xlabel( '时间/s' ); ylabel( '幅值' ); grid on; title( '频率为10Hz的载波信号byLDM' ); subplot( 2,1,2 ); plot (f( 1 :N/ 2 ),fft_ct( 1 :N/ 2 )); xlabel( '频率/Hz' ); ylabel( '幅值' ); grid on; title( '载波信号信号频谱' ); %% 调制PM sig_PM = cos (wc*t+beatFM* cos ( 2 * pi *wm*t)); %PM已调制信号 fft_sig_PM = abs (fft(sig_PM)); figure ( 2 ) subplot( 2,1,1 ); plot (t,sig_PM); xlabel( '时间/s' ); ylabel( '幅值' ); grid on; title( '调制信号tPM' ); subplot( 212 ); plot (f( 1 :N/ 2 ),fft_sig_PM( 1 :N/ 2 )); xlabel( '频率/Hz' ); ylabel( '幅值' ); grid on; title( '调制信号频谱' ); %% 加噪声 sig_noise_PM1 = awgn(sig_PM, 3 ); % 添加信噪比为6dB的高斯白噪声 fft_sig_noise_PM1 = abs (fft(sig_noise_PM1)); % 对AM解调信号进行傅里叶变换 sig_noise_PM2 = awgn(sig_PM, 20 ); % 添加信噪比为15dB的高斯白噪声 fft_sig_noise_PM2 = abs (fft(sig_noise_PM2)); % 对AM解调信号进行傅里叶变换 figure ( 3 ) subplot( 2 , 1 , 1 ); plot (t,sig_noise_PM1); xlabel( '时间/s' ); ylabel( '幅值' ); grid on; title( '添加信噪比为40dB噪声后的FM信号byLDM' ); subplot( 2 , 1 , 2 ); plot (f( 1 :N/ 2 ),fft_sig_noise_PM1( 1 :N/ 2 )); xlabel( '频率/Hz' ); ylabel( '幅值' ); grid; title( '添加信噪比为40dB噪声后的FM信号频谱' ); figure ( 4 ) subplot( 2 , 1 , 1 ); plot (t,sig_noise_PM2); xlabel( '时间/s' ); ylabel( '幅值' ); grid on; title( '添加信噪比为20dB噪声后的FM信号byLDM' ); subplot( 2 , 1 , 2 ); plot (f( 1 :N/ 2 ),fft_sig_noise_PM2( 1 :N/ 2 )); xlabel( '频率/Hz' ); ylabel( '幅值' ); grid; title( '添加信噪比为20dB噪声后的FM信号频谱' ); %% 解调 % 无噪声 diff_sig_PM = diff(sig_PM); % 调频信号经过一个微分器 d_m = abs(hilbert(diff_sig_PM)); % 用希尔伯特变换能够实现包络检波 fft_dm = abs(fft(d_m)); % 对AM解调信号进行傅里叶变换 % 3 dB高斯白噪声 diff_sig_PM1 = diff(sig_noise_PM1); % 调频信号经过一个微分器 d_m2 = abs(hilbert(diff_sig_PM1)); % 用希尔伯特变换能够实现包络检波 fft_dm2 = abs(fft(d_m2)); % 对AM解调信号进行傅里叶变换 % % 20 dB高斯白噪声 % diff_sig_fm2 = diff(sig_noise_fm2); % 调频信号经过一个微分器 % d_m3 = abs(hilbert(diff_sig_fm2)); % 用希尔伯特变换能够实现包络检波 % fft_dm3 = abs(fft(d_m3)); % 对AM解调信号进行傅里叶变换 % figure(5) subplot(2,1,1); plot(t(1:end-1),d_m); xlabel('时间/s'); ylabel('幅值'); grid on; title('无噪声条件下解调信号波形FMbyLDM'); subplot(2,1,2); plot(f(1:N/2),fft_dm(1:N/2)); xlabel('频率/Hz'); ylabel('幅值'); axis([0 100 0 100]); grid; title('FM解调信号频谱byCL'); figure(6) subplot(2,1,1); plot(t(1:end-1),d_m2); xlabel('时间/s'); ylabel('幅值'); grid on; title('信噪比为40dB条件下解调信号波形FMbyLDM'); subplot(2,1,2); plot(f(1:N/2),fft_dm2(1:N/2)); xlabel('频率/Hz'); ylabel('幅值'); axis([0 100 0 100]); grid; title('FM解调信号频谱byCL'); % figure(7) % subplot(2,1,1); % plot(t(1:end-1),d_m3); % xlabel('时间/s'); % ylabel('幅值'); % grid on; % title('信噪比为20dB条件下解调信号波形FMbyLDM'); % subplot(2,1,2); % plot(f(1:N/2),fft_dm3(1:N/2)); % xlabel('频率/Hz'); % ylabel('幅值'); % axis([0 100 0 100]); % grid; % title('FM解调信号频谱byCL'); % figure(5) % subplot(2,1,1); % plot(t(1:end-1),d_m); % xlabel('时间/s'); % ylabel('幅值'); % grid on; % title('无噪声条件下解调信号波形FMbyLDM'); % subplot(2,1,2); % plot(f(1:N/2),fft_dm(1:N/2)); % xlabel('频率/Hz'); % ylabel('幅值'); % axis([0 100 0 100]); % grid; % title('FM解调信号频谱byCL'); % % figure(6) % subplot(2,1,1); % plot(t(1:end-1),d_m2); % xlabel('时间/s'); % ylabel('幅值'); % grid on; % title('信噪比为40dB条件下解调信号波形FMbyLDM'); % subplot(2,1,2); % plot(f(1:N/2),fft_dm2(1:N/2)); % xlabel('频率/Hz'); % ylabel('幅值'); % axis([0 100 0 100]); % grid; % title('FM解调信号频谱byCL'); % % figure(7) % subplot(2,1,1); % plot(t(1:end-1),d_m3); % xlabel('时间/s'); % ylabel('幅值'); % grid on; % title('信噪比为20dB条件下解调信号波形FMbyLDM'); % subplot(2,1,2); % plot(f(1:N/2),fft_dm3(1:N/2)); % xlabel('频率/Hz'); % ylabel('幅值'); % axis([0 100 0 100]); % grid; % title('FM解调信号频谱byCL');