matlab习题结果解析

习题与解答

第1 章P30

1. Matlab 菜单栏的使用

1)利用菜单栏使Matlab 桌面只显示命令窗口和菜单栏；再恢复默认设置。

2)在我的文档文件夹中建立一个以自己学号为名的子文件夹，并把它添加至Matlab 搜索路

径的尾端。

解：

1)点击Desktop，在下拉菜单中将除Command Window 前的对号点击去除，则只显示命

令窗口和菜单栏；点击Desktop，Desktop Layout，Default 则恢复默认设置。

2) 依次点击File，Set Path，Add Folder，选中“我的文档”文件夹，点击新文件夹，在

窗口下部输入文件夹名，点击确定，最后点击“Move to bottom”，并保存结果。

2. Matlab 查找与帮助命令的使用

1) 利用help 命令学习使用ginput 函数的作用以及其用法，利用ginput 函数以50 个点描

绘你的左手形状，并将这些点的横纵坐标分别赋值于hand_x 和hand_y 变量中；

2) 将hand_x 和hand_y 两个变量保存到文件Hand.mat 中；

3) 在命令窗口输入clear 命令清除变量空间的变量，再使用load 命令，将Hand.mat 文件

中的命令调入变量空间中，采用plot 命令绘图。

解：

1) 采用help ginput 命令，可知ginput 函数的使用；运行[hand_x, hand_y]=ginput(50)命

令，待出现图形窗口后点击鼠标获取坐标值；

2) save hand hand_x hand_y

3) 命令依次为clear，load hand，plot(hand_x,hand_y)

3.标准条件下钢管中流动的空气压降可由下式给出：

1.84

1.24

0.03

1000

p L V

d

Δ = ?? ??

? ?

式中，L 为管长(单位：m)，V 为空气流速(单位：m/min)，d 为钢管直径(单位mm)。

求L＝3000m，d＝45mm，V＝1600m/min 时的Δp 值。计算结果按如下格式输出：

L=3000m d=45mm V=1600m/min

压降计算值为：

deltP= …

(注：…处为计算结果，保留6 位小数)

解：

L=3000;d=45;V=1600;

deltP=0.03*(V/1000)^1.84/d^1.24

disp('L=3000m d=45mm V=1600m/min')

disp('压降计算值为：')

fprintf('\tdeltP=%.6f\n',deltP)

deltP=0.000635

4. 如下图所示，圆形横截面开口管道流体的流量Q(单位m3/s)由下式给出：

3/ 2 5/ 2 3/ 2

5/ 2

2 ( 0.5sin(2 ))

8 sin (1 cos )

c D g

Q

θ θ

θ θ

?

=

?

其中，g=9.8m/s2 为重力加速度，Dc 由下式给出：

(1 cos )

c 2

D = d ? θ

1) 编写一个Script 文件，计算d＝2m，θ＝60°时，Q 的值；

2) 编写一个计算Q 的函数文件，使d 和θ作为输入变量，Q 作为输入变量；并利用此函数

文件计算d＝3m，θ＝50°时的Q 值。

解：

1)

Theta=60*pi/180;

d=2;g=9.8;

Dc=d*(1-cos(Theta))/2;

Q=2^(3/2)*Dc^(5/2)*sqrt(g)*(Theta-

0.5*sin(2*Theta))^(3/2)/(8*sqrt(sin(Theta))*(1-cos(Theta))^(5/2))

Q = 0.5725

2)

function Q=QCal(d,Theta)

Theta=Theta*pi/180;

g=9.8;

Dc=d*(1-cos(Theta))/2;

Q=2^(3/2)*Dc^(5/2)*sqrt(g)*(Theta-

0.5*sin(2*Theta))^(3/2)/(8*sqrt(sin(Theta))*(1-cos(Theta))^(5/2));

>> Q1=QCal(3,50)

Q1 = 0.8171

5. 采用plot 命令一次画出两条函数曲线sint 和e?t cost ，t ∈[0,3π ] 。其中，第一条曲线采

用红色实线，第二条曲线采用蓝色虚线。给图形加上图题，图例和坐标轴名。

解：

t=[0:0.1*pi:3*pi];

y1=sin(t);

y2=exp(-t).*cos(t);

plot(t,y1,'-r',t,y2,'.b')

title('The plots of sin(t) and exp(-t)*cos(t)')

legend('sin(t)','exp(-t)*cos(t)')

xla