%FM信號的產(chǎn)生
clear;clc;
fs=500000;%采樣頻率
dt=1/fs;
t=0:dt:0.002;
A_m=2;A_c=5;??%振幅
SNR=15;????%信噪比

fm=1e3;fc=2*1e4;
kf=2*pi*2.5*1e3;%調頻靈敏度
u_c_t=A_c*cos（2*pi*fc*t）;???%載波信號

%調制信號
m_t=A_m*cos（2*pi*fm*t）;????
%m_t=randn（1length（u_c_t））;

%求u_t的積分
int_m_t（1）=0;
for?i=1:length（t）
???int_m_t（i）=sum（m_t（1:i））*dt;
end

u_FM_t=A_c*cos（2*pi*fc*t+kf*int_m_t）;%FM信號