MATLAB二维三维画图仿真数组运算
1. 数学 (1) 数组运算; x=[1 2 3 4]; y=[3 4 5 6]; z=x+y %数组x与数组y相加得到数组z z = 4 6 8 10 z=x-y %数组x与数组y相减得到数组z z = -2 -2 -2 -2 z=x.*y %数组x与数组y相乘得到数组z z = 3 8 15 24 z=x./y %数组x与数组y相除得到数组z z = 0.3333 0.5000 0.6000 0.6667 (2)矩阵及其运算 a=[1 2 3 4;2 3 4 5;3 4 5 6] a= 1 2 3 4 2 3 4 5 3 4 5 6 a1=a %矩阵的转置 a1 = 1 2 3 2 3 4 3 4 5 4 5 6 b=[1 2 3;4 5 6;7 8 9]; c=[12 34 44;2 3 5;6 7 8]; d=b+c %矩阵的相加 d = 13 36 47 6 8 11 13 15 17 e=b-c %矩阵的相减 e = -11 -32 -41 2 2 1 1 1 1 A=b*c %矩阵的相乘 A = 34 61 78 94 193 249 154 325 420 A=c*b %矩阵的相乘 A = 456 546 636 49 59 69 90 111 132 A=[1 12 7;3 8 5;4 3 6]; B=[4;9;13]; C=A\B %矩阵的左除 C = 2.6442 -0.1731 0.4904 B=flipud(A) %矩阵B由矩阵A上下翻转得到 B = 4 3 6 3 8 5 1 12 7 C=A/B %矩阵的右除 C = 0 0 1 0 1 0 1 0 0 D=A^2 %矩阵A的二次方 D = 65 129 109 47 115 91 37 90 79 (3)函数及其运算 x=0:pi/100:2*pi; y=sin(x); plot(x,y) %正弦函数 (4)数值分析 %数值积分 d=pi/100; t=0:d:3*pi; nt=length(t); y=exp(-0.5*t).*sin(t+pi/6); z=trapz(y)*d z = 0.9008 sc=cumsum(y)*d; scf=sc(nt) scf = 0.9086 (5)数据分析 统计计算 randn( seed ,0); a=randn(100,5); a_max=max(max(a)) a_max = 3.2069 a_end=median(a) a_end = 0.1489 0.0301 0.0443 -0.0189 0.1375 a_mean=mean(a) a_mean = 0.0727 0.0264 0.0351 0.0424 0.0752 a_prod=prod(a) a_prod = 1.0e-020 * -0.0000 0.1268 -0.0000 -0.0000 0.0007 2. 图形技术 (1) 平面绘图 x=0:0.4:4*pi; y1=exp(-0.1*x).*sin(x); y2=exp(-0.1*x).*sin(x+1); plot(x,y1,x,y2) %绘制两条曲线 x=(0:pi/45:2*pi) ; k=0.3:0.1:1; y=cos(x)*k; plot(y) %绘制多条不同色彩的曲线 x=0:0.01:20; y1=200*exp(-0.5*x).*sin(x); y2=0.8*exp(-0.5*x).*sin(10*x); plotyy(x,y1,x,y2); %绘制双纵坐标图 (2)立体绘图(三维) [X,Y]=meshgrid(-10:0.5:10); B=sqrt(X.^2+Y.^2)+eps; Z=cos(B)./B; mesh(X,Y,Z); %绘制三维曲面网线 [X,Y]=meshgrid(-2:0.2:2); Z=X.*exp(-X.^2-Y.^2); plot3(X,Y,Z); %绘制多条三维曲线 [X,Y]=meshgrid(-10:0.5:10); B=sqrt(X.^2+Y.^2)+eps; Z=sin(B)./B; surf(X,Y,Z); %绘制三维曲面 (3)绘制子图 x=pi*(0:1000)/1000; y1=sin(x); y2=sin(10*x); y12=y1.*y2; subplot(221);plot(x,y1);axis([0,pi,-1,1]);%分割并绘制第一幅子图 subplot(222);plot(x,y2);axis([0,pi,-1,1]);%分割并绘制第二幅子图 subplot( position ,[0.2,0.05,0.6,0.45]);%分割第三幅子图的位置 plot(x,y12);axis([0,pi,-1,1]); %绘制第三幅子图 (4)图形标志 x=0:2*pi/90:2*pi; y1=sin(2*x); y2=exp(-15*x).*sin(2*x); plot(x,y1, r: ,x,y2, m* ); legend( 第一条 , 第二条 ) 3. 程序设计 (1) 控制流