I'm using a for loop to create a new random variable X with range [0 1], I want the forloop to ignore values above 1 and iterate untill all the values are in range.

Question

Fredrik Slungård 2015-2-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/180074-i-m-using-a-for-loop-to-create-a-new-random-variable-x-with-range-0-1-i-want-the-forloop-to-ignor

编辑： Sean de Wolski 2015-2-25

clc; clear all; close all;
n = 200; % Number of iteration
ro = 0.75; % Correlation factor
X = zeros(1,n);
x = zeros(1,n);
for ii = 1:n
    x = rand(1,ii);
    y = rand(1,ii);
    X = x*ro+y*sqrt(1-ro^2);
    if X > 1 % If X > 1 start iteration no. ii again until X is below 1 
       ii = ii - 1; 
    end
end
plot(x,X,'+')

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Niels 2015-2-25

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/180074-i-m-using-a-for-loop-to-create-a-new-random-variable-x-with-range-0-1-i-want-the-forloop-to-ignor#answer_169207

编辑：Niels 2015-2-25

在 MATLAB Online 中打开

Change your for-loop to a while-loop if you want to keep looping. I also believe that you stored your data incorrectly, as x, y and X were not stored at every time step, while your random number generation grows larger for every iteration ii...

 clc; clear all; close all;
 n = 200; % Number of iteration
 ro = 0.75; % Correlation factor
 X = zeros(1,n);
 x = zeros(1,n);
 y = zeros(1,n);
 ii = 1;
 while (ii <= n)
    x(ii) = rand();
    y(ii) = rand();
    X(ii) = x(ii)*ro+y(ii)*sqrt(1-ro^2);
    if (X(ii) <= 1)     % Don't increment if X > 1. 
       ii = ii + 1; 
    end
 end
 figure, plot(x,X,'+')

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Fredrik Slungård 2015-2-25

编辑：Fredrik Slungård 2015-2-25

Thank you for the explanation!

请先登录，再进行评论。

Answer 2

Joseph Cheng 2015-2-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/180074-i-m-using-a-for-loop-to-create-a-new-random-variable-x-with-range-0-1-i-want-the-forloop-to-ignor#answer_169094

在 MATLAB Online 中打开

you'd want to use "continue" to skip the iteration.

 for ind =1:5
    if ind==3
       continue
 else
    disp(ind)
 end

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Fredrik Slungård 2015-2-25

编辑：Fredrik Slungård 2015-2-25

在 MATLAB Online 中打开

Sorry for taking so long answering, as I'm relative new to Matlab I've never used the continue function. But I don't want to skip the iteration, only redo iteration no. i if X > 1, and I still get about 30 values where X > 1.

clc; clear all; close all;
n = 200;
ro = 0.75;
X = zeros(1,n);
x = zeros(1,n);
for ii = 1:n
    x = rand(1,ii);
    y = rand(1,ii);
    X = x*ro+y*sqrt(1-ro^2);
    if X > 1;
       continue
    else 
        disp(X)
    end
end
plot(x,X,'+')

请先登录，再进行评论。

Answer 3

Sean de Wolski 2015-2-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/180074-i-m-using-a-for-loop-to-create-a-new-random-variable-x-with-range-0-1-i-want-the-forloop-to-ignor#answer_169097

在 MATLAB Online 中打开

Why not just only generate random values that meet your criteria? Here's an example:

% Solve for x
syms x ro y 
X = x*ro+y*sqrt(1-ro^2);
xx = solve(X,x)
% Plug in values
yval = -0.1;
roval = 0.3;
xmax = double(subs(xx,{y, ro},{yval,roval}));
% Generate only random numbers that meet criteria
xrand = rand(100,1)./xmax;
% Create X
Xfinal = xrand.*roval+yval*sqrt(1-roval^2);
% Check all X < 1
assert(all(Xfinal < 1))

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Fredrik Slungård 2015-2-25

Thank you for quick answer! As this is only a small part of a greater script, I think it would be better suited using a forloop in this particular case.

请先登录，再进行评论。

Answer 4

Roger Stafford 2015-2-25

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/180074-i-m-using-a-for-loop-to-create-a-new-random-variable-x-with-range-0-1-i-want-the-forloop-to-ignor#answer_169210

编辑：Sean de Wolski 2015-2-25

在 MATLAB Online 中打开

I am guessing that what you actually want is the following. It will give you 200 values in X which are all less than or equal to 1.

   X = zeros(1,n);
   x = zeros(1,n);
   for ii = 1:n
     b = true;
     while b
       t1 = rand;
       t2 = t1*ro+rand*sqrt(1-ro^2);
       b = t2 > 1;
     end
     X(ii) = t2;
     x(ii) = t1;
   end
   plot(x,X,'+')