m0_51153866的博客

原创 xgrads库在ctl转nc中存在的不足及修改

xgrads库在ctl转nc中存在的不足及修改1、问题描述ctl文件dset C:\Users\马冠龙\Desktop\python\EOF\sstpx.grdtitle Coads SSTAundef -999.0xdef 18 linear 120 10ydef 12 linear -27.5 5zdef 1 linear 1000 1tdef 516 linear 1jan1948 1monthvars 1S 0 99 Coads SST anomaly interper

原创 第十三章 Emacs编辑器

第十三章 Emacs编辑器Emacs界面包括标题栏、菜单栏、工具栏、窗口区域、模式行、回显行模式行介绍-:--- buf pos L#(Pattern)冒号后面前两个字符表示缓冲区文本状态“–”表示缓冲区未被编辑“**”表示缓冲区文本已被修改“%%”表示只读缓冲区文本未被修改“%*”表示只读缓冲区文本已被修改buf文件名pos显示文本范围ALL显示了整个文件Top显示了开始部分Bot显示了末尾部分NN%窗口顶部以上部分显示缓冲区文件的百分比L#表示字母

原创 第十二章 VI编辑器

第十二章 VI编辑器一半编辑器有两种类型：行编辑器、全屏编辑器VI编辑器三种模式：命令模式、文本编辑模式、末行模式模式切换，任何模式下按Esc进入命令模式；命令模式按i、a、o插入命令进入文本编辑器模式；命令模式按：进入末行模式VI常用命令进入VI命令含义vi进入VIvi 文件名打开文件；如果文件不存在，则创建文件vi+n 文件名进入VI，光标停留在第n行开始处vi+ 文件名进入VI，光标停在最后一行开始处vi+/字符串 文件

原创 天气学诊断实习六 计算水汽流函数和水汽势函数

一、实习目的:熟悉水汽平流、水汽通量、水汽通量散度的实际编程计算。二、实习内容:编制计算水汽平流、水汽通量、水汽通量散度的程序，并且绘制出两个时次25日20时，26日20时的水汽平流、水汽通量、水汽通量散度分布（850hPa,500hPa)三、代码实现# -*- coding: utf-8 -*-import numpy as npimport pandas as pdimport mathimport matplotlib.pyplot as pltimport cartopy.crs

原创 第十一章 系统服务和日志

第十一章 系统服务和日志Linux系统服务也称守护程序，在Linux系统启动时自动加载，在Linux退出时自动停止linux系统常见服务及作用服务名称功能简介建议acpid电源管理接口。如果是笔记本电脑用户，则建议开启，可以监听内核层的相关电源事件开启anacron系统的定时任务程序。是 cron 的一个子系统，如果定时任务错过了执行时间，则可以通过 anacron 继续唤醒执行关闭alsasoundalsa 声卡驱动。如果使用 alsa 声卡，则开

原创 第十章 网络管理

第十章 网络管理联网操作是通过机器上的PCI设备，即网络接口卡（NIC）来实现的可以使用lspci命令查看PCI设备是否被检测Linux内核不允许用户将NIC组为文件进行访问，在/dev文件下没有NIC的设备节点，Linux 通过网络接口访问NIC，对于每一个识别出的NIC，内核都产生一个网络接口,常用的Linux接口名称名称类型eth0以太网lo虚拟回环设备ppp0使用PPP协议的串口设备tr0令牌环fddi0光纤用ifconf

原创 第九章 设备管理

第九章 设备管理Linux系统中，用户是通过文件系统与设备来连接的所有设备都被当做特别文件，每个设备都对应文件系统中的一个索引节点，都有一个文件名；每个设备都是在驱动控制下允许的，驱动程序则与/dev目录下的特殊文件联系在一起。Linux采用设备文件统一管理硬件设备。Linux系统中常见的的硬件设备名称如下：硬件设备文件名称IDE设备/dev/hd[a~d]SCSI/SATA/USB盘/dev/sd[a~p]软驱/dev/fd[0~1]打印机

原创 第八章 储存管理

第八章 储存管理机械硬盘主要由磁盘盘片、磁头、主轴、转动轴等组成；逻辑结构主要为磁道、扇区、柱面硬盘大小：磁头数×柱面数×扇区数×每个扇区的大小常见的机械硬盘接口：IDE、SATA、SCSI硬盘识别MBR主引导记录，占据0磁道的第一个扇区，在一个512字节的MBR中boot loader引导加载程序：前446个字节partition table分区表：16×4个字节，主要存储分区号、分区的起始柱面、分区柱面数量结束标志：最后2字节常见的分区号：0x5或0xf：可拓展分区0

原创 第七章 进程管理

第七章 进程管理进程是是正在执行的一个程序或命令，每个进程都是一个运行的实体，都有自己的地址空间，并占用一定的系统资源。程序是使用计算机语言编写可以实现特定目标或解决特定问题的代码集合。进程管理作用：判断服务器健康状态、查看系统中所有进程、“杀死”进行进程监控命令ps命令查看系统中所有进程的详细信息基本格式：ps [选项]-a：显示一个终端中的所用进程-u：显示进程的归属用户及内存使用情况-x：显示没有控制终端的进程-l：以长格式显示更加详细的信息-e：显示所有进

原创 第六章 文件管理

第六章 文件管理文件：文件系统中储存数据的一个命名对象，它是用户眼中储存文件的最小单元。目录：包含许多文件控制块项目的一类特殊文件，文件系统中每个文件都登记在一个或多个目录中子目录：背包含在另一目录中的目录，除根目录外所有目录都是子目录，有父目录存在，目录的主要作用是和文件名构成路径。Linux系统中所有文件都被组织成以根目录开始的倒置树状结构，文件系统按名访问绝对路径：从根目录写起的文件或目录名相对路径：相对于当前路径的文件或目录名文件命名：除了字符“/”之外，所有字符都

原创 天气学诊断实习五 计算水汽平流、水汽通量、水汽通量散度

一、实习目的:熟悉水汽平流、水汽通量、水汽通量散度的实际编程计算。二、实习内容:编制计算水汽平流、水汽通量、水汽通量散度的程序，并且绘制出两个时次25日20时，26日20时的水汽平流、水汽通量、水汽通量散度分布（850hPa,500hPa)三、代码实现# -*- coding: utf-8 -*-import numpy as npimport pandas as pdimport mathimport matplotlib.pyplot as pltimport cartopy.crs

原创 天气学诊断实习四 计算垂直速度

一、实习目的:熟悉垂直差分在气象中的应用，掌握垂直速度的实际编程计算。二、实习内容:编制计算垂直速度程序，并绘制500hPa垂直速度。用第二种修正方案，其中大气层顶的垂直速度可以直接采用0，也可以用绝热法。并且绘制出两个时次25日20时，26日20时的修正后的垂直速度分布（850hPa,500hPa)三、算法原理:垂直涡度计算：修正方案二：四、代码实现# -*- coding: utf-8 -*-import numpy as npimport pandas as pdimport

原创 天气学诊断实习三 计算涡度、散度、涡度平流和温度平流

一、实习目的:熟悉使用差分方法解决方程中的导数、偏导数，计算涡度、散度、涡度平流和温度平流的程序二、实习内容:编制计算涡度、散度、涡度平流和温度平流的程序，并且绘制出两个时次25日20时，26日20时的涡度、散度、涡度平流和温度平流分布（850hPa,500hPa)三、算法原理:涡度计算：散度计算：平流计算：差分方法：中间使用中央差分、边缘部分使用向前差分、向后差分四、代码实现:# -*- coding: utf-8 -*-import numpy as npimport pa

原创 气象数据处理

气象数据处理一、获取气象数据网站NOAA官网站点数据：https://www.ncei.noaa.gov/data/global-hourly卫星格点数据：https://psl.noaa.gov/data/gridded/index.html二、批量文件下载站点数据下载：for stid in `cat china_stid`do for i in `seq 1961 2021` do wget https://www.ncei.noaa.gov/data/

原创 天气学诊断实习二 位温和相当位温的计算

一、实习内容  编制计算位温和相当位温的程序，并且绘制出两个时次25日20时，26日20时的位温和相当位温分布（850hPa,500hPa)二、代码实现import pandas as pdimport numpy as npimport matplotlib.pyplot as pltimport cartopy.crs as ccrsimport cartopy.feature as cfeatureimport cartopy.mpl.ticker as ctickerimport

原创 天气学诊断实习一 露点和凝结抬升高度的计算

一、实习内容  已知南京站某一时次的相对湿度（92.8%）、气温（14.1ºC）、气压（1017.4, hPa）及站点高度（35.2m），请利用迭代法求出：（1） 该站的露点温度（2） 凝结高度和凝结温度二、算法原理露点温度计算：凝结高度，凝结温度计算：![在这里插入图片描述](https://img-blog.csdnimg.cn/9ec11368b40b478b9d6dca89a522a405.png?x-oss-process=image/watermark,type_ZHJvaWRzY

原创 第四章：Shell脚本编程基础

第四章：Shell脚本编程基础用户程序不能直接控制计算机硬件，只有操作系统内核才能真正控制计算机硬件shell本身是一个用C语言编写的程序,他是连接用户和linux/UNIX内核的桥梁。shell是一个交互式命令解释器，一个读入并解释用户命令的程序shell同时也是一个程序设计语言，是解释型语言，不需要编译。shell分为交互式和非交互式Linux中创建的文件默认是不可执行的，执行脚本首先要赋予文件权限shell程序第一行#！/bin/bash中#！说明该脚本是用哪

原创 第三章：Linux系统使用基础

第三章：Linux系统使用基础常见的Linux系统桌面环境：GNOME、KDE、Unity、Fluxbox[root@localhost ~]#root表示当前登录用户，localhost当前系统主机名，~用户当前所在目录，#表示超级用户家目录在/root，$表示普通用户家目录在/home/用户名。Linux命令基本格式：命令 [选项] [参数]四个简单指令：w和who、echo、data、passwdw和who查看已登录用户信息，w可以查看所有用户信息，who查看当前登录用

原创 第二章：Linux系统概述

第二章：Linux系统概述Linux是一类UNIX计算机操作系统的统称。Linux之父是Linus Torvalds内核是指一个提供设备驱动、文件系统、进程管理、网络通信等功能的系统软件。Linux内核版本主要有两种：稳定版和开发版Linux内核版本号由3组数字组成，“主版本.次版本.释出版本-修改版本”第一组数字：当前内核主版本号第二组数字：当前内核次版本号，偶数为稳定版，奇数为开发版第三组数字：当前内核更新次数和当前内核修补次数常见Linux系统发行版：Red Hat、Debi

原创 气象统计 第八次实习 EOF分析

气象统计 第八次实习一、实习题目二、实习程序这次实习本来打算先用如下命令合并CTL和GRD转化为nc再开始计算的，但是不知道什么原因，平常用的合并方法一直出错，代码如下：#ctl,grd数据转ncds = open_CtlDataset('SSTPX.ctl')ctl = CtlDescriptor(file='SSTPX.ctl')ds.attrs['pdef' ] = 'None'ds.to_netcdf('SSTPX.nc')那位大佬知道怎么错的可以给我讲讲，所以这次直接再lin

原创 气象统计 第七次实习 滑动平均累计距平

气象统计 第七次实习一、实习题目二、实习程序import numpy as npimport pandas as pd import matplotlib.pyplot as pltfrom pylab import *mpl.rcParams['font.sans-serif']=['SimHei']##滑动平均函数，也可以使用datapj= y.rolling(11).mean()def hdpj(data,k): n=int(len(data)-k) t=np.ze

原创 气象统计 第六次实习 降水线性倾向率

气象统计 第六次实习一、实习题目二、实习程序# -*- coding: utf-8 -*-import numpy as npimport pandas as pdimport matplotlib.pyplot as pltimport cartopy.crs as ccrsimport cartopy.feature as cfeatureimport cartopy.mpl.ticker as ctickerfrom cartopy.io import shapereader a

原创 气象统计 第五次实习 多元线性回归

气象统计 第五次实习一、实习题目二、实习程序import randomimport numpy as np##多元线性相关函数，这里是三个变量的，其他变量可以略加修改def dyxxxg(x1,x2,x3,y): n=len(x1) x0=np.ones((n,1)) x=np.array([x0[:,0],x1,x2,x3]).T y=y.reshape((n,1)) t1=np.linalg.inv(np.dot(x.T,x)) t2=np.

原创 气象统计 第四次实习 一元线性回归

气象统计 第四次实习一、实习题目二、实习程序import numpy as npimport pandas as pdimport randomimport matplotlib.pyplot as plt##读文件Circulation_Index=np.array(pd.read_table('Circulation_Index.txt',header=None))Temperature=np.array(pd.read_table('Temperature.txt',header=

原创 气象统计 第二次实习 简单相关系数和自相关系数

气象统计 第二次实习一、实习题目二、实习程序import numpy as npimport math import operatorimport random###自相关系数函数def zxgxs(a): n=len(a) N=int(n/2) ave=np.mean(a) r=np.zeros(int(n/2)) for i in range(0,N): sum=0 for j in range(0,n-i-1):

原创 气象统计 第三次实习 落后交叉相关系数和偏相关系数

气象统计 第三次实习一、实习题目二、实习程序import pandas as pdimport numpy as np##读文件data1=np.array(pd.read_table('mon1.txt',header=None))data2=np.array(pd.read_table('mon2.txt',header=None))data12=np.array(pd.read_table('mon12.txt',header=None))##简单相关系数,注意x,y长度相等d

原创 气象统计 第一次实习 气候场、距平场、均方差场

气象统计 第一次实习一、实习题目二、实习程序import pandas as pdimport numpy as npimport mathimport matplotlib.pyplot as pltimport cartopyimport cartopy.crs as ccrsimport cartopy.feature as cfeatureimport cartopy.mpl.ticker as ctickerfrom pylab import *

原创 全国计算机等级考试三级Linux考纲2021

原创 n皇后问题

一、问题描述  在nXn的方格棋盘上放置n个皇后要求每个皇后不同行、不同列、不同对角线，求解放置方法二、算法  利用一维数组a[n]来解决问题,建立横坐标为n,纵坐标为a[n]的棋盘，皇后位置用（i，a[i]）表示。题目要求每个皇后不同行、不同列、不同对角线.若同列则a[k]==j，若为对角线位置，则会构成等腰三角形，则|a[k]-j|==|i-k|所以以下位置不可放置皇后(a[k]==j)||(abs(a[k]-j)==abs(i-k))搜索时逐列寻找，在第i行上试探每一个列j,找到与

原创 冒泡排序算法

一、问题描述  对于一个给定的含有n个元素的数组，利用冒泡排序按元素值递增排序。二、算法  简单选择排序是将原数组分为有序区和无序区两个部分，有序区中的所有元素都不大与无序区中的元素，初始时有序区为空。经过n-1次排序，每次在无序区通过相邻元素两两比较，若相邻元素反序则两两调换，若正序则继续比较下一组相邻元素，最终将无序区中最小的元素移动到无序区开头（有序区末尾）。三、代码#include<stdio.h>//交换x、y的值 void fun1(double &x,doub

原创 简单选择排序

一、问题描述  对于一个给定的含有n个元素的数组，利用简单选择排序按元素值递增排序。二、算法  简单选择排序是将原数组分为有序区和无序区两个部分，有序区中的所有元素都不大与无序区中的元素，初始时有序区为空。经过n-1次排序，每次都将无序区中最小的元素移动到无序区开头（有序区末尾）。有序区无序区排序前a[0]、a[1]⋯a[i−1]a[0]、a[1]\cdots a[i-1]a[0]、a[1]⋯a[i−1]a[i]、a[i+1]⋯a[n−1]a[i]、a[i+1]\cdots

原创 主成分分析PAC

一、模型介绍  维度灾难，常是指在涉及到向量的计算的问题中，随着变量维数的增加，计算量呈指数倍增长的一种现象。变量过多导致原本简单的问题复杂化，甚至出现了无法解决的情况。在维度灾难的背景下，主成分分析法（Principal components analysis）孕育而生，由于各变量间存在一定的相关关系，因此有可以用较少的综合指标来反映多维数据里的信息，利用这一思路就可以通过数学手段，把多个变量化为少数变量，实现数据降维。二、符号说明符号说明XXX原始数据矩阵CCC协方差

原创 层次分析APH

这里写自定义目录标题欢迎使用Markdown编辑器新的改变功能快捷键合理的创建标题，有助于目录的生成如何改变文本的样式插入链接与图片如何插入一段漂亮的代码片生成一个适合你的列表创建一个表格设定内容居中、居左、居右SmartyPants创建一个自定义列表如何创建一个注脚注释也是必不可少的KaTeX数学公式新的甘特图功能，丰富你的文章UML 图表FLowchart流程图导出与导入导出导入欢迎使用Markdown编辑器你好！ 这是你第一次使用 Markdown编辑器 所展示的欢迎页。如果你想学习如何使用Mar

原创 趣味数学建模

以下文章转载自公众号：数学模型趣味建模用十分有趣的例子结合了背后的理论知识对生活中的问题给出了合理的答案，是非常适合闲暇之余消遣的文章。如何根据一张照片确定男朋友的心？流行天王迈克尔·杰克逊 45 度反重力前倾很难么？如果没有国家管控，北京鼠疫将如何漫延？澳大利亚山火蔓延的数学模型你突然变成一个“2B”的概率有多大“科比球”背后的数学情人节，数学教你如何科学地哄女朋友？情人节，数学告诉你什么没有女朋友?如何抢别人怀里的女孩子?疫情防控，哪些省市才是真硬核？“待我长发及腰”，女生说这

原创 mathtype导致无法粘贴解决方法

mathtype导致无法粘贴第一步关闭office第二步在MathType安装目录下点击Office Support\64找到 MathType Commands 2016.dotm复制；第三步打开C:\Program Files\Microsoft Office\root\Office16\STARTUP覆盖粘贴MathType Commands 2016.dotm重启office即可...

原创 WFL的安装

WFL的安装一、下载WFl在气象家园下载WFL库：气象家园或直接下载我上传的文件：二、安装WFI解压安装包后，运行.exe安装过程中不要修改安装路径，负责可能导致无法使用三、环境配置在vs2019中点击：工具→\rightarrow→选项→\rightarrow→打开选项卡后依次点击：Intel Compilers and Libraries→\rightarrow→IFX Intel Fortran →\rightarrow→Compilers分别Win32和Win6