小弟想计算一个空间曲柄摇杆机构输出件的角位移和时间的关系，得到了下面的公式,n为输出件的角位移取值范围0-90°，
（50*cos（n）-30*sin（360*t））^2+(30*cos(360*t)+45)^2=150^2-(165-50sin(n))^2
用ezplot画出了n和t的曲线，现在想求n对t的一节导数和t的关系即角速度和时间的关系，以及n对t的二阶导数和时间的关系，即角加速度和时间的曲线，请各位大神帮帮忙，小弟很急，第一次用matlab，请高手帮帮忙！