教育资源为主的文档平台

当前位置: 查字典文档网> 所有文档分类> 工程科技> 信息与通信> 数字仿真

数字仿真

上传者:李大成
|
上传时间:2015-05-10
|
次下载

数字仿真

二、仿真题

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