【端午礼物】随风摇曳的她——美蕨（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

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