随风摇曳的她——美蕨（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:REPSoldleft = newleft;newleft = -wave+(rand-.5)/20;right = +wave-(rand-.5)/20;delta = (right-newleft)/20;t1 = 0; t2 = 0; % costsfor w = [oldleft:delta:right, right:-delta:newleft]tstart = tic; % start time computex = .5; y = .5;for j = 1:nxs(j) = x;ys(j) = y;r = rand;if r < p1t = A11*x + w*y; % rotate and scaley = -w*x + A22*y + b1;x = t;elseif r < p2b1 = 1.5;t = B11*x + B12*y;y = B21*x + B22*y + b2;x = t;elseif r < p3t = C11*x + C12*y;y = C21*x + C22*y + b3;x = t;else % the stemy = D22*y;x = 0;endendt1 = t1 + toc(tstart); % end time computetstart = tic(); % start time plotset(h,'xdata', xs, 'ydata', ys);drawnow;pause(rand/100); % slow to human speedt2 = t2 + toc(tstart); % end time plotend%fprintf('rep=%d compute=%.1g secs plot=%.1g secs\n', reps, t1, t2);
end