MATLA上机作.docVIP

MATLAB上机作业1对以下问题,编写M文件：(1) 用起泡法对10个数由小到大排序。即将相邻两个数比较,将小的调到前头。function f=qipaofa(x)for j=9:-1:1for i=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矩阵，编程求出其最大值及其所处的位置。function f=zuidazhi(x)a=1;b=1;c=x(1,1);for i=1:4 for j=1:5 if x(i,j)c a=i;b=j;c=x(i,j); end endendf=[c,a,b]x=rand(4,5)zuidazhi(x)(3) 编程求。function f=qiuhe(x)b=0;for i=1:x a=prod(1:i); b=b+a;endf=bqiuhe(20)(4)一球从100米高度自由落下,每次落地后反跳回原高度的一半，再落下。 求它在第10次落地时，共经过多少米?第10次反弹有多高?function f=gao(x)b=x;for i=2:10 x=x/2; a=x*2; b=b+a;end f=[b x/2]gao(100)(5)有一函数，写一程序，输入自变量的值,输出函数值。Function f=fun(x)f=x(1)^2+sin(x(1)*x(2))+2*x(2)MATLAB上机作业21. 求和 。syms ks=4^k;S=symsum(s,k,0,40)2. 求函数的极值，并作图。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. 设，求y=exp(x)*sin(x)-7*cos(x)+5*x^2;dy=diff(y)d2y=diff(y,2)4. 求积分 。syms xy=sqrt(1-2*sin(2*x));F=int(y,0,pi/2)5. 求积分 ，其中D为直线围成的区域。syms x yf=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的整数矩阵（元素可在之间），求该随机阵的行最简形，秩，行列式，逆矩阵，列向量组的一个极大线性无关组，特征值和特征向量。a=round(-20+40*rand(6))rref(a)rank(a)det(a)inv(a)[V,D]=eig(a)7. 求解方程组的通解 A=[1 2 4 6 -3 2;2 4 -4 5 1 -5;3 6 2 0 5 -9;2 3 0 4 0 1;0 -4 -5 2 1 4 ;5 5 -3 6 6 -4 ];b=[4 3 -1 8 -5 2];B=[A b];C=rref(B)A=[1 2 4 6 -3 2;2 4 -4 5 1 -5;3 6 2 0 5 -9;2 3 0 4 0 1;0 -4 -5 2 1 4 ;5 5 -3 6 6 -4 ];b=[4 3 -1 8 -5 2];B=[A b];n=6;R_A=rank(A)R_B=rank(B)format ratif R_A==R_BR_A==nX=A\belseif R_A==R_BR_AnX=A\bC=null(A,r)else X=equation has no solvesendsyms k1 k2X=k1*C+Xpretty(X)8. 求符号微分方程 的通解和当y(0)=2的特解。eq=Dy+y*tan(x)=cos(x);y1=dsolve(eq,x)y2=dsolve(eq,y(0)=2,x)MATLAB上机作业31、某校60名学生的一次考试成绩如下:93 75 83 93 91 85 84 82 77 76 77 95 94 89 91 88 86 83 96 81 79 97 78 75 67 69 68 84 83 81 75 66 85 70 94 84 83 82 80 78 74 73 76 70 86 76 90 89 71 66 86 73 80 94 79 78 77 63 53 551)计算均值、标准差、极差、偏度、峰度，画出直方图；2)检验分布的正态性；3)若检验符合正态分布，估计正态分布的参数并检验参数。x=[93 75 83 93 91 85 84 82 77 76 77 95 94 89 91 88