matlab期中程序

1、测试(第2,3,4题选做，其余题必做)1.解答题(1) 写出执行以下代码后，MATLAB命令窗口上显示的x矩阵的值 x=0,2,0,2,0,3,0,4;for k=1:8if x(k)=0x(k)=k;elsex(k)=2*k+1;endend dis p(x);结果:39513717(2)写出执行以下代码后C, D , E的值A=1,2,3;4,5,6;7,8,9;C=A;10,17,12D=C(1:3,2 3)E=C(2,1 2)结果: C1012 D E2.MATLAB7.0的符号计算功能计算下列题目x(1 )求 exsin x 在3,18上的定积分;(2) 求解非线性方程组sin(x)

2、 2cos(y) 02cos(x) 3sin(y) 0 (3 )计算级数501S；n 1 n(4 )设 f Xyey 2t,求复合函数f (g(z)；(5)计算zxy32e xy 12x的二.x3.利用MATLAB7.0的符号计算功能计算下列题目(1 )求定积分 ：e dx ；求解非线性方程组x 0.7sin(x) 0.2cos(y)y 0.7cos(x) 0.2sin(y)计算级数S nS ；求极限 lim tanx 3sinxX 0sin X223Z(5)计算Z yx 3y x 3y的导数一。y4.利用matlab7.0 计算下列题目2小时的室外温度(C )(1 )某观测站测得某日6:00

3、时至18:00时之间每隔t=15,19,24,28,34,32,29用3次样条插值求该日室外6:30至17:30之间每隔2小时各点的近似温度(C )；(2)采用8阶多项式对函数 y x3 5sin(x),x0,10进行拟合;(3)求f(x) x3 2x 5在0,10内的最小值点。5.利用MATLAB7.0 计算下列题目:建立一个6阶魔方矩阵A,(1 )求A的逆;(2)求A的行列式;(3 )求A的秩；(4)求A的迹；(5)求A的特征值。解:6阶魔方矩阵A: A=magic(6)351626192433272123253192222720828331710153053412141643629131

4、8111.60840.0000-1.6084-1.60840.00001.60843.21690.0000-3.2169-3.21690.00003.21693.21690.0000-3.2169-3.21690.00003.21691.60840.00001.60841.60840.0000-1.6084(1 ) A 的逆:inv(A)结果: inv(A)ans =1.0e+014 *-3.2169-0.00003.21693.2169-0.0000-3.2169-3.2169-0.00003.21693.2169-0.0000-3.2169(2 ) A 的行列式：d=det(A) d=de

5、-t (A结果: d=det(A)0(3) A 的秩：rank(A)结果: ran k(A)ans =(4) A 的迹：trace(A)结果: trace(A)ans =111（5） A的特征值：eig（A）结果: eig(A)ans =111.000027.0000-27.00009.7980-0.0000-9.79806.利用matlab7.0 计算下列题目（1 ）某实验对一根长10米的钢轨进行热源的温度传播测试。用X表示测量点025:10（米），T=95,14,0,0,0;88,48,32,12,6;67,64,54,48,42( 秒、钢轨每隔1米处的温度TI;用h表示测量时间0:30:

6、60（秒），用T表示测试所得各点的温度C ）。试用线性插值求出在一分钟内每隔20（2 ）随机产生一向量 X,求向量x的最大值;解:(1)命令如下:x=0:2.5:10;h=0:30:60;T=95,14,0,0,0;88,48,32,12,4;47,44,54,48,41;xi=0:10;hi=0:20:60;TI=i nterp2(x,h,T,xi,hi)结果:TI =Columns 1 through 795.000062.600030.200011.20005.60000090.333368.866747.400033.600027.466721.333316.000074.333363

7、.266752.200045.200042.266739.333333.200047.000045.800044.600046.000050.000054.000051.6000Colu mns 8 through 11000010.66676.93334.80002.666727.066722.466719.400016.333349.200046.600043.800041.0000(2 ) x=rand(); m=max(x)0.4001运用函数rand生成了一组10个在01之间的随机数X中最大值：max(x)结果: 最大值为0.9830(3) 命令如下: syms x fun=in l

8、i ne(x -0x + 2 ,x) x,f=fsolve(fu n,0.55)结果:0.3758-9.2419e-0117.利用matlab7.0 绘制以下图形解:(1)( 1)在同一直角坐标系内画出y xsinx和y xsin1在区间1,3 的图形，添加xx轴和y轴的简单注解，对图形进行文本标注，并尝试不同的线形、颜色;命令如下x=li nsp ace(1,3* pi,1000);y1=x.*s in( x);y2=x.*si n( 1./x);plot(x,y1,-r,x,y2,-g)xlabel(x);ylabel(y=xs in x,y=xs in (1/x);(2)分别在同一图形窗

9、口的不同子图绘制y1 si nt, y2 cost, y3 si n 2t, y4 si nt在区间0,3 上的图像。命令如下:t=0: pi/10:3* pi;sub plot(2,2,1)plot(t,s in (t)sub plot(2,2,2)p lot(t,cos(t) sub plot(2,2,3) p lot(t,si n( 2*t) sub plot(2,2,4)plot(t,abs(si n( t)8.利用matlab7.0 绘制以下图形:解:(1)在同一窗口分别用红色实线和绿色虚线绘制y1 sinx和y2 cosx在区间0,2上的曲线，并用星号*标出两条曲线的交点以及建立图例命令如下:x=li nsp ace(0,2* pi,1000);y1=s in( x);y2=cos(x);k=fi nd(abs(y1-y2)=90dis p(A);elseif p=80dis p(B);elseif p=70dis p(C);else dis p(D);end结果:请输入分埶BEA请输入分數怕rB请输入分數M目C请输人分埶曲HD6=1+2+3 ，所以6是完数。10若一个数等于它的各个真因子之和，则称该数为完数，如求1,1000之间的全部完数。解:function fun cti on5()result=;for n=1:1