%《MATLAB在電子信息課程中的應用》第六章例6.7程序q607
% 方波表為多階正弦波之和
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

t = 0:.01:2*pi;					% 設定一個時間數組，有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，重編程。
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, figure(1),plot(t,y(1:9,: )),grid
line([0,pi+0.5],[pi/4,pi/4])			% 加上方波幅度線及標注
text(pi+0.5,pi/4,'pi/4')
%      將各半波形繪成三維網格圖，看出增加諧波階次對方波逼近程度的影響
halft=ceil(length(t)/2);
pause, figure(2),mesh(t(1:halft),[1:10],y(:,1:halft)), pause	%只用正半周波形
clc