<P>头都大了呵呵<BR>syms a t %主程序<BR>t=linspace(0,2*pi,361);<BR>a=[0:1:89];<BR>for a=[0:1:89]<BR>uA=source(t,a);<BR>uB=shift(uA,t);<BR>uC=shift(uB,t);<BR>uAO=(2*uA-uB-uC)/3;<BR><BR>function y=shift(u,x) %120度移相<BR>x1=(x/2/pi)*360;<BR>y1=u(x1<=120);<BR>y2=u(x1>120&x1<=360);<BR>y=[y2 y1];<BR><BR><BR>function y =source(x,a) %产生A相波形,是个360度内的分段函数<BR>x1=(t/2/pi-floor(t/pi/2))*360;<BR>uA(x1<=a)=0;<BR>uA(x1>a&x1<(180-a))=1;<BR>uA(x1>=(180-a)&x<=180)=0;<BR>uA(x1>180&x1<(180+a))=0;<BR>uA(x1>(180+a)&x1<(360-a))=-1;<BR>uA(x1>=(360-a)&x<=360)=0;<BR><BR><BR><BR>如何求uAO的有效值?</P>