精通MATLAB语言，有MATLAB编程问题的朋友，请直接联系我！
联系方式：
QQ：993878382
Email：993878382@https://www.doczj.com/doc/ed7594153.html,
MATLAB博客：https://www.doczj.com/doc/ed7594153.html,/46036333_d.html
更多学习资料和电子书，请直接访问我的博客！
如果我QQ不在线，请将问题发到我的邮箱，或者直接给我QQ留言，第一时间答复你！简单问题请直接留言，复杂问题可以提供解决思路！同探讨，共进步！
还可以为公司和科研单位设计各类算法，欢迎前来咨询！
郑重声明：本资源来源于网络，仅限交流使用，请勿用于商业用途！
自定义HTML载入中...
[和讯博客]
[发表文章]
[博客设置] [文章管理]
搜索
载入中
快速链接如何在MATLAB里绘制同心圆 源程序 [转贴 2010-01-11 00:01:20]
标签: MATLAB 源程序 绘制 同心圆 何在 .
分类: MATLAB编程 %绘制同心圆的函数
function circle(R)
theta=0:0.01:2*pi;
x=R*sin(theta);
y=R*cos(theta); plot(x,y) axis equal hold on end %主程序调用 clear;clc for R=1:2:20
circle(R)
end
下一篇: MATLAB waitbar 设置取消按钮 源程序
上一篇: 无穷个解的线性方程的基本解的计算方法 MATLAB源程序
发表评论
大 名:
(不网 址:
(您标 题:
内 容:
请根据下图中的字符输入验证码：
点这里显示验证码。
博客链
和讯博友免费美国之旅全记录高手教你暴跌中何时抢反弹
拉抬石化双雄用意何在
中国为何在G20之前提高准备金率（转载.
11.9收盘：休整的意义何在
短线持续突破的依托何在？
高校逼学生迁农村户口目的何在？
[微博]IPAD该不该被征高额关税 [微博]小贩对菜价上涨的超强分南方金砖四国现正发售财富人生 平安相伴
Re: 如何在MATLAB里绘制同心圆 源程序
字数上限为2000字
提 交
(您的评论将有可能审核后才能发表)
和讯个人门户 v1.0 | 和讯家园 | 意见反馈
%MATLAB程序：已知三个位置设计平面四杆机构求解程序（位移矩阵法） clear;clc; %凡是变量名前带v的为数值变量，不带的是符号变量 vxp1=0; vyp1=0; vsita1=0*pi/180; vxp2=-2; vyp2=6; vsita2=40*pi/180; vxp3=-10; vyp3=8; vsita3=90*pi/180; %精确位置P1，P2，P3及各角度 vsita12=vsita2-vsita1; vsita13=vsita3-vsita1; vxa=-10; vya=-2; vxd=-5; vyd=-2; %选定A,D点 %所有数值均在此确定，更改此处即可解出不同数值的四杆机构位移矩阵方程 syms xp1 yp1 xp2 yp2 xp3 yp3 sita12 sita13; syms xa ya xb1 yb1 xb2 yb2 xb3 yb3; f1='(xb2-xa)^2+(yb2-ya)^2=(xb1-xa)^2+(yb1-ya)^2'; f2='(xb3-xa)^2+(yb3-ya)^2=(xb1-xa)^2+(yb1-ya)^2'; %前两个机构方程 f3='xb2=cos(sita12)*xb1-sin(sita12)*yb1+xp2-xp1*cos(sita12)+yp1*sin(sita12)'; f4='yb2=sin(sita12)*xb1+cos(sita12)*yb1+yp2-xp1*sin(sita12)-yp1*cos(sita12)'; %由第一个位移矩阵方程得出 f5='xb3=cos(sita13)*xb1-sin(sita13)*yb1+xp3-xp1*cos(sita13)+yp1*sin(sita13)'; f6='yb3=sin(sita13)*xb1+cos(sita13)*yb1+yp3-xp1*sin(sita13)-yp1*cos(sita13)'; %由第二个位移矩阵方程得出 f1=subs(f1,{xa,ya},{vxa,vya}); f2=subs(f2,{xa,ya},{vxa,vya}); f3=subs(f3,{xp1,xp2,yp1,sita12},{vxp1,vxp2,vyp1,vsita12}); f4=subs(f4,{xp1,yp1,yp2,sita12},{vxp1,vyp1,vyp2,vsita12}); f5=subs(f5,{xp1,xp3,yp1,sita13},{vxp1,vxp3,vyp1,vsita13}); f6=subs(f6,{xp1,yp1,yp3,sita13},{vxp1,vyp1,vyp3,vsita13}); %代入具体数值 [xb1,xb2,xb3,yb1,yb2,yb3]=solve(f1,f2,f3,f4,f5,f6); %解方程 vxb1=vpa(xb1); vyb1=vpa(yb1); vxb2=vpa(xb2); vyb2=vpa(yb2); vxb3=vpa(xb3); vyb3=vpa(yb3); (vxb1-vxa)^2+(vyb1-vya)^2; (vxb2-vxa)^2+(vyb2-vya)^2; (vxb3-vxa)^2+(vyb3-vya)^2; %去掉这三行分号可验证B点三个位置是否距离A点相等 syms xd yd xc1 yc1 xc2 yc2 xc3 yc3;
Matlab 特殊图形和高维可视化 2009-10-20 01:06 7.4 特殊图形和高维可视化 7.4.1 特殊图形指令例示 7.4.1.1 面域图area 【* 例7.4.1 .1-1 】面域图指令area 。该指令的特点是：在图上绘制多条曲线时，每条曲线（除第一条外）都是把“前”条曲线作基线，再取值绘制而成。因此，该指令所画的图形，能醒目地反映各因素对最终结果的贡献份额。注意：（ 1 ）area 的第一输入宗量是单调变化的自变量。第二输入宗量是“各因素”的函数值矩阵，且每个“因素”的数据取列向量形式排放。第三输入宗量是绘图的基准线值，只能取标量。当基准值为0 （即以x 轴为基准线）时，第三输入宗量可以缺省。（2 ）本例第<4> 条指令书写格式x' , Y' ，强调沿列方向画各条曲线的事实。 clf;x=-2:2 % 注意：自变量要单调变化 Y=[3,5,2,4,1;3,4,5,2,1;5,4,3,2,5] % 各因素的相对贡献份额 Cum_Sum=cumsum(Y) % 各曲线在图上的绝对坐标 area(x',Y',0) %<4> legend(' 因素A',' 因素B',' 因素C'),grid on,colormap(spring) x = -2 -1 0 1 2 Y = 3 5 2 4 1 3 4 5 2 1 5 4 3 2 5 Cum_Sum = 3 5 2 4 1 6 9 7 6 2 11 13 10 8 7
图 7.4.1 .1-1 面域图表现各分量的贡献 7.4.1.2 各种直方图bar, barh, bar3, bar3h 【 * 例 7.4.1 .2-1 】二维直方图有两种图型：垂直直方图和水平直方图。而每种图型又有两种表现模式：累计式：分组式。本例选其两种加以表现。 x=-2:2; % 注意：自变量要单调变化 Y=[3,5,2,4,1;3,4,5,2,1;5,4,3,2,5]; % 各因素的相对贡献份额 subplot(1,2,1),bar(x',Y','stacked') % “累计式”直方图 xlabel('x'),ylabel('\Sigma y'),colormap(cool)% 控制直方图的用色legend(' 因素 A',' 因素 B',' 因素 C') subplot(1,2,2),barh(x',Y','grouped') % “分组式”水平直方图 xlabel('y'),ylabel('x') 图 7.4.1 .2-1 二维直方图 clf;x=-2:2; % 注意：自变量要单调变化 Y=[3,5,2,4,1;3,4,5,2,1;5,4,3,2,5]; % 各因素的相对贡献份额 subplot(1,2,1),bar3(x',Y',1) % “队列式”直方图 xlabel(' 因素 ABC'),ylabel('x'),zlabel('y') colormap(summer) % 控制直方图的用色 subplot(1,2,2),bar3h(x',Y','grouped') % “分组式”水平直方图 ylabel('y'),zlabel('x')
基于maｔlａｂ的连杆机构设计
————————————————————————————————作者: ————————————————————————————————日期：
目录 １平面连杆机构的运动分析 (1) 1．2 机构的工作原理 (1) 1.３机构的数学模型的建立 (1) 1.３.1建立机构的闭环矢量位置方程．.....．.......．...．....．...．..．．．..．....．...．．．...．１ 1.3.２求解方法..．..．..．.．..................．．...．........．....．...．．．..．．.．..．..．．．2 ２基于MATLAB程序设计 (4) 2.1 程序流程图 (4) 2．2 M文件编写 (6) ２.３程序运行结果输出 (7) 3 基于MATLAB图形界面设计 (11) 3.1界面设计……………………………………………………………………………………………１1 3.2代码设计……………………………………………………………………………………………１2
4 小结 (17) 参考文献 (18) 1平面连杆机构的运动分析 １.1 机构运动分析的任务、目的和方法 曲柄摇杆机构是平面连杆机构中最基本的由转动副组成的四杆机构，它可以用来实现转动和摆动之间运动形式的转换或传递动力。 对四杆机构进行运动分析的意义是:在机构尺寸参数已知的情况下,假定主动件（曲柄)做匀速转动，撇开力的作用,仅从运动几何关系上分析从动件（连杆、摇杆)的角位移、角速度、角加速度等运动参数的变化情况。还可以根据机构闭环矢量方程计算从动件的位移偏差。上述这些内容,无论是设计新的机械，还是为了了解现有机械的运动性能,都是十分必要的,而且它还是研究机械运动性能和动力性能提供必要的依据。 机构运动分析的方法很多,主要有图解法和解析法。当需要简捷直观地了解机构的某个或某几个位置的运动特性时，采用图解法比较方便，而且精度也能满足实际问题的要求。而当需要精确地知道或要了解机构在整个运动循环过程中的运动特性时，采用解析法并借助计算机，不仅可获得很高的计算精度及一系列位置的分析结果，并能绘制机构相应的运动线图，同时还可以把机构分析和机构综合问题联系起来,以便于机构的优化设计。 1.2 机构的工作原理 在平面四杆机构中,其具有曲柄的条件为: a.各杆的长度应满足杆长条件,即: 最短杆长度+最长杆长度≤其余两杆长度之和。 b.组成该周转副的两杆中必有一杆为最短杆，且其最短杆为连架杆或机架（当最短杆为连架杆时,四杆机构为曲柄摇杆机构;当最短杆为机架时，则为双曲柄机构）。 在如下图１所示的曲柄摇杆机构中，构件ＡB为曲柄,则B点应能通过曲柄与连杆两次共线的位置。
如何用M a t l a b绘制曲 线图 Document serial number【UU89WT-UU98YT-UU8CB-UUUT-UUT108】
各位同学：在写论文和报告时，为了很好地表达你研究和开发的结果，不仅要用文字详细地描述你方法、步骤和结果，还必须配以各种图来说明问题。下面是我们实验室张媛媛老师申请博士学位论文中的部分曲线图、硬件框图、软件流程图和实验装置原理框图。她将在部分曲线图下面给出绘制图形的Matlab程序和相关步骤，供大家学习和参考。 例一： 图2-3-6 动态线性环节的输入输出信号图2-3-7 模型输出和消噪后实验时数据比较1,输入信号u(k)；2,输出信号y(k) 1,实验数据；2,模型输出 绘图程序如下： figure(1) plot(t,y,'k',t,x,'k','LineWidth', xlabel('Time(s)','fontname','宋体','Fontsize',9);％绘制横坐标 ylabel('Voltage(v)','fontname','宋体','Fontsize',9); ％绘制纵坐标 %xlabel('时间(s)','fontname','宋体','Fontsize',9); %ylabel('电压(v)','fontname','宋体','Fontsize',9); ％设置合适的图框大小.可将下面四句变为子程序，以便调用。 set(gcf,'color',[1,1,1]); set(gca,'xcolor',[0,0,0],'ycolor',[0,0,0]); set(gcf,'units','centimeters','position',[5,10,,]); set(gca,'box','on','fontname','宋体','Fontsize',9);
Matlab图形绘制经典案例 1、 三维曲线 >> t=0:pi/50:10*pi; >> plot3(sin(2*t),cos(2*t),t) >> axis square >> grid on
2、一窗口多图形>> t=-2*pi:0.01:2*pi; >> subplot(3,2,1)
>> plot(t,sin(t)) >> subplot(3,2,2) >> plot(t,cos(t)) >> subplot(3,2,3) >> plot(t,tan(t)) >> axis([-pi pi -100 100]) >> subplot(3,2,4) >> plot(t,cot(t)) >> axis([-pi pi -100 100]) >> subplot(3,2,5) >> plot(t,atan(t)) >> subplot(3,2,6) >> plot(t,acot(t))
3、图形样式、标注、题字(也可以利用菜单直接Insert) >> x=0:pi/20:2*pi;
>> plot(x,sin(x),'b-.') >> hold on >> plot(x,cos(x),'r--') >> hold on >> plot(x,sin(x)-1,'g:') >> hold on >> plot(x,cos(x)-1) >> xlabel('x'); >> xlabel('x轴'); >> ylabel('y轴'); >> title('图形样式、标注等'); >> text(pi,sin(pi),'x=\pi'); >> legend('sin(x)','cos(x)','sin(x)-1','cos(x)-1'); >> [x1,y1]=ginput(1) %利用鼠标定位查找线上某点的值x1 = 2.0893 y1 = -0.5000 >> gtext('x=2.5') %鼠标定位放置所需的值在线上
MATLAB受到了广大理工科学生和学者青睐，除了Matlab强大的矩阵计算功能和功能齐全的toolbox以外，一个重要原因是因为它提供了方便的绘图功能。下面我们将详细介绍2维图形对象的生成函数及图形控制函数的使用方法以及一些图形的修饰与标注函数及操作和控制MATLAB各种图形对象的方法. 一、图形窗口与坐标系； A.图形窗口 1.MATLAB在图形窗口中绘制或输出图形,因此图形窗口就像一张绘图纸. 2.在MATLAB下,每一个图形窗口有唯一的一个序号h,称为该图形窗口的句 柄.MATLAB通过管理图形窗口的句柄来管理图形窗口; 3.当前窗口句柄可以由MATLAB函数gcf获得; 4.在任何时刻,只有唯一的一个窗口是当前的图形窗口(活跃窗口); figure(h)----将句柄为h的窗口设置为当前窗口; 5.打开图形窗口的方法有三种: 1)调用绘图函数时自动打开; 2)用File---New---Figure新建; 3)figure命令打开,close命令关闭. 在运行绘图程序前若已打开图形窗口,则绘图函数不再打开,而直接利用已打开的图形窗口;若运行程序前已存在多个图形窗口,并且没有指定哪个窗口为当前窗口时,则以最后使用过的窗口为当前窗口输出图形. 6.窗口中的图形打印:用图形窗口的File菜单中的Print项. 7.可以在图形窗口中设置图形对象的参数.具体方法是在图形窗口的Edit菜单中选择Properties项,打开图形对象的参数设置窗口,可以设置对象的属性. B.坐标系； 1.一个图形必须有其定位系统,即坐标系; 2.在一个图形窗口中可以有多个坐标系,但只有一个当前的坐标系; 3.每个坐标系都有唯一的标识符,即句柄值; 4.当前坐标系句柄可以由MATLAB函数gca获得; 5.使某个句柄标识的坐标系成为当前坐标系,可用如下函数:axes(h) h为指定坐标系句柄值.
首先创建函数FoutBarPosition，函数fsolve通过他确定。 function t=fourbarposition(th,th2,L2,L3,L4,L1) t=[L2*cos(th2)+L3*cos(th(1))-L4*cos(th(2))-L1;… L2*sin(th2)+L3*sin(th(1))-L4*sin(th(2))]; 主程序如下： disp ' * * * * * * 平面四杆机构的运动分析* * * * * *' L1=304.8;L2=101.6;L3=254.0;L4=177.8; %给定已知量，各杆长L1,L2,L3,L4 th2=[0:1/6:2]*pi; %曲柄输入角度从0至360度，步长为pi/6 th34=zeros(length(th2),2); %建立一个N行2列的零矩阵，第一列存放options=optimset('display','off'); %θ_3,第二列存放θ_3 for m=1:length(th2) %建立for循环，求解θ_3，θ_4 th34(m,:)=fsolve('fourbarposition',[1 1],…%调用fsove函数求解关于θ_3，θ_4 options,th2(m),L2,L3,L4,L1); %的非线性超越方程，结果保存在th34中 end y=L2*sin(th2)+L3*sin(th34(:,1)'); %连杆3的D端点Y坐标值 x=L2*cos(th2)+L3*cos(th34(:,1)'); %连杆3的D端点X坐标值 xx=[L2*cos(th2)]; %连杆3的C端点X坐标值 yy=[L2*sin(th2)]; %连杆3的C端点Y坐标值 figure(1) plot([x;xx],[y;yy],'k',[0 L1],[0 0],…%绘制连杆3的几个位置点 'k--^',x,y,'ko',xx,yy,'ks') title('连杆3的几个位置点') xlabel('水平方向') ylabel('垂直方向') axis equal %XY坐标均衡 th2=[0:2/72:2]*pi; %重新细分曲柄输入角度θ_2，步长为5度 th34=zeros(length(th2),2); options=optimset('display','off'); for m=1:length(th2)
用MATLAB画曲线族 （y-c)^2-2/3*(x-c)^3=0的包络线 1 求包络线的方程 syms x y c; f = (y-c)^2-2/3*(x-c)^3 dfc = diff(f, c) S = solve(f,dfc) S1x = S.x S1y = S.y 计算结果： 该曲线族有两条包络线： ① x1 = c1 ; y1 = c1 ; ② x1 = c1 + 2/3; y1 = c1 + 4/9; 2 画线 close all clear,clc warning('off') figure % 曲线族 hold on for c = -10:0.5:10 x = -10:0.1:10; y = (2/3)^0.5.*(x-c).^1.5 + c; plot(x,y) end % 包络线 c1 = -10:0.1:10; x1 = c1 ; y1 = c1 ; plot(x1,y1,'r','LineWidth',2)
figure % 曲线族 hold on for c = -10:0.5:10 x = -10:0.1:10; y = -(2/3)^0.5.*(x-c).^1.5 + c; plot(x,y) end % 包络线 c1 = -10-2/3:0.1:10-2/3; x1 = c1 + 2/3; y1 = c1 + 4/9; plot(x1,y1,'r','LineWidth',2)
............................ 包络线 跳转到：导航, 搜索
在几何学，某个曲线族的包络线（Envelope），是跟该曲线族的每条线都有至少一点相切的一条曲线。（曲线族即一些曲线的无穷集，它们有一些特定的关系。） 设一个曲线族的每条曲线C s可表示为 ，其中s是曲线族的参数，t是特定曲线的参数。若包络线存在，它是由 得出，其中h(s)以以下的方程求得：
实验二matlab图形绘制 一、实验目的 1、学习MATLAB图形绘制的基本方法； 2、熟悉和了解MATLAB图形绘制程序编辑的基本指令； 3、熟悉掌握利用MATLAB图形编辑窗口编辑和修改图形界面，并添加图形的各种标注； 二、实验原理 1．二维数据曲线图 （1）绘制单根二维曲线plot(x,y); （2）绘制多根二维曲线plot(x,y) 当x是向量，y是有一维与x同维的矩阵时，则绘制多根不同颜色的曲线。当x，y是同维矩阵时，则以x，y对应列 元素为横、纵坐标分别绘制曲线，曲线条数等于矩阵的列数。 （3）含有多个输入参数的plot函数plot(x1,y1,x2,y2,…,xn,yn) （4）具有两个纵坐标标度的图形plotyy(x1,y1,x2,y2) 2．图形标注与坐标控制 1）title （图形名称） 2）xlabel（x轴说明） 3）ylabel（y轴说明） 4）text（x，y图形说明） 5）legend（图例1，图例2，…） 6）axis （[xmin xmax ymin ymax zmin zmax]） 3．图形窗口的分割 subplot（m,n,p） 4．三维曲线 plot3（x1,y1,z1,选项1，x2,y2，选项2，…,xn,yn,zn,选项n）
5．三维曲面 mesh(x,y,z,c) 与surf(x,y,z,c)。一般情况下，x ，y ，z 是维数相同的矩阵。X ，y 是网格坐标矩阵，z 是网格点上的高度矩阵，c 用于指定在不同高度下的颜色范围。 三、实验内容及步骤 1．绘制下列曲线： (1) 2 1100 x y += x=0:0.02:10; y=100./(1+x.^2); plot(x,y) title('my first plot'); xlabel('x'); ylabel('y'); grid on 截图：
Matlab中Bode图的绘制技巧学术收藏2010-06-04 21:21:48 阅读54 评论0 字号：大中小订阅我们经常会遇到使用Matlab画伯德图的情况，可能我们我们都知道bode这个函数是用来画bode图的，这个函数是Matlab内部提供的一个函数，我们可以很方便的用它来画伯德图，但是对于初学者来说，可能用起来就没有那么方便了。譬如我们要画出下面这个传递函数的伯德图： 1.576e010 s^2 H(s= ------------------------------------------------------------------------------------------ s^4 + 1.775e005 s^3 + 1.579e010 s^2 + 2.804e012 s + 2.494e014 (这是一个用butter函数产生的2阶的，频率范围为[20 20K]HZ的带通滤波器。我们可以用下面的语句：num=[1.576e010 0 0]; den=[1 1.775e005 1.579e010 2.804e012 2.494e014]; H=tf(num,den; bode(H 这样，我们就可以得到以下的伯德图： 可能我们会对这个图很不满意，第一，它的横坐标是rad/s，而我们一般希望横坐标是HZ；第二，横坐标的范围让我们看起来很不爽；第三，网格没有打开（这点当然我们可以通过在后面加上grid on解决）。下面，我们来看看如何定制我们自己的伯德图风格：在命令窗口中输入：bodeoptions 我们可以看到以下
内容：ans = Title: [1x1 struct] XLabel: [1x1 struct] YLabel: [1x1 struct]TickLabel: [1x1 struct]Grid: 'off' XLim: {[1 10]}XLimMode: {'auto'}YLim: {[1 10]} YLimMode: {'auto'}IOGrouping: 'none'InputLabels: [1x1 struct]OutputLabels: [1x1 struct]InputVisible: {'on'} OutputVisible: {'on'}FreqUnits: 'rad/sec'FreqScale: 'log' MagUnits: 'dB' MagScale: 'linear'MagVisible: 'on' MagLowerLimMode: 'auto'MagLowerLim: 0PhaseUnits: 'deg'PhaseVisible: 'on'PhaseWrapping: 'off' PhaseMatching: 'off'PhaseMatchingFreq: 0 PhaseMatchingValue: 0我们可以通过修改上面的每一 项修改伯德图的风格，比如我们使用下面的语句画我 们的伯德图：P=bodeoptions;P.Grid='on'; P.XLim={[10 40000]};P.XLimMode={'manual'};P.FreqUnits='HZ'; num=[1.576e010 0 0];den=[1 1.775e005 1.579e010 2.804e012 2.494e014];H=tf(num,den; bode(H,P 这时，我们将会看到以下的伯德图： 上面这张图相对就比较好了，它的横坐标单位 是HZ，范围是[10 40K]HZ，而且打开了网格，便于我 们观察-3DB处的频率值。当然，你也可以改变bodeoptions中的其它参数，做出符合你的风格的伯
目录 1平面连杆机构的运动分析 (1) 1.2 机构的工作原理 (1) 1.3 机构的数学模型的建立 (1) 1.3.1建立机构的闭环矢量位置方程 (1) 1.3.2求解方法................................................................... ..2 2 基于MATLAB程序设计 (4) 2.1 程序流程图 (4) 2.2 M文件编写 (6) 2.3 程序运行结果输出 (7) 3 基于MATLAB图形界面设计 (11) 3.1界面设计 (11) 3.2代码设计 (12)
4 小结 (17) 参考文献 (18) 1平面连杆机构的运动分析 1.1 机构运动分析的任务、目的和方法 曲柄摇杆机构是平面连杆机构中最基本的由转动副组成的四杆机构，它可以用来实现转动和摆动之间运动形式的转换或传递动力。 对四杆机构进行运动分析的意义是：在机构尺寸参数已知的情况下，假定主动件（曲柄）做匀速转动，撇开力的作用，仅从运动几何关系上分析从动件（连杆、摇杆）的角位移、角速度、角加速度等运动参数的变化情况。还可以根据机构闭环矢量方程计算从动件的位移偏差。上述这些内容，无论是设计新的机械，还是为了了解现有机械的运动性能，都是十分必要的，而且它还是研究机械运动性能和动力性能提供必要的依据。 机构运动分析的方法很多，主要有图解法和解析法。当需要简捷直观地了解机构的某个或某几个位置的运动特性时，采用图解法比较方便，而且精度也能满足实际问题的要求。而当需要精确地知道或要了解机构在整个运动循环过程中的运动特性时，采用解析法并借助计算机，不仅可获得很高的计算精度及一系列位置的分析结果，并能绘制机构相应的运动线图，同时还可以把机构分析和机构综合问题联系起来，以便于机构的优化设计。 1.2 机构的工作原理 在平面四杆机构中，其具有曲柄的条件为： a.各杆的长度应满足杆长条件，即： 最短杆长度+最长杆长度≤其余两杆长度之和。 b.组成该周转副的两杆中必有一杆为最短杆，且其最短杆为连架杆或机架（当最短杆为连架杆时，四杆机构为曲柄摇杆机构；当最短杆为机架时，则为双曲柄机构）。 在如下图1所示的曲柄摇杆机构中，构件AB为曲柄，则B点应能通过曲柄与连杆两次共线的位置。
MATLAB画图入门篇--各种基本图形绘制的函数与实例【来自网络】 一．二维图形(Two dimensional plotting) 1.基本绘图函数(Basic plotting function)：Plot,semilogx,semilogy,loglog,polar,plotyy (1).单矢量绘图(single vector plotting)：plot(y),矢量y的元素与y元素下标之间在线性坐标下的关系曲线。 例1：单矢量绘图 y=[00.62.358.311.71517.719.420];plot(y) 可以在图形中加标注和网格， 例2：给例1的图形加网格和标注。 y=[00.62.358.311.71517.719.420];plot(y) title('简单绘图举例');xlabel('单元下标');ylabel('给定的矢量');grid (2).双矢量绘图(Double vector plotting)：如x和y是同样长度的矢量,plot(x,y)命令将绘制y元素对应于x元素的xy曲线图。 例：双矢量绘图。 x=0:0.05:4*pi;y=sin(x);plot(x,y) (3).对数坐标绘图(ploting in logarithm coordinate)：x轴对数semilogx,y轴对数semilogy,双对数loglog, 例：绘制数组y的线性坐标图和三种对数坐标图。 y=[00.62.358.311.71517.719.420]; subplot(2,2,1);plot(y);subplot(2,2,2);semilogx(y) subplot(2,2,3);semilogy(y);subplot(2,2,4);loglog(y) （4）极坐标绘图(Plotting in polar coordinate)： polar(theta,rho)theta—角度，rho—半径 例：建立简单的极坐标图形。 t=0:.01:2*pi;polar(t,sin(2*t).*cos(2*t)) 2.多重曲线绘图(Multiple curve plotting) （1）一组变量绘图(A group variable plotting) plot(x,y) (a)x为矢量，y为矩阵时plot(x,y)用不同的颜色绘制y矩阵中各行或列对应于x的曲线。 例1： x=0:pi/50:2*pi;y(1,:)=sin(x);y(2,:)=0.6*sin(x);y(3,:)=0.3*sin(x);plot(x,y) (b)x为矩阵，y为矢量时绘图规则与（a）的类似，只是将x中的每一行或列对应于y进行绘图。。 例2： x(1,:)=0:pi/50:2*pi;x(2,:)=pi/4:pi/50:2*pi+pi/4;x(3,:)=pi/2:pi/50:2*pi+pi/2; y=sin(x(1,:));plot(x,y) (c)x和y是同样大小的矩阵时,plot(x,y)绘制y矩阵中各列对应于x各列的图形。 例3： x(:,1)=[0:pi/50:2*pi]';x(:,2)=[pi/4:pi/50:2*pi+pi/4]';x(:,3)=[pi/2:pi/50:2*pi+pi/2]'; y(:,1)=sin(x(:,1));y(:,2)=0.6*sin(x(:,1));y(:,3)=0.3*sin(x(:,1)); plot(x,y) 这里x和y的尺寸都是101×3，所以画出每条都是101点组成的三条曲线。如行列转置后就会画出101条曲线，每条线
本科生毕业设计 基于MATLAB的双摇杆机构运动分析与仿真 Based on the MATLAB double rocker organization movement analysis and simulation
基于MATLAB/SIMULINK的双摇杆机构运动学分析与仿 真 邹凯旋 云南农业大学工程技术学院，昆明黑龙潭650201 摘要 平面连杆机构的应用十分广泛，它的分析及设计一直是机构学研究的一个重要课题。MATLAB的Simulink是一个对动态系统建模和仿真分析的软件包，为信号与系统仿真实验提供了很好的平台。借助其强大的模拟仿真分析功能可以方便的实现机构性能分析和动态仿真，降低分析的难度，有效提高设计工作效率、产品开发质量、降低开发成本。本设计课题以MATLAB的simulink\simMechanics 动态模拟仿真工具为平台，对双摇杆机构进行运动分析。结果表明该仿真方法能方便、准确的得到机构的运动、动力数据，能为机构的选择、优化设计提供参考依据。应用此工具可很好地对机械系统的各种运动进行分析，构造出平面连杆机构的数学模型。通过对此数学模型分析，分离出可独立求解的机构模型，并用相应的机构分析方法对它进行求解，建立了平面连杆机构运动学分析专家系统。系统可完成部分平面连杆机构的运动学分析及动画仿真，从而为机械系统的建模仿真提供一个强大而方便的工具。 关键词：连杆机构；动态仿真；SimMechanics；数学模型
Based on the MATLAB double rocker organization movement analysis and simulation Zou kaixuan Faculty of Engineering and Technology Yunan Agricultural University，Heilongtan Kunming 650201 ABSTRACT Planar linkage mechanism used widely, its analysis and design of the study of institutions has been an important subject. MATLAB Simulink is a dynamic system modeling and simulation software package, for signal and system simulation results provide a good platform. With its powerful simulation analysis function is realized the performance analysis and the dynamic simulation institutions, reduce the difficulties of analysis, effectively improve the design work efficiency and product development quality, reduce development costs. This design task to MATLAB simulink \ simMechanics dynamic simulation tools as the platform, on the double rocker organization motion analysis. The results show that the simulation method can conveniently, accurately to get the kinematic and dynamic data organization, for the choice of institutions, optimum design to provide the reference. This tool can application is mechanical system analysis of all kinds of sports, constructed the mathematical model of the planar linkage mechanism. Through mathematical model to analysis, separating out can be independent of solving mechanism model, and the corresponding institutions analysis method to solve it, a planar linkage mechanism kinematic analysis of the expert system. System can finish part of planar linkage mechanism kinematic analysis and animated simulation, thus for mechanical system modeling simulation provide a strong and convenient tool. Key words: linkage；Dynamic Simulation；SimMechanics；mathematical model
绘图功能是Matlab的一个强大的功能。 subplot是MATLAB中常用的的函数。在绘图过程中经常要在一个页面中绘制几张图 它的使用格式：subplot（m,n,p）或者subplot（m n p）。 函数subplot是将多个图画到一个平面上的工具。括号中的m表示是图排成m行，n 表示图排成n列，也就是整个figure中有n个图是排成一列的，一共m行，如果m=3就是表示3行图。p表示图所在的位置，p=3表示从左到右从上到下的第3个位置。 以下是对它的一些应用，还用到了其它的一些函数 程序的代码如下 x=0:0.01:10; y1=sin(x); subplot(3,3,1); plot(x,y1); xlabel('x'); ylabel('y1'); title('y1=sin(x)');axis([0 pi*2 -1 1]); y2=cos(x+2); subplot(3,3,2); plot(x,y2); xlabel('x'); ylabel('y2');title('y2=cos(x+2)');axis([0 pi*2 -1 1]); y3=sin(x)+y2;subplot(3,3,3); plot(x,y3); xlabel('x'); ylabel('y3');title('y3=sin(x)+y2'); y4=sin(x).^3+cos(x);subplot(3,3,4);axis([0 pi*2 -1 2]); plot(x,y4); xlabel('x'); ylabel('y4');title('y4=sin(x).^3+cos(x)');axis([0 pi*2 -1 2]); y5=9*x.^5+3*x.^4+x.^3+2*x.^2;
Matlab绘制频散曲线 程序代码 -CAL-FENGHAI-(2020YEAR-YI function disper %绘制平板频散曲线 %tic
clc;clear; cl=5790;%材料纵波波速(钢板) cs=3200;%材料横波波速(钢板) dfd=*le3; fdO=:dfd/le3:2O)*le3;%频厚积(MHz*mm) d_Q235二6; cps_mi n二2700; cpa_min=100; cp_max=10000; mode=3;%绘制的模式数 precision=le-8; cpa=zeros(length(fdO),mode); cps=zeros(le ng th(fdO),mode); for i=l:length(fdO) fd=fdO(i); [cpl2 n]=ss(cps_min/cp_max/fd/cl,cs,mode); for j=l:n cpl=cpl2(j,l); cp2=cpl2(j,2); cps(i,j)=serfe n(cpl,cp2,fctcl￡S'precisi on); end [cpl2 n]=aa(cpa_min,cp_max/fd/cl/cs,mode); for j=l:n cpl=cpl2(j,l); cp2=cpl2(j,2); cpa(ij)=aerfe n(cpbcp2,fd￡l‘cs,precisi on); end end h=zeros(mode,2); %相速度 figure(l) for j=l:2 ifj==l cp=cps; color=,b,; else cp=cpa; color二T; end for i=l:mode cpp=cp(:,i); in d=fi nd(cpp==0); if ^isempty(ind) h(i/j)=plot((fdO(ind(end)+l:end))/d_Q235/cpp(ind(end)+l:end),color); else h(i/j)=plot(fdO/d_Q235,cpp/color); end hold on end ifj==2 xlabel('f/(KHz)') ylabel('C_{p}/(km-sA{-l})')
基于matlab 的平面四连杆机构设计以及该机构的运动仿真分析 摘要 四连杆机构因其结构方便灵活,能够传递动力并实现多种运动形式而被广泛应用于各个领域,因此对其进行运动分析具有重要的意 义。传统的分析方法主要应用几何综合法和解析综合法,几何综合法简单直观,但是精确度较低;解析法精确度较高,但是计算工作量大。随着计算机辅助数值解法的发展,特别是MATLAB 软件的引入,解析法已经得到了广泛的应用。对于四连杆的运动分析，若应用MATLAB 则需要大量的编程，因此我们引入proe 软件，我们不仅可以在此软件中建立实物图，而且还可以对其进行运动仿真并对其运动分析。 在设计四连杆时，我们利用解析综合法建立数学模型，再根据数学模型在MATLAB 中编程可以求得其他杆件的长度。 针对范例中所求得的各连杆的长度，我们在proe 软件中画出其三维图（如图4）并在proe软件中进行仿真分析得出B,C的角加速度的变化，从而得到B,C两接触处所受到的力是成周期性变化的，可以看出B,C两点处极易疲劳断裂，针对B,C两点处的疲劳断裂，我们提出了在设计四连杆中的一些建议。 关键字：解析法MATLAB软件proe软件运动仿真
建立用解析法设计平面四杆机构模型 对于问题中所给出的连架杆AB的三个位置与连架杆CD的三个位置相对应，即三组对应位置为： f 1」2卜2,「3卜3，其中他们对应的值分别为：135 ,112 ,90 ,82 ,45 ,52，为了便于写代数式，可作出AB与CD对应的关系，其图如下： 图一2 AB与CD三个位置对应的关系 通过上图我们可以通过建立平面直角坐标系并利用解析法来求解，其直角坐标系图如下： 图一3平面机构直角坐标系 通过建立直角坐标系OXY,如上图所示，其中:0与°为AB杆与
首先创建函数FｏuｔBaｒPoｓｉｔioｎ，函数ｆsolvｅ通过他确定。 fｕnｃtion t＝fｏurｂarｐoｓitｉon(tｈ，th2,L2,L3,L4,L1) t=［L2*coｓ(th２）+L3*cos(th（1））-Ｌ4*cos(tｈ(2))-L1;… Ｌ2*ｓin(ｔh2)＋Ｌ3＊sin(ｔh（1))-Ｌ4*sin（th(2）)]; 主程序如下： ｄisp ＇* ＊* ＊＊* 平面四杆机构的运动分析*＊* **＊＇ L1=3０4.8；L２=1０1.６;L3＝254.0;Ｌ4=１77.8；%给定已知量,各杆长L1,L2,L３,L4 th2=［0:1/6：2］＊ｐi; %曲柄输入角度从0至360度,步长为pｉ/６ tｈ34=ｚeroｓ（ｌｅnｇth(tｈ2），2）; %建立一个Ｎ行2列的零矩阵，第一列存放opｔｉons=oｐtiｍsｅt（＇disｐlay＇,'ｏff'）; %θ_３,第二列存放θ＿３ for m=１:length(th2) %建立for循环,求解θ_3，θ_４ th3４（m,:)=fｓolｖe(＇ｆｏurｂarpｏｓｉtion＇,[1 １]，…%调用ｆｓｏｖｅ函数求解关于θ_3,θ＿4 opｔions,tｈ2(m),L２,L3，L4,Ｌ１）; %的非线性超越方程，结果保存在th３4中 end y=L2*sin(th2)+L３*sin(th34(:，１)')；%连杆3的D端点Y坐标值 x=Ｌ２＊ｃos(th2)+Ｌ3*cｏｓ(ｔｈ３4(：,1)')；%连杆3的D端点X坐标值 xｘ=［L２*ｃｏs(ｔh2)]; %连杆３的C端点X坐标值 yy=［L2*ｓin(th2)]; %连杆３的C端点Y坐标值 figure（1) pｌot(［x;xｘ],[y;yy］,'k',［0 L1]，[0 0]，…%绘制连杆3的几个位置点 'k-－＾',x,y，'ｋo',xｘ,yy,'ks'） titlｅ（'连杆3的几个位置点') ｘｌabｅｌ（'水平方向') ylabel('垂直方向'） ａｘｉs equal ％XY坐标均衡 th2＝[0:２/72:２]*pｉ; ％重新细分曲柄输入角度θ＿2,步长为5度 ｔh３4=zeros（ｌength(ｔh2）,2)； ｏｐｔｉons=ｏｐtiｍｓeｔ('diｓｐlay＇，＇off'); ｆoｒm＝1:leｎgth(ｔｈ2）