数字仿真
二、仿真题
1. 记第K个工地的位置为(xk,yk),中心的位置为(x0,y0),则目标函数为
内容需要下载文档才能查看
由上市可进行编程确定搅拌中心的位置使费用最少。
2.系统中,实体是水,属性是水的体积、含盐量和含盐率,活动是水的注入与流出,由于注入和流出活动的作用,池中水的体积与含盐量、含盐率均随时间变化,初始时刻含盐率为0.001kg/m3,以后每分钟注入含盐率为0.5 kg/m3的水6 m3,流出混合均匀的盐水4m3,当池中水的含盐率达到0.2kg/m3时,结束仿真过程。 h=1; s0=2; w0=2000; r0=s0/w0;
s(1)=s0+0.5*6*h-4*h*r0;
w(1)=w0+2*h;
r(1)=s(1)/w(1); t(1)=h;
y(1)=(2000000+3000000*h+3000*h^2+h^3)/(1000+h)^2;
for i=2:200
t(i)=i*h;
s(i)=s(i-1)+0.5*6*h-4*h*r(i-1);
w(i)=w(i-1)+2*h;
r(i)=s(i)/w(i);
y(i)=(2000000+3000000*t(i)+3000*t(i)^2+t(i)^3)/(1000+t(i))^2;
m=floor(i/10);
if i/10-m<0.1
tm(m)=m;
wm(m)= w(i),
sm(m)=s(i);
rm(m)=r(i);
end
if r(i)>0.2
t02=i*h;
r02=r(i);
break
end
end
[t02,r02]
[10*tm',sm',rm']
plot(t,s,'blue',t,y,'red')
3.
%**************************线性同余法求随机数******************************%
clc;
close all;
clear all;
n_randn=3000;
a=1220703125;b=35;m=2^b-1;x=2;
U=zeros(1,n_randn);
n=1;
while n<=n_randn
n=n+1;
x=rem(a*x,m);
r=x/m;
U(1,n-1)=r;
end
U
%*************************数字特征检验其均匀性*****************************%
sum1=0;sum2=0;sum3=0;
for i=1:n_randn
sum1=sum1+U(i);
sum2=sum2+(U(i)-1/2)^2;
sum3=sum3+U(i)^2;
end
y=sum1/n_randn;
s=sum2/(n_randn-1);
Y=sum3/n_randn;
z_std=2.58;
z(1)=sqrt(12*n_randn)*(y-1/2);
z(2)=sqrt(180*n_randn)*(s-1/12);
z(3)=sqrt(45*n_randn)*(Y-1/3);
for i=1:3
if z(i)>z_std&&z(i)<(-z_std)
R_ind(i)=1;
else
R_ind(i)=0;
end
end
R_ind %若其输出全为0,则表明此随机序列具有均匀性,否则,不具有均匀性
%******************************游程检验***********************************%
Q1=0;
Q2=0;
Q3=0;
Q4=0;
Q5=0;
Q6=0;
m=1;
for i=1:n_randn-1
j=i+1;
k=U(i);
l=U(j);
if k<=l
m=m+1;
else switch(m)
case 1
Q1=Q1+1;
case 2
Q2=Q2+1;
case 3
Q3=Q3+1;
case 4
Q4=Q4+1;
case 5
Q5=Q5+1;
otherwise
Q6=Q6+1;
end
m=1;
end
if i==n_randn-1
if k>l
Q1=Q1+1;
end
end
end
q=[Q1,Q2,Q3,Q4,Q5,Q6];
b=[1/6,5/24,11/120,19/720,29/5040,1/840];
a=[4529.4,9044.9,13568,18091,22615,27892;9044.9,18097,27139,36187,45234,55789;13568,27139,40721,54281,67852,83685;18091,36187,54281,72414,90470,111580;22615,45234,67852,90470,113262,139476;27892,55789,83685,111580,139476,172860];
z=0;
for i=1:6
for j=1:6
x=a(i,j)*(q(i)-n_randn*b(i))*(q(j)-n_randn*b(j));
z=z+x;
end
end
x=z/n_randn;
P_ind1= 14.4;
P_ind2= 1.24;
if x>P_ind1&&x<(P_ind2)
R_ind4=1;
else
R_ind4=0;
end
R_ind4 %若R输出为0时,则表明此随机序列具有独立性,否则,不具有独立性
%************************matlab自带函数产生随机数**************************%
clc;
close all;
clear all;
n_randn=3000;
U=rand(1,n_randn);
U
%*************************数字特征检验其均匀性*****************************%
sum1=0;sum2=0;sum3=0;
for i=1:n_randn
sum1=sum1+U(i);
sum2=sum2+(U(i)-1/2)^2;
sum3=sum3+U(i)^2;
end
y=sum1/n_randn;
s=sum2/(n_randn-1);
Y=sum3/n_randn;
z_std=2.58;
z(1)=sqrt(12*n_randn)*(y-1/2);
z(2)=sqrt(180*n_randn)*(s-1/12);
z(3)=sqrt(45*n_randn)*(Y-1/3);
for i=1:3
if z(i)>z_std&&z(i)<(-z_std)
R_ind(i)=1;
else
R_ind(i)=0;
end
end
R_ind %若其输出全为0,则表明此随机序列具有均匀性,否则,不具有均匀性
%******************************游程检验***********************************%
Q1=0;
Q2=0;
Q3=0;
Q4=0;
Q5=0;
Q6=0;
m=1;
for i=1:n_randn-1
j=i+1;
k=U(i);
l=U(j);
if k<=l
m=m+1;
else switch(m)
case 1
Q1=Q1+1;
case 2
下载文档
热门试卷
- 2016年四川省内江市中考化学试卷
- 广西钦州市高新区2017届高三11月月考政治试卷
- 浙江省湖州市2016-2017学年高一上学期期中考试政治试卷
- 浙江省湖州市2016-2017学年高二上学期期中考试政治试卷
- 辽宁省铁岭市协作体2017届高三上学期第三次联考政治试卷
- 广西钦州市钦州港区2016-2017学年高二11月月考政治试卷
- 广西钦州市钦州港区2017届高三11月月考政治试卷
- 广西钦州市钦州港区2016-2017学年高一11月月考政治试卷
- 广西钦州市高新区2016-2017学年高二11月月考政治试卷
- 广西钦州市高新区2016-2017学年高一11月月考政治试卷
- 山东省滨州市三校2017届第一学期阶段测试初三英语试题
- 四川省成都七中2017届高三一诊模拟考试文科综合试卷
- 2017届普通高等学校招生全国统一考试模拟试题(附答案)
- 重庆市永川中学高2017级上期12月月考语文试题
- 江西宜春三中2017届高三第一学期第二次月考文科综合试题
- 内蒙古赤峰二中2017届高三上学期第三次月考英语试题
- 2017年六年级(上)数学期末考试卷
- 2017人教版小学英语三年级上期末笔试题
- 江苏省常州西藏民族中学2016-2017学年九年级思想品德第一学期第二次阶段测试试卷
- 重庆市九龙坡区七校2016-2017学年上期八年级素质测查(二)语文学科试题卷
- 江苏省无锡市钱桥中学2016年12月八年级语文阶段性测试卷
- 江苏省无锡市钱桥中学2016-2017学年七年级英语12月阶段检测试卷
- 山东省邹城市第八中学2016-2017学年八年级12月物理第4章试题(无答案)
- 【人教版】河北省2015-2016学年度九年级上期末语文试题卷(附答案)
- 四川省简阳市阳安中学2016年12月高二月考英语试卷
- 四川省成都龙泉中学高三上学期2016年12月月考试题文科综合能力测试
- 安徽省滁州中学2016—2017学年度第一学期12月月考高三英语试卷
- 山东省武城县第二中学2016.12高一年级上学期第二次月考历史试题(必修一第四、五单元)
- 福建省四地六校联考2016-2017学年上学期第三次月考高三化学试卷
- 甘肃省武威第二十三中学2016—2017学年度八年级第一学期12月月考生物试卷
网友关注
- 2010年安徽公务员考试面试真题及解析(三)
- 2012年安徽公务员考试申论(B)作答思路
- 2007年安徽省公务员考试面试真题解析
- 2012年安徽省公务员面试真题(6月11日)
- 2012年安徽省公务员面试真题解析(6月10日)
- 2010年安徽公务员考试面试真题及解析(一)
- 2011年安徽省公务员考试《申论》冲刺预测试卷
- 历年安徽省公务员考试行测题型题量分析
- 2013年413省考行测试卷冲刺卷
- 2012年安徽公务员考试申论(A)作答思路
- 乡镇公务员公共基础模拟试卷(二)及参考答案
- 2010年安徽公务员考试面试真题及解析(四)
- 2011安徽省考申论B卷深度解读
- 2010年安徽省公务员考试面试真题
- 2013年安徽省公务员录用考试《申论》A卷
- 2011年安徽省从优秀村干部中考试录用乡镇机关公务员面试真题
- 2012年安徽省公务员面试真题(6月9日)
- 乡镇公务员公共基础模拟试卷(四)及参考答案
- 2011年安徽省公务员考试行测真题
- 2007年安徽省公务员考试面试真题
- 2012年安徽省公务员面试真题(6月10日)
- 2013年413多省公务员考试申论冲刺试卷(四) 参考答案
- 2010年安徽公务员考试面试真题及解析(二)
- 2011年安徽省国税系统招录公务员面试真题(3月1日)
- 乡镇公务员公共基础模拟试卷(一)及参考答案
- 2013年413省考行测冲刺试卷答案
- 2011年安徽省公务员录用考试行测冲刺预测试卷
- 2010年8月9日安徽省公务员考试面试真题解析
- 2009年安徽省公务员考试面试真题解析
- 乡镇公务员公共基础模拟试卷(三)及参考答案
网友关注视频
- 二年级下册数学第三课 搭一搭⚖⚖
- 沪教版牛津小学英语(深圳用) 四年级下册 Unit 7
- 七年级英语下册 上海牛津版 Unit3
- 沪教版八年级下册数学练习册21.3(3)分式方程P17
- 8 随形想象_第一课时(二等奖)(沪教版二年级上册)_T3786594
- 苏科版八年级数学下册7.2《统计图的选用》
- 外研版八年级英语下学期 Module3
- 七年级下册外研版英语M8U2reading
- 外研版英语七年级下册module3 unit1第二课时
- 《小学数学二年级下册》第二单元测试题讲解
- 冀教版小学数学二年级下册第二单元《有余数除法的简单应用》
- 小学英语单词
- 沪教版八年级下册数学练习册20.4(2)一次函数的应用2P8
- 第19课 我喜欢的鸟_第一课时(二等奖)(人美杨永善版二年级下册)_T644386
- 第12章 圆锥曲线_12.7 抛物线的标准方程_第一课时(特等奖)(沪教版高二下册)_T274713
- 三年级英语单词记忆下册(沪教版)第一二单元复习
- 北师大版小学数学四年级下册第15课小数乘小数一
- 19 爱护鸟类_第一课时(二等奖)(桂美版二年级下册)_T502436
- 冀教版小学数学二年级下册第二周第2课时《我们的测量》宝丰街小学庞志荣
- 外研版英语三起5年级下册(14版)Module3 Unit1
- 【部编】人教版语文七年级下册《泊秦淮》优质课教学视频+PPT课件+教案,广东省
- 【部编】人教版语文七年级下册《逢入京使》优质课教学视频+PPT课件+教案,安徽省
- 沪教版牛津小学英语(深圳用) 五年级下册 Unit 10
- 北师大版八年级物理下册 第六章 常见的光学仪器(二)探究凸透镜成像的规律
- 30.3 由不共线三点的坐标确定二次函数_第一课时(市一等奖)(冀教版九年级下册)_T144342
- 外研版英语三起5年级下册(14版)Module3 Unit2
- 冀教版小学英语五年级下册lesson2教学视频(2)
- 冀教版英语五年级下册第二课课程解读
- 七年级英语下册 上海牛津版 Unit9
- 第8课 对称剪纸_第一课时(二等奖)(沪书画版二年级上册)_T3784187
精品推荐
- 2016-2017学年高一语文人教版必修一+模块学业水平检测试题(含答案)
- 广西钦州市高新区2017届高三11月月考政治试卷
- 浙江省湖州市2016-2017学年高一上学期期中考试政治试卷
- 浙江省湖州市2016-2017学年高二上学期期中考试政治试卷
- 辽宁省铁岭市协作体2017届高三上学期第三次联考政治试卷
- 广西钦州市钦州港区2016-2017学年高二11月月考政治试卷
- 广西钦州市钦州港区2017届高三11月月考政治试卷
- 广西钦州市钦州港区2016-2017学年高一11月月考政治试卷
- 广西钦州市高新区2016-2017学年高二11月月考政治试卷
- 广西钦州市高新区2016-2017学年高一11月月考政治试卷
分类导航
- 互联网
- 电脑基础知识
- 计算机软件及应用
- 计算机硬件及网络
- 计算机应用/办公自动化
- .NET
- 数据结构与算法
- Java
- SEO
- C/C++资料
- linux/Unix相关
- 手机开发
- UML理论/建模
- 并行计算/云计算
- 嵌入式开发
- windows相关
- 软件工程
- 管理信息系统
- 开发文档
- 图形图像
- 网络与通信
- 网络信息安全
- 电子支付
- Labview
- matlab
- 网络资源
- Python
- Delphi/Perl
- 评测
- Flash/Flex
- CSS/Script
- 计算机原理
- PHP资料
- 数据挖掘与模式识别
- Web服务
- 数据库
- Visual Basic
- 电子商务
- 服务器
- 搜索引擎优化
- 存储
- 架构
- 行业软件
- 人工智能
- 计算机辅助设计
- 多媒体
- 软件测试
- 计算机硬件与维护
- 网站策划/UE
- 网页设计/UI
- 网吧管理