2011上机练习0203MATLAB编程基础训练.docVIP

上机二 运动分析的矩阵方法 一、上机目的 根据平面机构运动分析的矩阵方法，应用matlab编程求解给定四杆机构的位移、速度和加速度线图。通过上机调试，初步奠定学生应用matlab进行机构运动分析的基础。 二、上机地点及要求 1、上机地点: 教A楼网络中心3机房。 2、上机前必须编好程序，上机调试通过以后完成指定的计算实例。 3、编程语言限Matlab。 4、上交源程序打印稿，并附程序使用说明和实际算例的计算结果抓图。 三、编程要求 1.编程进行图示四杆机构的运动分析，已知：。 ? 要求：取的步长为1°，编程求解对应的杆位移、速度及加速度的值，并应用绘图指令绘制出相应的线图。 选做（如完成有适当的加分）： 应用matlab文件IO指令，将求得的数据按格式写入文件中。 绘制出P点的运行轨迹。 编制机构的运动仿真程序（需要自学相关内容） 2，部分参考代码： function y=fsolve_ex2() clc % clear the screen L=[30 60 90 100];%给定一组四杆机构的杆长 omig1=20;%omig1=20rad/s x0 = [d2r(30); d2r(120)]; % 猜测一个初始点，可以借助图解法 options=optimset(Display,off); % Option to display output [x,fval] = fsolve(@myfun,x0,options,0,L) % 循环前调用一次，用于获得较好的猜测点 for i=1:361 seit1(i)=d2r(i-1);%将seit1转换成弧度 x0=x; [x,fval] = fsolve(@myfun,x0,options,seit1(i),L); seit2(i)=x(1); seit3(i)=x(2); end h=plot(r2d(seit1),r2d(seit2)) set(h,LineWidth,2,color,red,LineStyle,--); hold on h=plot(r2d(seit1),r2d(seit3)); set(h,LineWidth,2,color,cyan,LineStyle,-); grid on legend(\itseit3--seit1,\itseit3--seit1) % 根据结果求解线性方程组，实现速度分析 for i=1:361 % 根据已知的seit1，写出对应的系数矩阵和常数矩阵，求线性方程组 % 并将对应的连杆2、3的角速度保存在相应的变量数组中 A=[-L(2)*sin(seit2(i)) L(3)*sin(seit3(i)) L(2)*cos(seit2(i)) -L(3)*cos(seit3(i))]; b=[omig1*L(1)*sin(seit1(i));-omig1*L(1)*cos(seit1(i))]; w=A\b; omig2(i)=w(1); omig3(i)=w(2); end % 绘制对应的函数曲线 figure %打开一个新的图形窗口 h=plot(r2d(seit1),omig2) set(h,LineWidth,2,color,blue,LineStyle,-.); hold on h=plot(r2d(seit1),omig3); set(h,LineWidth,2,color,black,LineStyle,-); grid on legend(\itomig2--seit1,\itomig3--seit1) % 根据结果求解线性方程组，实现加速度分析并绘图，由同学们上机实验时自 %行完成 %求p点的轨迹，并绘图表示，选做，建议进行研究并完成。 function F = myfun(x,seit1,L) % seit1代表每一个时刻对应曲柄的转角θ1 % L(1),L(2),L(3),L(4)分别代表各杆长l1,l2,l3,l4 % x(1)代表θ2，x(2)代表θ3 F(1)=L(2)*cos(x(1))-L(3)*cos(x(2))+L(1)*cos(seit1)-L(4); F(2)=L(2)*sin(x(1))-L(3)*sin(x(2))+L(1)*sin(seit1); function y=d2r(x) y=x*pi/180; function y=r2d(x) y=x*180/pi; 3.参考的运行结果： 上机三 优化工具箱使用训练 一、上机目的 熟悉matlab优化工具箱的使用方法，会根据给定的优化模型，设定适当的参数调