MATLAB程序设计与应用(刘卫国编)课后实验答案

1、实验一 MAT LAB运算基础1、先求下列表达式得值，然后显示MAT LAB工作空间得使用 情况并保存全部变量.(2)，其中0.3a0.3ae e0.3 a(3)Z3 sin(a 0.3) In, a 3.0, 2.9, 2.9, 3.022(4)，其中 t=0： 0、5 : 2、5解：M文件：z 1 = 2 * s in (85 * p i /1 8 0)/ (1+exp ( 2 )x= 21 +2* i ; -、4 5 5;z2=1/2*log (x + sqrt(1+x A 2 )a= 3、0:0、1: 3、0;z 3= (ex p ( 0、3、*a) -exp( 0、3、* a)、/2

2、、* s i n(a+0、3) + 1 og( (0、3+a)、/2 )t=0: 0、5: 2、5;z4 =( t =0&t=1& t 2)、* ( t、 A 21) +( t=2 &t3) 、*(t、八 2 -2 *t+1)运算结果：z1=2*sin( 8 5 pi/1 8 0 )/ ( 1 + e xp ( 2 )x=2 1+2*i ; 、4 5 5;z2=1 /2 * l o g（x+sqrt（ 1 +xA2）a = 3、0: 0、1 : 3、0 ;z3= （exp（0、3、*a） -ex p （0、3、*a）、/2、*sin（a+0、3） + log（ （。、3+a）、/2 ）t =0

3、:0、5:2、5 ;z4=（t=0&t =1& t =2&t 1 541 0、0044iz3 =Co 1 u m n s 1 t hro ugh 40、7388 +3、14 16i0、7696 + 3、1416i0、787 1 + 3、14 1 6 i 0、7 913 + 3、141 6 iColumns 5 t h rough 88、I 980 %LlOO、 乙 q6no jqi 62 suuu n |og9PZ0、氏一! 9L 打 L + L 86、乙一 ！ 9L 打 L 、 + 乙/乙、乙-!9L 打 L 、+9 % 、I 8 2 M nojq 】g乙 s uuun|oo! 9 M、 +

4、 L，1-!9 l 9、 + T OZ 3、l-! 9WV + 260 1 9Vn、 + 9 ，8、0 PZ q6noq 】 suuunio + S9Z 9 、0!9L 打 I 、 +EH，、0一 !9M、+ 99，、0一 !96、 + %C、0OS qBnoj q 】Z.L s u uuniog! 9配 T、 + VLLQ、0 !9 I h + 2 6,0、0 9Vn、+ 0891，、09Vn + 9Z Z 2、09 i6noj q 1 L s u uun| o 9!96、 + 082 、0!9M、+ 8 89。、0!9 T Pi、 + 96,9、0! 9M +9 6L9、0口 qBnojq

5、i 6 suuunio !9降 I、 + P819、0 !9LVL、+ P9Z L、0! 9L 打 I、 + 乙092、0! 9L正、 + ZZQL、0P1P9 、009、09 q q6n oqi f 9sum n i og91，89q6 n oj999 、06 Z9 、0suuun i og992、0gt? qBnojqi 处 suinniog乙620、09920、0 LZV、。-PPq 6 n oq , SUlLin|096L82 、0 1,91,9、0 Z999 、。-0t7 q6 n o j q 4 氏 suuun o 9808、0 6、0-9 qBnoj qi suuun|ooIL

6、68Z69、I -Columns 5 7 throu g h 6 00、5 7770、53270、4 7740、 4126Column 6 10、 3388z 4 =00、2 5 0001、25001、0 0 002、25002 、 已知：求下列表达式得值：(1) A+6*B= A& htf按 IE RT0 SIaa-(E,Jti-EtS30.3,2110ATMkigHLr. ID C*l 6输出结果：S =1002A =1、0 0 00000、5 3830、442701 、 000000、99 6 10、1067001、00000、0 7820、96190001 、0000000002、00

7、00a1、0 00 0001 、0 7 6 71、3 28001、 0 0 0 001 、992 30、32001 00000、15 642 、 88570000an s =0000 000000000000000000000001 00 0 0004 、00 0 0由an s,所以2、产生5阶希尔伯特矩阵 H与5阶帕斯卡矩阵 P,且求其行列式得值 H h与Hp以及它们得条件数Th与 Tp,判断哪个矩阵性能更好。为什么 ？解：M文件如下：输出结果：H =1 、 0 0000、5 0 000、33330、 250 00、 20 0 00、 50000、 33 3 30、 25000、2 0 00

8、0、16 6 70、3 3 3 30、 25000、 20000、 16670、 1 4 290、 1 2 500、 2 0000、1 6 670 14 2 90、1 2 5 00、1 1 11P = 11111123451361 01 50、 2 5 000、2 0 0 00、 1 42 90、 166 72 03 51410151535703、7493e012H p =1Th =4、7661e+0 0 5T p =8、 5 175e+00 3因为它们得条件数 T h Tp ,所以pascal矩阵性能更好3、建立一个5X5矩阵，求它得行列式值、迹、秩与范数解： M文件如下：输出结果为：A =

9、1724181 52357141 6461320221 01219213d =1 11 82 5295070 000t =6 5cl =6、 8 500c2 =5、 4618cinf =6、 8 50 04、 已知求A得特征值及特征向量，并分析其数学意义。解：M文件如图：输出结果为0、7 1 3 00、2 803-0、6084 0、78670、 27330、 87250、 34870、 55010、4 0 5 0-2 5、3 16 9010、 5 18 2016、 83 5 1数学意义：V得3个列向量就是 A得特征向量，D得主对角线上3个就是 A得特征值，特别得,V得3个列向量分别就是 D得3

10、个特征值得特征向量5、下面就是一个线性方程组:(1)求方程得解.(2)将方程右边向量元素b 3改为0、53再求解，并比较b 3得变化与解得相对变化。(3) 计算系数矩阵A得条件数并分析结论。解：M文件如下： Edit or - Untitled?*匚 I XEdiP ext 最 C 艺 115 Debug： 口 exlct 口 p tind 口 vr Hwl 口 .乂日。|3一 0可自，|* ，强IIB 帕 国 希回忑 :唱喟-卜.,。力 *峨党jq1 clcEA=L/2 3 4;345;4 5 6；3b=a+95;0.6ir(D. 52,4X=uw(A)*b% (D间的方程解-5gb2=0.

11、96;D4 S7;0.53;TX2=iny(A)*b*(2)向b3改变后的方程解一Q9C=cond(A)肮十算矩阵2-A的条件颊-10scriptlit 9 Cal 39输出结果：X =1、2 0 000、 6 00 00、 6 00 0X2 =1、2 0 000、60 0 00、6 0 00C =1、 3 533e +003由结果,X与X 2得值一样，这表示b得微小变化对方程解也影响较小，而A得条件数算得较小，所以数值稳定性较好，A就是较好得矩阵。6、 建立A矩阵，试比较 sqrt m ( A)与s q rt (A)，分析它们得区别。解：M文件如下：运行结果有:A =166182051 29

12、853、 8891-0、 11023、 2 1 0 33、 2 9 172、1 4 360、3 6980、38 5 52、07601 73 0 54、0 0 002、 44954、 24264、 47212 236 13、4 6413、00002、 8 2842、 2 36116、000 06、00001 8、000 020、 0 0 005、00001 2、0 0 009、 0 0 008、0 0 005 000 0bb 1b2分析结果知：s qr t m (A)就是类似A得数值平方根(这可由bl * b 1 = A得结果瞧出)而sqrt (A)则就是对A中得每个元素开根号 ，两则区别就在于

13、此。实验三选择结构程序设计、实验目得1、掌握建立与执行M文件得方法。2、掌握利用i f语句实现选择结构得方法 .3、掌握利用s wi t ch语句实现多分支选择结构得方法4、掌握try语句得使用。、实验内容1、求分段函数得值。用 if 语句实现，分别输出 x =-5、0, -3、0, 1、0,2、0,2、5,3、0,5、0 时得 y (1.解:M文件如下：运算结果有：f( 5)y =14 f ( 3)y =11 f (1)y =2 fy =1 f (2、5)y =0、2500 f (3)y =5 f (5) y =192、 输入一个百分制成绩，要求输出成绩等级A、B、C、D、E。其中90分10

14、0分为A,80分 8 9分为 B, 7 9分7 9分为C,60分69分为D , 6 0分以下为 E.要求：(1) 分别用if语句与s wit c h语句实现。(2)输入百分制成绩后要判断该成绩得合理性，对不合理得成绩应输出出错信息.解:M文件如下试算结果：score= 8 8grade = Bscore= 1 23错误:输入得成绩不就是百分制成绩3、 硅谷公司员工得工资计算方法如下:1 1)工作时数超过120小时者，超过部分加发15%。2 2)工作时数低于60小时者，扣发700元。(3)其余按每小时84元计发。试编程按输入得工号与该号员工得工时数，计算应发工资。 解：M文件下4、设计程序，完成

15、两位数得加、减、乘、除四则运算，即产生两个两位随机整数，再输入一个运算符号，做相应得运算，并显示相应得结果。解：M文件如下;运算结果例：a =3 8b =3 3输入一个运算符：Ac =falsea =9 2b =4 0输入一个运算符：+c =1 3 25、 建立5X6矩阵,要求输出矩阵第 n行元素。当n值超过矩阵得行数时，自动转为输出 矩阵最后一行元素，并给出出错信息。解：M文件如下：1 1 E di tor 一 Uni tie d8*Fil* Edit Text 典 Qtll Tooli 山如工 Baiktop Window Hlpn | - X：Q d* fej 3i HiQ日 国|1.Q

16、+ 占卜1x 噱森| Q1elc2Af input C输入一个5行矩阵A=);3 n= input C 输入1正、数;4 if nD5 disp(Atrij :),Gelse7 disp(A(5, l) ;|8 end9 lasterr1011121314151617scriptLa 7 Cel 18运算结果如下:输入一个5行6列矩阵A= 1 23 4 5 6 7输入一正整数 n=41123输入一个5行6列矩阵A=13; 23 45 6 7输入一正整数n=62342 3 4 5 5;2972 3 45 5 ;23 4 5 7 6;2 2 2 2 2 3; 11 2 397 3;33 4 5 7

17、 6; 222 2 2 3;11 2 3 9 7a n s =Error u s i ng = = dis pToo many input argume n ts、实验四循环结构程序设计一、实验目得1、掌握利用f。r语句实现循环结构得方法.2、 掌握利用wh 1 le语句实现循环结构得方法。3、熟悉利用向量运算来代替循环操作得方法。二、实验内容1、根据，求兀得近似值。当 n分别取10 0、1000、100 0 0时,结果就是多少？要求:分别用循环结构与向量运算(使用s um函数)来实现。解：M文件如下：卷-UntitledSr匚：叵|冈工件如)编辑 文本9 跳转&)单元 力十，X :na哈力目

18、，回*回工i3 q56T8 g10111213;驾唱1不1+JT不燃病级环结构计算口 i值|匚y=o：input ( n= )；f l=1;nk y+l/i/i;号口 dpi=s(irt (6*y)先向量方法i十算P i值n=input ( n=)；i= 1. / (1 ： n)2;B=5unt(i);pi=sqrt脚本17 13 列13运行结果如下：K %循环结构计算p i值y=0;n= i npu t (n=);for i=1 : ny = y +1/i/i ;endp i =sq r t ( 6*y)n= 100p i =3、 132 1n= 100 0p i =3、 1 4 06n=

19、1000 0p i =3、1 4 15%向量方法计算P i值n=input( n=1 = 1、/ (1 ： n)、A2;s=sum (i);pi= s qrt (6 * s)n=1 0 0p i =3 、 1 3 21n=1 0 00 p i =3、 140 6n=1 0 000pi =3、 14152、根据，求：(1 ) y y=0 ; n=0 ;while y 3n =n+1;y = y+1/( 2 * n-1); endy n if y 3n=n 1;e ndny =3、00 3 3n =575 63、考虑以下迭代公式:其中a、b为正得学数。（1 ）编写程序求迭代得结果，迭代得终止条件为

20、I Xn+i xn | w 1 0 -5,迭代初值xo=1、 0,迭代次数不超过 5 0 0次.（2）如果迭代过程收敛于r ，那么r得准确值就是，当（a, b）得值取（1 , 1）、（8, 3）、 （10,0、1）时，分别对迭代结果与准确值进行比较。解：M文件如下：运算结果如下；请输入正数a=1请输入正数b=1X =0、6180r =0、61804、7 016r =0、6 1 80-1、6180s =-0、 0 0 0 02、2 3 6 1请输入正数a = 8请输入正数b=3x =1、7 0 16r =1、70161、61 8 0r =1、7 01 64、7 0 16s =0.0-6、4 0

21、3 1请输入正数a = 1 0请输入正数b=0、1x =3、 1 1 27r =3、1 1 274、7 0 16r =3、112 7 3、2 12 7s =-0、00006、32544、 已知求f 1f 100中:(1 )最大值、最小值、各数之与。(2)正数、零、负数得个数。解:M-文件以下就是运算结果： m ax (f)=437 7 6 3 2 8 2 635 min ( f )= 89941211 3 52 8 sum ( f) = - 7 427 4 5601 9 51 cl = 49 c2= 2 c 3 = 495若两个连续自然数得乘积减1就是素数，则称这两个边疆自然数就是亲密数对,该

22、素数就是亲密素数.例如，2X3-1 =5,由于5就是素数，所以2与3就是亲密数，5就是 亲密素数。求2, 5 0区间内：(1 )亲密数对得对数.(2) 与上述亲密数对对应得所有亲密素数之与。解：M文件：运算结果为：j =29s = 2361 5实验五函数文件一、实验目得1、 理解函数文件得概念.2、掌握定义与调用 MATLAB函数得方法。二、实验内容1、定义一个函数文件，求给定复数得指数、对数、正弦与余弦，并在命令文件中调用 该函数文件。解：M文件如下: 函数f u shu、M文件：function e,l,s,c= fu s h u(z)%fush u 复数得指数，对数，正弦，余弦得计算%e

23、复数得指数函数值% 1复数得对数函数值% s复数得正弦函数值%c复数得余弦函数值e= c x p(z);l= 1 og (z)；s= s in(z);c=cos(z );命令文件M：z =input (请输入一个复数z=/ );a, b,c,d = fushu (z)运算结果如下：z= i npu t (请输入一个复数 z=);a, b ,c,d= f ushu (z)请输入一个复数z = 1+ia =1、46 8 7 + 2、2 87 4 ib =0、3466 + 0、7854ic =1、2985 +0 635 0 id =0、8337 0、9 8 8 9 i2、 一物理系统可用下列方程组来

24、表示：从键盘输入 m1、m2与。得值，求a 1、a2、Ni与N2得值.其中g取9、8,输入。时以角度为 单位。要求：定义一个求解线性方程组AX =B得函数文件，然后在命令文件中调用该函数文件。解：M文件函数fc、M文件：func t io n X= f c ( A, B )% fcfc就是求解线性方程得函数% AA就是未知矩阵得系数矩阵X=AB;命令M文件：c 1 c;ml = input(/ 输入m 1=);m2=in p ut(/ 输入m 2=z );theta = inpu t (/ 输入 the t a =/ );x= t het a * p i/18 0 ;g=9、8;A=m1*co

25、s (x) m1 - s in (x)0m 1 * sin (x) 0 cos (x) 00 m2 -sin (x)00 0 cos ( x ) 1 ;B= 0;m1*g ; 0; m2 * g;X= f c (A,B)isisn输入m1= 1输入m2=1输入 t h eta=30X =7、840 03、 39486、 789615、 68003、一个自然数就是素数，且它得数字位置经过任意对换后仍为素数。例如13就是绝对素数。试求所有两位绝对素数。要求：定义一个判断素数得函数文件。解：M文件：函数prime、m文件fun c tion p = prime ( p )% 输入p得范围，找出其中得

26、素数m=p (length (p);for i=2 : sqrt (m )n=find( rem (p,i) = =0 &p1=i);P (n 户口 ;%将中能被i整除，而却不等于i得元素，即下标为 n得元素剔除，其余得即为素数 e ndP ;命令文件：clc;p=10: 99;p= p r i me ( p ) ;%找出1。到99内得所有素数p = 1 0*rem (p, 10 )+(p - rem ( p,10) ) /1 0 ;%将p素数矩阵每个元素个位十位调换顺序p=prime ( p)%再对对换后得素数矩阵找出所有得素数运算结果:p =1 1317 11 37 3173797794、

27、设，编写一个 MATLA B函数文件fx、m,使得调用f(x)时，x可用矩阵代入，得 出得f ( x )为同阶矩阵.解：函数fx、m文件：f u nctio n f= fX(x)%fxfx求算x矩阵下得f(x)得函数值A=0、1+ (x-2)、人 2;B=0、0 1+ ( x 3)、A4;f=1、/A +1、/B;命令文件：c 1 c;x= input (输入矩阵 x=/ );f=fx ( x )运算结果: x= i np u t C 输入矩阵 x= /);f=fx(x)输入矩阵x=7 2 ; 1 2 5f =0、 043710、 99 0 10、01010、17 2 45、 已知(1) 当

28、f (n) = n + 10ln(n 2+5 )时，求 y 得值.(2) 当 f(n) =1X2+ 2X3+3X4 +、+nX (n+1)时，求 y 得值.解：(1)函数f、m文件：fu n c tio n f=f ( x )f=x+10*1 o g ( xA2+ 5 );命令文件：c 1 c;nl = input( ，n1= / );n 2 = input ( n2=z )；n3=input ( n 3 =);y1=f (n1 )；y 2 =f ( n2 );y 3 =f(n3)；y = y1/(y2+y3)运算结果如下：n1 =40n 2 = 3 0n3=2 0y =0、6 39 0(2)

29、、函数g、m文件f u n c tion s = g(n ) for i = 1: ng (i ) = i * ( i+1 );ends =s u m (g);命令文件：c lc ;n 1 =i n put( / n 1= / ); n2= input (n2=);n 3=input ( n3 =);y 1 =g ( n 1 );y 2= g ( n2);y 3 = g( n 3);y=y1/(y2 +y3)运算结果如下：n 1 = 40n 2 =3 0n3=2 0y =1、 76 6 2实验六高层绘图操作、实验目得1、 掌握绘制二维图形得常用函数。2、掌握绘制三维图形得常用函数。3、掌握绘制

30、图形得辅助操作。、实验内容1、设，在x=02兀区间取101点，绘制函数得曲线。解：M文件如下：clc ；x=linsp a ce( 0, 2 * pi,101)；y= (0、5 + 3* s i n(x )、/( 1 + x、人 2); plot ( x , y )运行结果有：2、 已知 y1=x2, y2 = cos (2x),y 3 =y1 x y2 ,完成下列操作:(1)在同一坐标系下用不同得颜色与线型绘制三条曲线。(2)以子图形式绘制三条曲线。(3)分别用条形图、阶梯图、杆图与填充图绘制三条曲线.解：(1 ) M文件：c 1 c ;x= - p i :pi/10 0 :pi;y1 =

31、x、A2;y 2=cos (2 * x);y 3 = y 1、*y2;plot (x, y1, b x, y2,r:, x, y3, k)运行结果:(2) M文件：c lc;x= pi:pi/10 0: pi;y 1 =x、人 2 ;y2 = CoS (2*x)；y3=y1 、 * y2 ；subplot(1 , 3 , 1 )； pl o t (x , y1, b -);tit 1 e ( y1=x 人 2)；subplot (1, 3, 2)； p lot ( x ,y 2 , r :);t i t le ( y2= cos(2x ); s u b p lot(1 ,3,3);plot(x

32、 , y3 , k )； t itle ( y3=y 1 * y2 / )；、运行结果:(3) M文件：c lc ;x =-pi: pi/ 1 00 ： pi ；y 1 =x、人 2 ;y2= c o s ( 2 * x);y 3 =y1、*y 2 ;sub plot(2 , 2 , 1 );plot (x,y 1 , b / ,x , y2 , r : ,x, y 3 , k-； su b p 10t(2,2,2);b ar (x, y 1 , b)；t itl e ( y1 = xA 2 );s u bp 1 o t ( 2 ,2, 3 );b a r(x, y 2, J );ti t l

33、e (y2= cos (2x )；subpl o t ( 2, 2, 4 );bar (x , y3 , / k);t it 1 e( y3 = y1 * y2)；由上面得M文件，只要依次将“ bar”改为stairs、“stem”、“fill”,再适当更改区间取得点 数，运行程序即可，即有下面得结果：3、 已知 在一5x01: 5;y=(x + sqrt ( p i )/( e xp (2) 、* ( x = 0 ) + 0、5*log ( x+ s qrt (1 + x、人 2)、* ( X 0 );p 1 ot(x,y)运行结果:Fil Edit Vi ew Instrt. Tools

34、Dtskt&p W i ndow Kelp in m tJ以* fe ，等货4口国由图可瞧出，函数在零点不连续.4、 绘制极坐标曲线p =asin（b+n。），并分析参数a、b、n对曲线形状得影响解：M文件如下：clc;theta=0:pi /100: 2 *pi;a = inp u t （输入 a=）;b=input （/ 输入 b=/）;n= i n p u t（输入 n =）;r ho = a * sin （b+n * t h eta）;pol a r（t h eta , rho , m）采用控制变量法得办法，固定两个参数，变动第三个参数观察输出图象得变化。fu 13 ImtI IiiU

35、 lukl- Lidn ItliOd d % CSkX-|9 00 口目口!； Ul f X，l Iida 1阳.1，。1U9 J d j b、，$X/”。目口皿！j|匕1 1， $- 3n 17 1似” mU Pb,1CU； Q -S)、10 号口”上一我, g 里必 giUir 1，2d j “ ?09，、/；|石。口| 口J U d i k ?) k B B分析结果：由这8个图知道，当a, n固定时，图形得形状也就固定了，b只影响图形得旋转得角度；当a小固定时，n只影响图形得扇形数，特别地，当n就是奇数时，扇叶数就就是n,当就是偶数时，扇叶数则就是2n个；当b, n固定时，a影响得就是图

36、形大小，特别地，当a就是整数时，图形半径大小就就是 a。5、绘制函数得曲线图与等高线。其中x得2 1个值均匀分布5,5 范围，y得31个值均匀分布在0,10,要求使用subp 1 o t (2,1 , 1)与s u bplo t (2, 1 ,2)将产生得曲面图与等高线图画在同一个窗口上。解：M文件：c 1 c;x =linspace(-5 , 5,21);y =l i n s pac e (0, 10, 3 1);x, y = meshgr id (x, y);z =cos (x)、*cos (y)、*e xp (-sqrt (x、人 2 + y、A2)/4);su b p 1 ot(2 ,

37、 1, 1);su r f ( x , y,z);title (/ 曲面图)；s ub p lot(2 , 1 , 2);s u r fc(x , y,z);t itle(等高线图/ );运行结果:6、绘制曲面图形，并进行插值着色处理。解:M文件：cl c ;S =0: p i/1 0 0:p i / 2 ;t=0: p i/100 : 3* p i / 2;s ,t =me s hgrid(s, t);x= cos( s )、* cos (t);y= cos( s)、* s in(t);z=sin(s);subplot(2,2, 1);me s h ( x,y , z);t i tle(未着

38、色得图形);s u bplot( 2 , 2, 2);s ur f (x ,y,z);title ( shading f a ceted(缺省)；subplo t (2, 2, 3 );surf (x, y,z)； shading f 1 at；title( hadi n g flat);su bp lot(2 ,2,4);surf (x, y, z); shading i nterp; t it 1 e(/ s hading i n te r p ;)运行结果有:实验七低层绘图操作二、实验内容1、建立一个图形窗口，使之背景颜色为红色，并在窗口上保留原有得菜单项，而且在按下鼠标器得左键之后显示

39、出Left Bu tton Pressed字样.解:M文件如下：clc;hf= f igu r c ( col or, 1 0 0,、Wind o w B utto n Dow n Fcn, / dis p ( eft But ton Press e dz );运行结果:左击鼠标后:2、先利用默认属性绘制曲线 y = x2e2x,然后通过图形句柄操作来改变曲线得颜色、线型与线宽，并利用文件对象给曲线添加文字标注解：M文件:cl c ;x=- 2: 0、01 : 2;y=x、A2、* e x p ( 2*x);h= p lo t(x ,y);set( h , color； 0、4 ,0、2,0、

40、5,linesty 1 e, z ,、/ lin e w idt h ,2);，,9;text (1、5,1、5A2*exp(2*l、 5)J1 e ftar row xA2exp(2x) z,z f ontsiz e3、利用曲面对象绘制曲面v (x,t)=1 0 e-0 01 x sin( 2 000 7tt- 0、2x+ 兀)。解:M文件：cl c ;x=0 : 0、1:2 * pi;x, t =m e s h gr i d(x);v = 10*ex p ( 0、01 * x 卜 *sin(2 0 0 0* p i * t -0、2* x +pi);axes (view -3 7, 30)

41、;hs= surf a ce( x , t , v,z f a cecol o r，,、0、 2,0、 3,0、 3 , e dg e color ,z flat);g r i d on ;xlabe 1 (x axi s ) ; y label ( y -ax i s);zlabel (/ z - axi s );titl e (/ mesh s u r V);pause%按任意键继续set (hs,兄 c eColor flat );text(0 , 0,0/ 曲面);运行结果:File Edi t fiew Insert Tools Desktop Window HlpEru按任意键继续

42、:4、以任意位置子图形式绘制出正弦、余弦、正切与余切函数曲线。5、生成一个圆柱体，并进行光照与材质处理。解：M文件：x,y,z =cylin d er (3,500); %c ylind e r 就是生成柱体得函数 s urf(x, y , z);t i tle (圆柱体得光照与材料处理)；X la b el (/ X-axis);丫 1 ab e 1 ( Y -axis7);Zla b el ( Zf x is);axis( 5, 5, -5,5, 0,1)gr i d off;lig h t (/ Colo r r Po s i tio n -4 ,0,0,s t y le, infini

43、te);shadi n g i n t e rp；material s h i ny;view(0, 10);li g hting ph o ng;axis off；运行结果:实验八数据处理与多项式计算一、实验目得1、掌握数据统计与分析得方法 .2、掌握数值插值与曲线拟合得方法及其应用。3、 掌握多项式得常用运算.二、实验内容1、利用MAT LAB提供得r and函数生成300 0 0个符合均匀分布得随机数，然后 检验随机数得性质：(1)均值与标准方差。(2)最大元素与最小元素。(3 )大于0、5得随机数个数占总数得百分比。解：M文件：c 1 c ;x= r an d (1 , 30000);mu=mea n ( x) %求这3 0 000个均匀分布随机数得平均值sig=std (x) %求其标准差(T iy= l e n g t h(find(x0、5) ; %找出大于 0、5 数得个数p = y/3 0O00%大于0、5得所占百分比运行结果:mu =0、499 4 885532 3 104 3sig =0、2 8 859 9 933559786P =0、4 9 9 42) 将1 0 0个学生5门功课得成绩存入矩阵 P中，