北京交通大学matlab上机实验2.docx
《北京交通大学matlab上机实验2.docx》由会员分享,可在线阅读,更多相关《北京交通大学matlab上机实验2.docx(44页珍藏版)》请在冰豆网上搜索。
北京交通大学matlab上机实验2
一,建立向量。
观察结果。
1,>>x=1:
5,x=(1:
5)'
x=
12345
x=
1
2
3
4
5
2,>>y=0:
pi/4:
pi
y=
00.78541.57082.35623.1416
3,>>x=(0:
0.2:
3)',y=exp(-x).*sin(x)
x=
0
0.2000
0.4000
0.6000
0.8000
1.0000
1.2000
1.4000
1.6000
1.8000
2.0000
2.2000
2.4000
2.6000
2.8000
3.0000
y=
0
0.1627
0.2610
0.3099
0.3223
0.3096
0.2807
0.2430
0.2018
0.1610
0.1231
0.0896
0.0613
0.0383
0.0204
0.0070
4,>>k=linspace(-pi,pi,5),k=logspace(-3,1,5)
k=
-3.1416-1.570801.57083.1416
k=
0.00100.01000.10001.000010.0000
二,已知x=[123],y=[456],试计算z=x.*y,x.\y和x./y。
>>x=[123],y=[456]
x=
123
y=
456
>>z=x.*y
z=
41018
>>z=x.\y
z=
4.00002.50002.0000
>>z=x./y
z=
0.25000.40000.5000
三,解线性方程。
>>A=[5,7,6,5,1;7,10,8,7,2;6,8,10,9,3;5,7,9,10,4;1,2,3,4,5]
A=
57651
710872
681093
579104
12345
>>B=[24,96;34,136;36,144;35,140;15,60]
B=
2496
34136
36144
35140
1560
>>X=A\B
X=
1.00004.0000
1.00004.0000
1.00004.0000
1.00004.0000
1.00004.0000
四,求顶点是A(2,5,6),B(11,3,8),C(5,1,11)的三角形各边长。
>>A=[256],B=[1138],C=[5,1,11]
A=
256
B=
1138
C=
5111
>>a=sqrt(dot(A-C,A-C))
a=
7.0711
>>b=sqrt(dot(A-B,A-B))
b=
9.4340
>>c=sqrt(dot(C-B,C-B))
c=
7
五,进行如下逻辑运算,观察结果
>>P=[100]
P=
100
>>~P
ans=
011
>>P|(~P)
ans=
111
>>P&(~P)
ans=
000
>>C=rem(P,2)
C=
100
>>C&P
ans=
100
>>C|P
ans=
100
>>(C-1)&P
ans=
000
>>any(P)
ans=
1
>>all(P)
ans=
0
>>all(P|(~P))
ans=
1
>>any(P&(~P))
ans=
0
六,进行如下关系运算,观测结果。
>>y=[4215306]
y=
4215306
>>i=find(y>3.0)
i=
147
>>t=1/0
t=
Inf
>>t=NaN
t=
NaN
>>isnan(t)
ans=
1
七,实现下面分段函数
x,|x|<=1
y=f(x)={1,x>1
-1,x<-1
>>x=input('x=');
ifx>1
y=1
elseifx<-1
y=-1
else
y=x
end
八,分别用for和while循环语句编写程序,求出
>>K=0
fori=0:
63
K=K+2.^i
end
K=
0
K=
1
K=
3
K=
7
K=
15
K=
31
K=
63
K=
127
K=
255
K=
511
K=
1023
K=
2047
K=
4095
K=
8191
K=
16383
K=
32767
K=
65535
K=
131071
K=
262143
K=
524287
K=
1048575
K=
2097151
K=
4194303
K=
8388607
K=
16777215
K=
33554431
K=
67108863
K=
134217727
K=
268435455
K=
536870911
K=
1.0737e+009
K=
2.1475e+009
K=
4.2950e+009
K=
8.5899e+009
K=
1.7180e+010
K=
3.4360e+010
K=
6.8719e+010
K=
1.3744e+011
K=
2.7488e+011
K=
5.4976e+011
K=
1.0995e+012
K=
2.1990e+012
K=
4.3980e+012
K=
8.7961e+012
K=
1.7592e+013
K=
3.5184e+013
K=
7.0369e+013
K=
1.4074e+014
K=
2.8147e+014
K=
5.6295e+014
K=
1.1259e+015
K=
2.2518e+015
K=
4.5036e+015
K=
9.0072e+015
K=
1.8014e+016
K=
3.6029e+016
K=
7.2058e+016
K=
1.4412e+017
K=
2.8823e+017
K=
5.7646e+017
K=
1.1529e+018
K=
2.3058e+018
K=
4.6117e+018
K=
9.2234e+018
K=
1.8447e+019
>>K=0
i=0
whilei<63
i=i+1
K=K+2.^i
end
K=
0
i=
0
i=
1
K=
2
i=
2
K=
6
i=
3
K=
14
i=
4
K=
30
i=
5
K=
62
i=
6
K=
126
i=
7
K=
254
i=
8
K=
510
i=
9
K=
1022
i=
10
K=
2046
i=
11
K=
4094
i=
12
K=
8190
i=
13
K=
16382
i=
14
K=
32766
i=
15
K=
65534
i=
16
K=
131070
i=
17
K=
262142
i=
18
K=
524286
i=
19
K=
1048574
i=
20
K=
2097150
i=
21
K=
4194302
i=
22
K=
8388606
i=
23
K=
16777214
i=
24
K=
33554430
i=
25
K=
67108862
i=
26
K=
134217726
i=
27
K=
268435454
i=
28
K=
536870910
i=
29
K=
1.0737e+009
i=
30
K=
2.1475e+009
i=
31
K=
4.2950e+009
i=
32
K=
8.5899e+009
i=
33
K=
1.7180e+010
i=
34
K=
3.4360e+010
i=
35
K=
6.8719e+010
i=
36
K=
1.3744e+011
i=
37
K=
2.7488e+011
i=
38
K=
5.4976e+011
i=
39
K=
1.0995e+012
i=
40
K=
2.1990e+012
i=
41
K=
4.3980e+012
i=
42
K=
8.7961e+012
i=
43
K=
1.7592e+013
i=
44
K=
3.5184e+013
i=
45
K=
7.0369e+013
i=
46
K=
1.4074e+014
i=
47
K=
2.8147e+014
i=
48
K=
5.6295e+014
i=
49
K=
1.1259e+015
i=
50
K=
2.2518e+015
i=
51
K=
4.5036e+015
i=
52
K=
9.0072e+015
i=
53
K=
1.8014e+016
i=
54
K=
3.6029e+016
i=
55
K=
7.2058e+016
i=
56
K=
1.4412e+017
i=
57
K=
2.8823e+017
i=
58
K=
5.7646e+017
i=
59
K=
1.1529e+018
i=
60
K=
2.3058e+018
i=
61
K=
4.6117e+018
i=
62
K=
9.2234e+018
i=
63
K=
1.8447e+019
>>K=1.*(1-2.^64)/(1-2)
K=
1.8447e+019
九,用对分发解超越方程
>>a=0
b=2
f=@(x)log(x)-cos(x)
Tol=1e-12
while(b-a)>Tol
c=(a+b)/2
iff(c)>0
b=c
elseiff(c)<0
a=c
else
disp(c)
break
end
end
disp((a+b)/2)
a=
0
b=
2
f=
@(x)log(x)-cos(x)
Tol=
1.0000e-012
c=
1
a=
1
c=
1.5000
b=
1.5000
c=
1.2500
a=
1.2500
c=
1.3750
b=
1.3750
c=
1.3125
b=
1.3125
c=
1.2813
a=
1.2813
c=
1.2969
a=
1.2969
c=
1.3047
b=
1.3047
c=
1.3008
a=
1.3008
c=
1.3027
a=
1.3027
c=
1.3037
b=
1.3037
c=
1.3032
b=
1.3032
c=
1.3030
b=
1.3030
c=
1.3029
a=
1.3029
c=
1.3029
a=
1.3029
c=
1.3029
a=
1.3029
c=
1.3030
a=
1.3030
c=
1.3030
b=
1.3030
c=
1.3030
b=
1.3030
c=
1.3030
b=
1.3030
c=
1.3030
b=
1.3030
c=
1.3030
a=
1.3030
c=
1.3030
a=
1.3030
c=
1.3030
b=
1.3030
c=
1.3030
b=
1.3030
c=
1.3030
b=
1.3030
c=
1.3030
a=
1.3030
c=
1.3030
a=
1.3030
c=
1.3030
a=
1.3030
c=
1.3030
a=
1.3030
c=
1.3030
a=
1.3030
c=
1.3030
b=
1.3030
c=
1.3030
a=
1.3030
c=
1.3030
b=
1.3030
c=
1.3030
b=
1.3030
c=
1.3030
b=
1.3030
c=
1.3030
a=
1.3030
c=
1.3030
b=
1.3030
c=
1.3030
a=
1.3030
c=
1.3030
a=
1.3030
c=
1.3030
b=
1.3030
1.3030
>>a=0
b=2
f=@(x)log(x)-cos(x+pi/8)
Tol=1e-12
while(b-a)>Tol
c=(a+b)/2
iff(c)>0
b=c
elseiff(c)<0
a=c
else
disp(c)
break
end
end
disp((a+b)/2)
a=
0
b=
2
f=
@(x)log(x)-cos(x+pi/8)
Tol=
1.0000e-012
c=
1
a=
1
c=
1.5000
b=
1.5000
c=
1.2500
b=
1.2500
c=
1.1250
b=
1.1250
c=
1.0625
a=
1.0625
c=
1.0938
b=
1.0938
c=
1.0781
a=
1.0781
c=
1.0859
a=
1.0859
c=
1.0898
a=
1.0898
c=
1.0918
b=
1.0918
c=
1.0908
a=
1.0908
c=
1.0913
b=
1.0913
c=
1.0911
b=
1.0911
c=
1.0909
a=
1.0909
c=
1.0910
b=
1.0910
c=
1.0910
b=
1.0910
c=
1.0910
b=
1.0910
c=
1.0910
b=
1.0910
c=
1.0909
b=
1.0909
c=
1.0909
b=
1.0909
c=
1.0909
a=
1.0909
c=
1.0909
b=
1.0909
c=
1.0909
a=
1.0909
c=
1.0909
a=
1.0909
c=
1.0909
a=
1.0909
c=
1.0909
b=
1.0909
c=
1.0909
a=
1.0909
c=
1.0909
a=
1.0909
c=
1.0909
a=
1.0909
c=
1.0909
b=
1.0909
c=
1.0909
b=
1.0909
c=
1.0909
b=
1.0909
c=
1.0909
b=
1.0909
c=
1.0909
a=
1.0909
c=
1.0909
b=
1.0909
c=
1.0909
b=
1.0909
c=
1.0909
b=
1.0909
c=
1.0909
a=
1.0909
c=
1.0909
a=
1.0909
c=
1.0909
a=
1.0909
c=
1.0909
a=
1.0909
1.0909
>>a=0
b=2
f=@(x)log(x)-cos(x+pi/4)
Tol=1e-12
while(b-a)>Tol
c=(a+b)/2
iff(c)>0
b=c
elseiff(c)<0
a=c
else
disp(c)
break
end
end
disp((a+b)/2)
a=
0
b=
2
f=
@(x)log(x)-cos(x+pi/4)
Tol=
1.0000e-012
c=
1
b=
1
c=
0.5000
a=
0.5000
c=
0.7500
a=
0.7500
c=
0.8750
a=
0.8750
c=
0.9375
b=
0.9375
c=
0.9063
b=
0.9063
c=
0.8906
a=
0.8906
c=
0.8984
b=
0.8984
c=
0.8945
a=
0.8945
c=
0.8965
b=
0.8965
c=
0.8955
a=
0.8955
c=
0.8960
b=
0.8960
c=
0.8958
b=
0.8958
c=
0.8956
a=
0.8956
c=
0.8957
a=
0.8957
c=
0.8957
a=
0.8957
c=
0.8957
b=
0.8957
c=
0.8957
a=
0.8957
c=
0.8957
a=
0.8957
c=
0.8957
b=
0.8957
c=
0.8957
a=
0.8957
c=
0.8957
b=
0.8957
c=
0.8957
b=
0.8957
c=
0.8957
b=
0.8957
c=
0.8957
b=
0.8957
c=
0.8957
b=
0.8957
c=
0.8957
a=
0.8957
c=
0.8957
b=
0.8957
c=
0.8957
b=
0.8957
c=
0.8957
b=
0.8957
c=
0.8957
a=
0.8957
c=
0.8957
a=
0.8957
c=
0.8957
a=
0.8957
c=
0.8957
a=
0.8957
c=
0.8957
a=
0.8957
c=
0.8957
a=
0.8957
c=
0.8957
a=
0.8957
c=
0.8957
a=
0.8957
c=
0.8957
b=
0.8957
c=
0.8957
a=
0.8957
c=
0.8957
b=
0.8957
0.8957
>>a=0
b=2
f=@(x)log(x)-cos(x+3.*pi/8)
Tol=1e-12
while(b-a)>Tol
c=(a+b)/2
iff(c)>0
b=c
elseiff(c)<0
a=c
else
disp(c)
break
end
end
disp((a+b)/2)
a=
0
b=
2
f=
@(x)log(x)-cos(x+3.*pi/8)
Tol=
1.0000e-012
c=
1
b=
1
c=
0.5000
a=
0.5000
c=
0.7500
b=
0.7500
c=
0.6250
a=
0.6250
c=
0.6875
a=
0.6875
c=
0.7188
a=
0.7188
c=
0.7344
b=
0.7344
c=
0.7266
b=
0.7266
c=
0.7227
a=
0.7227
c=
0.7246
b=
0.7246
c=
0.7236
b=
0.7236
c=
0.7231
b=
0.7231
c=
0.7229
a=
0.7229
c=
0.7230
b=
0.7230
c=
0.7230
a=
0.7230
c=
0.7230
a=
0.7230
c=
0.7230
a=
0.7230
c=
0.7230
b=
0.7230
c=
0.7230
b=
0.7230
c=
0.7230
b=
0.7230
c=
0.7230
b=
0.7230
c=
0.7230
a=
0.7230
c=
0.7230
b=
0.7230
c=
0.7230
a=
0.7230
c=
0.7230
a=
0.7230
c=
0.7230
b=
0.7230
c=
0.7230
a=
0.7230
c=
0.7230
a=
0.7230
c=
0.7230
b=
0.7230