FM信号调制信号和载波

Posted 2022-06-24 studyer_domi

1、内容简介

略
353-可以交流、咨询、答疑

2、内容说明

略

3、仿真分析

%已调信号为FM信号，该调频信号由1KHz的调制信号和10kHz的载波而成
clear all;            %初始化
clc;
close all;
f=5000;               %载波频率  5KHz
fs=1000000;            %采样频率
N=5000;               %采样点数            
Ts=1/fs;                %采样周期
t=(0:Ts:(N-1)*Ts);       %仿真时间
f1=500;               %调制信号频率  500Hz
msg=sin(2*pi*f1*t);               %调制信号
figure,plot(msg);
figure,plot(cumsum(msg));
kf=0.0716;                       %调频灵敏度
Signal1=exp(j*(2*pi*f*t));          %载波
Signal=exp(j*(2*pi*f*t+2*pi*kf*cumsum(msg))); 
figure,plot(imag(Signal1));
figure,plot(real(Signal1));
%已调信号，FM信号，cumsum为累加函数；

phi_hat(1)=10;        %已调信号（锁相环）输入初值 10
e(1)=0;               %频差初值
phd_output(1)=0;      %解调信号输出初值
vco(1)=0;             %压控振荡器输出初值
kp=0.125;              %环路滤波器参数
ki=0.125;               %环路滤波器参数

%锁相环


startplot=1;                %绘图开始位置
endplot=2000;              %绘图结束位置


figure(1);                  %在窗口1绘图
subplot(2,2,2);              %多个图绘到一个平面
plot(t(startplot:endplot),msg(startplot:endplot));         %绘图函数，横坐标为时间t
title('调制信号（频率1KHz）');                     %图标标题
xlabel('时间（seconds）');                          %横坐标注释
ylabel('幅度');                                     %纵坐标注释
grid;
figure(1);
subplot(2,2,1);
plot(t(startplot:endplot),real(Signal1(startplot:endplot)));
title('载波');
ylabel('幅度');
grid;

figure(1);
subplot(2,2,3);
plot(t(startplot:endplot),real(Signal(startplot:endplot)));
title('FM已调信号');
ylabel('幅度');
grid;


figure(1);
subplot(2,2,4);
plot(t(startplot:endplot),e(startplot:endplot));
title('解调信号');
ylabel('幅度');
grid;

[f_Signal,y_Signal]=change_fft(fs,Signal);
[f_msg,y_msg]=change_fft(fs,msg);
[f_e,y_e]=change_fft(fs,e);
figure
plot(f_Signal,y_Signal)
xlabel 频率/Hz
ylabel 幅值
title Signal
xlim([0 100000])
figure
plot(f_msg,y_msg)
xlabel 频率/Hz
ylabel 幅值
title msg
xlim([0 5000])
figure
plot(f_e,y_e)
xlabel 频率/Hz
ylabel 幅值
title e
xlim([0 5000])

[f_Signal_imag,y_Signal_imag]=change_fft(fs,imag(Signal));
figure
plot(f_Signal_imag,y_Signal_imag)
xlabel 频率/Hz
ylabel 幅值
title Signal虚部
xlim([0 100000])

 

 

 

 

4、参考论文

略