%《MATLAB及其在理工課程中的應用指南》第八章例8－3－1程序
% 可控硅點火角與電壓的關系
% 西安電子科技大學出版社出版  陳懷琛編著 1999年10月初版，2004年11月第二版
%
% 計算負載上得到的有效值電壓與點火角(1 到 179 度)的函數關系,并畫出曲線
wt = [1:180]*pi/180; 				            % 把半波分割為180份 
volts = 220 * sqrt(2) * sin(wt);		        % 完整的半波波形
for ii = 1:180				                    % 求不同點火角ii時的波形
    waveform = [zeros(1,ii-1),volts(ii:180)];   % 求不同點火角ii時的波形
    if ii == 45   waveform45 = waveform; 		% 記錄點火角為45度時的波形
    elseif ii == 90    waveform90 = waveform; 	% 記錄點火角為90度時的波形
    end
    temp = sum(waveform.^2); 		            % 計算各點波形的平方積分：
    rms(ii) = sqrt(temp/180); 	                % 計算積分的均方根：
end
% 畫負載上的有效值電壓與點火角關系曲線：
figure(1)
plot([1:180],rms,'Linewidth',2.0);
set(gcf,'color','w')                            % 置圖形背景色為白色
legend('\bf負載電壓有效值(V)');grid on;
gtext('\bf點火角（度)');pause

%畫點火角45度時負載上的電壓波形：
figure(2)
plot(wt,waveform45, 'Linewidth',2.0);
set(gcf,'color','w')                            % 置圖形背景色為白色
grid on;
gtext('\bf\omega\itt\rm (弧度)');