#|她很焦虑，是时候送一波更高级的玫瑰(Python&Matlab实现) python|matlab|开发语言|女朋友的

目录
1 为什么她不快乐了
2 生活成本提高了
3 最喜欢的一句话
4 Python实现
【#|她很焦虑，是时候送一波更高级的玫瑰(Python&Matlab实现)】5 Matlab实现

1 为什么她不快乐了

昨天她给我说自己很焦虑，玩着心里不踏实，我打算先从哲学的角度分析这个问题，然后然后给点实际的，三维的立体玫瑰花，我分别用Python和Matlab实现。
《华商报》统计，说现在把一个孩子从幼儿园养活到大学毕业，研究生毕业，你没有上百万元是很难实现的。所以我就很奇怪了，你说你今天物质生活丰富了，怎么叫丰富了?请注意，这不叫物质生活丰富了，这叫生存成本提高了。概念别搞错。

2 生活成本提高了

你今天的生存安全感还在一路丧失。玩着心里也不踏实。我举一个例子，古人是物物交换，早期从来没有通胀，到后来重金属作为货币，即使在欠收年间粮食也不过涨价两三倍。可我们今天使用的钱是什么?纸币。是假钱，叫钞。这种东西最大的特点是什么?通货膨胀。它在正常年份每年通货膨胀2%到3%，也就是每十五年购买力下降一半，更重要的是它动辄发生恶性通货膨胀。

1918年第一次世界大战结束德国战败，在魏玛共和国期间，在短短两三年时间里，德国马克通货膨胀八万倍。1948 年到1949年，民国政府行将崩溃，它的法币和金圆券在一年时间里通货膨胀三万倍，直到1991年前苏联解体，苏联卢布一夕之间通货膨胀一万两千倍，原来一卢布换三美元，短短几天时间变成四千卢布换一美元，我们就拿通货膨胀一万倍来计算，它是个什么含义：一个千万富翁一夜之间变成不到一千块钱的穷光蛋。
这就是为什么今天即使是百万富翁、千万富翁、亿万富翁，他都没有安全感，他得拿自己的钱去干什么?去投资，投到汇率、房产、做市、期货。这些东西是什么?是一个每分每秒都在不断震荡的曲线。
请大家想想，你把你的身心抛洒在这条每分每秒都让你心惊肉跳的曲线上，你还有什么幸福可言?这就是为什么古人从来用不着编制幸福指数，人类一定是缺什么才喊什么。你今天是因为丢失了那个幸福感，处在极度的焦虑之中和不安定之中，因此你才得编制幸福指数。

是因为文明已经长成一棵枝繁叶茂的大树。今天的人对于文明进程的生存压力，觉得是物质生活极大的丰富了，其实是生存成本提高了；我们可以说古人有电、有手机、有空调、有飞机吗......?
可是大家想想，我们今天能缺少这些东西吗？突然哪一天停电了，估计你的生活都会错乱。我们总说古人不开心不快乐，我们这叫情景评价，是因为你站在今天文明进程的角度上去评价他们。评价他们是否快乐，我们应该回到他那个时代：那时候的人，过得很安详，没有压力的悠闲着，过年一直要到腊月十五，今天的上班族，大年初二就要匆忙地去为生活奔波了。然后你又会问：古代战乱的人民呢？其实大家得知道，中国古代每两三百年才会有一次战乱，我们今天的新中国，也才73年，放在古代也正是一个朝代最繁荣的时候。所以我们看待一件事，要换位思考，不能只以自我为中心。这一系列的例子数不胜数，这里先不展开了，后面的文章会有阐述。
大家对于文明进程的理解。相当于盲人摸象、一叶障目。什么意思?你指着文明大树的一枝一叶说这就叫文明社会全貌，会对文明做出一些误判。其实文明刚发生的时候，像一棵文明幼苗。被老子看得很清楚，在他的《道德经》中有了深刻的认识。网上和书上有很多的阐释，我有自己的理解：《道德经》开头带来的感悟
后面有时间我会写更多自己关于《道德经》的理解。
所以我再强调一遍，这叫生活成本提高了，不叫物质生活丰富了。

3 最喜欢的一句话

我一直坚信勿忘初心，方得始终。只有自己能力提高了，面对文明进程中生活成本的提高才能从容应对。
我现在就想带着她一起慢慢变成优秀的大人。祝我俩未来的命运的旖旎风景绚丽多姿吧！

4 Python实现

文章图片

import numpy as np import matplotlib.pyplot as plt from mpl_toolkits.mplot3d import Axes3D fig = plt.figure(figsize=(8,8)) ax = fig.gca(projection='3d') [x, t] = np.meshgrid(np.array(range(25)) / 24.0, np.arange(0, 575.5, 0.5) / 575 * 30 * np.pi - 4*np.pi) p = (np.pi / 2) * np.exp(-t / (8 * np.pi)) change = np.sin(20*t)/50 u = 1 - (1 - np.mod(3.3 * t, 2 * np.pi) / np.pi) ** 4 / 2 + change y = 2 * (x ** 2 - x) ** 2 * np.sin(p) r = u * (x * np.sin(p) + y * np.cos(p)) * 1.5 h = u * (x * np.cos(p) - y * np.sin(p)) c= plt.get_cmap('magma') surf = ax.plot_surface(r * np.cos(t), r * np.sin(t), h, rstride=1, cstride=1, cmap= c, linewidth=0, antialiased=True)ax.set_xticks([]) ax.set_yticks([]) ax.set_zticks([])from matplotlib.font_manager import FontProperties font_set = FontProperties(fname=r"C:\Windows\Fonts\STFANGSO.TTF",size=20) plt.title('送给LZ！', fontproperties=font_set) plt.show()

5 Matlab实现

文章图片

function drawrose grid on [x,t]=meshgrid((0:24)./24,(0:0.5:575)./575.*20.*pi+4*pi); p=(pi/2)*exp(-t./(8*pi)); change=sin(15*t)/150; u=1-(1-mod(3.6*t,2*pi)./pi).^4./2+change; y=2*(x.^2-x).^2.*sin(p); r=u.*(x.*sin(p)+y.*cos(p)); h=u.*(x.*cos(p)-y.*sin(p)); map=[0.91760.94121.0000 0.90940.93410.9992 0.90110.92710.9984 0.89290.92000.9977 0.88470.91300.9969 0.87640.90590.9961 0.86820.89880.9953 0.86000.89180.9945 0.85180.88470.9938 0.84350.87770.9930 0.83530.87060.9922 0.83370.86980.9910 0.83220.86900.9898 0.83060.86820.9887 0.82900.86740.9875 0.82750.86670.9863 0.82590.86590.9851 0.82430.86510.9839 0.82270.86430.9828 0.82120.86350.9816 0.81960.86270.9804 0.80780.85210.9765 0.79610.84150.9726 0.78430.83100.9686 0.77260.82040.9647 0.76080.80980.9608 0.74900.79920.9569 0.73730.78860.9530 0.72550.77810.9490 0.71380.76750.9451 0.70200.75690.9412 0.68360.74000.9396 0.66510.72320.9381 0.64670.70630.9365 0.62820.68940.9349 0.60980.67250.9334 0.59140.65570.9318 0.57290.63880.9302 0.55450.62190.9286 0.53600.60510.9271 0.51760.58820.9255 0.50270.57760.9271 0.48780.56700.9286 0.47290.55650.9302 0.45800.54590.9318 0.44310.53530.9334 0.42820.52470.9349 0.41330.51410.9365 0.39840.50360.9381 0.38350.49300.9396 0.36860.48240.9412 0.36230.47420.9404 0.35610.46590.9396 0.34980.45770.9388 0.34350.44940.9380 0.33720.44120.9373 0.33100.43300.9365 0.32470.42470.9357 0.31840.41650.9349 0.31220.40820.9341 0.30590.40000.9333 0.29810.39180.9235 0.29020.38350.9137 0.28240.37530.9039 0.27450.36700.8941 0.26670.35880.8843 0.25890.35060.8745 0.25100.34230.8647 0.24320.33410.8549 0.23530.32580.8451 0.22750.31760.8353 0.21690.30860.8165 0.20630.29960.7977 0.19570.29060.7788 0.18510.28160.7600 0.17460.27260.7412 0.16400.26350.7224 0.15340.25450.7036 0.14280.24550.6847 0.13220.23650.6659 0.12160.22750.6471]; set(gca,'CameraPosition',[2 2 2]) hold on surface(r.*cos(t),r.*sin(t),h,'EdgeAlpha',0.1,... 'EdgeColor',[0 0 0],'FaceColor','interp') colormap(map)end