数学实验答案

1.（1）[1 2 3 4;0 2 -1 1;1 -1 2 5;]+(1/2).*([2 1 4 10;0 -1 2 0;0 2 3 -2]) 2.A=[3 0 1;-1 2 1;3 4 2],B=[1 0 2;-1 1 1;2 1 1]X=(B+2*A)/23.A=[-4 -2 0 2 4;-3 -1 1 3 5]abs(A)>34.A=[-2 3 2 4;1 -2 3 2;3 2 3 4;0 4 -2 5]det(A),eig(A),rank(A),inv(A)求计算机高手用matlab解决。

>> A=[-2,3,2,4;1,-2,3,2;3,2,3,4;0,4,-2,5]求|A|>> abs(A)ans =2 3 2 41 2 3 23 2 3 40 4 2 5求r（A）>> rank(A)ans =4求A-1>> A-1ans =-3 2 1 30 -3 2 12 1 2 3-1 3 -3 4求特征值、特征向量>> [V,D]=eig(A) %返回矩阵A的特征值矩阵D 与特征向量矩阵VV =- ++ - - +- +- + - +D =+ 0 0 0 0 - 0 0 0 0 + 0 0 0 0 -将A的第2行与第3列联成一行赋给b>> b=[A(2,:),A(:,3)']b =1 -23 2 2 3 3 -21．a=round(unifrnd(1,100))i=7;while i>=0i=i-1;b=input('请输入一个介于0到100的数字： '); if b==adisp('You won!');break;else if b>adisp('High');else if b<adisp('Low');endendendend结果a =82请输入一个介于0到100的数字： 50Low请输入一个介于0到100的数字： 75Low请输入一个介于0到100的数字： 85请输入一个介于0到100的数字： 82You won!2．clear all；clc;n=input('请输入数字 n=');n1=floor(n/100); %取出百位数字n1n2=mod(floor(n/10),10); %取出十位数字n2n3=mod(n,10) ; %取出个位数字n3if n1^3+n2^3+n3^3==nfprintf('%d是“水仙花数”', n) % 注意输出格式前须有%符号elsefprintf('%d不是“水仙花”', n) % 注意输出格式前须有%符号end结果请输入数字 n=234234不是“水仙花数”>>3．price=input('请输入商品价格');switch fix(price/100)case {0,1} %价格小于200rate=0;case {2,3,4} %价格大于等于200但小于500rate=3/100;case num2cell(5:9) %价格大于等于500但小于1000 rate=5/100;case num2cell(10:24) %价格大于等于1000但小于2500 rate=8/100;case num2cell(25:49) %价格大于等于2500但小于5000 rate=10/100;otherwise %价格大于等于5000rate=14/100;endprice=price*(1-rate) %输出商品实际销售价格结果请输入商品价格250price =Function f=myfun(x)x=input;s=pi*x*xl=pi*x^24、Function y=circle(r)s=pi*x*xl=pi*x^24.syms rs=pi*r*rl=2*pi*r5. function fibonacci(n,m) f(1)=1;f(2)=1;for i=3:max(n,m)f(i)=f(i-1)+f(i-2); endfprintf('第%d项',m)x=f(m)fprintf('前%d项',n)s=f(1:n)COMMAND WINDOW输入：fibonacci(20,50)1.绘制])4,0[)(3sin(3π∈=x x e y x 的图像，要求用蓝色的星号画图；并且画出器官包络线3x e y ±=的图像，用红色的点划线画图。

2.用fplot 和ezplot 命令绘出函数)21sin(32t ey t +=-在区间]10,0[上的图像。

3.在同一图像窗口画三个子图要求使用指令gtext,axis,legend,title,xlabel,和ylabel ： （3）]8,1[,sin 1∈=x x e y x1.x=0:pi/25:4*pi;y1=exp(x/3).*sin(3*x);y2=exp(x/3);y3=-exp(x/3);plot(x,y1,'b*',x,y2,'r-.',x,y3,'r-.')2.t=1::10y=exp(-2*t/3).*sin(1+2*t);plot(t,y);figurefplot('exp(-2*t/3).*sin(1+2*t)',[1,10])ezplot('exp(-2*t/3).*sin(1+2*t)',[1,10])3.x=1:1/50:8;y=exp(1./x).*sin(x);subplot(1,3,3),plot(x,y,'b-');legend('y=exp(1/x)sinx'),grid on ,title('y=exp(1/x)sinx'),xlabel('x 轴'),ylabel('y 轴')gtext('真棒'),axis([x1 x8 y1 y10])4.x=0:pi/50:2*pi;y1=sin(x);y2=cos(x);y3=sin(2*x);plot(x,y1,'k*--',x,y2,'rs-',x,y3,'bo--'),grid ontitle('曲线y1=sin(x),y2=cos(x)与y3=sin(2*x)')xlabel('x 轴'),ylabel('y 轴')gtext('y1=sin(x)'),gtext('y2=cos(x)'),gtext('y3=sin(2*x)')legend('y1=sin(x)','y2=cos(x)','y3=sin(2*x)')5.绘制圆锥螺线的图像并加各种标注，圆锥螺线的参数方程为：⎪⎪⎪⎩⎪⎪⎪⎨⎧=≤≤==t z t t t t y t t t x 2)200(,6cos 6cos π 6.在同一图形窗口画半径为1的球面，柱面122=+y x 以及极坐标]2,0[,4sin 21πρ∈=t t7.用mesh 与surf 命令绘制三维曲面223y x z +=的图像，并使用不同的着色效果及光照效果。

8.绘制由函数14169222=++z y x 形成的立体图，并通过改变观测点获得该图形在各个坐标平面上的投影。

9.画三维曲面)2,2(522≤≤---=y x y x z 与平面3=z 的交线。

5.t=1:pi/50:20*pi;x=t.*cos(pi/6.*t);y=t.*sin(pi/6.*t);z=2*t;plot3(x,y,z);grid on ;title('圆锥螺线')xlabel('x 轴');ylabel('y 轴');zlabel('z 轴');axis square6. v=[-2 2 -2 2 -2 2];subplot(1,3,1),spheretitle('以半径为1的球面'),xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴'),axis(v)subplot(1,3,2),cylindertitle('柱面'),xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴')subplot(1,3,3),t=0:pi/100:2*pi;polar(t,(1/2)*sin(4*t))title('p=(1/2)*sin(4t)')7.[X,Y]=meshgrid(-8::8);Z=X.^2+3*Y.^2;subplot(1,2,1);mesh(X,Y,Z);shading interpsubplot(1,2,2);surf(X,Y,Z);shading flat8.[xx,yy,zz]=sphere(40);x=xx*3;y=yy*4;z=zz*2;surf(x,y,z)axis equal9.[X,Y]=meshgrid(-2::2);Z1=5-X.^2+Y.^2;subplot(1,3,1),mesh(X,Y,Z1),title('曲面');Z2=3*ones(size(X));% 创建一个和y矩阵相同大小的纯1矩阵subplot(1,3,2),mesh(X,Y,Z2),title('平面');r0=abs(Z1-Z2)<=1;ZZ=r0.*Z2;YY=r0.*Y;XX=r0.*X;subplot(1,3,3)plot3(XX(r0~=0),YY(r0~=0),ZZ(r0~=0),'*')title('交线')10.v=[-22 -22 -22];[x,y,z]=sphere(30);surf(4*x,4*y,4*z)title('半径为4的球面'),axis(v)。