第一章习题解答.pptVIP

* 第一章习题解答 一、实验一 1.实质是矩阵的拼接 E=eye(3) R=rand(3,2) O=zeros(2,3) S=[2,0,;0 3] A=[E,R;O,S] A=[E,R;O,S]; A^2 B=[E R+R*S;O,S^2] 2.求解线性齐次线性方程组的通解 A=[1 1 -1 1 1;4 2 -2 1 2;2 1 -1 -1 1] rref(A) ans = 1 0 0 0 0 0 1 -1 0 1 0 0 0 1 0 即线性方程组化简为同解方程组为 3.利润收入计算 A=[1 2 3 4 5 6 7 8 9;7.15 8.25 3.20 10.3 6.68 12.03 16.85 17.51 9.30;11.1 15 6 16.25 9.9 18.25 20.8 24.15 15.5;568 1205 753 580 395 2104 1538 810 694] B(1,:)=A(1,:); for i=1:2 B(i+1,:)=A(i+1,:).*A(4,:); %计算成本和收入,第2行为成本,第三行为收入 end B(4,:)=B(3,:)-B(2,:) %计算利润 for i=2:9 if B(4,1)B(4,i) B(:,[1,i])=B(:,[i,1]); %求出最小利润 end end min=B(:,1) for i=2:9 if B(4,1)B(4,i) B(:,[1,i])=B(:,[i,1]); %求出最大利润 end end max=B(:,1) for i=1:8 for j=i+1:9 if B(4,i)B(4,j) B(:,[j,i])=B(:,[i,j]); % 将所有商品及其收入按收入由小到大排序 end end end C=B % 用C记录排序后的矩阵 totalr=sum(B(3,:)) totalli=sum(B(4,:)) 实验二 1.用两种方法作图，并要求用两种方法加标注 法一）plot((x,y) s=-1.1:0.01:1.1; y1=x.^2;y2=x.^3;y3=x.^4;y4=x.^5; subplot(2,2,1);plot(x,y1);gtext(y=x^2) subplot(2,2,2);plot(x,y2);gtext(y=x^3) subplot(2,2,3);plot(x,y3);gtext(y=x^4) subplot(2,2,4);plot(x,y4);gtext(y=x^5) s=-1.1:0.01:1.1; y1=x.^2;y2=x.^3;y3=x.^4;y4=x.^5; plot(x,y1,y);hold on plot(x,y2,b);hold on plot(x,y3,r);hold on plot(x,y4);hold off 法二）用hold on命令 2。用subplot命令在不同的坐标系下作出四条曲线 x=-3:0.01:3; y1=exp(-x.^2/2); subplot(2,2,1);plot(x,y1); subplot(2,2,2);ezplot(‘3*t/(1+t^3)’,’3*t^3/(1+t^3’); y=-0.9:0.01:0.9;x=log(1+sqrt(1-y.^2)./y)-sqrt(1-y.^2); subplot(2,2,3);plot(x,y) Subplot(2,2,4);ezplot(sin(2*t)*cos(t),sin(2*t)*sin(t)) 3.作出曲面的3维图形 x=-10:0.5:10;y=x; [X,Y]=meshgrid(x,y);Z=sin(pi*sqrt(X.^2+Y.^2)); surf(X, Y, Z) *