1、Matlab 程序设计上机报告第 2 章1.求下列表达式的值。（1） )103425.(w6 w=sqrt(2)*(1+0.34245*10(-6)w =1.4142（2） ，其中acbex)tn(2a=3.6,b=5,c=-9.8 a=3.5;b=5;c=-9.8;x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2)/(tan(b+c)+a)x =0.9829（3） ,其483.0412ay中 =3.32， =-7.9。a=3.32;b=-7.9;y=2*pi*a2*(1-pi/4)*b-(0.8333-pi/4)*a)y =-128.4271（4） ，其中 t=221lntt
2、ez。65.031i t=2,1-3i;5,-0.65;z=0.5*exp(2*t)*log(t+sqrt(1+t*t)z =1.0e+04 *0.0057 - 0.0007i 0.0049 - 0.0027i1.9884 - 0.3696i 1.7706 - 1.0539i2.已知 023518,7613045BA求下列表达式的值（1）A+6B 和 A2-B+I（2）A*B、A.*B 和 B*A（3）A/B 及 B/A（4）A,B 和A(1,3,:);B2 A=-1,5,-4;0,7,8;3,61,7; B=8,3,-1;2,5,3;-3,2,0; A+6*Bans =47 23 -1012
3、 37 26-15 73 7 A2-B+eye(3)ans =-18 -217 1722 533 10921 867 526 A*Bans =14 14 16-10 51 21125 328 180 A.*Bans =-8 15 40 35 24-9 122 0 A.*Bans =-8 15 40 35 24-9 122 0 A/Bans =1.2234 -0.9255 2.9787-0.9468 2.3511 -0.95744.6170 3.8723 13.8936 B/Aans =-5.3173 -3.5661 0.8942-1.2768 -0.4743 0.24112.1229 1.31
4、73 -0.2924 A,Bans =-1 5 -4 8 3 -10 7 8 2 5 33 61 7 -3 2 0 A(1,3,:);B2ans =-1 5 -43 61 773 37 117 37 13-20 1 93.已知14.354.9620326178.A完成下列操作：（1）输出 A 在10，25 范围内的全部元素。（2）取出 A 前 3 行构成矩阵 B，前两列构成矩阵 C，右下角 32 只矩阵构成矩阵D，B 与 C 的乘积构成矩阵 E。（3）分别求表达式 E find(A=10 P=pascal(5); Hh=det(H)Hh =3.7493e-12 Hp=det(P)Hp =1 T
5、h=cond(H)Th =4.7661e+05 Tp=cond(P)Tp =8.5175e+03Tp 比 Th 好，更接近于 1 。5.已知： 582069A求A的特征值及特征向量，并分析其数学意义。 A=-29,6,18;20,5,12;-8,8,5;V,D=eig(A)V =0.7130 0.2803 0.2733-0.6084 -0.7867 0.87250.3487 0.5501 0.4050D =-25.3169 0 00 -10.5182 00 0 16.8351设A为n阶矩阵，若存在常数及非零的n维向量x，使得Ax=x，则称是矩阵A的特征值， x是A属于特征值的特征向量第 3 章
6、1.从键盘输入一个 3 位整数，将它反向输出。 a=input(输入一个 3 位整数:,s);fliplr(a)输入一个 3 位整数:345ans =5432.输入一个百分之成绩，要求输出等级A、B、C 、D、E。其中 90100 分为A，8089 分为 B,70 分79 分为 C,60 分69分为 D,60 分以下为 E。要求：（1）分别用 if 语句和 switch 语句实现。（2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。 m=input(请输入一个百分制的成绩: );if m=0 if m=90a(i)=input(s);endb=a(1);c=a(1);fo
7、r i=1:20if a(i)bb=a(i);endif a(i) max(a) min(a)5. ,当 a23.0ln).si(23.0. aeya 取-3.0、-2.9 、-2.8 、-2.8、2.8、2.9、3.0 时，求各点的函数值。要求分别用顺序结构和循环结构实现。 a=-3:0.1:3;b=(exp(0.3.*a)-exp(-0.3.*a)./2);c=sin(a+0.3);y=b.*c+log(0.3+a)./2) a=-3:0.1:3;for i=1:61y=(exp(0.3*a(i)-exp(-0.3*a(i)/2*sin(a(i)+0.3)+log(0.3+a(i)/2)E
8、nd5.当 n 分别取 100、1000、10000 时，求下列格式的值：要求分别用顺序结构和向量运算（使用 sum函数）来实现。（1） )6.(1.3222n A=0,0,0;n=100,1000,10000;for i=1:3for j=1:n(i)A(i)=1/(j2)+A(i);endendAA =1.6350 1.6439 1.6448 A=0,0,0;n=100,1000,10000;for i=1:3j=1:1:n(i);A(i)=sum(1./j.2);endAA =1.6350 1.6439 1.6448（2） 2.)1(2.756341 nA=1,1,1;n=100,100
9、0,10000;for i=1:3for j=1:n(i)A(i)=A(i)*(2*j)2)/(2*j-1)*(2*j+1);endend AA =1.5669 1.5704 1.57086.建立 56 矩阵，要求输出矩阵第 n 行元素。当n 值超过矩阵的行数时，自动转为输出矩阵最后一行元素，并给出出错信息。 A=rand(5,6)n=input(输入第几行？ );if n5n=5;disp(超出范围，输出最后一行的元素);A(n,:)elseA(n,:)end第 4 章1.绘制下列曲线（1） !3xy x=0:pi/100:2*pi; y=x-x.3/factorial(3); plot(x
10、,y)（2） 642yx ezplot(x.2+2*y.2-64,-8,8,-6,6)2,设 ,t,1tey在同一图形窗口采用子图的形式绘制条形图、阶梯图、杆图、和对数坐标图等不同图形，并对不同图形加标注说明 t=-pi:pi/100:pi; y=1./(1+exp(-t); subplot(2,2,1), stairs(t,y),title 阶梯图; t=-pi:pi/10:pi; y=1./(1+exp(-t); subplot(2,2,1), stairs(t,y),title 阶梯图; subplot(2,2,2), bar(t,y),title 条形图; subplot(2,2,3)
11、, stem(t,y),title 杆图; subplot(2,2,4), semilogx(y),grid on,title 对数坐标图3.绘制下列及坐标图（1） 4cos5 x=-2*pi:pi/20:pi; y=5.*cos(x)+4; polar(x,y)（2） 3,cos5in2 x=-pi/3:pi/30:pi/3; y=5.*(sin(x).2)./cos(x); polar(x,y)4.绘制下列三维图形。（1） 。20,sinco2ttzeyxtt t=0:pi/100:2*pi;x=exp(-t/20).*cos(t);y=exp(-t/20).*sin(t);z=t;plo
12、t3(x,y,z)（2）z=5,|x|5,|y|5。要求应用插值着色处理。 x,y=meshgrid(-10:0.1:10);z=x;z(:,:)=5;shading interp i=find(abs(x)5z1=z;z1(i)=NaN;surf(x,y,z1);shading interp;5.播放一个直径不断变化的球体。 n=30;x,y,z=sphere;m=moviein(n);for j=1:nsurf(j*x,j*y,j*z);m(:,j)=getframe; end movie(m,2)第 9 章1.建立一个图形窗口，使之背景颜色为红色，并在窗口上保留原有的菜单项，而且在按下鼠
13、标左键后显示出 Left Button Pressed 字样。 co=figure(Color,r,WindowButtonDownFcn,axis off;text(0.1,0.5,Left Button Pressed,fontsize,30)co =22.利用图形对象回执下列曲线，要求先利用默认属性绘制曲线，然后通过图形句柄操作来改变曲线的颜色、线性和线宽，并利用文本对象给曲线添加文字标注。(1) 21ln2xy x=0:0.1:20;y=log(x+sqrt(1+x.2)./2;fi=plot(x,y);set(fi,color,r,linestyle,:,linewidth,1)（2
14、） 25tyx t=0:0.1:20;x=t.2;y=5.*t.2;fi=plot3(t,x,y);set(fi,color,g,linestyle,-.,linewidth,2)4.利用图形对象绘制下列三维图形，要求对图形进行光照处理。（1） )sin(52xyxz x,y=meshgrid(0:0.1:10);z=x.*2+y.*2-5*sin(x.*y); surf(x,y,z);light(style,local)（2） 3zy x,y=meshgrid(0:0.1:10); z=y.3; surf(x,y,z);light(posi,0,0,1)4.以任意位置子图形绘制出正弦、余弦、
15、正切和余切函数曲线。 t=linspace(0,4*pi,50);x=sin(t);y=cos(t);z=tan(t);w=cot(t);axes(Position,0.2,0.6,0.3,0.3,GridLineStyle,-); plot(t,x);title(sint);axes(Position,0.6,0.6,0.3,0.3,GridLineStyle,-);plot(t,y);title(cost);axes(Position,0.2,0.2,0.3,0.3,GridLineStyle,-);plot(t,z);title(tant);axes(Position,0.6,0.2,0.3,0.3,GridLineStyle,-);plot(t,w);title(cott)5.用 patch 函数绘制一个填充渐变色的正五变形。 t=linspace(0,2*pi,6);k=t(1:length(t)-1);patch(sin(k),cos(k),rand(1,length(k),facecolor,interp);axis equal
