[求助]请帮我看下积分程序哪错了
<P>a=0;<br>u=invert2(t,a);<br>q = quad(@invert2,0,pi/6,1.0e-6,0,a)<br>plot(t,u)<br><br><br>function y=invert2(x,a)<br>x1=(x/2/pi)*360;</P><P>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;</P>
<P>uB1=uA(x1<120);<br>uB2=uA(x1>=120&x1<=360);<br>uB=;</P>
<P>uC1=uB(x1<120);<br>uC2=uB(x1>=120&x1<=360);<br>uC=;</P>
<P>uAO=(2*uA-uB-uC)/3;<br>y=uAO.^2;<br><br><br>没有出错<br>可是结果是<br>q=0<br><br></P>
[此贴子已经被aspen于2006-3-30 9:28:11编辑过]
图：<BR>http://classpic2.chinaren.com/album/pic/864/850/2917034545011850864/114355853950845.jpg <P>不好意思，漏了一行<BR>t=linspace(0,2*pi,361);<BR><BR>不过q=0<BR>还是让我不知道怎么得来的</P>
回复：（cycr1234）[求助]请帮我看下哪错了
<DIV class=quote><B>以下是引用<I>cycr1234</I>在2006-3-28 22:48:54的发言：</B><BR><P>a=0;<BR>u=invert2(t,a);<BR>q = quad(@invert2,0,pi/6,1.0e-6,0,a)<BR>plot(t,u)<BR><BR><BR>function y=invert2(x,a)<BR>x1=(x/2/pi)*360;</P>
<P>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;</P>
<P><FONT color=#ff0000>uB1=uA(x1<120);<BR>uB2=uA(x1>=120&x1<=360);<BR></FONT>uB=;</P>
<P><FONT color=#ff0000>uC1=uB(x1<120);<BR>uC2=uB(x1>=120&x1<=360);<BR></FONT>uC=;</P>
<P>uAO=(2*uA-uB-uC)/3;<BR>y=uAO.^2;<BR><BR><BR>没有出错<BR>可是结果是<BR>q=0<BR><BR></P></DIV>
<P>红色部分表示什么意思？</P>
回复：（cycr1234）[求助]请帮我看下哪错了
不行你就用if吧 <FONT color=#ff0000>uB1=uA(x1<120);<BR>uB2=uA(x1>=120&x1<=360);</FONT><BR><BR>是循环移位，移120度相位<BR><BR>明明画出的图没错，就是我想要的图<BR>可是对它求积分怎么可能是零嘛，<BR>回复：（cycr1234）uB1=uA(x1<120);uB2=uA(x1>=...
<DIV class=quote><B>以下是引用<I>cycr1234</I>在2006-3-30 10:21:52的发言：</B><BR><FONT color=#ff0000>uB1=uA(x1<120);<BR>uB2=uA(x1>=120&x1<=360);</FONT><BR><BR>是循环移位，移120度相位<BR><BR>明明画出的图没错，就是我想要的图<BR>可是对它求积分怎么可能是零嘛，<BR></DIV><P>这个语句我还是没搞懂，你可以随便输入一个数值看看，返回的肯定都是0<BR><BR>例如invert2(1,0)</P>
页:
[1]