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