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.

2015 2 24
4 个回答
回答已采纳
更新时间：2015 2 25
15 次查看（30 天）

0 个投票

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,'+')

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

 采纳的回答

Niels
Niels 2015-2-25
编辑：Niels 2015-2-25

0 个投票

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,'+')

更多回答（3 个）

Joseph Cheng
Joseph Cheng 2015-2-24

0 个投票

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
Fredrik Slungård 2015-2-25
编辑：Fredrik Slungård 2015-2-25
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,'+')

请先登录，再进行评论。

Sean de Wolski
Sean de Wolski 2015-2-24

0 个投票

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
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.

请先登录，再进行评论。

Roger Stafford
Roger Stafford 2015-2-25
编辑：Sean de Wolski 2015-2-25

0 个投票

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,'+')

类别

在 帮助中心 和 File Exchange 中查找有关 Get Started with Signal Processing Toolbox 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by