ch1(实验--Matlab的基本操作编程).docx

1设有分块矩阵，其中E,R,O,S分别为单位阵、随机阵、零阵和对角阵，试通过数值计算验证。 程序代码：clear;E=eye(3);R=rand(3,2);O=zeros(2,3);X=3,4;S=diag(X);A=E,R;O,S;% EQ1,EQ2分别表示待验证的等式左边和右边EQ1=A2;EQ2=E,R+R*S;O,S2;EROSAEQ1EQ2if EQ1=EQ2fprintf(左边=右边，该等式成立n)elsefprintf(左边不等于右边，该等式不成立n)end运行结果：E = 1 0 0 0 1 0 0 0 1R = 0.9218 0.4057 0.7382 0.9355 0.1763 0.9169O = 0 0 0 0 0 0S = 3 0 0 4A = 1.0000 0 0 0.9218 0.4057 0 1.0000 0 0.7382 0.9355 0 0 1.0000 0.1763 0.9169 0 0 0 3.0000 0 0 0 0 0 4.0000EQ1 = 1.0000 0 0 3.6873 2.0285 0 1.0000 0 2.9528 4.6773 0 0 1.0000 0.7051 4.5845 0 0 0 9.0000 0 0 0 0 0 16.0000EQ2 = 1.0000 0 0 3.6873 2.0285 0 1.0000 0 2.9528 4.6773 0 0 1.0000 0.7051 4.5845 0 0 0 9.0000 0 0 0 0 0 16.0000左边=右边，该等式成立2.某零售店有9种商品的单件进价（元）、售价（元）及一周的销量如下表，问哪种商品的利润最大，哪种商品的利润最小；按收入由小到大，列出所有商品及其收入；求这一周该9种商品的总收入和总利润。 货号123456789单件进价7.158.253.2010.306.6812.0316.8517.519.30单件售价11.1015.006.0016.259.9018.2520.8024.1515.50销量568120575358039521041538810694程序代码：clear;inPrice=7.15,8.25,3.20,10.30,6.68,12.03,16.85,17.51,9.30;%进价outPrice=11.10,15.00,6.00,16.25,9.90,18.25,20.80,24.15,15.50;%售价sold=568,1205,753,580,395,2104,1538,810,694;%销量inCome=outPrice.*sold;%收入profit=(outPrice-inPrice).*sold;%利润max1,pos1=max(profit);min1,pos2=min(profit);fprintf(n1.求利润最大、最小的商品n)fprintf(货号为%d的商品利润最大，为%f元n,pos1,max1)fprintf(货号为%d的商品利润最小，为%f元n,pos2,min1)fprintf(n2.按收入由小到大，列出所有商品及其收入n)A,pos=sort(inCome);i=1;while i=10e-5)x1=(x2+a/x2)/2;temp=x1;x1=x2;x2=temp;endx=x2;运行结果：A=sq(2),sq(3);sq(4),sq(5);sq(6),sq(7);sq(8),sq(9)A = 1.41421356237469 1.73205081001473 2 2.23606797749998 2.44948974278755 2.64575131111137 2.82842712504986 3.000000001396985 求函数的极限、导数或积分：程序代码：syms x y z n;fprintf(n1)求lim(x+3x)(1/x),当x-时)F1=(x+3x)(1/x);f1=limit(F1,x,inf)fprintf(n2)求lim(exp(x)*sin(x)-x*(x+1)/(x3),当x-0时)F2=(exp(x)*sin(x)-x*(x+1)/(x3);f2=limit(F2,x,0)fprintf(n3)f(x)=(x2+2*x-1)/(exp(-x)*sin(x)+1),求f(1)(x)F3=(x2+2*x-1)/(exp(-x)*sin(x)+1);f3=diff(F3,x)fprintf(n4)f(x)=x2/(1-x2),求f(n)(0)n)fprintf(nn=20,为本组成员学号后两位(04 17 39)的平均值)n=(4+17+39)/3;F4=x2/(1-x2);f4=subs(diff(F4,x,n),x,0)fprintf(n5)atan(y/x)=log(sqrt(x2+y2),求dx/dy)F5=atan(y/x)-log(sqrt(x2+y2);df_dx=diff(F5,x);df_dy=diff(F5,y);dx_dy=-(df_dy/df_dx)fprintf(n6)z=x*atan(y),求偏导dz/dx,偏导dz/dy)F6=x*atan(y);dzdx=diff(F6,x)dzdy=diff(F6,y)fprintf(n7)exp(2*x)/(exp(x)+2)dx,求此积分)F7=exp(2*x)/(exp(x)+2);f7=int(F7,x)运行结果：1)求lim(x+3x)(1/x),当x-时 f1 =32)求lim(exp(x)*sin(x)-x*(x+1)/(x3),当x-0时 f2 = 1/33)f(x)=(x2+2*x-1)/(exp(-x)*sin(x)+1),求f(1)(x) f3 =(2*x+2)/(exp(-x)*sin(x)+1)-(x2+2*x-1)/(exp(-x)*sin(x)+1)2*(-exp(-x)*sin(x)+exp(-x)*cos(x)4)f(x)=x2/(1-x2),求f(n)(0)n=20,为本组成员学号后两位(04 17 39)的平均值f4 = 2.4329e+0185)atan(y/x)=log(sqrt(x2+y2),求dx/dy dx_dy =-(1/x/(1+y2/x2)-1/(x2+y2)*y)/(-y