專欄轉(zhuǎn)視頻-北太天元學(xué)習(xí)21-隨機(jī)模型的最后一個(gè)圖用的腳本

%北太天元模擬隨機(jī)游走

clc;

clear all;

clf;

close all;

load_plugin("time") %

N = 100; % 每一次出發(fā)后走的總步數(shù)，每一步x方向和y方向上走+1或者-1的距離

總次數(shù) = 200; % 走的總次數(shù)

x_t(1) = 0;??%初始位置的x坐標(biāo)

y_t(1) = 0;?%初始位置的y坐標(biāo)

for m=1:總次數(shù)

???%plot(-25:25,-25:25)

???%hold on

?for n = 1:N %?隨機(jī)確定每一步走到哪個(gè)位置

????% rand生成(0,1)之間的均勻分布, rand-0.5 的符號(hào)為正負(fù)的概率都是0.5

????%根據(jù) rand-0.5的符號(hào)來(lái)確定是走+1或者-1,

????% rand-0.5 的符號(hào)可能取為0，這個(gè)概率是0，可以忽略。

???sx = sign(rand-0.5);

???x_t(n+1) = x_t(n) + sx;

???sy = sign(rand-0.5); % 和 sx 一樣

???y_t(n+1) = y_t(n) + sy;

???%plot(x_t(n:n+1), y_t(n:n+1),'LineWidth', 10);

???%str = sprintf("第%d次出發(fā), 第%d步",m,n)

???%title(str)

???%pause(0.2)

?end

?%???clf;

?%???plot(-25:25,-25:25)

?%???hold on

?plot(x_t, y_t,'LineWidth',10);

??% pause(1)

?%hold off

?hold on

end