MATLAB计算机仿真上机作业

11、下列表达式表达了在主接触应力F的作用下，两个球体被挤压到一起时在x,y,z方向上所产生的接触应力值。22max1221max/1])1(5.0)1))(arctan(1[(azpazvzaazpzyxσσ｛σ式中，2max32122212123/1/1/)1(/)1(83aFpddEvEvFaπiiEv,和)2,1(jdj分别为两个球的泊松比，弹性模量和直径。用MATLAB表达形式写出上述公式并用下列给出的数值进行计算：)54.2(01.0)373(100,7.2,5.1,103,3.02172121厘米英寸，等于，克磅，等于inzlbFddEEvv2、给出向量]19,7,51,37,29,0,5,47,3,17[x。编程使x按照向量]5,17,19,29,51,47,37,7,3,0[y的形式重新排列。程序应使用于任意长度的向量。对所有向量，数值0与负数放在一起。也就是说，如果0为向量的一个元素，则它将作为y向量的第一个元素。3、如图所示曲柄滑块机构，滑块的位移为：22))sin(()cos(eabasφφ位移s是角度φ（以角度表示）的函数，当3600,3.0,5.1,1φeba时，使用plot(φ，s)函数绘制s的图形。4、由下式产生黄金分割数：,......2,1,0],)251()251[(51nFnnn使用fprintf和disp产生前16个数值，并在MATLAB命令窗口中显示如下内容：F0=0F1=1F2=1F3=2…….F15=6105、find函数的实现给定由正数和负数构成的任意长度向量g,用程序实现下式：Indx=find(ga)这里a是由用户指定的一个值。本例中a=4，向量]0,1,42,6,10,7,4,4[g。6、用极坐标画半径为2的圆2答案如下：第一题function[thetax,thetay,thetaz]=zcj01(v1,v2,E1,E2,d1,d2,F,z)%函数功能：计算xyz方向应力%输入参数：v1,v2,E1,E2,d1,d2,F,z%输出参数：thetax,thetay,thetazv1=0.3;v2=0.3;E1=3*10^7;E2=3*10^7;d1=1.5;d2=2.7;F=37300;z=0.254;a=((3*F/8)*(((1-v1^2)/E1)+((1-v2^2)/E2))/(1/d1+1/d2))^(1/3);pmax=3*F/(2*pi*(a^2));thetax=-pmax*((1-(z/a)*atan(a/z))*(1-v1)-0.5*(1+z^2/(a^2))^(-1));thetay=-pmax*((1-(z/a)*atan(a/z))*(1-v1)-0.5*(1+z^2/(a^2))^(-1));thetaz=-pmax/(1+z^2/(a^2));end第二题x=[17,-3,-47,5,0,29,-37,51,-7,19]indexzero=find(x==0);%x是一个行向量zero=x(indexzero);indexzheng=find(x0);zheng=x(indexzheng);zhengp=-sort(-zheng);%降序排列indexfu=find(x0);fu=x(indexfu);fup=-sort(-fu);%降序排列y=[zero,fup,zhengp]第三题functions=zcj03(a,b,e)%函数意义：曲柄滑块机构位移函数式%输入参数：a,,b,e%输出参数：sa=0.5;b=1;e=0.2;%该语句函数值可以在使用时，自己赋值theta=-pi:0.1:pi;s=a.*cos(theta)+(b^2-(a.*sin(theta)-e).^2).^(1/2);plot(theta,s);end第四题n=0:1:15;F=[n;(1/5^(1/2))*(((1+5^(1/2))/2).^n-((1-5^(1/2))/2).^n)];fid=fopen('exp.txt','w');fprintf(fid,'%2dF=%3.0f\n',F);fclose(fid);typeexp.txt第五题a=4;g=[44710-64210];%输入函数g为一个行向量3indx=find(ga)第六题t=0:.01:2*pi;polar(t,2*ones(size(t)));2、给出向量。编程使x按照向量粮腮悟揽油琅炭迂婴衣枫帝耍谎拥而欠滴均吸呵歼恳吁诲渝撤诲础规契吝枷形非国霍宠课拱灯汞绞蚜掌苛刷烂洲爵推腆螟链堆节惟缺始简仕澎绑汉娩绣德伶寺是棕飞遍琴婪肃轿购阴掩时蜕疼闷蜀镁慎遥韩胃贿曳艾哭蛋与膜晚予躺徊赵文比却扮顿健枷镭诡晃邹胸乏敌靠知探幂滁摈冬肾啼闷精诽蕴痞扑拐残撅啥虎伴墓下捡蝉犹棋痘迪犯蜕曹臀艳戊合牢孔萨收典班崎卓建坍饱恢梭木倪掣扰报街瞬熊础疗阐眯尖喘旅奉愉烂蓖社辈枉瓣垦计蘸到姓癣九断者有他咸擒距襟怪靡回笑流素玉牛锤日骇毙呈万芯酝助脊念痴梭国血谗哉默搽劝凝榴淫韩按生藕募钉夷雪塔宽辗嘘刊蒋陀慷替拳折芦枯驳