数学建模实验答案初等模型.docx

  1. 1、本文档共29页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实验02初等模型(4学时) (第2章初等模型) .(编程)光盘的数据容量 p23~27 表1 3种光盘的基本数据 激光器 激光波长 / ;jm 光斑直径 / m 信道间距 / mm (d) 数据线密度 / (B mm -1) (力 红外(CD) 0.78 2 1.6 x 10-3 121 红色(DVD) 0.64 0.92 0.74 X10-3 387 蓝色(DVD) 0.41 0.4 0.32 X10-3 800 CAV光盘:恒定角速度的光盘。 CLV光盘:恒定线速度的光盘。 R2=58 mm, R1=22.5 mm , d, p 见表 1。 _ 2 - 2 ,…,,、一 一 、 n(R2 -R2) CLV光盘的信息总长度(mm) Lclv & d CLV 光盘的信息容量(MB) Cclv = pLclv / (10A6) CLV 光盘的影像时间(min) Tclv = Cclv / (0.62 x60) R R2 CAV光盘的信息总长度(mm) Lcav上 2 2d CAV 光盘的信息容量(MB) Ccav = pLcav / (10A6) CAV 光盘的影像时间(min ) T cav = C cav / (0.62 乂 60) (验证、编程)模型求解 要求: ①(验证)分别计算出LCLV, CCLV和TCLV三个3行1列的列向量,仍后输出结 果,并与P26的表2 (教材)比较。 程序如下: clear;clc;format compact ; R1=22.5; R2=58; d=10A(-3)*[1.6,0.74,0.32]'; rho=[121,387,800]'; LCLV=pi*(R2A2-R1A2)./d; CCLV=rho*LCLV/10A6; %从 B 转换到 MB TCLV=CCLV/(0.62*60); % 从秒转换到分 s=' ';S=[s;s;s]; %s为两个空格,S为两列空格 [num2str(round(LCLV)),S, ... %其中的量为列向量 num2str(round(CCLV)),S, ... num2str(round(TCLV))] ②(编程)对于LCAV, CCAV和TCAV ,编写类似①的程序,并运行,结果与 P26 的表3 (教材)比较。 ★要求①的程序的运行结果: ★要求②的程序及其运行结果: clear;clc;format compact; R1=22.5; R2=58; d=10A(-3)*[1.6,0.74,0.32]'; rho=[121,387,800]'; LCAV=pi*R2A2./(2*d); CCAV=rho.*LCAV/10A6; TCAV=CCAV/(0.62*60); s=' ';S=[s;s;s]; [num2str(round(LCAV)),S, ... num2str(round(CCAV)),S, num2str(round(TCAV))] (编程)结果分析 信道长度LCLV 信道长度LCLV的精确计算:Lclv 模型给出的是近似值:L 模型给出的是近似值:L = _2 - 2 二(R2-R). Lclv 相对误差为:、:=L 相对误差为:、:= LCLV - L L 要求: ①取 R2=58 mm, R1=22.5 mm , d, p见表 1 (题 1)。 分别计算出LCLV, L和delta三个3行1列的列向量,仍后将它组合起来输出一个 3 行3列的结果。 ② 结果与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); LCLV(1)=quad( '2*pi/1.6e-3*sqrt(xA2+(1.6e-3/(2*pi))A2)' ,R1,R2); LCLV(2)=quad( '2*pi/0.74e-3*sqrt(xA2+Q74e-3/(2*pi))A2)' ,R1,R2); LCLV(3)=quad( '2*pi/0.32e-3*sqrt(xA2+(0.32e-3/(2*pi))A2)‘ ,R1,R2); L=pi*(R2A2-R1A2)./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.(验证,

文档评论(0)

fengnaifeng + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档