matlab复习1

•7• 结束
显示： E=1 0 0 O=0 0 0 1 0 0 0 0 0 1 0 0 R = 0.9501 0.4860 0.2311 0.8913 0.6068 0.7621 S = 0.4565 0 0.9501 0.2311 0.6068 0.4565 0 0 0.0185 0.4860 0.8913 0.7621 0 0.0185
的3维图形。 维图形。
•20•
结束
（1）程序代码如下： clear x=[-5:.1:5]; y=[-5:.1:5]; [X,Y]=meshgrid(x,y); Z=sin(pi.*sqrt(X.^2+Y.^2)); mesh(X,Y,Z)
（2）程序代码如下： clear u=[-.1:.1:2*pi]; v=[-.1:.1:2*pi]; [U,V]=meshgrid(u,v); x=(1+cos(U)).*cos(V); y=(1+cos(U)).*sin(V); z=sin(U); mesh(x,y,z)
•21•
结束
（1）显示：
•22•
结束
（2）显示：
•23•
结束
4．自由发挥：自己提出问题，实验探索，广泛 自由发挥：自己提出问题，实验探索， 联想，发现规律，大胆猜想。比如函数sin(1/x) 联想，发现规律，大胆猜想。比如函数sin(1/x) x=0附近的振荡现象 有无规律可寻？ 附近的振荡现象， 在x=0附近的振荡现象，有无规律可寻？ 程序代码如下： clear x=[-100:.01:100]; y=x.*sin(x)-1; plot(x,y) grid on f=solve('x*sin(x)=1') n1=numeric(f) 显示： f=
•12• 结束
第一章习题（2）
1. 1 在同一个坐标下作出 y1=1+x, y2=1+x+x2/2, y3=1+x+x2/2!+x3/3!, y4=ex 这四条曲线的图形,观察、发现、联想、 猜想，给出验证及理论证明。 。
•13• 结束
程序代码如下： clear x=-3:0.01:3; y0=0*x; % x轴 y1=1+x; y2=y1+x.^2/2; y3=y2+x.^3/6; y4=exp(x); plot(x,y0,'g') hold on;plot(x,y1,'r') plot(x,y2,‘m’) % 紫色 plot(x,y3,'b') plot(x,y4,'k') grid on, hold off
•5• 结束
3．设有分块矩阵
E 3×3 A= O 2×3 R 3×2 S 2×2
其中E,R,O,S分别为单位阵、随机阵、 其中E,R,O,S分别为单位阵、随机阵、零阵 E,R,O,S分别为单位阵 和对角阵，试编写一个命令M文件， 和对角阵，试编写一个命令M文件，计算验 证
E R + RS A = 0 S2
x = ln 1± 1− y2 m y 1− y2
结束
•16•
程序代码如下： clear x=[-4:.01:4]; subplot(2,2,1);y1=exp(-x.^2);plot(x,y1) subplot(2,2,2);y2=sin(2*x);polar(x,y2) subplot(2,2,3);t=[-0.5:0.1:10 ,-10:0.1:-1.5]; x3=3*t./(1+t.^3);y3=3*t.^2./(1+t.^3); plot(x3,y3) subplot(2,2,4);y4=[0.01:.001:1]; %定义域(0,1] x41=log((1+sqrt(1-y4.^2))./y4)-sqrt(1-y4.^2); plot(y4,x41);hold on x42=log((1-sqrt(1-y4.^2))./y4)+sqrt(1-y4.^2); plot(y4,x42), hold off
•26•
s= 153 s= 370 s= 371 s= 407 n= 4
结束
5．建立一个命令M-文件：求所有的“水仙花数” 建立一个命令M 文件：求所有的“水仙花数” 所谓“水仙花数”是指一个三位数， ，所谓“水仙花数”是指一个三位数，其各位数字 的立方和等于该数本身。例如，153是一个水仙花 的立方和等于该数本身。例如，153是一个水仙花 因为153=1 数，因为153=13+53+33。 显示： 程序代 码如下： clear;n=0; for i=1:9;for j=0:9;for k=0:9 a=k+j*10+i*100; b=i^3+j^3+k^3; if a==b n=n+1;s=a end end;end;end n
•14• 结束
显示：
•15•
结束
2．用subplot分别在不同的坐标系下作出 subplot分别在不同的坐标系下作出 四条曲线： 四条曲线： 1）概率曲线 =sin2 2）四叶玫瑰线 ρ=sin2θ； 直角坐标方程， 3t x = 1 + t3 极坐标方程， 3）叶形线 2 参数方程 y = 3t 3 1+ t 各领风骚 4）曳物线
-1.1141571408719300873005251781692
n1 = -1.1142
•24•
结束
显示：
•25•
结束
5．建立一个命令M-文件：求所有的“水仙花数” 建立一个命令M 文件：求所有的“水仙花数” 所谓“水仙花数”是指一个三位数， ，所谓“水仙花数”是指一个三位数，其各位数字 的立方和等于该数本身。例如，153是一个水仙花 的立方和等于该数本身。例如，153是一个水仙花 因为153=1 数，因为153=13+53+33。 显示： 程序代 码如下： clear;n=0; for m=100:999 i=mod(m,10); j=mod((m-i)/10,10); k=floor(m/100); b=i^3+j^3+k^3; if m==b n=n+1;s=m end,end n
•11• 结束
显示： MinLi = 1.2719e+003 inmin = 5 MaxLi = 1.3087e+004 inmax = 6 D2 = 1.0e+004 * Columns 1 through 6 0.3911 0.4518 0.6305 0.9425 1.0757 1.8075 Columns 7 through 9 1.9562 3.1990 3.8398 inD2 = 5 3 1 4 9 2 8 7 6 ZongShouRu =1.4294e+005 ZongLiRu结束
A =1.0000 0 0 0 1.0000 0 0 0 1.0000 0 0 0 0 0 0
•8•
X =1.0000 0 0 1.3838 0.4950 0 1.0000 0 0.3366 0.9078 0 0 1.0000 0.8838 0.7762 0 0 0 0.2084 0 0 0 0 0 0.0003 Y = 1.0000 0 0 1.3838 0.4950 0 1.0000 0 0.3366 0.9078 0 0 1.0000 0.8838 0.7762 0 0 0 0.2084 0 0 0 0 0 0.0003 ans = 成立
2
•6•
结束
程序代码如下： clear E=eye(3); %三阶单位矩阵 R=rand(3,2); % 3行2列随机矩阵 O=zeros(2,3); % 2行3列零矩阵 v=rand(1,2);S=diag(v); % 以随机数为对角元素的对角阵 B1=[E,R];B2=[O,S];A=[B1;B2] C1=R+R*S;C2=S^2;X=A^2 % 等式左边 D1=[E,C1];D2=[O,C2];Y=[D1;D2] % 等式右边 if X= =Y '成立' else '不成立' end
结束
程序代码如下： clear A=[2,1,-1,1;4,2,-2,1;2,1,-1,-1];b=[1;2;1]; %求无穷个解中零最多的解 xz=A\b xn=inv(A)*b %求无穷个解中最小范数解 r1=rank(A) AA=[A,b] %增广矩阵 r2=rank(AA) if r1==r2 eA=rref(AA) %初等行变换 else c='无解' end
编写M 文件求最大利润、 编写M-文件求最大利润、最小利润和相应的商 品号；按收入由小到大， 品号；按收入由小到大，列出所有商品及其收 求这一周该10种商品的总收入和总利润。 10种商品的总收入和总利润 入；求这一周该10种商品的总收入和总利润。
•10• 结束
程序代码如下： clear A=[7.15,8.25,3.20,10.30,6.68,12.03,16.85,17.51,9.30]; B=[11.10,15.0,6.00,16.25,9.90,18.25,20.80,24.15,15.5]; C=[568,1205,753,580,395,2104,1538,810,694]; P=(B-A).*C; [MinLi, inmin]=min(P) %最小利润及其标号 %最小利润及其标号 [MaxLi, inmax]=max(P) %最大利润及其标号 最大利润及其标号 D1=B.*C; %各商品收入 收入 [D2，inD2]=sort(D1); %收入升序排列及其标号 收入升序排列及其标号 ZongShouRu=sum(D1) %总收入 总 ZongLiRun=sum(P) %总利润 总
•17• 结束
显示：
•18•
结束
显示：第4幅图有误，未考虑定义域