MATLAB上机作业.doc

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) 有一个45矩阵，编程求出其最大值及其所处的位置。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,bx=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/2gao(100)(5)有一函数，写一程序，输入自变量的值,输出函数值。Function f=fun(x)f=x(1)2+sin(x(1)*x(2)+2*x(2)MATLAB上机作业21. 求和 。syms ks=4k;S=symsum(s,k,0,40)2. 求函数的极值，并作图。y=2*x3-6*x2-18*x+7;y_=-2*x3+6*x2+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*x2;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=x2/y2;int(int(f,y,x/2,x*2),x,0,4)+int(int(f,y,x/2,12-x),x,4,8)6. 分别随机产生一个66的整数矩阵（元素可在之间），求该随机阵的行最简形，秩，行列式，逆矩阵，列向量组的一个极大线性无关组，特征值和特征向量。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_B&