Help on a problem

4 次查看（过去 30 天）
Gaëtan Poirier
Gaëtan Poirier 2017-9-22
评论： Walter Roberson 2017-9-22
I have this problem and I'm wondering if anyone can help me to figure it out. It goes as follows: write a program to iteratively generate points in two-dimensional space using the following rules:
(x_{n+1},y_{n+1}) =
/ (0.5,0.27*y_n), with 2% probability;
|
| (-0.139*x_n + 0.263*y_n + 0.57, 0.246*x_n + 0.224*y_n - 0.036), with 15% probability;
<
| (0.17*x_n - 0.215*y_n + 0.408, 0.222*x_n + 0.176*y_n + 0.0893), with 13% probability;
|
\ (0.781*x_n + 0.034*y_n + 0.1075, -0.032*x_n + 0.739*y_n + 0.27), with 70% probability.
Start from an initial point (x_1,y_1)=(0.5,0.0). Carry out the iteration at least 30,000 times

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

采纳的回答

Walter Roberson
Walter Roberson 2017-9-22
for iter = 1: 30000
r = rand() ;
if r < 0.02
....
elseif r < 0.02+0.15
.....
  3 个评论
显示 1更早的评论
Gaëtan Poirier
Gaëtan Poirier 2017-9-22
This is what I have so far, doesn't seem repeat 30000 times, only comes out with one varible in r, x, and y.
for iter = 1: 30000
r = rand();
x = 0.5;
y = 0;
if r < 0.02
y = y * 0.27*r;
elseif r < 0.02 + 0.15
x = x -0.139 * r + 0.246 * r;
y = y + 0.263 * r + 0.224 * r - 0.036;
elseif r < 0.02 + 0.15 + 0.13
x = x + 0.17 * r - 0.032 * r;
y = y + 0.034 * r + 0.739 * r + 0.27;
elseif r < 0.02 + 0.15 + 0.13 + 0.70
x = x + 0.781 * r - 0.032 * r;
y = y + 0.034 * r + 0.1075 + 0.739 * r + 0.27;
end
end
Walter Roberson
Walter Roberson 2017-9-22
x_ = 0.5;
y_ = 0;
for n = 1: 30000
r = rand();
if r < 0.02
y_(n+1) = y_(n) * 0.27*r;
elseif r < 0.02 + 0.15
x_(n+1) = x_(n) -0.139 * r + 0.246 * r;
y_(n+1) = y_(n) + 0.263 * r + 0.224 * r - 0.036;
etc
After all of that you are probably expected to do
scatter(x, y)

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Get Started with MATLAB 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by