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]) 
