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

1、vi,Ei和dj(j 1,2)分別為兩個(gè)球的泊松比，彈性模量和直徑。用 MATLAB 表達(dá)形式寫(xiě)出上述公式并用下列給出的數(shù)值進(jìn)行計(jì)算：V1V20.3, E1E23 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的形式重新排列。程序應(yīng)使用于任意長(zhǎng)度的向量。對(duì)所有向量，數(shù) 值 0 與負(fù)數(shù)放在一起。也就是說(shuō)，如果0 為向量的一個(gè)元素，則它將作為y 向量的第一個(gè)元素。a 1

2、,b1.5, e 0.3,0360時(shí)，使用11.5、n1.5、niFn()(),n 0,1,2,5 2 2使用 fprintf 和 disp 產(chǎn)生前 16 個(gè)數(shù)值，并在 MATLAB 命令窗口中顯示如下內(nèi)容:F0=0F1=1F2=1F3=21、下列表達(dá)式表達(dá)了在主接觸應(yīng)力F 的作用下，兩個(gè)球體被擠壓到一起時(shí)在x,y,z 方向上所產(chǎn)生的接觸應(yīng)力值。2M xQyPmax(1Zarctan(a)(1 w)0.5(1$)1aza（TPmaxz2/a2a式中，3F （1 v）/E1（1V；）/E21/d11/d2pmax3F2na23、 如圖所示曲柄滑塊機(jī)構(gòu)， 滑塊的位移為:sa cos)、b2(a s

3、i n() e)2位移 s 是角度以角度表示)的函數(shù)，當(dāng)4、由下式產(chǎn)生黃金分割數(shù)：F15=6105、find 函數(shù)的實(shí)現(xiàn)給定由正數(shù)和負(fù)數(shù)構(gòu)成的任意長(zhǎng)度向量g,用程序?qū)崿F(xiàn)下式:In dx=fi nd(ga)這里 a 是由用戶(hù)指定的一個(gè)值。本例中a=4,向量g 4,4,7,10, 6,42,1,0。6、 用極坐標(biāo)畫(huà)半徑為 2 的圓答案如下： 第一題function thetax,thetay,thetaz= zcj01(v1,v2,E1,E2,d1,d2,F,z)% 函數(shù)功能：計(jì)算 xyz 方向應(yīng)力% 輸入?yún)?shù)： v1,v2,E1,E2,d1,d2,F,z% 輸出參數(shù)： thetax,thetay

4、,thetaz 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 是一個(gè)行向量 zero=x(indexzero);indexzheng=find(x0); zheng=x(indexzheng);zhengp=-sort(-zheng);% 降序排列