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 1

2 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]); %绘制第三幅子图