matlab期中 程序.doc

测试（第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 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）求在内的最小值点。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 -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(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 ()。试用线性插值求出在一分钟内每隔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 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个在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); plot(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(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）分别在同一图形窗口的不同子图绘制在区间上的图像。命令如下：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('请输入分数:'); 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之间的全部完数。解：function function5()result=;for n=1:10000 divisors=1:(n/2); remainder=mod(n,divisors); sumDivisors=sum(divisors(remainder=0); if n=sumDivisors result=result,n; endendresult结果：result = 6 28 496 812810 果 , = )= ( = ) ( = ) : 00 ( 解数完的0数是以 =，完称和子个它个果     ;)(   0=   '(    = ;)(      =>分入 (   解 下0 分0 分-，分-中， 级出输，百个序列现 ) ( ,( )( (), ) ) ) ) , ; */:下像的间绘图口图同'*, ' ,- )   - ; ; ;0 * 下例立以交线两号用曲区制虚绿实红口解：形制绘 ) (, , ) ( ,), ) ) () ) ; 下像图上制子不图同别;/ ) ) ' -, ; ;)(*;) ( 如色颜的不尝本文形注的和加形的和内标直在（解图图绘 用 - . 果 ', + 0果 (： 机机 个了 0 0 ( ) 00 0. 0 0 . 0 . 0 . 0 0 00 0. .0 0 0 0. 0 0 00. 0 . 00 0 . 0. 0 0 0. 00 果 , , '0 0 0 , ;00,'0:0 下如解根的附值大 求 产机 的 每0每钟出求性) , , ;0, =温得测表)秒:0间表，(: 0示表试传的行钢0一验题题0 用 -00. 00.00 00 (果 特 果 ： ( 果 ( 秩 (=果 式的 - 00 . 00 . 0. .- 00 . 00 . - 00 . 00 . 0. . - 00- . 00 - 00 . 00 . 果 ( 的 0 阵解值征）；的 ； 式的求；阵方魔：下0 利点值内合行函项 ；温似各隔之: :室求插 度室的每时:至0:测观题列计 用数算限限数数组方线；分题下计计的0 用的计数函求数级组程性求分积在题题能计的 果,(, ( ,; ,;,; 的，码行写 果)( =) )( = : ;,0,值的的示令 ，代以 题解做题，题, （ ，题代 示,;: ) = 果 行， , ( 的题在性级函数 计下；组限数 观:室 求:之似 行值利0：；的 ）征阵 的( . 0 . - . . 0 0 .0 .0 . 0 = 果 ( 0.0 -用0验一行示:表0表温;, 性钟0的产 值解下:', 0', 果 . 0 0 0 00 0 0 0 . . 0 00 机 果 , 果 . 绘（直内和注本不色( *; - ) 别同制像 ) () ,) ,( 制：红制区两以 ; ; -  -',*图绘:/, ) ,)( )列输级，-   分>  = '  (;  果个称=以的 0 )= ,