%《MATLAB及其在理工課程中的應(yīng)用指南》(十一五規(guī)劃版)第八章例8－3－1程序
% 可控硅點火角與電壓的關(guān)系
% 西安電子科技大學出版社出版  陳懷琛編著 2007年7月
%
% 計算負載上得到的有效值電壓與點火角(1 到 179 度)的函數(shù)關(guān)系,并畫出曲線
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
% 畫負載上的有效值電壓與點火角關(guān)系曲線：
figure(1)
plot([1:180],rms,'Linewidth',2.0);
legend('\bf負載電壓有效值(V)');grid on;
% gtext('\bf點火角（度)');pause

%畫點火角45度時負載上的電壓波形：
figure(2)
plot(wt,waveform45, 'Linewidth',2.0);
grid on;set(gcf,'color','w')
% 下一條語句中的\omega和\itt說明如何同時標注希臘字母和漢字
gtext('\bf\omega\itt\rm (弧度)');