%《MATLAB及其在理工課程中的應用指南》第九章例9－1－6程序
% 方波表為多階正弦波之和
% 西安電子科技大學出版社出版  陳懷琛編著 1999年10月初版，2004年11月第二版
%
t = 0:.1:10;						        % 設定一個時間數組，有101個點
y = sin(t);	plot(t,y),figure(gcf),pause     % 頻率為w=1(f=1/2π)的正弦基波
y = sin(t) + sin(3*t)/3; plot(t,y), pause 	% 疊加三次諧波
%  用1，3，5，7，9次諧波疊加
y = sin(t) + sin(3*t)/3 + sin(5*t)/5 + sin(7*t)/7 + sin(9*t)/9;plot(t,y)
% 為了繪制三維曲面，要把各次波形數據存為一個三維數組，因此必須
% 重新定義y，重編程.由于擬求至19次諧波，把點取密一些。
t = 0:.031:3.14;	
y = zeros(10,max(size(t))); x = zeros(size(t));
for k=1:2:19
    x = x + sin(k*t)/k; y((k+1)/2,: ) = x;
end                                     
%  將各波形迭合繪出
pause, plot(t,y(1:9,: ))
set(gcf,'color','w')                        % 置圖形背景色為白色
%  將各波形繪成三維網格圖，看出增加諧波階次對方波逼近程度的影響
pause, mesh(t,[1:10],y), pause
clc