随風搖曳的她——美蕨（matlab實作）

目錄

​​1 随風搖曳的她​​

​​2  搖曳帶來的哲思​​

​​3 Matlab代碼實作​​

1 随風搖曳的她

夢幻的場景、浪漫的氣息，帶上心愛的人，擁抱在這片花海之下，便有了電影男女主角的氛圍感；

就算閱盡了世間風貌，也抵不上和她在一起時錦短情長，一束束光交織纏繞來來往往皆是有緣之人；

欠她的浪漫，用往後餘生一點點彌補，蕨很美，想極了随風搖曳的她。

2  搖曳帶來的哲思

最早的古文明古埃及，它今天發展狀況紊亂，社會發展頭緒不清，為什麼?對自己的傳統文化沒有清理。中國人由于上個世紀初葉，不斷地檢討和清理自己的傳統文化，進而在某種程度上掃除了社會轉型的思想文化障礙，由此構成中國現代改革開放的突進，這都是大家要注意了解的。是以當我們講傳統文化價值極高的時候，請大家記住，它在人類社會演進發展上是一個時間範疇,而不是一個永恒價值。

同時我們永遠不能抛棄它的基底層參考系。大家知道我們近代讨論中國傳統文化不斷地用兩個詞彙，叫糟粕和精華。可沒有哪個學者能說清，中國傳統文化哪個部分是精華?哪個東西是糟粕?

3 Matlab代碼實作

function wulizi(n, REPS)
%% ====歡迎關注公衆号：電力系統與算法之美=======
if nargin < 1, n = 50000; end
if n <= 2000, dotsize = 6; else dotsize = 1; end

%% 使曲率矩陣變平，以使JIT可用
p1  =  .85;  p2  =  .92;  p3  = .99;

A11 =  .85;  A12 =  .04;              A22 = A11;
B11 =  .20;  B12 = -.26;  B21 = .23;  B22 = .22;
C11 = -.15;  C12 =  .28;  C21 = .26;  C22 = .24;
D22 =  .16;

b1 = 1.4;  b2 = 1.6; b3 = 0.44;

wave    = A12;
newleft = wave;
xs = zeros(1,n);
ys = zeros(1,n);

%% 可視化
figure(gcf)                                       % on top
set(gcf, 'color','white','menubar','none', 'numbertitle','off',...
    'name', '随風搖曳的清香', 'doublebuff', 'on');
darkgreen = [0 2/3 0];                            % fern-ish
h = plot(xs, ys, '.','markersize',dotsize,'color',darkgreen);
axis([-4 4 0 10])                                 % set size
axis off                                          % make pretty

%% 開始示範
if nargin < 2; REPS = 10; end
for rp = 1:REPS
    oldleft = newleft;
    newleft = -wave+(rand-.5)/20;
    right   = +wave-(rand-.5)/20;
    delta   =  (right-newleft)/20;
    t1 = 0; t2 = 0;                                 % costs
    for w = [oldleft:delta:right, right:-delta:newleft]
        tstart = tic;                                 % start time compute
        x  = .5;  y = .5;
        for j = 1:n
            xs(j) = x;
            ys(j) = y;
            r = rand;
            if r < p1
                t = A11*x + w*y;                          % rotate and scale
                y = -w*x + A22*y + b1;
                x = t;
            elseif r < p2
                b1 = 1.5;
                t = B11*x + B12*y;
                y = B21*x + B22*y + b2;
                x = t;
            elseif r < p3
                t = C11*x + C12*y;
                y = C21*x + C22*y + b3;
                x = t;
            else                                        % the stem
                y = D22*y;
                x = 0;
            end
        end
        t1 = t1 + toc(tstart);                        % end time compute
        tstart = tic();                               % start time plot
        set(h,'xdata', xs, 'ydata', ys);
        drawnow;
        pause(rand/100);                              % slow to human speed
        t2 = t2 + toc(tstart);                        % end time plot
    end
    %fprintf('rep=%d compute=%.1g secs plot=%.1g secs\n', reps, t1, t2);
end