|
s****1
发帖数: 425 | 2 我有一个200行的matlab code，由于数据量比较大，运行时间也比较长（几天），我
32G 内存已经不够用了，不知道转成C或者java是不是有帮助？谢谢
maxepoch=1500;
cmax=1500000;
for nb=1:3
for m=1:30
Qw(:,:,m,nb)=ones(cmax,1)*constant; % 1500000*4*30*3
mu(:,:,m,nb)=
sigma(:,:,m,nb)=
end
end
for g=1:2000
for nb=1:3
x=batchdata(:,:,nb);%20000*800*3
.....
for m=1:30
for i=1:numwits(:,nb)%1000000*3
w(i,m,nb)= normrnd(mu(i,j,m,nb),sigma(i,j,m,nb));% j=one of
[1 2 3 ... 16]
... 阅读全帖 |
|
s****1
发帖数: 425 | 3 我有一个200行的matlab code，由于数据量比较大，有几个百万量级的循环，运行时间
也比较长（几天），我
32G内存已经不够用了，不知道转成C或者java是不是有帮助？目前没有HPC。谢谢
为了简单，我只保留了结构。
maxepoch=1500;
cmax=1500000;
for nb=1:3
for m=1:30
Qw(:,:,m,nb)=ones(cmax,1)*constant; % 1500000*4*30*3
mu(:,:,m,nb)=
sigma(:,:,m,nb)=
end
end
for g=1:2000
for nb=1:3
x=batchdata(:,:,nb);%20000*800*3
.....
for m=1:30
for i=1:numwits(:,nb)%1000000*3
w(i,m,nb)= normrnd(mu(i,j,m,nb),sigma(i,j,m,n... 阅读全帖 |
|
b***i
发帖数: 3043 | 4 为什么要清除wx wxx wxxx wxxxx h hh hhh hhhh yx yxx yxxx yxxxx？他们的结果被
谁用了？
ii loop能否不用循环？
Qw， mu, sigma在哪里被用？
%Initialization
for nb=1:3
for m=1:30
Qw(:,:,m,nb)=ones(cmax,1)*constant; % 1500000*4*30*3
mu(:,:,m,nb)=
sigma(:,:,m,nb)=
end
end
%g loop
for g=1:2000
%nb loop
for nb=1:3
x=batchdata(:,:,nb);%20000*800*3
.....
% m loop
for m=1:30
for i=1:numwits(:,nb)%1000000*3
w(i,m,nb)= normrnd(mu(i,j,m,nb),sigma(i,j,m,nb)... 阅读全帖 |
|
s****1
发帖数: 425 | 5 在此谢谢所有回复，给了我很多启发。我发现matlab自带的函数bi2de 和normrnd比较
耗时，所以以后不要太相信它的函数了 |
|
s****1
发帖数: 425 | 6 我有一个200行的matlab code，由于数据量比较大，有几个百万量级的循环，能不用循
环的我都用矢量替代了。需要迭代上千次，运行时间也比较长（几天），我
32G 内存已经不够用了，不知道转成C或者java是不是有帮助？谢谢
maxepoch=1500;
cmax=1500000;
for nb=1:3
for m=1:30
Qw(:,:,m,nb)=ones(cmax,1)*constant; % 1500000*4*30*3
mu(:,:,m,nb)=
sigma(:,:,m,nb)=
end
end
for g=1:2000
for nb=1:3
x=batchdata(:,:,nb);%20000*800*3
.....
for m=1:30
for i=1:numwits(:,nb)%1000000*3
w(i,m,nb)= normrnd(mu(i,j,m,nb),sigma(i,j,m... 阅读全帖 |
|
s****1
发帖数: 425 | 7 我有一个200行的matlab code，由于数据量比较大，有几个百万量级的循环，能不用循
环的我都用矢量替代了。需要迭代上千次，运行时间也比较长（几天），我
32G 内存已经不够用了，不知道转成C或者java是不是有帮助？谢谢
maxepoch=1500;
cmax=1500000;
for nb=1:3
for m=1:30
Qw(:,:,m,nb)=ones(cmax,1)*constant; % 1500000*4*30*3
mu(:,:,m,nb)=
sigma(:,:,m,nb)=
end
end
for g=1:2000
for nb=1:3
x=batchdata(:,:,nb);%20000*800*3
.....
for m=1:30
for i=1:numwits(:,nb)%1000000*3
w(i,m,nb)= normrnd(mu(i,j,m,nb),sigma(i,j,m... 阅读全帖 |
|
s*******g
发帖数: 37 | 8 u made a stupid mistake in your code.
U didnot generate the norm-distributed data, since the variable [enm_bf_mcnx
] in your code could not come from standard norm.
U should try the following:
enm_bf_mcnx=normrnd(0,1,1000,1);
[mu,sigma]=normfit(enm_bf_mcnx) |
|
k**********g
发帖数: 33 | 9 n=20;
sigma=1;
A=normrnd(0,sigma,n,n);
Let V(i) = ( norm( A(:,i) ) )^2; i=1:n
what's the expected value of K=min {V(1),...,V(n)}
thanks a lot;
Simulation result shows that E(K)equals about 0.45 for all n>=2;
======================================================
if [Q,R]=qr(A);
Let S(i) = ( R(i,i))^2; i=1:n
T=min {S(1),...,S(n)}
Simulation result shows that E(K)=E(S)?
=====================================================
one lemma for reference:
2mR(i,i)^2 are chi-square random variables with 2( |
|
