%《MATLAB在電子信息課程中的應用》第五章例5.9程序q509
%  非正弦信號的頻譜分析和濾波
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

clear, format compact
L=0.4;C=10e-6;R=200;					% 輸入元件參數
Um=100;w1=100*pi;
% 輸入只有偶次諧波，故只分析偶次諧
N=input('需分析的諧波次數2N=（鍵入偶數）  ');
n=1:N/2;w=[eps,2*n*w1];					% 設定頻率數組
Us = 4*Um/pi*[0.5,-1./(4*n.^2-1)];	% 輸入信號頻譜數組
z1=j*w*L;z2=1./(j*w*C);z3=R;		% 計算阻抗數組
z23=z2.*z3./(z2+z3);			% 注意為何用元素群運算符
Ur=Us.*z23./(z1+z23)			% 求Ur
disp('  諧波次數   諧波幅度  諧波相移(度)')
disp([2*[0,n]',abs(Ur)',angle(Ur)'*180/pi])