haukuai.m

机械原理的上机作业,电气学院的,应该有用的啊

function huakuai(A,B,D,R0,DR,S0,N,E)
x=0;
while x==0
dr=DR/N;
R=zeros(1,N+1);S=zeros(1,N+1);Q=zeros(1,N+1);V=zeros(1,N+1);
for i=0:N
    R(i+1)=R0+i*dr;
    S(i+1)=A*(sin(R(i+1)*pi/180)-sin(R0*pi/180))+B*(sqrt(1-((A*cos(R0*pi/180)-D)/B)^2)-sqrt(1-((A*cos(R(i+1)*pi/180)-D)/B)^2));
    Q(i+1)=i*S0/N;
    V(i+1)=abs(S(i+1)-Q(i+1))/S0;
    if abs(V(i+1))>E
        x=0;
        disp('错误，超出线性范围，请重新输入!');
        break
    else 
        x=1;
    end
end  
end
hold on;
plot(R,S,'m--');
hold on;
plot(R,Q);
legend('实际图','理想图');