Matlab程序设计上机报告.doc

Matlab程序设计上机报告第2章1. 求下列表达式的值。（1） w=sqrt(2)*(1+0.34245*10(-6)w =1.4142（2），其中a=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）,其中=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=。 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.0027i 1.9884 - 0.3696i 1.7706 - 1.0539i2. 已知求下列表达式的值（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 -10 12 37 26 -15 73 7 A2-B+eye(3)ans = -18 -217 17 22 533 10921 867 526 A*Bans = 14 14 16 -10 51 21 125 328 180 A.*Bans = -8 15 4 0 35 24-9 122 0 A.*Bans = -8 15 4 0 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.3173 -0.2924 A,Bans = -1 5 -4 8 3 -1 0 7 8 2 5 3 3 61 7 -3 2 0 A(1,3,:);B2ans = -1 5 -4 3 61 7 73 37 1 17 37 13 -20 1 93. 已知完成下列操作：（1） 输出A在10，25范围内的全部元素。（2） 取出A前3行构成矩阵B，前两列构成矩阵C，右下角32只矩阵构成矩阵D，B与C的乘积构成矩阵E。（3） 分别求表达式E find(A=10&A B=A(1:3,:)B = 23.0000 10.0000 -0.7780 0 41.0000 -45.0000 65.0000 5.0000 32.0000 5.0000 0 32.0000 C=A(:,1:2)C = 23.0000 10.0000 41.0000 -45.0000 32.0000 5.00006.0000 -9.5400 D=A(2:4,3:4)D = 65.0000 5.0000 0 32.0000 54.0000 3.1400 E=B*CE = 1.0e+03 * 0.9141 -0.2239 1.2080 2.71231.1330 -0.2103 E E E&Dans = 1 1 0 1 1 1 E|Dans = 1 1 1 1 1 1 E|Dans = 0 0 1 0 0 04. 产生5 阶希尔伯特矩阵H和5 阶帕斯卡矩阵P，且求其行列式的值Hh和Hp以及它们的条件数Th和Tp，判断哪个矩阵性能更好，为什么？ H=hilb(5); P=pascal(5); Hh=det(H)Hh = 3.7493e-12 Hp=det(P)Hp = 1 Th=cond(H)Th = 4.7661e+05 Tp=cond(P)Tp = 8.5175e+03Tp比Th好，更接近于1 。5. 已知：求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.8725 0.3487 0.5501 0.4050D = -25.3169 0 0 0 -10.5182 0 0 0 16.8351设A为n阶矩阵，若存在常数及非零的n维向量x，使得Ax=x，则称是矩阵A的特征值，x是A属于特征值的特征向量第3章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 & m m=input(请输入一个百分制的成绩： );if m=90&m=80&m=70&m=60&m=0&mfor i=1:20s=strcat(输入第,num2str(i),个数字);a(i)=input(s);endb=a(1);c=a(1);for i=1:20if a(i)bb=a(i);endif a(i) max(a) min(a)5. ,当a取-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)End5. 当n分别取100、1000、10000时，求下列格式的值：要求分别用顺序结构和向量运算（使用sum函数）来实现。（1） 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） A=1,1,1;n=100,1000,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） x=0:pi/100:2*pi; y=x-x.3/factorial(3); plot(x,y)（2） ezplot(x.2+2*y.2-64,-8,8,-6,6)2,设在同一图形窗口采用子图的形式绘制条形图、阶梯图、杆图、和对数坐标图等不同图形，并对不同图形加标注说明 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), stem(t,y),title 杆图; subplot(2,2,4), semilogx(y),grid on,title 对数坐标图3. 绘制下列及坐标图（1） x=-2*pi:pi/20:pi; y=5.*cos(x)+4; polar(x,y)（2） x=-pi/3:pi/30:pi/3; y=5.*(sin(x).2)./cos(x); polar(x,y)4. 绘制下列三维图形。（1） t=0:pi/100:2*pi;x=exp(-t/20).*cos(t);y=exp(-t/20).*sin(t);z=t;plot3(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)5&abs(y)5);z1=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. 建立一个图形窗口，使之背景颜色为红色，并在窗口上保留原有的菜单项，而且在按下鼠标左键后显示出Left Button Pressed字样。 co=figure(Color,r,WindowButtonDownFcn,axis off;text(0.1,0.5,Left Button Pressed,fontsize,30)co = 22. 利用图形对象回执下列曲线，要求先利用默认属性绘制曲线，然后通过图形句柄操作来改变曲线的颜色、线性和线宽，并利用文本对象给曲线添加文字标注。(1) 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） 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） x,y=meshgrid(0:0.1:10);z=x.*2+y.*2-5*sin(x.*y); surf(x,y,z);light(style,local)（2） x,y=meshgrid(0:0.1:10); z=y.3; surf(x,y,z);light(posi,0,0,1)4. 以任意位置子图形绘制出正弦、余弦、正切和余切函数曲线。 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