LCLV(1)=quad('2*pi/1.6e-3*sqrt(x.^2+(1.6e-3/(2*pi))^2)',R1,R2);
LCLV(2)=quad('2*pi/0.74e-3*sqrt(x.^2+(0.74e-3/(2*pi))^2)',R1,R2);
LCLV(3)=quad('2*pi/0.32e-3*sqrt(x.^2+(0.32e-3/(2*pi))^2)',R1,R2);
0
1
2
3
4
5
6
7
流量f(m3·h-1)
150.12
115.56
84.96
66.60
68.04
71.64
82.08
132.84
时间t(h)
8
9
10
11
12
13
14
15
流量f(m3·h-1)
185.04
226.80
246.60
250.92
261.00
271.44
273.96
279.00
时间t(h)
16
17ቤተ መጻሕፍቲ ባይዱ
程序如下：
clear;clc;formatcompact;
R1=22.5; R2=58;
d=10^(-3)*[1.6,0.74,0.32]';
rho=[121,387,800]';
LCLV=pi*(R2^2-R1^2)./d;
CCLV=rho.*LCLV/10^6;%从B转换到MB
TCLV=CCLV/(0.62*60);%从秒转换到分
[num2str(round(LCAV)),S,...
num2str(round(CCAV)),S,...
num2str(round(TCAV))]
1.2
信道长度LCLV的精确计算：
模型给出的是近似值：
相对误差为：
要求：
①取R2=58 mm, R1=22.5 mm，d,ρ见表1（题1）。
分别计算出LCLV, L和delta三个3行1列的列向量，仍后将它组合起来输出一个3行3列的结果。
参考图形结果：
要求：
①运行问题(1)中的程序。
②编程解决问题(2)：实际值与计算值比较（数据比较和和拟合图形）。
★（验证）
模型为：
★（编程）
程序：
n=[1 2 4 8]';
t=[7.21 6.88 6.32 5.84]';
logt=log(t); logn=log(n);
p=polyfit(logn,logt,1);
s=' ';S=[s;s;s];%s为两个空格，S为两列空格
[num2str(round(LCLV)),S,...%其中的量为列向量
num2str(round(CCLV)),S,...
num2str(round(TCLV))]
②（编程）对于LCAV, CCAV和TCAV，编写类似①的程序，并运行，结果与P26的表3（教材）比较。
L=pi*(R2^2-R1^2)./d;
delta=abs(LCLV-L)./abs(L);
s=' ';S=[s;s;s];
[num2str(round(LCLV)),S,...
num2str(round(L)),S,...
num2str(round(1000*delta)/100)]
运行结果：
2.
800
CAV光盘：恒定角速度的光盘。
CLV光盘：恒定线速度的光盘。
R2=58 mm, R1=22.5 mm，d,ρ见表1。
CLV光盘的信息总长度(mm) LCLV
CLV光盘的信息容量(MB) CCLV=ρLCLV/ (10^6)
CLV光盘的影像时间(min) TCLV= CCLV/ (0.62×60)
beta=p(1);
alfa=exp(p(2));
t2=alfa*n.^beta;
[n,t,t2]
a=0:0.01:10;
t3=alfa*a.^beta;
plot(n,t,'x',a,t3);
数值结果：
图形结果：
3.
表2 (p35)社区一天以小时为单位间隔的生活污水流量（单位：m3/h）
时间t(h)
logt=log(t);logn=log(n);
p=polyfit(logn,logt,1);% polyfit函数使用格式见提示
beta=p(1)
alfa=exp(p(2))
(2)实际值与计算值比较（数据比较和和拟合图形）
参考数据结果：
第1列为桨手人数，第2列为实际比赛平均成绩，第3列为计算比赛平均成绩。
★
★
clear;clc;formatcompact;
R1=22.5; R2=58;
d=10^(-3)*[1.6,0.74,0.32]';
rho=[121,387,800]';
LCAV=pi*R2^2./(2*d);
CCAV=rho.*LCAV/10^6;
TCAV=CCAV/(0.62*60);
s=' ';S=[s;s;s];
模型：t=αnβ
其中，t为比赛成绩（时间），n为桨手人数，α和β为参数。
为适合数据拟合，将模型改为：logt=logα+βlogn
桨手人数n
比赛平均成绩t
1
7.21
2
6.88
4
6.32
8
5.84
(1)参数α和β估计
程序如下：
clear;clc;
n=[1 2 4 8];%桨手人数
t=[ 7.21 6.88 6.32 5.84];%比赛平均成绩
CAV光盘的信息总长度(mm) LCAV
CAV光盘的信息容量(MB) CCAV=ρLCAV/ (10^6)
CAV光盘的影像时间(min ) TCAV= CCAV/ (0.62×60)
1.1
要求：
①（验证）分别计算出LCLV, CCLV和TCLV三个3行1列的列向量，仍后输出结果，并与P26的表2（教材）比较。
18
19
20
21
22
23
流量f(m3·h-1)
291.60
302.04
310.68
290.52
281.16
248.40
210.24
②结果与P26的表2和P27（教材）的结果比较。
[提示]
定积分计算用quad、quadl或trapz函数，注意要分别取d的元素来计算。要用数组d参与计算，可用quadv（用help查看其用法）。
★
程序：
R1=22.5;R2=58;
d=[1.6e-3,0.74e-3,0.32e-3]';
LCLV=zeros(3,1);
实验02初等模型（4学时）
（第2章 初等模型）
1.
表1 3种光盘的基本数据
激光器
激光波长/μm
光斑直径/μm
信道间距/mm
(d)
数据线密度/
(B·mm-1)
(ρ)
红外(CD)
0.78
2
1.6×10-3
121
红色(DVD)
0.64
0.92
0.74×10-3
387
蓝色(DVD)
0.41
0.4
0.32×10-3