fujinhai
2009-03-01, 16:53
以下是我的一个方程:
t=[0:0.01:1.5];y=acos(((55*t/3-55/(2*pi)*sin(2*pi*t/3)+27.5).^2+(2*t/3-(1/(2*pi))*sin(4*pi*t/3)+21).^2+48.75)/(22*sqrt((55*t/3-55/(2*pi)*sin(2*pi*t/3)+27.5).^2+(2*t/3-(1/(2*pi))*sin(4*pi*t/3)+21).^2)))+atan((55*t/3-(55/(2*pi))*sin(2*pi*t/3)+27.5)/(2*t/3-(1/(2*pi))*sin(4*pi*t/3)+21));plot(t,y);
画出图是一条多点的直线,代入数之后发现acos中的值是大于一的,所以会出现这种情况。按道理来说这应该是错误的,因为余弦值最大才为1,这种情况下的matlab是怎样计算出来的数值,虚部实部又分别代表什么,希望高手回答,谢谢。
t=[0:0.01:1.5];y=acos(((55*t/3-55/(2*pi)*sin(2*pi*t/3)+27.5).^2+(2*t/3-(1/(2*pi))*sin(4*pi*t/3)+21).^2+48.75)/(22*sqrt((55*t/3-55/(2*pi)*sin(2*pi*t/3)+27.5).^2+(2*t/3-(1/(2*pi))*sin(4*pi*t/3)+21).^2)))+atan((55*t/3-(55/(2*pi))*sin(2*pi*t/3)+27.5)/(2*t/3-(1/(2*pi))*sin(4*pi*t/3)+21));plot(t,y);
画出图是一条多点的直线,代入数之后发现acos中的值是大于一的,所以会出现这种情况。按道理来说这应该是错误的,因为余弦值最大才为1,这种情况下的matlab是怎样计算出来的数值,虚部实部又分别代表什么,希望高手回答,谢谢。