MatlabAssignment2.docVIP

Matlab Assignment 2 学号：05101030318 班级：机05A-3 姓名：任海艳 问题一 利用Matlab解决物理学中的速度加速度问题——求导和画图象问题 1．具体问题：某质点的运动方程为x=6t3-9t2+52t-3试求：该质点的速度和加速度表达式以及在t=3s时的速度和加速度大小 2． 建模过程 （1）由质点的运动方程为：x=6t3-9t2+52t-3，求速度和加速度。应用物理学原理分别对时间t求一阶导数和二阶导数； 即v=x′=dx/dt a=x′′=d2x/dt2 （2）t=3s，5s，6s时刻的速度和加速度直接代入即可 （3）利用Matlab进行求导计算 %已知支点运动方程为x=6t3-9t2+52t-3，求其对时间t的一阶导数和二阶导数 syms t %声明变量 f=sym(6*t^3-9*t^2+52*t-3) %定义一个函数f v=diff(f) %对时间求得一阶导数即是速度 a=diff(f,2) %对时间求得二阶导数即是加速度 subs(v,t,3) %代入t=3，求得速度值 subs(a,t,3) %代入t=3，求得加速度值 f = 6*t^3-9*t^2+52*t-3 v = 18*t^2-18*t+52 a = 36*t-18 ans = 160 ans = 90 3. 结论： （1）该质点的速度和加速度表达式为：v =18*t^2-18*t+52 a =36*t-18 在t=3s时，该质点的速度和加速度的大小分别为160m/s 和 90m/s2 （2）做题感受及遇到问题： 这道题是微分问题，总体来说不难，只是容易与积分等问题弄混淆，我开始时用inline定义函数，但diff执行不了，这是一个问题；另外，我想引用function，这样可以求得任意时刻质点的速度和加速度，但在Word里如何调用，我还不清楚，由于时间比较紧，我没有解决这个问题 问题二 利用Matlab求解四连杆的角加速度——解线性方程 1．具体问题：四连杆如图（1）所示,图中未知量如下表（2）所示,设输入连杆2的角速度为100rad/s，试求图示位置时连杆3、4的角速度。 图（1） 表（2） 四连杆机构的长度和转角 连杆 长度（R）/cm 转角（th）/（o） 1 12.0 0 2 4.0 45 3 10.0 24.65 4 7.0 90.68 2．建摸过程 （1）由图（1）画出四连杆机构的闭环矢量图，如图（3）所示 （2）利用矢量夹角的正弦和余弦函数可以得闭环矢量方程及其两个分量式如下： R2+R3=R1+R4 r2cos(th2)+r3cos(th3)=r1cos(th1)+r4cos(th4) r2sin(th2)+r3sin(th3)=r2sin(th2)+r4sin(th4) （3）分别对时间求一阶导数得（其中th1=0） -w2 r2sin(th2)-w3 r3sin(th3)=-w4 r4sin(th4) w2 r2cos(th2)+w3 r3cos(th3)=w4 r4cos(th4) （4）将上式矩阵形式 [-r3sin(th3) r4sin(th4)] （5）利用Matlab编程求解w3 和w4，编程如下 %线性代数：根据已知方程及系数计算矩阵的解 %输入已知量的值 r2=4; r3=10; r4=7; %将角度转换为弧度 th2=45*pi/180;