版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、实验02 初等模型(4学时)(第2章 初等模型)1.(编程)光盘的数据容量p2327表1 3种光盘的基本数据激光器激光波长/m光斑直径/m信道间距/mm(d)数据线密度/(B·mm-1)()红外(CD)0.7821.6×10-3121红色(DVD)0.640.920.74×10-3387蓝色(DVD)0.410.40.32×10-3800CAV光盘:恒定角速度的光盘。CLV光盘:恒定线速度的光盘。R2=58 mm, R1=22.5 mm,d, 见表1。CLV光盘的信息总长度(mm) LCLV CLV光盘的信息容量(MB) CCLV = LCLV / (1
2、06)CLV光盘的影像时间(min) TCLV = CCLV / (0.62×60)CAV光盘的信息总长度(mm) LCAVCAV光盘的信息容量(MB) CCAV = LCAV / (106)CAV光盘的影像时间(min ) TCAV = CCAV / (0.62×60)1.1(验证、编程)模型求解要求:(验证)分别计算出LCLV, CCLV和TCLV三个3行1列的列向量,仍后输出结果,并与P26的表2(教材)比较。程序如下:clear;clc;format compact;R1=22.5; R2=58;d=10(-3)*1.6,0.74,0.32'rho=121,
3、387,800' LCLV=pi*(R22-R12)./d;CCLV=rho.*LCLV/106;%从B转换到MBTCLV=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;R
4、1=22.5; R2=58;d=10(-3)*1.6,0.74,0.32'rho=121,387,800'LCAV=pi*R22./(2*d);CCAV=rho.*LCAV/106;TCAV=CCAV/(0.62*60);s=' 'S=s;s;s;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和de
5、lta三个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(x.2+(1.6e-3/(2*pi)2)',R1,R2);LCLV(2)=quad('2*pi/0.74e
6、-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);L=pi*(R22-R12)./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.(验证,编程)划艇比赛的成绩p2931模型:t=n其中,t为比赛成绩(时间),n为
7、桨手人数,和为参数。为适合数据拟合,将模型改为:log t=log + log n桨手人数n比赛平均成绩t17.2126.8846.3285.84(1) 参数和估计程序如下:clear; clc;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);% polyfit函数使用格式见提示beta=p(1)alfa=exp(p(2)(2) 实际值与计算值比较(数据比较和和拟合图形)参考数据结果:第1列为桨手人数,第2列为实际比赛平均成绩,第3列为计算比赛平均成绩。
8、参考图形结果:要求: 运行问题(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);beta=p(1);alfa=exp(p(2);t2=alfa*n.beta;n,t,t2a=0:0.01:10;t3=alfa*a.beta;plot(n,t,'x'
9、;,a,t3);数值结果:图形结果:3.(编程,验证)污水均流池的设计p3437均流池净化设备g 恒速 成 f 变速 成 表2 (p35) 社区一天以小时为单位间隔的生活污水流量(单位:m3/h)时间t (h)01234567流量f (m3·h-1)150.12115.5684.9666.6068.0471.6482.08132.84时间t (h)89101112131415流量f (m3·h-1)185.04226.80246.60250.92261.00271.44273.96279.00时间t (h)1617181920212223流量f (m3·h-1)2
10、91.60302.04310.68290.52281.16248.40210.24186.843.1(编程)均流池的恒定流出量和最大容量模型(离散)每小时污水流入均流池的流量为f (t), t=0, 1, 2, , 23。一天的平均流量 均流池中污水的空量 c(t), t=0, 1, 2, , 23。c(t+1)=c(t)+f(t)-g, t=0, 1, 2, , 22 (模型)要求: 求g,画f(t)和g的图形(与P35图1比较)。 求c(t), t=0, 1, 2, , 23, c(0)=0,并求其中的最小值M(与P36表3比较)。 求c(t), t=0, 1, 2, , 23, c(0)
11、=-M(与P36表4比较)。画c(t)分别当c(0)和c(-M)时的图形(与P37图2比较)。 要求的程序和运行结果:程序:t=0:23;f=150.12 115.56 84.96 66.60 68.04 71.64,. 82.08 132.84 185.04 226.80 246.60 250.92,. 261.00 271.44 273.96 279.00 291.60 302.04,. 310.68 290.52 281.16 248.40 210.24 186.84;s=0;for i=1:24 s=s+f(i);endg=s/24t2=0:0.01:23;plot(t,f,t2,g,
12、'r-');text(10,213.67,'g=203.67');grid on命令窗口的结果:图形窗口的结果: 要求的程序和运行结果:程序:t=0:23;f=150.12 115.56 84.96 66.60 68.04 71.64,. 82.08 132.84 185.04 226.80 246.60 250.92,. 261.00 271.44 273.96 279.00 291.60 302.04,. 310.68 290.52 281.16 248.40 210.24 186.84;s=0;for i=1:24 s=s+f(i);endg=s/24;c
13、1(1)=0;for j=1:23 c1(j+1)=c1(j)+f(j)-g;endc1M=min(c1)c2(1)=-M;for k=1:23 c2(k+1)=c2(k)+f(k)-g;endc2plot(t,c1,t,c2);text(7,160,'c(0)=876.15');text(7,-760,'c(0)=0');xlabel('t');ylabel('c');grid on命令窗口的结果:图形窗口的结果:3.2(验证)均流池的恒定流出量和最大容量模型(连续)p56习题3每小时污水流入均流池的流量为f (t), t=0,
14、 1, 2, , 23。用3次样条插值得到连续函数f(t), 0t23。(仍用f(t)表示)一天的平均流量 均流池中污水的容量 c(t) , 0t23。c(t+t)-c(t)=(f(t)-g) t (模型)(1) 求g,画f(t)和g的图形(与P35图1比较)。程序:function y=f(t)tt=0:23;ft=150.12 115.56 84.96 66.60 68.04 71.64,. 82.08 132.84 185.04 226.80 246.60 250.92,. 261.00 271.44 273.96 279.00 291.60 302.04,. 310.68 290.52
15、 281.16 248.40 210.24 186.84;y=interp1(tt,ft,t,'spline'); %3次样条插值clear; clc;g=1/(23-0)*quad('f',0,23)t=0:0.0001:23;plot(t,f(t),0,25,g,g);text(10+0.7,g+13,'g=',num2str(g);xlabel('itt');ylabel('itf');grid on;(2) 求c(t), 0t23, c(0)=0时的最小值M。画c(t)初值条件分别为c(0)=0和c(0)=
16、-M时的图形(与P37图2比较)。程序:function y=c(t,c0,g) %c0, g将作为参数y=quad('f',0,t)-t*g+c0;clear; clc;t=0:0.1:23;g=1/(23-0)*quad('f',0,23);ct0=zeros(size(t);for i=1:length(t) ct0(i)=c(t(i),0,g);endtt,m=fminbnd(t)c(t,0,g),0,23) %求最小值,注意函数c的参数格式ctm=zeros(size(t);for i=1:length(t) ctm(i)=c(t(i),-m,g);e
17、ndplot(t,ct0,t,ctm);text(7,c(8,0,g)+100,'c(0)=',num2str(ct0(1);text(7,c(8,-m,g)+100,'c(0)=',num2str(ctm(1);xlabel('itt');ylabel('itf');grid on;要求 运行(1)中的程序,结果与P35图1比较。 运行(2)中的程序,结果与P37图2比较。 阅读并理解程序。 要求的运行结果:命令窗口的结果:图形窗口的结果: 要求的运行结果:命令窗口的结果:图形窗口的结果:4.(编程)天气预报的评价p495431
18、天4种(AD)预报方法的有雨预报(%)及实际观测结果function M=tab()%日期ABCD有雨=1/无雨=0M=1903090601;2403050801;3603080701;4603090701;560300200;6303010501;7803010400;8703020300;9803040300;10603060400;11803020801;124030304001450306020016603050801180300500;19903060400;2070301000;2120300300;22
19、403020300;23403010100;24803050400;2530300200;26303010300;2730302000;2803060401;2960300200;30203010100;31803050100;4.1(编程求解)计数模型p5052若预报有雨概率>50%,则认为明天有雨,<50%则认为无雨,且依照明天是否有雨的实际观测,规定预报是否正确,从而统计预报的正确率。求出4种预报的结果计数矩阵: 预报和实测都有雨的天数 预报有雨而实测无雨的天数 预报无雨而实测有雨的天数 预报和实测都无雨的天数 预报的正确率:对角线数字之和/全部数之和。要求: 编写程序求出4
20、种预报的结果计数(天数),并分别计算出它们的预报正确率(取2位小数)。 结果与p51中的结果比较。 程序和运行结果:程序:function M=tab()M=1 90 30 90 60 1; 2 40 30 50 80 1; 3 60 30 80 70 1; 4 60 30 90 70 1; 5 60 30 0 20 0; 6 30 30 10 50 1; 7 80 30 10 40 0; 8 70 30 20 30 0; 9 80 30 40 30 0; 10 60 30 60 40 0; 11 80 30 20 80 1; 12 40 30 30 40 0; 13 90 30 90 40
21、1; 14 50 30 60 20 0; 15 10 30 20 10 0; 16 60 30 50 80 1; 17 20 30 10 30 0; 18 0 30 0 50 0; 19 90 30 60 40 0; 20 70 30 10 0 0; 21 20 30 0 30 0; 22 40 30 20 30 0; 23 40 30 10 10 0; 24 80 30 50 40 0; 25 30 30 0 20 0; 26 30 30 10 30 0; 27 30 30 20 0 0; 28 0 30 60 40 1; 29 60 30 0 20 0; 30 20 30 10 10 0;
22、 31 80 30 50 10 0;for j=2:5 a=0;b=0;c=0;d=0; for i=1:31 if M(i,j)>50 && M(i,6)=1 a=a+1; end if M(i,j)>50 && M(i,6)=0 b=b+1; end if M(i,j)<50 && M(i,6)=1 c=c+1; end if M(i,j)<50 && M(i,6)=0 d=d+1; end end x=a,b;c,d p=vpa(a+d)/(a+b+c+d),2)end运行结果:4.2(编程求解)记分
23、模型p5253将预报有雨概率的大小与实测结果(有雨或无雨)比较,给予记分。注意:要将M中的预报概率值转换为小数。模型1记第k天某种预报有雨概率为pk,第k天实测有雨为vk=1,无雨为vk=0,令第k天的某种预报得分为将sk对k求和得到某预报的分数S1(越大越好)。模型2sk = | pk - vk |将sk对k求和得到某预报的分数S2(越小越好)。模型3sk = ( pk - vk )2将sk对k求和得到某预报的分数S3(越小越好)。要求: 编程求4种预报在模型1、2、3下的相应分数S1、S2、S3。 运行结果与p52的结果比较。 程序和运行结果:function M=tab()M=1 90
24、30 90 60 1; 2 40 30 50 80 1; 3 60 30 80 70 1; 4 60 30 90 70 1; 5 60 30 0 20 0; 6 30 30 10 50 1; 7 80 30 10 40 0; 8 70 30 20 30 0; 9 80 30 40 30 0; 10 60 30 60 40 0; 11 80 30 20 80 1; 12 40 30 30 40 0; 13 90 30 90 40 1; 14 50 30 60 20 0; 15 10 30 20 10 0; 16 60 30 50 80 1; 17 20 30 10 30 0; 18 0 30 0
25、 50 0; 19 90 30 60 40 0; 20 70 30 10 0 0; 21 20 30 0 30 0; 22 40 30 20 30 0; 23 40 30 10 10 0; 24 80 30 50 40 0; 25 30 30 0 20 0; 26 30 30 10 30 0; 27 30 30 20 0 0; 28 0 30 60 40 1; 29 60 30 0 20 0; 30 20 30 10 10 0; 31 80 30 50 10 0;M(:,2:5)=M(:,2:5)/100;for i=2:5 S1=0; S2=0; S3=0; for j=1:31 s1=(-
26、1)M(j,6)*(0.5-M(j,i); S1=S1+s1; s2=abs(M(j,i)-M(j,6); S2=S2+s2; s3=(M(j,i)-M(j,6)2; S3=S3+s3; end i-1 S1 S2 S3endans = 1S1 = 1S2 = 14.5S3 = 8.95ans = 2S1 = 2.6S2 = 12.9S3 = 6.39ans = 3S1 = 7S2 = 8.5S3 = 4.23ans = 4S1 = 6.7S2 = 8.8S3 =3.224.3(部分编程求解)图形模型模型1p53以预报有雨概率p(值为小数)为横轴,实测值v(值为0或1)为纵轴,奖表tab的数据
27、在图上用符号*标出,其中*上面的数字是坐标在*的天数。预报A的程序:clear; clc;M=tab();M(:,2:5)=M(:,2:5)/100;%概率值都改为小数K=logical(M(:,6);%将M(:,6)中非0值的位置赋逻辑值true,MATLAB建议用logical代替findpp=M(K,2);%K中为true的对应位置的M元素保留,其中有重复概率值pp=sort(pp);%按升序排序p,I,=unique(pp,'last');%去掉重复值,p(i)值在pp中最后出现的行下标为I(i)c=I-0;I(1:end-1);%求p(i)在pp中重复个数c(i)pl
28、ot(0,1,1,1,':',p,ones(size(p),'b*');%画v=1的水平线和概率值p处的*for i=1:length(c)%标注重复数T=text(p(i),1+0.05,num2str(c(i);set(T,'color','b');end%自己完成v=0时的程序,把该图与上面的图合在一起xlabel('itp');ylabel('itv');text(0.5,0.5,'预报A');运行结果示例:要求: 自己完成上面未完整的程序并运行。 修改预报A的程序,分别用于
29、B、C、D,并运行。 运行结果与p53中的结果比较。 预报A的完整程序:M=1 90 30 90 60 1; 2 40 30 50 80 1; 3 60 30 80 70 1; 4 60 30 90 70 1; 5 60 30 0 20 0; 6 30 30 10 50 1; 7 80 30 10 40 0; 8 70 30 20 30 0; 9 80 30 40 30 0; 10 60 30 60 40 0; 11 80 30 20 80 1; 12 40 30 30 40 0; 13 90 30 90 40 1; 14 50 30 60 20 0; 15 10 30 20 10 0; 16
30、 60 30 50 80 1; 17 20 30 10 30 0; 18 0 30 0 50 0; 19 90 30 60 40 0; 20 70 30 10 0 0; 21 20 30 0 30 0; 22 40 30 20 30 0; 23 40 30 10 10 0; 24 80 30 50 40 0; 25 30 30 0 20 0; 26 30 30 10 30 0; 27 30 30 20 0 0; 28 0 30 60 40 1; 29 60 30 0 20 0; 30 20 30 10 10 0; 31 80 30 50 10 0;M(:,2:5)=M(:,2:5)/100;K
31、=logical(M(:,6);pp=M(K,2);pp=sort(pp);p,I,=unique(pp,'last');c=I-0;I(1:end-1);plot(0,1,1,1,':',p,ones(size(p),'b*');for i=1:length(c) T=text(p(i),1+0.05,num2str(c(i); set(T,'color','b');endK=logical(M(:,6)-1);pp=M(K,2);pp=sort(pp);p,I,=unique(pp,'last');c=I-0;I(1:end-1);hold onplot(p,zeros(size(p),'r*');hold offfor i=1:l
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2024教练车转让合同范文
- 2024临时工聘用合同规范
- 2024购销板材合同范本范文
- 颈椎病健康科普知识宣教
- 2024市内运输合同范本
- 苏州科技大学天平学院《展示设计》2022-2023学年第一学期期末试卷
- 2024脚手架搭设合同
- 2024办公文档范本土地使用权出让转让合同公证书
- 农业科学与农村公共设施考核试卷
- 家用纺织品的品牌知名度与市场认知度考核试卷
- 衢州江山市交投人力资源有限公司招聘笔试题库2024
- 职业资格证书遗失补发申请表
- 九年级化学上册 第1单元 走进化学世界教案 (新版)新人教版
- 18慈母情深 第二课时公开课一等奖创新教学设计
- 2024年部编新改版语文小学三年级上册第五单元测试题及答案
- 2024年全国数据应用大赛“数字安全赛”备赛试题库(含答案)
- 上市公司并购重组对股价的影响分析开题报告(修改)
- (必会)企业人力资源管理师(三级)近年考试真题题库(含答案解析)
- 2024年职业院校技能大赛数字化装备制造技术赛项理论题库(备赛500题)
- 《中外历史纲要》(下)
- 电话销售人员聘用合同
评论
0/150
提交评论