matlab习题第三章.docx

matlab习题第三章.docx

《matlab习题第三章.docx》由会员分享，可在线阅读，更多相关《matlab习题第三章.docx（17页珍藏版）》请在冰豆网上搜索。

matlab习题第三章

第三章

实验指导

1、

m=input（‘请输入一个三位数：

’）；

m1=fix（m/100）;

m2=rem（fix（m/10）,10）;

m3=rem（m,10）;

m=m1+m2*10+m3*100;

disp（m）;

请输入一个三位数:

639

936

2、

m=input（'请输入一个百分制的成绩:

'）;

ifm>=0&m<=100

switchfix（m/10）

case{9,10}

disp（'这个成绩所属的成绩等级：

A'）;

case{8}

disp（'这个成绩所属的成绩等级：

B'）;

case{7}

disp（'这个成绩所属的成绩等级：

C'）;

case{6}

disp（'这个成绩所属的成绩等级：

D'）;

otherwise

disp（'这个成绩所属的成绩等级：

E'）;

end

else

disp（'错误：

输入的成绩不合理'）;

end

3、

a=[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20]

b=max（a）

c=min（a）

b=0

c=0

fori=1:

20

ifa（i）>b

b=a（i）;

end

ifa（i）

c=a（i）;

end

end

4、

a=[-3:

0.1:

3];

b=（（exp（0.3.*a）-exp（-0.3.*a））./2）

c=sin（a+0.3）

y=b.*c+log（（0.3+a）./2）

fori=1:

61

y=（（exp（0.3.*a（i））-exp（-0.3.*a（i）））./2）*sin（a（i）+0.3）+log（（0.3+a（i））/2）

end

b=

Columns1through10

-1.0265-0.9840-0.9423-0.9015-0.8615-0.8223-0.7838-0.7461-0.7090-0.6725

Columns11through20

-0.6367-0.6014-0.5666-0.5324-0.4986-0.4653-0.4325-0.4000-0.3678-0.3360

Columns21through30

-0.3045-0.2733-0.2423-0.2115-0.1810-0.1506-0.1203-0.0901-0.0600-0.0300

Columns31through40

00.03000.06000.09010.12030.15060.18100.21150.24230.2733

Columns41through50

0.30450.33600.36780.40000.43250.46530.49860.53240.56660.6014

Columns51through60

0.63670.67250.70900.74610.78380.82230.86150.90150.94230.9840

Column61

1.0265

c=

Columns1through10

-0.4274-0.5155-0.5985-0.6755-0.7457-0.8085-0.8632-0.9093-0.9463-0.9738

Columns11through20

-0.9917-0.9996-0.9975-0.9854-0.9636-0.9320-0.8912-0.8415-0.7833-0.7174

Columns21through30

-0.6442-0.5646-0.4794-0.3894-0.2955-0.1987-0.09980.00000.09980.1987

Columns31through40

0.29550.38940.47940.56460.64420.71740.78330.84150.89120.9320

Columns41through50

0.96360.98540.99750.99960.99170.97380.94630.90930.86320.8085

Columns51through60

0.74570.67550.59850.51550.42740.33500.23920.14110.0416-0.0584

Column61

-0.1577

y=

Columns1through5

0.7388+3.1416i0.7696+3.1416i0.7871+3.1416i0.7913+3.1416i0.7822+3.1416i

Columns6through10

0.7602+3.1416i0.7254+3.1416i0.6784+3.1416i0.6196+3.1416i0.5496+3.1416i

Columns11through15

0.4688+3.1416i0.3780+3.1416i0.2775+3.1416i0.1680+3.1416i0.0497+3.1416i

Columns16through20

-0.0771+3.1416i-0.2124+3.1416i-0.3566+3.1416i-0.5104+3.1416i-0.6752+3.1416i

Columns21through25

-0.8536+3.1416i-1.0497+3.1416i-1.2701+3.1416i-1.5271+3.1416i-1.8436+3.1416i

Columns26through30

-2.2727+3.1416i-2.9837+3.1416i-37.0245-3.0017-2.3085

Columns31through35

-1.8971-1.5978-1.3575-1.1531-0.9723

Columns36through40

-0.8083-0.6567-0.5151-0.3819-0.2561

Columns41through45

-0.1374-0.02550.07920.17660.2663

Columns46through50

0.34780.42060.48410.53790.5815

Columns51through55

0.61450.63660.64740.64700.6351

Columns56through60

0.61190.57770.53270.47740.4126

Column61

0.3388

y=

0.7388+3.1416i

y=

0.7696+3.1416i

y=

0.7871+3.1416i

y=

0.7913+3.1416i

y=

0.7822+3.1416i

y=

0.7602+3.1416i

y=

0.7254+3.1416i

y=

0.6784+3.1416i

y=

0.6196+3.1416i

y=

0.5496+3.1416i

y=

0.4688+3.1416i

y=

0.3780+3.1416i

y=

0.2775+3.1416i

y=

0.1680+3.1416i

y=

0.0497+3.1416i

y=

-0.0771+3.1416i

y=

-0.2124+3.1416i

y=

-0.3566+3.1416i

y=

-0.5104+3.1416i

y=

-0.6752+3.1416i

y=

-0.8536+3.1416i

y=

-1.0497+3.1416i

y=

-1.2701+3.1416i

y=

-1.5271+3.1416i

y=

-1.8436+3.1416i

y=

-2.2727+3.1416i

y=

-2.9837+3.1416i

y=

-37.0245

y=

-3.0017

y=

-2.3085

y=

-1.8971

y=

-1.5978

y=

-1.3575

y=

-1.1531

y=

-0.9723

y=

-0.8083

y=

-0.6567

y=

-0.5151

y=

-0.3819

y=

-0.2561

y=

-0.1374

y=

-0.0255

y=

0.0792

y=

0.1766

y=

0.2663

y=

0.3478

y=

0.4206

y=

0.4841

y=

0.5379

y=

0.5815

y=

0.6145

y=

0.6366

y=

0.6474

y=

0.6470

y=

0.6351

y=

0.6119

y=

0.5777

y=

0.5327

y=

0.4774

y=

0.4126

y=

0.3388

5、

（1）n=input（'输入一个数：

'）;

a=0;

fori=i:

n

a=a+1/（i^2）;

end

disp（a）

j=1:

n;

b=sum（1./（j.^2））;

disp（b）

输入一个数：

100

1.6350

1.6350

输入一个数：

1000

1.6439

1.6439

输入一个数：

10000

1.6439

1.6439

（2）

n=input（'输入一个数：

'）;

a=1;

fori=1:

n

a=（（2*i*i*2）/（（2*i-1）*（2*i+1）））*a;

end

disp（a）;

clear

i=1:

n;

b=prod（（2.*i）.*（2.*i）./（（2.*i-1）.*（2.*i+1）））

disp（b）

输入一个数：

100

1.5669

输入一个数：

1000

1.5704

输入一个数：

10000

1.5708

6、

a=rand（5,6）;

n=input（'请输入一个数：

'）;

ifn>5

n=5

a（n,:

）

disp（'输入超出范围'）;

else

a（n,:

）

end

请输入一个数：

5

ans=

0.89130.44470.17630.89360.13890.1988

请输入一个数：

6

n=

5

ans=

0.46600.67210.83180.18970.15090.5936

输入超出范围

7、

f1.m代码如下

functionm=f1（n）

m=n+10*log（n^2+5）;

f2.m代码如下:

functionm=f2（n）

m=0;

fori=1:

n

m=m+i*（i+1）;

end

在command界面中输入：

y1=f1（40）/（f1（30）+f1（20））

y1=

0.6390

>>y2=f2（40）/（f2（30）+f2（20））

y2=

1.7662

8、

函数文件factor1.m代码如下：

functionf=factor1（k）

ifk==1

f=1;

else

f=factor1（k-1）+k;

end

函数文件factor2.m代码如下:

functionf=factor2（k）

ifk==1

f=1;

else

f=factor2（k-1）^2+k^2;

end

函数文件factor3.m代码如下：

functionf=factor3（k）

ifk==1

f=1;

else

f=factor31/（k-1）+1/k;

end

思考练习

1、用MATLAB语言编写的程序称为M文件。

书P57下部。

2、

N=[1,2,3,4,5]

N=

12345

>>A=N.*2

A=

246810

>>B=N./4

B=

0.25000.50000.75001.00001.2500

>>C=1./N

C=

1.00000.50000.33330.25000.2000

>>D=1./（N./2）./（N./2）/4

D=

1.00000.25000.11110.06250.0400

3、M文件

B=0;

fori=1:

20

A（i）=fix（（99-10+1）*rand

（1）+10）;

B=A（i）+B;

end

B=B/20;

fori=1:

20

ifA（i）

A（i）

end

end

输出结果：

ans=

66

ans=

18

ans=

24

ans=

22

4、

（1）M文件

A=1;

n=input（'输入一个数'）;

fori=2:

100

A=（-1）^（i+1）*1/（2*n-1）+A;

end

A

输出结果：

输入一个数100

A=

0.9950

>>factor3

输入一个数1000

A=

0.9995

>>factor3

输入一个数10000

A=

0.9999

（2）M文件

A=0;

n=input（‘输入一个数'）;

fori=1:

1000

A=1/（（4）^i）+A;

end

A

输出结果：

输入一个数10000

A=

0.9999

>>factor3

输入一个数100

A=

0.3333

>>factor3

输入一个数1000

A=

0.3333

>>factor3

输入一个数10000

A=

0.3333

5、

M文件

function[m,n]=chengji（A,B）

m=A.*B;

n=A*B;

输出程序及结果：

>>a=[123;456;789];

>>b=[-101;1-10;011]

[m,n]=chengji（a,b）;

disp（'a.*b='）;

disp（m）;

disp（'a*b='）;

disp（n）;

b=

-101

1-10

011

a.*b=

-103

4-50

089

a*b=

114

1110

1116