频率调制制 matlab 实现

t=0:0.000001:0.004;
f=2000;
fc=20000;
fs=1000000;
m1=cos（2*pi*f*t）;
w1=0;w2=0;
for m=1:length（t）
    w1=m1（m）+w2;
    w2=m1（m）+w1;
    fi（m）=w1/（2*fs）;
end
fi=fi*2*pi/max（abs（fi））;
Kfm=1;
I=cos（Kfm*fi）;
Q=sin（Kfm*fi）;
A=1;
y0=A*cos（2*pi*fc*t）.*I-A*sin（2*pi*fc*t）.*Q;
subplot（321）
plot（ty0）;
xlim（[0.0010.004]）;
xlabel（‘t/s‘）;ylabel（‘FM时域波形‘）;

Y0=fft（y0）;
f1=（0:4000）*fs/4001-fs/2;
subplot（322）
plot（f1abs（fftshift（Y0）））;
xlim（[-4000040000]）;
xlabel（‘f/Hz‘）;ylabel（‘FM频谱‘）;

snr=30;
y1=awgn（y0snr）;
subplot（323）;
plot（ty1）;
hold on 
subplot（323）;
plot（ty0）;
xlim（[0.0010.004]）;
xlabel（‘t/s‘）;ylabel（‘FM加噪后信号‘）;
Y1=fft（y1）;
f1=（0:4000）*fs/4001-fs/2;
subplot（324）;
plot（f1fftshift（abs（Y1）））;
xlim（[-4000040000]）;
xlabel（‘f/HZ‘）;ylabel（‘F