【C语言】计算实时太阳角度（高度角、方位角），以及使用stm32单片机实时获取时间戳

最新推荐文章于 2024-08-13 15:52:32 发布
最新推荐文章于 2024-08-13 15:52:32 发布
阅读量1.7k 收藏 11
点赞数 12
文章标签： c语言 开发语言
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_70042862/article/details/134554703
版权

整体计算方法

在编写该代码的过程中寻找了多篇博文和论文，综合所有文章且按网上的以0时的方位角的0°，且随时间累加累加至360度。我修改了博文和论文的一些角度的计算方法。得到一下代码与网站计算的方位角相互验证过，误差不超过1°
验证网站

太阳赤纬角

在这里插入图片描述
n为该天在当年的序号
另一种方法没尝试成功：
在这里插入图片描述
代码如下读者可以自行尝试，我没尝试成功

	int n, n0;
	double t, theta, delta;

	// 计算n和n0
	n = month * 30 + day;
	n0 = 79.6764 + 0.2422 * (year - 1985) - ((year - 1985) / 4);

	// 计算t
	t = n ;

	// 计算θ（日角）
	theta = 2 * M_PI * t / 365.2422;

	// 计算δ（太阳赤纬）
	delta = 0.3723 + 23.2567 * sin(theta) + 0.1149 * sin(2 * theta) - 0.1712 * sin(3 * theta)
	        - 0.758 * cos(theta) + 0.3656 * cos(2 * theta) + 0.0201 * cos(3 * theta);

太阳时角

在这里插入图片描述
时角这里的话中午时时角为0下午时角为负，可以简化为时角=15*(t-12)
不过这里计算的是时区的时间，不是当地时间所以存在误差，我怀疑代码计算的误差大部分都是出自这里
当地时角位于东7区的换算方法: (120度 – 当地经度) * 4，1个小时15度，平均下来，1度4分钟。
既时角=15*(t-12)+(120-经度)

太阳高度角

在这里插入图片描述
通过上述两个方程，即可求出高度角的各个参数。

太阳方位角

在这里插入图片描述
h为太阳高度角
Φ为纬度数据
σ为太阳的赤纬角

代码部分

C语言三角函数运算

因为角度计算涉及大量三角函数，所以C语言的三角函数计算不可避免需要使用。

C语言的三角函数计算函数

需要包含头文件#include <math.h>

  double __cdecl sin(double _X);
  double __cdecl cos(double _X);
  double __cdecl tan(double _X);
  double __cdecl asin(double _X);
  double __cdecl acos(double _X);
  double __cdecl atan(double _X);
  double __cdecl atan2(double _Y,double _X);

其中sin、cos、tan 函数的参数是一个用弧度表示的角度，分别返回这个角度的正弦、余弦、正切值。度和弧度都是衡量角的大小的单位。度用 ° 来表示，弧度用rad表示。

1度＝π/180 弧度 ( ≈0.017453弧度 )

asin()、acos、atan 函数分别返回它们的参数的反正弦、反余弦、反正切值。功能刚好和sin、cos、tan 函数相反。

赤纬角计算

通过公式：
在这里插入图片描述

/********************************************************************************************************************************************
*@*名称 ： Declination_Angle_Calculation
*@*功能 ： 计算太阳赤纬角
*@*形参 ： YearData年数据	DayData日数据
*@*返回值 ： DeclinationAngle太阳赤纬角 
*********************************************************************************************************************************************/
static double Declination_Angle_Calculation(double Daydata,short YearData)
{
	double DeclinationAngle;
	DeclinationAngle = 23.45 * sin((Daydata-80)/370.0 * 360*PI/180);		//赤纬角公式 
	
	return DeclinationAngle;												//返回赤纬角 
}

时角计算

通过公式：
时角=15*(t-12)

/********************************************************************************************************************************************
*@*名称 ： Hour_Angle_Calculation
*@*功能 ： 计算太阳时角
*@*形参 ： HourData小时数据 
*@*返回值 ： HourAngle太阳时角 
*********************************************************************************************************************************************/
static double Hour_Angle_Calculation(double HourData)
{
	double HourAngle;
	
	HourAngle = 15*(HourData-12);
	
	return HourAngle;
	
}

太阳高度角计算

通过公式
在这里插入图片描述

/********************************************************************************************************************************************
*@*名称 ： Height_Angle_Calculation
*@*功能 ： 计算太阳高度角
*@*形参 ： LatitudeData纬度数据		HourData小时数据 	DeclinationAngle赤纬角数据 
*@*返回值 ： ElevatingAngle太阳高度角 
*********************************************************************************************************************************************/
static double Height_Angle_Calculation(double LatitudeData,double HourAngle,double DeclinationAngle)
{
	double ElevatingAngle;													//高度角
	double ElevatingAngle_sin;												//高度角正弦值
	
	double DeclinationAngle_sin = sin(DeclinationAngle * PI / 180);			//赤纬角的正弦值 
	double DeclinationAngle_cos = cos(DeclinationAngle * PI / 180);			//赤纬角余弦值 
	double LatitudeData_sin = sin(LatitudeData * PI / 180);					//纬度正弦值 
	double LatitudeData_cos = cos(LatitudeData * PI / 180);					//纬度余弦值 
	double HourAngle_cos = cos(HourAngle * PI/180);
		
	ElevatingAngle_sin = (DeclinationAngle_sin*LatitudeData_sin) + (DeclinationAngle_cos*LatitudeData_cos*HourAngle_cos);
	ElevatingAngle = asin(ElevatingAngle_sin)*180/PI;
	return ElevatingAngle;
}

太阳方位角计算

通过公式
在这里插入图片描述

/********************************************************************************************************************************************
*@*名称 ： Azimuth_Angle_Calculation
*@*功能 ： 计算太阳方位角
*@*形参 ： ElevatingAngle太阳高度角		HourAngle太阳时角	DeclinationAngle赤纬角数据	HourData小时数据	LatitudeData纬度数据
*@*返回值 ： 太阳方位角（0点开始360度） 
*********************************************************************************************************************************************/
static double  Azimuth_Angle_Calculation(double ElevatingAngle,double HourAngle,double DeclinationAngle,double HourData,double LatitudeData)
{
	double AzimuthAngle,AzimuthAngle_cos;
	
 	double ElevatingAngle_cos = cos(ElevatingAngle * PI / 180);				//高度角余弦值 
 	double ElevatingAngle_sin = sin(ElevatingAngle * PI / 180);				//高度角正弦值 
	double LatitudeData_sin = sin(LatitudeData * PI / 180);					//纬度正弦值 
	double LatitudeData_cos = cos(LatitudeData * PI / 180);					//纬度余弦值  	
 	double HourAngle_sin = sin(HourAngle * PI / 180); 						//时角正弦值 
 	double DeclinationAngle_cos = cos(DeclinationAngle * PI / 180);			//赤纬角余弦值 
 	double DeclinationAngle_sin = sin(DeclinationAngle * PI / 180);			//赤纬角正弦值
 	
 	AzimuthAngle_cos = (ElevatingAngle_sin*LatitudeData_sin-DeclinationAngle_sin)/(ElevatingAngle_cos*LatitudeData_cos);		//初步计算太阳方位角弧度值的cos值 
 	AzimuthAngle = acos(AzimuthAngle_cos)*180/PI;																				//计算方位角的角度 
 	
	if (HourData > 12) 														//不判断的话12点为180度最高，然后两边的角度都是180度递减，所以需要一个360度的角度反馈需要将方位角判断一下。 
	{
		return AzimuthAngle + 180;
	}
	else 
	{
		return -AzimuthAngle + 180;
	}
	
}

其中函数尾部的返回时180度加减计算得出的方位角，是为了对齐网页计算得出的方位角，以0时的方位角的0°，且随时间累加累加至360度。

角度计算到此就结束了

时间戳转换

下文内容多数引用自：博客

因为我打算用stm32f103mcu进行太阳方位角以及高度角的计算，所以需要用stm32单片机获取时间戳信息

我采用的方法是使用esp8266-01s模块来通过联网获取时间戳的方式得到时间。

配置 ESP8266 入网

AT+CWMODE=1                     // station模式，设置成3也行
AT+RST                          // 重启模块
AT+CWJAP="ssid","password"      // 连接AP

连接授时服务器

AT+CIPMUX=0                                // 设置WiFi模块为单路链接
AT+CIPSTART="UDP","1.cn.pool.ntp.org",123  
// 链接到 NTP 服务器，时间服务器默认的端口号是 123，协议为UDP.AT+CIPSEND=48                              
// 发送48byte的数据若ESP8266返回'>'则继续注意，严格遵照以下数据且使用16进制发送
//切勿发送新行E3 0 6 EC 0 0 0 0 0 0 0 0 31 4E 31 34 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
//用16进制接收数据
E3 0 6 EC 0 0 0 0 0 0 0 0 31 4E 31 34 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

不知道为什么我好像需要发送两遍下边的16进制数据，然后接收两帧数据

接收数据

  • 执行完上面两步，此时如果顺利的话会受到一大堆16进制数据
  • 诸如：41 54 2B 43 49 50 53 45 4E 44 3D 34 38 0D 0D 0A 0D 0A 4F 4B 0D 0A 3E 20 0D 0A 52 65 63 76 20 34 38 20 62 79 74 65 73 0D 0A 0D 0A 53 45 4E 44 20 4F 4B 0D 0A 0D 0A 2B 49 50 44 2C 34 38 3A 24 02 06 E9 00 00 00 52 00 00 08 F2 0A 45 02 22 E0 4F 57 31 79 AF 20 82 00 00 00 00 00 00 00 00 E0 4F 59 9C 2F 4A 40 07 E0 4F 59 9C 2F 4C 22 B9
  • 此时需要找到数据中的关键帧，也就是服务器返回的48位时间数据，它们通常具有以下格式：（如果找不到的话请重试2中的3、4步）
    • ① 它前面的数据为"2B 49 50 44 2C 34 38 3A"(+IPD,48:）,他后面没有东西;
    • ② 它通常以 24 开头.3.对数据提取后编号，其中[40][41][42][43]为我们需要的时间数据，在本例中为"E0 4F 5C 00".

处理数据

①将”E04F5C00“转为十进制为3763297280，此即为1900距今的秒数。
②把这个时间减去2208988800（70年），得到标准的Unix时间戳。[3763297280-2208988800=1554308480]
③[1554308480/31556736+1970=2019]，得到当下年份。
④计算月份和日期需要考虑平年闰年和每月天数，太麻烦我也没搞懂。
⑤[(1554308480/3600)%24+8=24]，得到北京时间-时(别忘了满24归零)
⑥[(1554308480/60)%60=21]，得到北京时间-分
⑦[1554308480%60=20]，得到北京时间-秒，既当下时间是 24：21：20
⑧将时间戳(1554308480)输入 时间戳生成网页 与计算结果进行比较

时间戳解码并加入代码

因为我们计算太阳角度的代码，只需要获取年份，当天在一年中的序号，以及小时分钟信息。
所以有以下计算

/********************************************************************************************************************************************
*@*名称 ： Timestamp_Calculation
*@*功能 ： 计算时间戳
*@*形参 ： Timestamp时间戳数据 	Year_data年数据存放的地址	DayData日数据存放的地址		 Hour_data小时数据存放的地址 
*@*返回值 ： 无
*********************************************************************************************************************************************/
static void  Timestamp_Calculation(unsigned int Timestamp,short *Year_data,short *DayData,double *Hour_data)
{
	Timestamp  = Timestamp-0x83AA7E80;						//将网络获取的时间戳减去70年得到标准时间戳 
	*Year_data = Timestamp/31556736+1970;					//用时间戳获取年份 
	*DayData = (short)((Timestamp%31556736)/86400)+1;		//获取日期 
	*Hour_data = (Timestamp/3600)%24+8;						//获取小时数 
	if(*Hour_data==24)
	{
		*Hour_data = 0;										//当时间为24点时归零 
	} 
	*Hour_data += ((Timestamp/60)%60)/60.0;					//获取分钟信息并且转为小数形式加入小时数中 
}

步骤：

  1. 将网络获取的时间戳减70年得到标准时间戳
  2. 然后用标准时间戳除一年的秒数，就得到距离时间戳开始运行的1970年过去了多少年
  3. 标准时间戳除于一天的秒数，为今年过去了多少天，还要加上没过完的今天，即加1。
  4. 小时数为将时间戳除一小时的秒数，得到至今运行了多少小时，再进行24小时取余，余今天过了多少小时
  5. 小时数在加上分钟数，得到带小数点的精确时间。（方法类似步骤5）
房东的哈士奇
关注
定量遥感：计算地方时和太阳高度角（C++代码）
「 刘一哥与GIS的故事」
06-16 1577
在定量遥感中，通常需要计算地方时和太阳高度角，本文采用C++语言实现。 #include <cmath> #include <iostream> #include <fstream> using namespace std; void main() { int JD,NF,Y,R,s[5],F[5]; float JF,WD; float N0; ifstream data1("d:\\result\\data-input.txt"); .
STM32开发 -- UTC、UNIX时间戳、北京时间之间的转换
热门推荐
不积跬步，无以至千里
02-18 1万+
我碰到的问题，GPS上UTC时间转北京时间和STM32上UNIX时间戳转北京时间。 这部分之前讲RTC和GPS的时候有涉及到一部分。 具体的RTC如何UNIX时间戳，和GNRMC如何解析得到UTC时间可以参看一下。 参看：STM32开发 – RTC详解 参看：STM32开发 – GPS模块开发详解 这里主要看三者转换方法： 一、UTC时间转换为北京时间 参看:UTC时间转换为北京时间 时间类型结构...
太阳方位角高度角计算C代码
07-17
有些项目中需要计算太阳的位置，如太阳能板追日功能等，代码使用NOAA（美国国家海洋和大气管理局）的算法，解算太阳的位置，包括太阳方位角和高度角，需要的输入信息：经纬度、日期、时间、时区。太阳方位角以正北为参考零度角，南北东西半球统一适用
STM32，根据时间和经纬度计算太阳位置（高度角方位角）
m0_73001789的博客
03-25 693
首先安利一个国外的网站，上面根据上述信息算位置网站说测量精度为0.01°咱们也不知道他说的是不是准的，但是博主参与的追光项目（向日葵）反馈的一天的角度变化与网站完全重合（咱们姑且认为他是准的）网站分享了一篇来自国外的论文，上面讲述了该网站的测量方法（博主学术不精，照着写出来的代码时乱码......）如果有路过的大佬像攻克一下我会将PDF放到博客后面供大家查看。
STM32—RTC实时时钟
最新发布
尘的博客
08-13 1563
BKP(Backup Registers)备份寄存器BKP可用于存储用户应用程序数据。当VDD(2.0~3.6V)电源被切断，他们仍然由VBAT(1.8~3.6V)维持供电。当系统在待机模式下被唤醒，或系统复位或电源复位时，他们也不会被复位TAMPER引脚产生的侵入事件将所有备份寄存器内容清除RTC引脚输出RTC校准时钟、RTC闹钟时钟脉冲或者秒脉冲存储RTC时钟校准寄存器用户数据存储容量：20字节（中容量和小容量）/84字节（大容量和互联型）RTC(Real Time Clock)实时时钟。
zzut太阳高度角方位角计算C语言代码+配套excel
03-24
C程序为根据经纬度、时钟时间计算太阳位置即太阳高度角、方位角等相关数据。 采用的是【1985年bourges提出的赤纬角度】计算方案。。 配套excel文件公式与C语言版本计算方法一致，vba可生成太阳相关系列数据。 用途： 可用于太阳能电池板自动对准太阳，经测试C程序可简单移植到arduino
太阳高度角与方位角计算（c语言）
qq_64777806的博客
03-21 337
原理不多说，直接给出代码
【数学建模】利用C语言来实现 太阳赤纬 太阳高度角 太阳方位角 计算和求解分析 树木树冠阴影面积与种植间距的编程计算分析研究
shaozheng0503的博客
07-17 5154
函数中，根据传入的经度（latitude）、太阳赤纬（solar_declination）和时角（hour_angle），计算太阳的高度角。根据给定的公式，利用经度、太阳赤纬和时角的三角函数计算出太阳高度角的正弦值sin_elevation，然后使用反正弦函数asin将其转换为太阳高度角的度数形式，并返回结果。函数中，根据传入的参数alpha（太阳高度角）、delta（太阳赤纬）、phi（观测地点纬度）、hour（当前时间，以小时为单位）和lambda（观测地点经度），计算太阳方位角。
stm32_时间戳
m0_56082077的博客
07-31 1455
1.RTC（Real Time Clock）实时时钟2.RTC是一个独立的定时器，可为系统提供时钟和日历的功能3.RTC和时钟配置系统处于后备区域，系统复位时数据不清零，VDD（2.0~3.6V）断电后可借助VBAT（1.8~3.6V）供电继续走时4.32位的可编程计数器，可对应Unix时间戳的秒计数器5.20位的可编程预分频器，可适配不同频率的输入时钟。
基于stm32的c语言版（HAL库函数）获取rtc时间和时间戳之间的转换例程
Robben.Han
06-04 8659
代码使用的是HAL库 const u8 Days[12] = {31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31}; const u16 monDays[12] = {0,31, 59, 90, 120, 151, 181, 212, 243, 273, 304, 334}; u32 time2Stamp(RTC_DateTypeDef date,...
时间戳功能（HAL库）
03-26
STM32的时间戳功能主要通过实时时钟(RTC, Real-Time Clock)来实现，这是一个能够独立于微控制器主时钟运行的低功耗时钟模块，可以用来追踪和记录特定事件的发生时间。在本项目中，我们将利用STM32的RTC时间戳功能，...
STM32 RTC（一）Unix时间戳
2301_76874885的博客
07-24 888
时间戳是不考虑闰秒的，那么闰秒又是什么？Unix时间戳（Unix TimeStamp），实质上是一个向上计数的秒计数器，不分年月日时分秒，它的起始点为伦敦标准时间的1970年1月1日0时0秒，截止到2023年，时间戳的计数已经有16亿多，一般的计算机采用的都是32位的整形数据，计数最高为21亿多，经人计算，2038年会出现计算机危机，也就是计数溢出，STM32的RTC采用的为无符号32位整形数据，2100年才会有这个问题不必担心，但随着时代发展，目前大多数使用的是64位。mktime是将日期转换成时间戳。
C语言实现字符转unix时间戳的简单实例
01-20
C语言实现字符转unix时间戳，需要先转成tm类型，再得到它的Unix时间戳。附上实现代码： #include #include int strtotime(char datetime[]) { struct tm tm_time; int unixtime; strptime(datetime, %Y-%m-%d %H...
STM32——RTC实时时钟（时间戳）
hehe的博客
11-06 4740
Unix 时间戳（Unix Timestamp）定义为从UTC/GMT的1970年1月1日 0时0分0秒开始所经过的秒数，不考虑闰秒。
太阳能最佳倾角c语言编程,用C语言编程计算太阳的高度角，用于控制太阳能电板最佳的倾角...
weixin_36472718的博客
05-21 734
满意答案supercatdo2014.01.11采纳率：45%等级：6已帮助：508人我们用h表示这个角度，它在数值上等于太阳在天球地平坐标系中的地平高度。太阳高度角随着地方时和太阳的赤纬的变化而变化。太阳直射纬度以δ表示，观测地地理纬度用φ表示，地方时(时角)以t表示，有太阳高度角的计算公式：sin h＝sin φ sin δ＋sin φ cos δ cos t日升日落，同一地点一天内太...
基于STM32的F1的RTC实时时钟+CUBEMX实时获取时间戳实现断电跨天同步时间戳
昊月光华的博客
10-08 3644
实现断电跨日期同步时间戳
[c++] 计算太阳高度角
XIEXINGHUA2010的博客
01-15 1571
/* 输入参数： Longitude - 经度（单位"度"） Latitude - 纬度（单位"度"） Year - 年 Month - 月 Day - 日 Hour - 时 Minute - 分 Second - 秒 输出参数： Height ...
太阳时角、太阳高度角、天顶角、太阳方位角和剖面角计算的Python程序
Fo*(Bi)的博客
01-15 4815
基本概念 经度，Longitude L，经度(longitude)是地球上一个地点离一根被称为本初子午线的南北方向走线以东或以西的度数。本初子午线的经度是0°，地球上其它地点的经度是向东到180°或向西到180°。在本初子午线以东的经度叫东经，在本初子午线以西的叫西经。东经用"E"表示，西经用"W"表示。 维度，Latitude φ(phi)，维度 ， 赤道以北或以南的角位置，北正；−90°≤φ≤90°。 赤纬角，Declination δ，赤纬，太阳正午（即，当太阳位于本初子午线时）太阳相对于赤道平面
基于STM32单片机采集数据并上云
weixin_43517147的博客
06-15 4191
基于STM32单片机采集数据并上云一、实验目的及要求二、实验软硬件环境三、实验过程3.1用STM32CudeMx工具创建HAL库工程模板3.2使用Keil5在HAL工程模板的基础上进行开发四、实验结果 一、实验目的及要求 利用STM32单片机通过连接DHT11温湿度传感器测出温湿度数据，并通过连接Esp8266模块把数据上传到云服务器中 二、实验软硬件环境 1、硬件 STM32F103C8T6开发板、DHT11温湿度传感器、ESP8266 2、软件 MDK5（Keil5）、STM32CudeMx 3、接线
嵌入式单片机基础篇（三十九）之CAN通讯以及STM32F1单片机程序
lostlll的博客
08-26 1万+
CAN通讯以及STM32F1单片机程序 1、CAN通讯基本概念： （1）定义： CAN是Controller Area Network 的缩写（以下称为CAN），是ISO国际标准化的串行通信协议。由德国电气商博世公司在1986 年率先提出。此后，CAN 通过ISO11898 及ISO11519 进行了标准化。现在在欧洲已是汽车网络的标准协议。 CAN协议经过ISO标准化后有两个标准：ISO11898标准和ISO11519-2标准。 其中ISO11898是针对通信速率为125Kbps~1Mbps的高速通信标准
评论 3
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值