tulun.m

来自「机械原理的上机作业,电气学院的,应该有用的啊」· M 代码 · 共 55 行

function tulun3(R,H)
syms S V A J
for j=0:0.01:pi*2
   if j<pi/2
        S=H/2*(1-cos(pi/(pi/2)*J));
        V=diff(S,J);
        A=diff(V,J);
        SS=subs(S,J,j);
        VV=subs(V,J,j);
        AA=subs(A,J,j);
   elseif j<pi
        S=H;
        V=0;
        A=0; 
        V=diff(S,J);
        A=diff(V,J);
        SS=subs(S,J,j);
        VV=subs(V,J,j);
        AA=subs(A,J,j);
   elseif j<3*pi/2
        S=H-H/2*(1-cos(pi/(pi/2)*J));
        V=diff(S,J);
        A=diff(V,J);
        SS=subs(S,J,j);
        VV=subs(V,J,j);
        AA=subs(A,J,j);
   else 
        S=0;
        V=0;
        A=0;
        V=diff(S,J);
        A=diff(V,J);
        SS=subs(S,J,j);
        VV=subs(V,J,j);
        AA=subs(A,J,j);
        
    end

grid on;
subplot(2,2,1);
hold on
polar(j,SS+R);
hold on
polar(j,R);
axis equal
subplot(2,2,2);
hold on
plot(20*j,VV,'c');
subplot(2,2,3);
hold on
plot(20*j,AA,'g');
subplot(2,2,4);
hold on
plot(20*j,SS,'m');
end