這是西安電子科技大學《matlab及其在理工課程中的應用》第174頁的程式：

clf,formatcompact,t=0:.001:1;

A0=10;w0=100;phi0=0;

x=A0*sin(w0*t-phi0);subplot(3,1,1),plot(t,x),hold

xlabel("t"),ylabel("載頻訊號")

dA=input("振幅增量=");dw=input("頻偏=");dphi=input("相移=");

y=(10+dA).*sin((w0+dw).*t+dphi);

subplot(3,1,2),

plot(t,y),ylabel("調製波形")

subplot(3,1,3);

Y=fft(y);plot(abs(Y)),grid

axis([0,50,0,5000])

1.原先的w0=100,要修改；

2.程式執行後，在命令提示處輸入5sin(200t)；

3.如果不明白，下載那本書吧。