儿子说了几次要爬泰山，我决定这次暑假带他爬一次。媳妇一放假就给我们买了去泰安的票。人挺多的，沈阳直达泰安的票都没了。只能到北京中转，一路急急忙忙从朝阳站到北京南站，终于赶上了去泰安的高铁。 车上我们就开始商量明早如何爬山，车厢座位对面有一对年轻情侣，他们是唐山的学生，从北京上车到泰安，准备下了车当晚夜爬泰山，儿子也想夜爬，但天气预报说明天阴雨天，同时因晚上还需要干点儿小活儿，我没同意。 下车后我们住到了具体泰安高铁站不远的快捷酒店。 第二天早起用过早餐过后，打出租车直奔红门。我们的计划是这样，走最经典的路线，全程不做缆车，从红门爬到中天门，在从中天门排到南天门，最后到玉皇顶。争取不落下每一处景点。儿子爬山的情绪非常高，看起来他比我还有劲儿。 到红门才8：09，时间挺早，但人已经不少，这时下山的人比上山的人多。我们随机采访下山的人是否看到今早的日出，他们说是，我们一方面感叹天气不准确，另一方面些许遗憾没有夜爬看日出，不过美妙的景色很快就让人留恋忘返，我们走了好久才走到验票口。我和儿子笑谈说，泰山的序幕好长，就像电影一样，看了好久才出题目，怪不得有人说泰山是免费爬的，原来他从来都没有爬到收费的地方。 红门过后，右侧水声潺潺，路上树木茂密，人在树荫下行走，所有烦恼均会忘去。 有一岔道——经石峪，岔出去还得原路返回，儿子要实践计划，不落每处美景，那就去吧。这真是个人迹罕至的好去处，从主路上叉出一条道，先下了一个大坡，然后又爬好长的一个大坡，绕过一个山岗，就到了经石峪。在一篇巨大的山石上，留着一层薄水，走进看时，石头上还刻着经文，真是巧夺天工。一回首，才发现左侧的涯壁上刻着高山流水的石刻。右侧还有石桌、石凳，桌上放着茶壶茶杯，有沏好的茶，能在这里赏景、听风、品茗真的人生乐事啊，过一会儿我才弄清楚，这是景点管理所的人员的茶壶，真是让人羡慕的工作。 儿子有点走不动了，让让他吃点儿喝点，我在经石峪下做了个倒立，绕着经石峪转了一圈，返回主道，继续爬山。返回主道的过程中，我们发现了一直爬泰山的蜗牛，我就拍了个视频，儿子配流浪地球中太空电梯升空的音乐，这只蜗牛也太励志了，哈哈。   一直爬，路遇泰山的跳山工，他们挑着重物每天上上下下10多趟，真是辛苦。因为中天门通车，看到挑山工从往下挑水，我们就知道距离中天门不远了。继续加油，马上就到了中天门。   中天门，稍作休息，给儿子买了一碗开水泡面（20元），我吃了一晚从沈阳一路背过去的泡面，加开水5元。刚吃完，漂泊大雨就来了，这雨太大了，有伞也不好使，只好躲在酒店的雨棚下面躲雨。雨稍微变小些的时候，工作人员就催促游客下山，并且说山门已经封山了。   我和儿子商议后，决定继续爬上，中天门后的一段平路，人很少，又爬了一段时间，人都不知道从哪里冒了出来，渐渐多了起来。穿着雨衣也有些热。经过步云桥，桥从瀑布前跨过，刚下过雨，瀑布颇为壮观，植物和石头都是湿的，石刻红色字迹更加清晰和醒目。   经过飞来石和五大夫松。来到升仙坊的时候，一阵凉风刮过，雾气从我们身后的山下扑掩而来，身后的山下很快就看不见了。   终于到达了南天门。在山顶上还是比较冷的，淅淅沥沥地下着小雨。我们在店里要了些吃的，稍微休息。然后，等雨稍微小些了，打算在山顶做个环游。山顶的景色是非常辽阔的，极目远眺，整个泰安尽收眼底，颇有一览众山小之感。我领着我们家的学生去拜访了孔子庙，但因为是假期，孔老师也放假了。但在庙门口的景色非常好。近景有苍松翠柏，远景也极为开阔。   游览的著名的五岳独尊处、无字碑、泰山极顶1545m。14:00我们开始下山。下山时，天已不再下雨，石头已没有上山时湿滑了。我们感觉下山时，人非常少，特别是过了中天门以后，有几段走了好久都不见一个人，好像整个景区都是为我们开的。等到出来景区，才发现是因为下雨，门口一直没有放人进山，而且，大多数游客都选择从中天坐车下山。因此，游览的体验感极佳。 […]

同样的数据，用二次函数拟合，Matlab的结果和R的回归参数为什么不一样，但R2都是一样的？   %% 河北粮食主产区水分生产率和灌溉量二次函数关系x=H.I_mm;y=H.WUE_2;% 二次函数拟合p = polyfit(x, y, 2); % 二次拟合，p(1)是最高项系数，p(2)是一阶项系数，p(3)是直流项（即y轴的截距）% 计算拟合值x_nihe=[min(x):0.1:max(x)]’yfit = polyval(p,x_nihe );% 绘制原始数据、拟合线和拟合公式figure;plot(x, y, ‘o’, x_nihe, yfit, […]

% 张旭东更新20240724 % DLc——平均日照时数[hour] % Rn——冠层表面净辐射[MJ/m^2/day] % Rns——净短波辐射[MJ/m^2/day] % Rnl——净长波辐射[MJ/m^2/day] % Rs——入射的太阳辐射[MJ/m^2/day] % Ra24——碧空太阳总辐射[MJ/m^2/day] % latitude——地理纬度转换为[度] % elevation——海拔高度，[m] % RHmean——平均相对湿度，% […]