【寒月、寒雪、凄景】（matlab实现）

1 水云间，俏无言

《江城子·斜风细雨作春寒》

斜风细雨作春寒。

对尊前，忆前欢，

曾把梨花，寂寞泪阑干。

芳草断烟南浦路，

和别泪，看青山。

昨宵结得梦夤缘。

水云间，俏无言，

争奈醒来，愁恨又依然。

展转衾裯空懊恼，

天易见，见伊难。

时当春日，词人独对孤樽，或许欲以解闷而已。不想这“斜风细雨”、这料峭“春寒”，却勾起了她对许多“前欢”往事的回顾。

这首词写失恋的悲愁，充满作者心灵深处的凄厉哀鸣。朱淑真在少女时期曾有过一段自由婚恋的幸福，可是后来由父母主婚，强嫁一俗吏，志趣难合，遂愤然离去。这棒打鸳鸯散的忧伤，这琼枝错插、忍遭摧损的隐痛，萦盘郁结于心，使她在恨、愁、悲、病、酒五字生涯中凄然以终。

2 心里跌宕的芥蒂

老子是中国思想文化之元祖。这句话什么意思呢?中国先秦诸子百家，其实不仅包括文人、学者、思想家，而且包括政治家、军事家、外交家当时也叫纵横家等等。从时间排序上讲，诸子百家第一子很有可能是春秋五霸之首齐桓公的宰相管仲。

我们之所以说老子是中国思想文化之元祖，是因为老子作为文人、思想家，是中国先秦诸子百家第一子，《老子》一书是中国古代典籍中最深奥最难读的一本书。老子思想高深，相较于其后中国所有的文人学者，他所达到的高度，后人不可企及。

中国文化思想史有一个很奇怪的现象。东西方文化总是相反的，它们不仅在内容上形式上经常相反，它们甚至连发展方式，运行脉络都大相径庭。西方文化从简单到复杂，从低级到高级，历经神学阶段、哲学阶段、科学阶段，但中国文化反着来，它是由高级向低级下流。

中国的整个国学思想奠基，全部发生在先秦时代，也就是秦始皇统一中国之前，此后中国在思想文化的发展上，从此几乎没有任何重大建树，只有一件事情例外，那就是东汉中期，印度的佛教传入中国。即使放在先秦时代看，在先秦时代，思想最高点仍然在先秦思想第一元祖老子那去。所以中国文化发展都跟西方的思想发展相反，不是从低级向高级发展，反而是从高级向低级下流。

3 Matlab实现凄凉的场景

function QiLiang

%% 欢迎关注公众号：电力系统与算法之美

MainFig=figure('units','pixels','position',[300 80 500 500],...
    'Numbertitle','off','menubar','none','resize','off',...
    'name','凄凉的一个夜晚');
axes('parent',MainFig,'position',[0 0 1 1],...
    'XLim', [0 500],...
    'YLim', [0 500],...
    'NextPlot','add',...
    'layer','bottom',...
    'Visible','on',...
    'Color',[0 59 129]./255,...
    'XTick',[], ...
    'YTick',[]);
hold on
%%
layerBEPos=[150;140;100;90];
layerColor=[146 210 245;101 183 231;0 101 181;0 59 129];
excursion=15;
interval=10;
xSep=25;
layerPos=zeros(size(layerBEPos,1),length(xSep:xSep:500));
layerPos(:,1)=layerBEPos(:,1);
for i=1:size(layerBEPos,1)
    for j=xSep:xSep:500
        listPos=round(j./xSep+1);
        tempRandi=randi([-excursion,excursion]);
        yPos=tempRandi+layerPos(i,listPos-1);
        if i>1&&yPos>=layerPos(i-1,listPos)-5
            yPos=layerPos(i-1,listPos)-interval;
        end
        yPos(yPos<0)=0;
        layerPos(i,listPos)=yPos;
    end   
end
for i=1:size(layerBEPos,1)
    XData=0:xSep:500;
    YData=layerPos(i,:);
    Yq=interp1(XData,YData,0:1:500,'spline');
    Xq=[0,0:1:500,500];
    Yq=[0 Yq 0];
    fill(Xq,Yq,layerColor(i,:)./255,'EdgeColor','none')
end
%%
XRandiS=randi([50 450],[1,15]);
YRandiS=randi([300 460],[1,15]);
scatter(XRandiS,YRandiS,10,'o','filled','CData',[252 241 0]./255)
XRandiB=randi([50 450],[1,15]);
YRandiB=randi([300 460],[1,15]);
scatter(XRandiB,YRandiB,18,'o','filled','CData',[252 241 0]./255)
%%
snowXpos=[randi([0,30],[1,6]),randi([470,500],[1,6]),randi([0,500],[1,12])];
snowYpos=[randi([220,470],[1,6]),randi([220,470],[1,6]),randi([470,500],[1,12])];
snowSize=randi([25,50],[1,24]);
snowAngle=2*pi*rand([1,24]);
snowWidth=1+0.6*rand([1,24]);
for i=1:24
    drawSnow(snowXpos(i),snowYpos(i),snowSize(i),snowAngle(i),snowWidth(i))
end
%%
moonPos=[320 300]+[randi([0 50]),randi([0 80])];
moonSize=randi([28 32]);
blankPos=randi([floor(0.3*moonSize) ceil(0.5*moonSize)],[1,2]);
drawMoon(moonPos(1),moonPos(2),moonSize,moonPos(1)-blankPos(1),moonPos(2)+blankPos(2))
%%
    function drawMoon(x,y,R,bx,by)
        t=0:pi/50:2*pi;
        X=x+cos(t).*R;
        Y=y+sin(t).*R;
        BX=bx+cos(t).*R;
        BY=by+sin(t).*R;
        fill(X,Y,[255 251 219]./255,'EdgeColor','none')
        fill(BX,BY,[0 59 129]./255,'EdgeColor','none')
    end
    function drawSnow(x,y,len,angle,width)
        for theta=0:pi/3:2*pi-pi/3
            xTail=cos(theta+angle)*len+x;
            yTail=sin(theta+angle)*len+y;
            plot([x,xTail],[y,yTail],'color',[0 134 207]./255,'lineWidth',width)
            for branchRatio=[0.4 0.55 0.7]
                BX=x+cos(theta+angle)*len*branchRatio;
                BY=y+sin(theta+angle)*len*branchRatio;
                LX=BX+cos(theta+angle+pi/3).*len.*0.4;
                LY=BY+sin(theta+angle+pi/3).*len.*0.4;
                RX=BX+cos(theta+angle-pi/3).*len.*0.4;
                RY=BY+sin(theta+angle-pi/3).*len.*0.4;
                plot([BX,LX],[BY,LY],'color',[0 134 207]./255,'lineWidth',width.*0.8)
                plot([BX,RX],[BY,RY],'color',[0 134 207]./255,'lineWidth',width.*0.8)
            end
        end
    end
end