MATLAB第三次实习

MATLAB

上机实验报告

086092班

魏旭阳

20091000567

一、实验内容

实验六循环结构程序设计

（1）编写程序实现求两矩阵A*B,A.*B.

A*B

A=input('输入矩阵A:');

B=input('输入矩阵B:');

C=zeros(size(A,1),size(B,2));

for ii=1:size(A,1)

for jj=1:size(B,2)

for k=1:size(A,2)

C(ii,jj)=C(ii,jj)+A(ii,k)*B(k,jj);

end

end

end

C

输入矩阵A:[1 2;3 4]

输入矩阵B:[4 3;2 1]

C =8 5

20 13

A.*B

A=input('请输入矩阵A:');

B=input('请输入矩阵B:');

if size(A)==size(B)

[m,n]=size(A);

for i=1:m

for j=1:n

C(i,j)=A(i,j)*B(i,j);

end

end

disp('输出A.*B:')

C

else

disp('无效')

end

请输入矩阵A:[1 2;3 4]

请输入矩阵B:[4 3;2 1]

输出A.*B:

C =4 6 6 4

（2）编写程序实现求矩阵的转置。

A=input('请输入矩阵A:'); [m,n]=size(A); B=zeros(n,m) for i=1:m for j=1:n

B(i,j)=A(j,i) end end disp(B)

请输入矩阵A:[5 6;7 8] B =5 7 6 8

实验八 二维曲线的绘制 （1）绘制下列曲线：

①3

3

x

x y -

=

x=-10:1:10;y=x-(x.^3)/3;plot(x,y)

-10

-8-6-4-20246810

-400-300-200-1000100200300

400

②2

2

21x

e

y π

=

x=-1:0.001:1;y=1/2*pi*exp(1/2*x.^2);plot(x,y)

-1

-0.8-0.6-0.4-0.200.20.40.60.81

1.41.6

1.8

2

2.2

2.4

2.6

2.8

③64222=+y x ezplot('x^2+2*y^2-64',[-8,8])

-8

-6-4-2

02468

-8-6-4-202468x

y

x 2+2 y 2-64 = 0

④⎩⎨

⎧==t

r y t r x cos sin

t=0:2*pi/100:2*pi;x=sin(t);y=cos(t);plot(x,y)（r 视为常量）

-1

-0.8

-0.6

-0.4

-0.2

0.2

0.4

0.6

0.8

1

-1-0.8-0.6-0.4-0.200.20.40.60.81

（2）通过用plot 和fplot 函数绘制x

y 1sin

的曲线，并分析其区别。

x=-1:0.01:1;y=sin(1./x);plot(x,y)

-1

-0.8

-0.6

-0.4

-0.2

00.20.40.60.81

-1-0.8-0.6-0.4-0.200.20.40.60.81

fplot('sin(1/x)',[-1,1],1e-4)

-1

-0.8

-0.6

-0.4

-0.2

0.2

0.4

0.6

0.8

1

-1-0.8-0.6-0.4-0.200.20.40.60.8

1

实验九 三维曲线和曲面的绘制 （1）绘制下列三维曲线

1) ⎪⎩

⎪⎨⎧=≤≤==--t z t t e

y t e x t t π20,sin cos 20

/20/ t=0:2*pi/100:2*pi;x=exp(-t/20).*cos(t);y=exp(-t/20).*sin(t),z=t;plot3(x,y,z);title('Line in

3-D Space');

xlabel('X');ylabel('Y');zlabel('Z'); grid on;

-1

-0.5

0.5

1

-1

-0.50

0.5

1

02

4

6

8

X

Line in 3-D Space

Y

Z

2) ⎪⎩⎪

⎨⎧=≤≤==3210,t

z t t y t x

t=0:0.001:1;x=t;y=t.^2;z=t.^3;plot3(x,y,z);title('Line in 3-D Space'); xlabel('X');ylabel('Y');zlabel('Z'); grid on;

0.2

0.40.6

0.8

1

0.5

1

00.20.40.60.81X

Line in 3-D Space

Y

Z

（2）绘制下列曲面图

1) z=5,|x|≤5,|y|≤5

[x,y]=meshgrid(-5:1:5);z=0*x+0*y+5;mesh(x,y,z)

-5

5

-5

5

44.5

5

5.5

6

2)3,3,15),(2

2

≤≤++-

=y x y

x y x f

[x,y]=meshgrid(-3:1:3);z=-5./(1+x.^2+y.^2);mesh(x,y,z)