MATLAB計算機仿真上機作業(yè)

1、1、下列表達式表達了在主接觸應力F的作用下，兩個球體被擠壓到一起時在x,y,z方向上所產(chǎn)生的接觸應力值。2M x Q y Pmax(1Zarctan(a)(1 w)0.5(1$) 1aza（TPmaxz2 /a2a式中，3F （1 v）/E1（1 v；）/E21/d11/d2p max3F2n a23、如圖所示曲柄滑塊機構(gòu)，滑塊的位移為: s a cos)、b2 (a si n() e)2位移s是角度以角度表示)的函數(shù)，當4、由下式產(chǎn)生黃金分割數(shù)：vi,Ei和dj(j 1,2)分別為兩個球的泊松比，彈性模量和直徑。用MATLAB表達形式寫出上述公式并用下列給出的數(shù)值進行計算：V1 V2 0.3

2、, E1E2 3 107,d11.5,d22.7,F 100lb(磅，等于 373克),z 0.01in(英寸，等于2.54厘米)2 、 給 出向量 x 17, 3, 47,5,0,29, 37,51, 7,19。 編程使 x 按 照向量 y 0, 3, 7, 37, 47,51,29,19,17,5的形式重新排列。程序應使用于任意長度的向量。對所有向量，數(shù) 值0與負數(shù)放在一起。也就是說，如果0為向量的一個元素，則它將作為y向量的第一個元素。a 1,b1.5, e 0.3,0360 時，使用11.5、n 1.5、niFn()(),n 0,1,2,5 2 2使用fprintf和disp產(chǎn)生前16

3、個數(shù)值，并在 MATLAB命令窗口中顯示如下內(nèi)容:F0=0F1=1F2=1F3=2F15=6105、find函數(shù)的實現(xiàn)給定由正數(shù)和負數(shù)構(gòu)成的任意長度向量g,用程序?qū)崿F(xiàn)下式:In dx=fi nd(ga)這里a是由用戶指定的一個值。本例中a=4,向量g 4,4,7,10, 6,42,1,0。6、 用極坐標畫半徑為2 的圓答案如下： 第一題function thetax,thetay,thetaz= zcj01(v1,v2,E1,E2,d1,d2,F,z)% 函數(shù)功能：計算 xyz 方向應力% 輸入?yún)?shù)： v1,v2,E1,E2,d1,d2,F,z% 輸出參數(shù)： thetax,thetay,the

4、taz v1=0.3;v2=0.3;E1=3*10A7;E2=3*10A7;d1=1.5;d2=2.7;F=37300;z=0.254;a=(3*F/8)*(1-v1A2)/E1)+(1-v2A2)/E2)/(1/d1+1/d2)A(1/3); pmax=3*F/(2*pi*(aA2);thetax=-pmax*(1-(z/a)*atan(a/z)*(1-v1)-0.5*(1+zA2/(aA2)A(-1);thetay=-pmax*(1-(z/a)*atan(a/z)*(1-v1)-0.5*(1+zA2/(aA2)A(-1); thetaz=-pmax/(1+zA2/(aA2);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);% 降序排列