目录
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
