資源描述：
《MATLAB結(jié)課作業(yè)》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫(kù)。
1、計(jì)算應(yīng)用軟件MATLAB程序設(shè)計(jì)及繪圖班級(jí)：學(xué)號(hào)：姓名：2016.12計(jì)算應(yīng)用軟件.工程力學(xué)專業(yè)題目1：分別利用For循環(huán)和While循環(huán)求下值：。For循環(huán)代碼：k=0;fori=0:1:63k=k+2^i;endkWhile循環(huán)代碼:k=0;i=0;while(i<64)k=k+2^i;i=i+1;endk結(jié)果：k=1.8447e+019題目2：已知方程組：，求方程組的解。代碼：A=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10];B=[4;-3;9;-8];X=AB11/11計(jì)算應(yīng)用軟件.工程力學(xué)專業(yè)結(jié)果：X=-1.4841-0.68
2、160.5337-1.2429題目3：某試驗(yàn)試件在受火后，其彈性模量發(fā)生變化，數(shù)據(jù)如下：受火溫度T=[20,100,200,300,400,500,600,700,800,900];受火后彈性模E=[21.95,1.85,1.76,1.70,1.55,1.50,1.42,1.35,1.3,1.32]*1011;用合適的多項(xiàng)式擬合彈性模量隨受火溫度變化的規(guī)律。代碼：T=[20,100,200,300,400,500,600,700,800,900];E=[21.95,1.85,1.76,1.70,1.55,1.50,1.42,1.35,1.3,1.28]*10^11;plot(
3、T,E,'LineWidth',4)結(jié)果：11/11計(jì)算應(yīng)用軟件.工程力學(xué)專業(yè)題目4：（1）編寫一個(gè)函數(shù)式文件，使其產(chǎn)生如下函數(shù)段：（2）編寫命令式文件，繪制x在[0,10]范圍內(nèi)f(x).f(x+2)的曲線。（1）代碼：functionf=Untitled(x)ifx<=2f=2*x;elseifx>6f=x-4;elsef=5-0.5*x;end（2）代碼：i=0;forx=0:0.1:10i=i+1;y1(1,i)=Untitled(x);endi=0forx=0:0.1:10i=i+1;y2(1,i)=Untitled(x+2);endy3=y1.*y2;x=0:0.
4、1:10;plot(x,y3,'LineWidth',5)11/11計(jì)算應(yīng)用軟件.工程力學(xué)專業(yè)結(jié)果：題目5：設(shè)f(x)=x5-4x4+3x2-2x+6；（1）取x=[-2,8]之間函數(shù)的值（取100個(gè)點(diǎn)），畫出曲線，看它有幾個(gè)零點(diǎn)。（提示：用polyval函數(shù)）（2）用roots函數(shù)求此多項(xiàng)式的根。（1）代碼：p=[1,-4,3,-2,6];x=-2:0.1:8;y=polyval(p,x);plot(x,y,'LineWidth',5)11/11計(jì)算應(yīng)用軟件.工程力學(xué)專業(yè)結(jié)果：（2）代碼：r=roots(p)n=length(r);fori=1:nifr(i)>=-2&&r
5、(i)<=8if(isreal(r(i)))r(i)endendend結(jié)果：ans=3ans=1.695611/11計(jì)算應(yīng)用軟件.工程力學(xué)專業(yè)題目6：在x=[-10,10]，y=[-10,10]范圍內(nèi)畫出函數(shù)的三維圖形。代碼：x=-10:0.01:10;y=-10:0.01:10;[X,Y]=meshgrid(x,y);Z=sin(sqrt(X.^2+Y.^2))./sqrt(X.^2+Y.^2);mesh(X,Y,Z)結(jié)果：題目7：由100*rand(30)作為全班成績(jī)S(i)，用餅圖分別表示、、、、所占比例。代碼：S=100*rand(30);k1=0;k2=0;k3=0
6、;k4=0;k5=0;fori=1:900ifS(i)<60k1=k1+1;11/11計(jì)算應(yīng)用軟件.工程力學(xué)專業(yè)elseifS(i)>=60&S(i)<70k2=k2+1;elseifS(i)>=70&S(i)<80k3=k3+1;elseifS(i)>=80&S(i)<90k4=k4+1;elsek5=k5+1;endendx=[k1k2k3k4k5];y=[11111];pie(x,y,{'Si<60','60<=Si<70','70<=Si<80','80<=Si<90','90<=Si'});結(jié)果：11/11計(jì)算應(yīng)用軟件.工程力學(xué)專業(yè)題目8：列舉“MATLAB程序設(shè)計(jì)
7、課程”在本專業(yè)的一個(gè)相關(guān)應(yīng)用實(shí)例。求解任意長(zhǎng)度為L(zhǎng)的簡(jiǎn)支梁在任意大小的集中F的作用下的彎矩圖和剪力圖，最大剪力值:F1max和最大彎矩值Mmax：。代碼：clcclear;disp('輸入集中力到A點(diǎn)的距離a,梁的長(zhǎng)度l,集中力F')a=input('a（/m）值：');L=input('L（/m）值：');F=input('F（/KN）值：');M=[];n=0;forx=0:0.1:L;ifx>=0&&x<=an=n+1;F1(n)=(L-a)/L*F;M(n)=(L-a)/F*x;elsen=n+
