matlab期中 程序

1、 测试（第2,3,4题选做，其余题必做）1. 解答题 (1) 写出执行以下代码后，MATLAB命令窗口上显示的x矩阵的值 x=0,2,0,2,0,3,0,4; for k=1:8 if x(k)=0 x(k)=k; else x(k)=2*k+1; endenddisp(x);结果：1 5 3 9 5 13 7 17(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)结果：>> CC = 1 2 3 4 5 6 7 8 9 10 17 12>> DD = 2 3 5 6 8

2、9>> EE = 4 52.MATLAB7.0的符号计算功能计算下列题目（1）求在上的定积分；（2）求解非线性方程组；（3）计算级数；（4）设 求复合函数；（5）计算的.3. 利用MATLAB7.0的符号计算功能计算下列题目（1）求定积分；（2）求解非线性方程组；（3）计算级数；（4）求极限；（5）计算的导数。4.利用matlab7.0计算下列题目（1）某观测站测得某日6:00时至18:00时之间每隔2小时的室外温度()t=15,19,24,28,34,32,29，用3次样条插值求该日室外6:30至17:30之间每隔2小时各点的近似温度()；（2）采用8阶多项式对函数进行拟合；（3

3、）求在内的最小值点。5. 利用MATLAB7.0计算下列题目： 建立一个6阶魔方矩阵，（1）求的逆；（2）求的行列式；（3）求的秩；（4）求的迹；（5）求的特征值。解：6阶魔方矩阵A:>> A=magic(6)A = 35 1 6 26 19 24 3 32 7 21 23 25 31 9 2 22 27 20 8 28 33 17 10 15 30 5 34 12 14 16 4 36 29 13 18 11（1）A的逆:inv(A)结果：>> inv(A)ans = 1.0e+014 * -3.2169 -0.0000 3.2169 3.2169 -0.0000 -

4、3.2169 -3.2169 -0.0000 3.2169 3.2169 -0.0000 -3.2169 1.6084 0.0000 -1.6084 -1.6084 0.0000 1.6084 3.2169 0.0000 -3.2169 -3.2169 0.0000 3.2169 3.2169 0.0000 -3.2169 -3.2169 0.0000 3.2169 -1.6084 0.0000 1.6084 1.6084 0.0000 -1.6084（2）A的行列式：d=det(A)结果：>> d=det(A)d = 0（3）A的秩：rank(A)结果：>> rank

5、(A)ans = 5（4）A的迹：trace(A)结果：>> trace(A)ans = 111（5）A的特征值：eig(A)结果：>> eig(A)ans = 111.0000 27.0000 -27.0000 9.7980 -0.0000 -9.79806. 利用matlab7.0计算下列题目（1）某实验对一根长10米的钢轨进行热源的温度传播测试。用x表示测量点0:2.5:10(米)，用h表示测量时间0:30:60(秒)，用T表示测试所得各点的温度 T=95,14,0,0,0;88,48,32,12,6;67,64,54,48,42 ()。试用线性插值求出在一分钟内

6、每隔20秒、钢轨每隔1米处的温度TI；（2）随机产生一向量x，求向量x的最大值；（3）求在附近的根。解：（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=interp2(x,h,T,xi,hi)结果：TI = Columns 1 through 7 95.0000 62.6000 30.2000 11.2000 5.6000 0 0 90.3333 68.8667 47.4000 33.6000 27.4667 21.3333 16.0000

7、 74.3333 63.2667 52.2000 45.2000 42.2667 39.3333 33.2000 47.0000 45.8000 44.6000 46.0000 50.0000 54.0000 51.6000 Columns 8 through 11 0 0 0 0 10.6667 6.9333 4.8000 2.6667 27.0667 22.4667 19.4000 16.3333 49.2000 46.6000 43.8000 41.0000（2）>> x=rand();>> m=max(x)m = 0.4001 运用函数rand生成了一组10个在

8、01之间的随机数 X中最大值：max(x)结果：最大值为0.9830（3）命令如下：>> syms x>> fun=inline(x10x + 2,'x')>> x,f=fsolve(fun,0.55)结果：x = 0.3758f = -9.2419e-0117. 利用matlab7.0绘制以下图形解：（1）（1）在同一直角坐标系内画出和在区间的图形，添加x轴和y轴的简单注解，对图形进行文本标注，并尝试不同的线形、颜色；命令如下x=linspace(1,3*pi,1000); y1=x.*sin(x); y2=x.*sin(1./x); pl

9、ot(x,y1,'-r',x,y2,'-g') xlabel('x');ylabel('y=xsinx,y=xsin(1/x)'); （2）分别在同一图形窗口的不同子图绘制在区间上的图像。命令如下：t=0:pi/10:3*pi; subplot(2,2,1) plot(t,sin(t) subplot(2,2,2) plot(t,cos(t) subplot(2,2,3) plot(t,sin(2*t) subplot(2,2,4)plot(t,abs

10、(sin(t)8. 利用matlab7.0绘制以下图形： 解：（1） 在同一窗口分别用红色实线和绿色虚线绘制和在区间上的曲线，并用星号*标出两条曲线的交点以及建立图例.命令如下：x=linspace(0,2*pi,1000); y1=sin(x); y2=cos(x); k=find(abs(y1-y2)<1e-2); x1=x(k); y3= sin(x1);plot(x,y1,'-r',x,y2,'-g',x1,y3,'*')（2）分别在同一图形窗口的不同子图绘制在区间上的图像

11、。命令如下：t=0:pi/10:4*pi; subplot(2,2,1) plot(t,sin(t) subplot(2,2,2) plot(t,cos(t) subplot(2,2,3) plot(t,sin(2*t) subplot(2,2,4)plot(t,abs(sin(t)9 用if语句实现下列程序：输入一个百分制成绩，要求输出成绩等级为A、B、C、D，其中90-100分为A，80-90分为B，70-79分为C，70分以下为D。解：function function4()p=input('请输入分

12、数:'); if p>=90                disp('A'); elseif p>=80         disp('B');   elseif p>=70         disp('C'); else  disp('D');         end结果：10 若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。求1,1000之间的全部完数