首页 > 其他分享 >2024-3-13 气象中的常用代码 154828

2024-3-13 气象中的常用代码 154828

时间:2024-10-03 21:24:27浏览次数:8  
标签:index 13 plt 2024 Tm time date 154828 unique

目录

画图

横坐标添加月份

Python

import matplotlib.pyplot as plt
import pandas as pd
import numpy as np

# 准备时间和温度数据
start_date = pd.to_datetime('1996-12-01') # the next date
end_date = pd.to_datetime('1998-12-01') # the current data, so 36
time_index = pd.date_range(start_date, end_date, freq='M')
temperature_data = np.random.uniform(0, 40, len(time_index))  # 这里使用随机温度数据,您应该替换为您的实际数据
print("shape:",temperature_data.shape)
# 创建时间序列图
plt.figure(figsize=(12, 6))
plt.plot(time_index, temperature_data, marker='o', linestyle='-')
plt.title('Temperature-Time-Series')
plt.xlabel('Time')
plt.ylabel('Temperature')
plt.xticks(rotation=45)  # 旋转x轴标签以包含月份
plt.grid(True)
plt.show()

Matlab

startYear = 2023;startMonth=1;endYear=2023;endMonth=12;
time_index = datetime(startYear,startMonth, 1):calmonths(1):datetime(endYear,endMonth, 1);
plot(time_index,rand(length(time_index),1));

空间映射

标准化

To = 1:10;        % 实际观测数据
Tm = 30:-1:21;    % 模型数据
% 标准化To和Tm序列
standardized_To = (To - mean(To)) / std(To);
standardized_Tm = (Tm - mean(Tm)) / std(Tm);

% 获取原始To中最后一个点在标准化后的To中的位置
relative_position = (To(5) - mean(To)) / std(To);

% 根据相对位置在标准化后的Tm中找到对应的值
% 使用unique确保样本点的唯一性
[unique_standardized_Tm, unique_idx] = unique(standardized_Tm);
corresponding_Tm_value = interp1(unique_standardized_Tm, Tm(unique_idx), relative_position);

disp('To序列:');
disp(To);
disp('对应的Tm序列:');
disp(corresponding_Tm_value);

标签:index,13,plt,2024,Tm,time,date,154828,unique
From: https://www.cnblogs.com/memokeerbisi/p/18446023

相关文章

  • 2024-1-16 三年总结 192623
    三年总结因果有关于感情,似乎不知从何说起,但是可以肯定的是有因果安排。这份因果一则还债,二则让我领悟一些道理。其中除了需要学会判断之外,还需要明白色即是空,也即放下执着,但并不仅限于爱情。对于一般人而言,必须要找一个活下去的意义,在世界上如果没有自己存在的价值,或者说优越感......
  • 2023-12-15 博士挑战--达成 113823
    目录现状改进未来现状这个世界最神奇的,莫过于永远都有意外,第二步以和平的地方式提前达成了。1.意外一:迫于各种压力(我的态度、项目进度、其他学者进展、外界评价),导师已于上周开始给师姐改论文并尽快投出去。这结局是最好的,只是有点过早。2.意外二:没想到师妹论文写得太不勤快了,......
  • 20241001
    桌游制造我们可以对于每种图案记录拥有这种图案的有那些圆片,然后我们枚举每一个圆片,枚举这个圆片上面的图案,枚举拥有这种图案的圆片还有哪些,然后分别打上标记,如果有一个圆片明明已经有标记了,然而又要被打一次标记,那么我们可以直接输出\(NO\)如果标记都已经打完了,可还是......
  • 2024.10.[2, 3]训练记录
    10.2上午noip模拟比赛是8:00开始的,人是8:40起床的。T1猜了结论,秒了。结论是,一开始按照倒序排,连续是\(1\)的段\(reverse\)成正序。这样逆序对最多。感觉做法太简单\(O(n\logn)\)肯定不放。于是想了\(O(n)\)做法。最开始有\(\dfrac{n*(n-1)}{2}\)个逆序对,按段考虑......
  • 20241003
    缩进优化我们可以枚举\(i\)的所有倍数,我们让每一块中的数除以\(i\)相等,显然这是调和集数#include<bits/stdc++.h>usingnamespacestd;#defineintlonglongconstintN=1e7+5,INF=1e18;intn,a[N],sum[N],ans=INF,cnt[N];signedmain(){cin......
  • 2024/10/03
    \(100+20+0+55=175\),T4数组开小挂了\(45\),T3暴力写挂挂了\(20\)#A.旋律的总数这真的是提高组的题吗不考虑同构有\(m^n\)种排法,一种同构的排法可以偏移\(m\)次,直接相除得到答案\(m^{n-1}\)//BLuemoon_#include<bits/stdc++.h>usingnamespacestd;usingLL=l......
  • CSP-S 2024 第八次
    忘记写了,补一下A依次加入每个\(a_i\),拿个大根堆维护当前以\(i\)结尾的和最大子段,和超过\(s\)了就弹堆顶直到和不超过\(s\)。不过好像出现了一些语文事故,先不管了。B倍增预处理出\(f_i\)表示\(i\)上方第一个大于\(a_i\)的点,询问\(u,v,c\)时,先倍增找到\(u\)上......
  • 『模拟赛』多校A层冲刺NOIP2024模拟赛01
    Rank打得还可以总A.构造字符串签,但是挂了40pts。发现判条件只有相等和不相等,于是想到并查集维护连通块,将强制相同的两个位置的连通块合并,强制不同的先记下,最后统一判断。重点在细节处理,合并连通块时要将位置靠后的合并到靠前的上,注意\(LCP(x,y)=z\)在\(x+z,y+z\le......
  • 多校A层冲刺NOIP2024模拟赛【衡中】
    多校A层冲刺NOIP2024模拟赛01构造字符串咕咕咕寻宝咕咕咕点击查看代码#include<bits/stdc++.h>usingnamespacestd;constintmaxn=50009;intn,m,k,q,tot,cnt,vis[32767];inta[4]={1,-1,0,0};intb[4]={0,0,-1,1};map<int,short>mp[maxn];queue<pair<int,int>>......
  • [39] (多校联训) A层冲刺NOIP2024模拟赛01
    你们不感觉最近机房网越来越慢了吗,现在下个10M的东西要用三分钟,而且期间访问不了网站整个机房分1000Mbps的带宽为啥只能分这么一点,huge拿我电脑挖矿了?本来以为多校就是多校的,结果是真的多校,一百一十多个人在一块考huge:参加的都是咱们北方这几个强校你说得对,但是广东......