[求助]请帮我看下积分程序哪错了

cycr1234 · 发表于 2006-3-28 22:48

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

a=0; u=invert2(t,a); q = quad(@invert2,0,pi/6,1.0e-6,0,a) plot(t,u) function y=invert2(x,a) x1=(x/2/pi)*360;
uA(x1<=a)=0; uA(x1>a&x1<(180-a))=1; uA(x1>=(180-a)&x<=180)=0; uA(x1>180&x1<(180+a))=0; uA(x1>(180+a)&x1<(360-a))=-1; uA(x1>=(360-a)&x<=360)=0;
uB1=uA(x1<120); uB2=uA(x1>=120&x1<=360); uB=[uB2 uB1];
uC1=uB(x1<120); uC2=uB(x1>=120&x1<=360); uC=[uC2 uC1];
uAO=(2*uA-uB-uC)/3; y=uAO.^2; 没有出错 可是结果是 q=0

[此贴子已经被aspen于2006-3-30 9:28:11编辑过]

cycr1234 · 发表于 2006-3-28 23:10

图：

登录/注册后可看大图

cycr1234 · 发表于 2006-3-29 14:31

不好意思，漏了一行 t=linspace(0,2*pi,361); 不过q=0 还是让我不知道怎么得来的

happy · 发表于 2006-3-29 20:52

<DIV class=quote>以下是引用cycr1234在2006-3-28 22:48:54的发言： 
a=0; u=invert2(t,a); q = quad(@invert2,0,pi/6,1.0e-6,0,a) plot(t,u) function y=invert2(x,a) x1=(x/2/pi)*360;
uA(x1<=a)=0; uA(x1>a&x1<(180-a))=1; uA(x1>=(180-a)&x<=180)=0; uA(x1>180&x1<(180+a))=0; uA(x1>(180+a)&x1<(360-a))=-1; uA(x1>=(360-a)&x<=360)=0;
uB1=uA(x1<120); uB2=uA(x1>=120&x1<=360); uB=[uB2 uB1];
uC1=uB(x1<120); uC2=uB(x1>=120&x1<=360); uC=[uC2 uC1];
uAO=(2*uA-uB-uC)/3; y=uAO.^2; 没有出错 可是结果是 q=0 </DIV>
红色部分表示什么意思？

happy · 发表于 2006-3-29 21:02

不行你就用if吧

cycr1234 · 发表于 2006-3-30 10:21

uB1=uA(x1<120); uB2=uA(x1>=120&x1<=360); 是循环移位，移120度相位 明明画出的图没错，就是我想要的图 可是对它求积分怎么可能是零嘛，

happy · 发表于 2006-3-30 18:47

<DIV class=quote>以下是引用cycr1234在2006-3-30 10:21:52的发言： uB1=uA(x1<120); uB2=uA(x1>=120&x1<=360); 是循环移位，移120度相位 明明画出的图没错，就是我想要的图 可是对它求积分怎么可能是零嘛， </DIV>
这个语句我还是没搞懂，你可以随便输入一个数值看看，返回的肯定都是0 例如invert2(1,0)

账号		自动登录	找回密码
密码			我要加入

[求助]请帮我看下积分程序哪错了

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（cycr1234）[求助]请帮我看下哪错了

回复：（cycr1234）[求助]请帮我看下哪错了

回复：（cycr1234）uB1=uA(x1<120);uB2=uA(x1>=...