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

MATLAB計(jì)算機(jī)仿真上機(jī)作業(yè)MATLAB計(jì)算機(jī)仿真上機(jī)作業(yè)MATLAB計(jì)算機(jī)仿真上機(jī)作業(yè)資料僅供參考文件編號：2022年4月MATLAB計(jì)算機(jī)仿真上機(jī)作業(yè)版本號：A修改號：1頁次：1.0審核：批準(zhǔn):發(fā)布日期：1、下列表達(dá)式表達(dá)了在主接觸應(yīng)力F的作用下，兩個(gè)球體被擠壓到一起時(shí)在x,y,z方向上所產(chǎn)生的接觸應(yīng)力值。式中，和分別為兩個(gè)球的泊松比，彈性模量和直徑。用MATLAB表達(dá)形式寫出上述公式并用下列給出的數(shù)值進(jìn)行計(jì)算：2、給出向量。編程使x按照向量的形式重新排列。程序應(yīng)使用于任意長度的向量。對所有向量，數(shù)值0與負(fù)數(shù)放在一起。也就是說，如果0為向量的一個(gè)元素，則它將作為y向量的第一個(gè)元素。3、如圖所示曲柄滑塊機(jī)構(gòu)，滑塊的位移為：位移s是角度φ（以角度表示）的函數(shù)，當(dāng)時(shí)，使用plot(φ，s)函數(shù)繪制s的圖形。4、由下式產(chǎn)生黃金分割數(shù)：使用fprintf和disp產(chǎn)生前16個(gè)數(shù)值，并在MATLAB命令窗口中顯示如下內(nèi)容：F0=0F1=1F2=1F3=2…….F15=6105、find函數(shù)的實(shí)現(xiàn)給定由正數(shù)和負(fù)數(shù)構(gòu)成的任意長度向量g,用程序?qū)崿F(xiàn)下式：Indx=find(g>a)這里a是由用戶指定的一個(gè)值。本例中a=4，向量。6、用極坐標(biāo)畫半徑為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,thetazv1=;v2=;E1=3*10^7;E2=3*10^7;d1=;d2=;F=37300;z=;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)*(1+z^2/(a^2))^(-1));thetay=-pmax*((1-(z/a)*atan(a/z))*(1-v1)*(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是一個(gè)行向量zero=x(indexzero);indexzheng=find(x>0);zheng=x(indexzheng);zhengp=-sort(-zheng);%降序排列indexfu=find(x<0);fu=x(indexfu);fup=-sort(-fu);%降序排列y=[zero,fup,zhengp]第三題functions=zcj03(a,b,e)%函數(shù)意義：曲柄滑塊機(jī)構(gòu)位移函數(shù)式%輸入?yún)?shù)：a,,b,e%輸出參數(shù)：sa=;b=1;e=;%該語句函數(shù)值可以在使用時(shí)，自己賦值theta=-pi::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('','w');fprintf(fid,'%2dF=%\n',F);fclose(fid);type第五題a=4;g=[44710-64210