资源描述：
《MATLAB结课作业》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、计算应用软件MATLAB程序设计及绘图班级：学号：姓名：2016.12计算应用软件.工程力学专业题目1：分别利用For循环和While循环求下值：。For循环代码：k=0;fori=0:1:63k=k+2^i;endkWhile循环代码:k=0;i=0;while(i<64)k=k+2^i;i=i+1;endk结果：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计算应用软件.工程力学专业结果：X=-1.4841-0.68
2、160.5337-1.2429题目3：某试验试件在受火后，其弹性模量发生变化，数据如下：受火温度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;用合适的多项式拟合弹性模量随受火温度变化的规律。代码：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)结果：11/11计算应用软件.工程力学专业题目4：（1）编写一个函数式文件，使其产生如下函数段：（2）编写命令式文件，绘制x在[0,10]范围内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计算应用软件.工程力学专业结果：题目5：设f(x)=x5-4x4+3x2-2x+6；（1）取x=[-2,8]之间函数的值（取100个点），画出曲线，看它有几个零点。（提示：用polyval函数）（2）用roots函数求此多项式的根。（1）代码：p=[1,-4,3,-2,6];x=-2:0.1:8;y=polyval(p,x);plot(x,y,'LineWidth',5)11/11计算应用软件.工程力学专业结果：（2）代码：r=roots(p)n=length(r);fori=1:nifr(i)>=-2&&r
5、(i)<=8if(isreal(r(i)))r(i)endendend结果：ans=3ans=1.695611/11计算应用软件.工程力学专业题目6：在x=[-10,10]，y=[-10,10]范围内画出函数的三维图形。代码：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)结果：题目7：由100*rand(30)作为全班成绩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计算应用软件.工程力学专业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'});结果：11/11计算应用软件.工程力学专业题目8：列举“MATLAB程序设计
7、课程”在本专业的一个相关应用实例。求解任意长度为L的简支梁在任意大小的集中F的作用下的弯矩图和剪力图，最大剪力值:F1max和最大弯矩值Mmax：。代码：clcclear;disp('输入集中力到A点的距离a,梁的长度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+
