%《MATLAB及其在理工課程中的應用指南》第八章例8－4－3程序
% 調幅、調頻、調相波形及頻譜
% 西安電子科技大學出版社出版  陳懷琛編著 1999年10月初版，2004年11月第二版
%
clf,format compact,t=0:0.001:1;
A0=10,w0=100,phi0=0
x=A0.*cos(w0*t+phi0);                   % 原始信號表達式
subplot(3,1,1),plot(t,x),hold on	    % 未調制波形的繪制
set(gcf,'color','w')                    % 置圖形背景色為白色
xlabel('t'),ylabel('載頻信號')
dA=input('振幅增量=(三次運行,書上分別取 5*sin(0.1*w0*t), 0, 0) ');
dw=input('頻偏=(三次運行,書上分別取 0, 5*sin(0.1*w0*t, 0) ');
dphi=input('相移=(三次運行,書上分別取 0, 0, 5*sin(0.1*w0*t)) ');
y=(10+dA).*cos((w0+dw).*t+dphi);        % 經綜合調制后的信號表達式
subplot(3,1,2),            
plot(t,y),ylabel('調制波形')    		% 調制后波形的繪制
subplot(3,1,3),
plot(t,y-x),ylabel('y-x')    		    % 調制前后波形之差的繪制
pause,subplot(3,1,1),
plot(t,y,'r')                   		% 把調制后波形疊加在原來波形上