通信原理实践(二)——幅度调制

一、幅度调制,并画出时域和频域波形

1、代码如下:

function [ p_n ] = AM_func( N,fs,fm,Am,fc,Ac,Ma )
%UNTITLED 此处显示有关此函数的摘要
% 此处显示详细说明
% N = 3000; % 长度
% fs = 48e3; % 采样率
% fm = 300; % 调制信号频率
% Am = 0.5; % 调制信号幅度
% fc = 8e3; % 载波频率
% Ac = 1; % 载波幅度
% A0 %直流分量 A0 = Am / Ma; %%
%时域波形,第一小题
t = (0:N-1)'/fs; % 时间t
m_t = Am*sin(2*pi*fm*t); % 产生调制信号
c_t = Ac*cos(2*pi*fc*t); % 产生载波信号
s_t = (m_t + A0) .* c_t; % 调制 figure()
subplot(1,3,1);
plot(t, s_t , 'b'); % 绘波形
xlabel('time');
ylabel('amplitude');
title('时域波形'); %%
%绘制功率谱,第二小题
L = length(s_t); % 取得序列长度
u = fftshift(fft(s_t )); % 离散傅里叶变换,求频谱
u_pow = pow2db(abs(u).^2); % 幅度转为dB
w = (0:L-1)'*fs/L - 1/2*fs; % 横坐标-频率 subplot(1,3,2);
plot(w, u_pow);
grid on;
xlabel('frequency(Hz)');
ylabel('magnitude(dB)');
title('功率谱');

2、结果分析

由于已经封装成一个函数,按照函数说明进行参数设置即可完成幅度调制的功能。波形图如下:

AM_func(3000,48e3,300,0.5,8e3,1,0.3)

通信原理实践(二)——幅度调制

二、频率谱分析——计算AM调制效率

1、代码展示:

代码需要放在上面代码的后面,或者基于上者的数据

u_pow_decsend = sort( abs(u),'descend');
p_all = sum( abs(u).^2 );
p_main = p_all - max(abs(u)).^2 * 2;
p_n = p_main / p_all * 100;

2、结果分析:

一般来说AM调制最大调制效率为33%,且是在Ma=1的情况下取得。

三、非相干解调——包络检波

1、代码展示

代码需要放在上面代码的后面,或者基于上者的数据

b = fir1(80, 2000/(fs/2));               % 设计截止频率为2KHz的FIR低通滤波器,如何截止频率设为300,则无法显示600Hz的波形,即可以检验出过调制的波形
s_t = s_t .* (s_t > 0);
filtered_t = filter(b,1,s_t); % 滤波
demod_t = filtered_t - sum(filtered_t)/length(filtered_t); subplot(1,3,3);
plot(t, demod_t , 'b'); % 绘制解调波形
xlabel('time');
ylabel('amplitude');
title('包络检波之后的波形');

2、结果分析

这边的滤波器需要设计好,两个参数的后者的范围是0~1,表示频率对于fs/2的归一化频率。

通信原理实践(二)——幅度调制

四、相干解调

1、代码展示

已封装成函数,可直接调用

function [] = jt_func( N,fs,fm,Am,fc,Ac,Ma,fc_c,phi_c )
%UNTITLED5 此处显示有关此函数的摘要
% 此处显示详细说明
% %发送端参数声明
% N = 3000; % 长度
% fs = 48e3; % 采样率
% fm = 300; % 调制信号频率
% Am = 0.5; % 调制信号幅度
% fc = 8e3; % 载波频率
% Ac = 1; % 载波幅度
% Ma = 10; %
% %接收端参数声明
% fc_c = 8e3; % 本地振荡器频率
% phi_c = 0; % 本地振荡器初相 %%
%发送端产生波形
A0 = Am / Ma;
t = (0:N-1)'/fs; % 时间t
m_t = Am*sin(2*pi*fm*t); % 产生调制信号
c_t = Ac*cos(2*pi*fc*t); % 产生载波信号
s_t = (m_t + A0) .* c_t; % 调制 %%
%接收端解调波形
r_t = s_t;
figure()
subplot(1,3,1);
plot(t, r_t , 'b'); % 绘波形
xlabel('time');
ylabel('amplitude');
title('时域波形'); c_t = cos(2*pi*fc_c*t+phi_c); % 本地振荡器波形
mix_t = r_t .* c_t; % 混频
subplot(1,3,2);
plot(t, mix_t , 'b'); % 绘制混频波形
xlabel('time');
ylabel('amplitude');
title('混频后的波形'); b = fir1(20, fc/(fs/2)); % 设计截止频率为fcFIR低通滤波器
filtered_t = filter(b,1,mix_t); % 滤波
demod_t = filtered_t - sum(filtered_t)/length(filtered_t);
% 去除DC分量
subplot(1,3,3);
plot(t, demod_t ,'b' ); % 绘制解调波形
xlabel('time');
ylabel('amplitude');
title('滤波后的波形'); end

2、结果分析

通信原理实践(二)——幅度调制

五、模拟幅度调制、解调过程——音频

1、代码如下:

主要是播放声音这个函数使用audioplayer,可以停止声音播放,否则。。。。。

clc,clear;

[m_t,fs] = audioread('MayDay.mp3');   % 得到音频信号以及采样率

N = length(m_t);                      % 调制信号长度
Am = max(m_t); % 音频信号最大幅度
fc = 200e3; % 载波频率
Ac = 1; % 载波幅度
Ma = 0.3;
A0 = Am / Ma; t = (0:N-1)'/fs; % 时间t
c_t = Ac*cos( 2*pi*fc*t ); % 产生载波信号
s_t = (m_t + A0) .* c_t; % 调制 r_t = s_t;
figure()
subplot(1,3,1);
plot(t, r_t , 'b'); % 绘波形
xlabel('time');
ylabel('amplitude');
title('时域波形'); c_t = cos(2*pi*fc*t); % 本地振荡器波形
mix_t = r_t .* c_t; % 混频
subplot(1,3,2);
plot(t, mix_t , 'b'); % 绘制混频波形
xlabel('time');
ylabel('amplitude');
title('混频后的波形'); b = fir1(20,fs/4/(fs/2)); % 设计截止频率为fcFIR低通滤波器
filtered_t = filter(b,1,mix_t); % 滤波
demod_t = filtered_t - sum(filtered_t)/length(filtered_t);
% 去除DC分量
subplot(1,3,3);
plot(t, demod_t ,'b' ); % 绘制解调波形
xlabel('time');
ylabel('amplitude');
title('滤波后的波形'); player = audioplayer(demod_t, fs);
play(player);

2、结果展示:

通信原理实践(二)——幅度调制

上一篇:自定义javascript log方法


下一篇:在Asp.Net的Global.asax中Application_Error跳转到自定义错误页无效的解决办法