matlab 点名程序

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

华中科技大学数学与统计学院
一、随机点名clear;
disp('统计0801随机点名');
n=input('输入人数');
A=1:30;
A=A(randperm(30));
for x=1:n
switch A(x)
case {1}
disp('程传颖')
case {2}
disp('方涵')
case {3}
disp('朱秋红')
case {4}
disp('邹伦')
case {5}
disp('付路')
case {6}
disp('周亚航')
case {7}
disp('赵玉')
case {8}
disp('査若晨')
case {9}
disp('吴新骍')
case {10}
disp('李丹青')
case {11}
disp('吴雪惠')
case {12}
disp('邓军')
case {13}
disp('冯昌')
case {14}
disp('晏丽')
case {15}
disp('裴爽')
case {16}
disp('张矜矜')
case {17}
disp('杨益新')
case {18}
disp('朱一杰')
case {19}
disp('邱林丹')
case {20}
disp('代竣')
case {21}
disp('钟新宇')
case {22}
disp('肖志鸿')
case {23}
disp('刘垚')
case {24}
disp('沈兴')
case {25}
disp('余栋')
case {27}
disp('孟文洁')
case {28}
disp('许飞')
case{29}
disp('彭杨可')
case{30}
disp('程凌珺')
end
end
二、热轧机
clf,clear, m=20;l=3;mu=3.5;sigma=0.3; mu1=3.5;sigma1=0.9;
a=0; b=max([mu+4*sigma,mu1+4*sigma1]); %设定坐标的范围。

subplot(2,2,1),
for k=1:m
rand('seed',1), x=normrnd(mu,sigma);
plot([0,x],[k,k],'linewidth',5),hold on, axis([0 mu+4*0.6 -2 m+5])
end
plot([l,l],[-2,m+5],'r-');hold on,axis([a b -2 m+5])
subplot(2,2,3),
t=linspace(a,b,50);f=normpdf(t,mu,sigma);
y0=max(f)+0.1;plot(t,f),hold on,axis([a b 0 y0]),
plot([l,l],[0,y0],'r-',[mu,mu],[0,y0],'r:');hold on
t=linspace(a,b,30);f=normpdf(t,mu,sigma);plot(t,f)
subplot(2,2,2)
for k=1:m
rand('seed',1),x=normrnd(mu1,sigma1);
plot([0,x],[k,k],'linewidth',5),hold on,axis([a b -2 m+5])
end
plot([l,l],[-2,m+5],'r-');hold on
subplot(2,2,4),
t=linspace(a,b,50);f=normpdf(t,mu1,sigma1);y0=max(f)+0.1;
plot(t,f),hold on,axis([a b 0 y0])
plot([l,l],[0,y0],'r-',[mu1,mu1],[0,y0],'r:');hold on
clf,clear
l=2;sigma=0.2;
n=10000;m=50;a=2.2;b=3;
mu=linspace(a,b,m);
for k=1:m
x=normrnd(mu(k),sigma,1,n);
k_chenpin=find(x>=l);
k_feipin=find(x<l);
w_chenpin=x(k_chenpin)-l;
w_feipin=x(k_feipin);
if length(k_chenpin)==0
waste(k)=NaN;
else
waste(k)=(sum(w_chenpin)+sum(w_feipin))/length(k_chenpin); end
end
[wmin,i]=min(waste);
[mu(i) wmin]
plot(mu,waste,'.-',mu(i),wmin,'ro'),set(gca,'fontsize',15)
text(mu(i),wmin,['\bullet\leftarrow The Min Value is ',num2str(wmin),' at \it{mu}= ',num2str(mu(i))]);。

相关文档
最新文档