您好,欢迎访问三七文档
当前位置:首页 > 高等教育 > 理学 > MATLAB上机作业
MATLAB上机作业1对以下问题,编写M文件:(1)用起泡法对10个数由小到大排序。即将相邻两个数比较,将小的调到前头。functionf=qipaofa(x)forj=9:-1:1fori=1:jif(x(i)x(i+1))t=x(i);x(i)=x(i+1);x(i+1)=t;endendendf=xx=round(10*rand(1,10))qipaofa(x);(2)有一个4×5矩阵,编程求出其最大值及其所处的位置。functionf=zuidazhi(x)a=1;b=1;c=x(1,1);fori=1:4forj=1:5ifx(i,j)ca=i;b=j;c=x(i,j);endendendf=[c,a,b]x=rand(4,5)zuidazhi(x)(3)编程求201!nn。functionf=qiuhe(x)b=0;fori=1:xa=prod(1:i);b=b+a;endf=bqiuhe(20)(4)一球从100米高度自由落下,每次落地后反跳回原高度的一半,再落下。求它在第10次落地时,共经过多少米?第10次反弹有多高?functionf=gao(x)b=x;fori=2:10x=x/2;a=x*2;b=b+a;endf=[bx/2]gao(100)(5)有一函数yxyxyxf2sin),(2,写一程序,输入自变量的值,输出函数值。Functionf=fun(x)f=x(1)^2+sin(x(1)*x(2))+2*x(2)MATLAB上机作业21.求和4024441Y。symsks=4^k;S=symsum(s,k,0,40)2.求函数71862)(23xxxxf的极值,并作图。y='2*x^3-6*x^2-18*x+7';y_='-2*x^3+6*x^2+18*x-7';[x_min,y_min]=fminbnd(y,-7,7)[x_max,y_max]=fminbnd(y_,-7,7)ezplot(y)3.设25cos7sinxxxeyx,求??,22dxyddxdyy=exp(x)*sin(x)-7*cos(x)+5*x^2;dy=diff(y)d2y=diff(y,2)4.求积分dxxI2/012sin21。symsxy=sqrt(1-2*sin(2*x));F=int(y,0,pi/2)5.求积分DdxdyyxI222,其中D为直线xyxyxy12,2/,2围成的区域。symsxyf=x^2/y^2;int(int(f,y,x/2,x*2),x,0,4)+int(int(f,y,x/2,12-x),x,4,8)6.分别随机产生一个6×6的整数矩阵(元素可在]20,20[之间),求该随机阵的行最简形,秩,行列式,逆矩阵,列向量组的一个极大线性无关组,特征值和特征向量。a=round(-20+40*rand(6))rref(a)rank(a)det(a)inv(a)[V,D]=eig(a)7.求解方程组的通解246635554254843219526335544242364265432165432642165321654321654321xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxA=[1246-32;24-451-5;36205-9;230401;0-4-5214;55-366-4];b=[43-18-52]';B=[Ab];C=rref(B)A=[1246-32;24-451-5;36205-9;230401;0-4-5214;55-366-4];b=[43-18-52]';B=[Ab];n=6;R_A=rank(A)R_B=rank(B)formatratifR_A==R_B&R_A==nX=A\belseifR_A==R_B&R_AnX=A\bC=null(A,'r')elseX='equationhasnosolves'endsymsk1k2X=k1*C+Xpretty(X)8.求符号微分方程xxydxdycostan的通解和当y(0)=2的特解。eq='Dy+y*tan(x)=cos(x)';y1=dsolve(eq,'x')y2=dsolve(eq,'y(0)=2','x')MATLAB上机作业31、某校60名学生的一次考试成绩如下:9375839391858482777677959489918886839681799778756769688483817566857094848382807874737670867690897166867380947978776353551)计算均值、标准差、极差、偏度、峰度,画出直方图;2)检验分布的正态性;3)若检验符合正态分布,估计正态分布的参数并检验参数。x=[937583939185848277767795948991888683968179977875676968848381756685709484838280787473767086769089716686738094797877635355];mean(x)std(x)range(x)skewness(x)kurtosis(x)hist(x)h=normplot(x)[muhat,sigmahat,muci,sigmaci]=normfit(x)[H,sig,ci]=ttest(x,80.1)2、据说某地汽油的价格是每加仑115美分,为了验证这种说法,一位学者开车随机选择了一些加油站,得到某年一月和二月的数据如下:一月:119117115116112121115122116118109112119112117113114109109118二月:1181191151221181211201221281161201231211191171191281261181251)分别用两个月的数据验证这种说法的可靠性;2)分别给出1月和2月汽油价格的置信区间;3)给出1月和2月汽油价格差的置信区间。x=[119117115116112121115122116118109112119112117113114109109118];y=[118119115122118121120122128116120123121119117119128126118125];[H,sig,ci]=ttest(x,115)[H,sig,ci]=ttest(y,115)[h,sig,ci]=ttest2(x,y)1、在同一平面中的两个窗口分别画出心形线和马鞍面。要求:1、在图形上加格栅、图例和标注2、定制坐标3、以不同角度观察马鞍面a=3;b=1;theta=linspace(-2*pi,2*pi);rho=a*(1-cos(theta));subplot(1,2,1);polar(theta,rho)gtext('rho');gtext('theta');title('心形线');axis([-1010-1010]);gridon;x=-a:0.1:a;y=-b:0.1:b;[x,y]=meshgrid(x,y);z=x.^2/(a^2*2)-y.^2/(b^2*2);subplot(1,2,2);surf(z)xlabel('x轴');ylabel('y轴');zlabel('z轴');title('马鞍面');axis([0100030-11]);gridon2、以不同的视角观察球面x2+y2+z2=r2和圆柱面x2+y2=rx所围区域。[x,y,z]=sphere(200);i=find(x.^2+y.^2-x=0);z(i)=NaN;mesh(x,y,z)holdonezplot('x^2+y^2=x')
本文标题:MATLAB上机作业
链接地址:https://www.777doc.com/doc-2406129 .html