t=0:0.05:8;
r=37;L=140.5;beta1=pi/36;beta2=pi/18;beta3=pi/9;
e1=((1/cos(beta1))-1)*r/2;
e2=((1/cos(beta2))-1)*r/2;
e3=((1/cos(beta3))-1)*r/2;
B1=(r-e1-2*e1*cos(2*pi*t/2))*cos(pi*t/2)*sin(beta1)+e1*cos(3*pi*t/2)*sin(beta1);
B2=(r-e2-2*e2*cos(2*pi*t/2))*cos(pi*t/2)*sin(beta2)+e2*cos(3*pi*t/2)*sin(beta2);
B3=(r-e3-2*e3*cos(2*pi*t/2))*cos(pi*t/2)*sin(beta3)+e3*cos(3*pi*t/2)*sin(beta3);
plot(t,B1,'k',t,B2,'b',t,B3,'r');
结果出现如下错误:
??? Error using ==> mtimes
Inner matrix dimensions must agree. 后来我知道是t=0:0.05:8;
>> C=cos(pi*t)*cos(pi*t/2);这个表达式matlab算不出来,是不是需要改编程,望高手指点,能将上面的B1,B2等计算出来,先在此谢过 |