clear all
clc
x=1e7;
chou=1e5;
R_o = normrnd(200,0.5^2,[1 1e7]);
L_o = normrnd(400.921,0.056^2,[1 1e7]);
%%
k_R = unidrnd(x,[1,chou]);
k_L = unidrnd(x,[1,chou]);
R_k = R_o(k_R);
L_k = L_o(k_L);
R_mean = mean(R_k);
L_mean = mean(L_k);

y_mean=L_mean-R_mean
y=L_k-R_k;
plot([1:chou],y)

2楼: Originally posted by acalephs at 2015-01-16 22:47:23
是这个意思吗？
y_mean是做差结果，4.009209989639458e+02
clear all
clc
x=1e7;
chou=1e5;
R_o = normrnd(200,0.5^2,);
L_o = normrnd(400.921,0.056^2,);
%%
k_R = unidrnd(x,);
k_L = unidrnd(x,);
R ...

3楼: Originally posted by 涛涛屈 at 2015-01-18 18:40:02
谢谢了，我先看看，不懂得再请教你。祝你事事顺心！...

clear all
clc
x=1e7;
chou=1e5;
R_o = normrnd(200,0.5,[1 1e7]);           %生成1e7个随机数，满足（200，.5的正态分布）
L_o = normrnd(400.921,0.056,[1 1e7]);  %生成1e7个随机数，满足（400.921,0.056^2）的正态分布
%%
k_R = unidrnd(x,[1,chou]);                     %语句的字面意思：从1到1e7中按均匀分布产生1行*1e5列的随机整数；语句的目的：从1e7个随机数抽取1e5个R的index，或者理解为数组下标
k_L = unidrnd(x,[1,chou]);                     %从1e7个随机数抽取1e5个L的index，或者理解为数组下标
R_k = R_o(k_R);                                   %从R_o按k_R下标抽1e5个数
L_k = L_o(k_L);                                   %从R_o按k_L下标抽1e5个数
R_mean = mean(R_k);
L_mean = mean(L_k);

y_mean=L_mean-R_mean
y=L_k-R_k;
plot([1:chou],y)

4楼: Originally posted by acalephs at 2015-01-19 09:58:42
不好意思，上一个程序不应该带平方。。。
clear all
clc
x=1e7;
chou=1e5;
R_o = normrnd(200,0.5,);           %生成1e7个随机数，满足（200，.5的正态分布）
L_o = normrnd(400.921,0.056,);  %生成1e7个随 ...

6楼: Originally posted by acalephs at 2015-01-22 10:45:36
把for循环删掉，找大于200的点有多少用这个命令。
length(find(y>200))

另外提问的时候发帖子，不要站内信，别人也可能有这样的问题，大家一起进步。

6楼: Originally posted by acalephs at 2015-01-22 10:45:36
把for循环删掉，找大于200的点有多少用这个命令。
length(find(y>200))

另外提问的时候发帖子，不要站内信，别人也可能有这样的问题，大家一起进步。