北斗编码实现

最新推荐文章于 2024-06-15 10:56:09 发布
最新推荐文章于 2024-06-15 10:56:09 发布
阅读量630 收藏 8
点赞数 16
文章标签： c++ 空间计算
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/chijingjing/article/details/135968625
版权

本文根据北斗格网编码标注编写的测试代码, 北斗国标描述网址 http://c.gb688.cn/bzgk/gb/showGb?type=online&hcno=77B7EA113926D3247F9688324D4A91C8

我将北斗编码整理成一张图, 如下(代码在文末附上):

#include "stdafx.h"
#include <vector>
#include <algorithm>
#include <iosfwd>
#include <string>
#include <map>
#include <algorithm>
#include <iostream>
#include <sstream>
#include <string>
#include <stdio.h>
#include <math.h>

std::string GBBeiDouCode(int level, double x, double y, double z)
{
	//if(y >=88)
	std::vector<std::string>  code;
	std::vector<std::string>  hcode;

	LocType ltype = XYLocType(x, y);
	std::pair<int, int> sym = LocTypeSym(ltype);

	int pre_r, pre_c;
	double posLon = -180.00000000, posLat = 0.00000000;
	if (x >= 180.00000)
		x = x - D3600_2_8888 / 100;
	if (x <= -180.000000)
		x = x - D3600_2_8888 / -100;
	//行列号按照左闭右开,最后180按照闭算到59
	for (int i = 1; i <= level; i++)
	{
		if (i == 1) {
			if (y >= 0)
				code.push_back("N");
			else
				code.push_back("S");
			int nr = floor((x - (posLon)) / 6.00000000); //0-59
			int nc = floor(fabs(y - posLat) / 4.00000000);//0-21
			std::string nrc;
			if (nr < 10)
			{
				nrc = nrc.append("0");
				int tmpnr = nr == 0 ? 1 : nr;
				nrc = nrc.append(std::to_string(tmpnr));
			}
			else
			{
				nrc = std::to_string(nr > 59 ? 60 : nr + 1);
			}
			if (nc > 21)
				nc = 21;
			char nccc = (char)(nc + 65);
			char f[2] = { nccc,0 };
			code.emplace_back();
			code[code.size() - 1].append(nrc);
			code[code.size() - 1].append(f);
			code[code.size() - 1].resize(3);
			nr = nr > 59 ? 59 : nr;
			nc = nc > 21 ? 21 : nc;
			pre_r = nr;
			pre_c = nc;

			continue;
		}
		else if (i == 2)
		{
			if (x > 0)
				posLon = posLon + (pre_r) * 6;
			else
			{
				/*	if (x <= -180.00000000000)
						posLon = posLon + (pre_r + 1) * 6;
					else*/
				posLon = posLon + (pre_r + 1) * 6;
			}
			if (y > 0)
				posLat = posLat + (pre_c) * 4;
			else
				posLat = posLat - (pre_c) * 4;

			int nc = (int)floor(fabs(y - posLat) / 0.5);//0-11
			int nr = (int)floor(fabs(x - posLon) / 0.5);//0-7
			nr = nr > 11 ? 11 : nr;
			nc = nc > 7 ? 7 : nc;
			char nrc = (char)(nr > 9 ? nr + 55 : nr + 48);
			char ncc = (char)(nc + 48);
			pre_r = nr;
			pre_c = nc;
			code.emplace_back();

			code[code.size() - 1].resize(2);
			code[code.size() - 1][0] = nrc;
			code[code.size() - 1][1] = ncc;
			continue;
		}
		else if (i == 3)
		{

			if (x > 0)
				posLon = posLon + (pre_r) * 0.5;
			else
				posLon = posLon - pre_r * 0.5;
			if (y > 0)
				posLat = posLat + (pre_c) * 0.5;
			else
				posLat = posLat - (pre_c) * 0.5;

			int nc = (int)floor(fabs(y - posLat) / D60_10);//0-2 
			int nr = (int)floor(fabs(x - posLon) / D60_15); //0-1
			nr = nr > 1 ? 1 : nr;
			nc = nc > 2 ? 2 : nc;
			pre_r = nr;
			pre_c = nc;
			int ncc = nc * 2 + nr;
			char nccc = (char)(ncc + 48);
			code.emplace_back();
			code[code.size() - 1].resize(1);
			code[code.size() - 1][0] = nccc;
			continue;

		}
		else if (i == 4)
		{
			if (x > 0)
				posLon = posLon + (pre_r)* D60_15;
			else
				posLon = posLon - (pre_r)* D60_15;
			if (y > 0)
				posLat = posLat + (pre_c)* D60_10;
			else
				posLat = posLat - (pre_c)* D60_10;

			int nc = (int)floor(fabs(y - posLat) / D60_1);
			int nr = (int)floor(fabs(x - posLon) / D60_1);
			nr = nr > 14 ? 14 : nr;
			nc = nc > 9 ? 9 : nc;
			pre_r = nr;
			pre_c = nc;
			char nrc = (char)(nr > 9 ? nr + 55 : nr + 48);
			char ncc = (char)(nc + 48);
			code.emplace_back();
			code[code.size() - 1].resize(2);
			code[code.size() - 1][0] = nrc;
			code[code.size() - 1][1] = ncc;
			continue;
		}
		else if (i == 5)
		{
			if (x > 0)
				posLon = posLon + (pre_r) * (D60_1);
			else
				posLon = posLon - (pre_r) * (D60_1);
			if (y > 0)
				posLat = posLat + (pre_c) * (D60_1);
			else
				posLat = posLat - (pre_c) * (D60_1);


			int nc = (int)floor(fabs(y - posLat) / D3600_4);//0-14
			int nr = (int)floor(fabs(x - posLon) / D3600_4);//0-14
			nr = nr > 14 ? 14 : nr;
			nc = nc > 14 ? 14 : nc;
			pre_r = nr;
			pre_c = nc;
			char nrc = (char)(nr > 9 ? nr + 55 : nr + 48);
			char ncc = (char)(nc > 9 ? nc + 55 : nc + 48);
			code.emplace_back();
			code[code.size() - 1].resize(2);
			code[code.size() - 1][0] = nrc;
			code[code.size() - 1][1] = ncc;

			continue;
		}
		else if (i == 6)
		{
			if (x > 0)
				posLon = posLon + (pre_r)* D3600_4;
			else
				posLon = posLon - (pre_r)* D3600_4;
			if (y > 0)
				posLat = posLat + (pre_c)* D3600_4;
			else
				posLat = posLat - (pre_c)* D3600_4;

			int nc = (int)floor(fabs(y - posLat) / D3600_2);//0-1
			int nr = (int)floor(fabs(x - posLon) / D3600_2);//0-1
			nr = nr > 1 ? 1 : nr;
			nc = nc > 1 ? 1 : nc;
			pre_r = nr;
			pre_c = nc;

			int ncc = nc * 2 + nr;
			char nccc = (char)(ncc + 48);
			code.emplace_back();
			code[code.size() - 1].resize(1);
			code[code.size() - 1][0] = nccc;

			continue;
		}
		else if (i >= 7 && i <= 10)
		{
			int d = i - 7 + 1;
			double dd2 = D3600_2 / pow(8, d);
			double dd1 = D3600_2 / pow(8, d - 1);
			if (x > 0)
				posLon = posLon + (pre_r)* dd1;
			else
				posLon = posLon - (pre_r)* dd1;
			if (y > 0)
				posLat = posLat + (pre_c)* dd1;
			else
				posLat = posLat - (pre_c)* dd1;

			int nc = (int)floor(fabs(y - posLat) / dd2);
			int nr = (int)floor(fabs(x - posLon) / dd2);
			nr = nr > 7 ? 7 : nr;
			nc = nc > 7 ? 7 : nc;
			pre_r = nr;
			pre_c = nc;
			char nrc = (char)(nr + 48);
			char ncc = (char)(nc + 48);
			code.emplace_back();
			code[code.size() - 1].resize(2);
			code[code.size() - 1][0] = nrc;
			code[code.size() - 1][1] = ncc;

			continue;
		}
	}

	//n0
	double r0 = 6378137.000000000;



	// a=1+sta0 b=(h+r0)/r0
	//n=sta0/sta * logb /log a b //a 为底, b为指

	double PI = 3.1415926535897932;
	double sta0 =  PI / 180.00000000;
	double sta = g_zdel[level-1]  * 0.017453292519943295;
	double rn = (r0 + z) / r0;
	double fg = (sta0 / sta) * (log(rn) / log(1 + sta0));
	int ffg = floor(fg);

	std::vector<std::pair<int, int>>g_levelStartEndLocation =
	{ { 32,32 },		{ 26, 31 },		{ 23, 25 },		{ 22, 22 },		{ 18, 21 },
	{ 14, 17 },		{ 13, 13 },		{ 10, 12 },		{ 7, 9 } ,		{ 4, 6 } ,		{ 1, 3 } };
	// 1 6 3 1 4 4 1 3 3 3 3 
	int left = 0;
	for (int j = 0; j <= level; ++j)
	{
		int left = 32 - g_levelStartEndLocation[j].second;
		int right = g_levelStartEndLocation[j].first + 32 - g_levelStartEndLocation[j].second - 1;
		int nVal = (((unsigned int)ffg << left) >> right);
		if (j == 3 || j == 6) {
			if (nVal > 1)
				nVal = 1;
		}
		else if (j == 4 || j == 5) {
			if (nVal > 14)
				nVal = 14;
		}
		else if (j == 2 || j >= 7)
		{
			if (nVal > 7)
				nVal = 7;
		}
		else if (j == 1)
		{
			if (nVal > 63)
				nVal = 63;
		}

		if (j == 4 || j == 5) {
			std::string s;
			char nrc = (char)(nVal > 9 ? nVal + 55 : nVal + 48);
			s.resize(1);
			s[0] = nrc;
			hcode.push_back(s);
		}
		else if (j == 1 && nVal < 10) {
			std::string f = "0";
			f.append(std::to_string(nVal));
			hcode.push_back(f.c_str());
		}
		else {
			hcode.push_back(std::to_string(nVal).c_str());
		}
	}

	std::stringbuf buf;
	for (int i = 0; i < code.size(); i++)
	{
		buf.sputn(code[i].c_str(), code[i].size());
		buf.sputn(hcode[i].c_str(), hcode[i].size());
	}

	return buf.str();

}

反算代码后面再写,

cejutue
关注
  • 16
    点赞
  • 踩
  • 8
    收藏
    觉得还不错? 一键收藏
  • 3
    评论
北斗得与失（1-二次编码的缺点）
LoseOfLock2的博客
07-30 3891
北斗得与失（1-二次编码的缺点） 摘要：北斗B1I信号在D1导航电文中引入现代化GPS和Galileo系统常用的二次编码调制，但却没有深入理解二次编码的应用场景。在现阶段，北斗系统不具备应用二次编码条件的情况下，这“画虎不成反类猫”的“借鉴”，降低了系统的捕获灵敏度。
bch编码matlab,BCH码的编码方法.doc
weixin_36084613的博客
03-28 3088
一、实验目的1、掌握循环码的编码原理2、掌握BCH码的编码方法3、了解编码与对误码性能的改善二、实验内容1、自行设置BCH码的参数，给出生成的BCH码；2、利用encode库函数实现编码；3、搭建一个通信仿真模块，并给出运行结果，分析BCH码对通信性能的影响；3、整理好所有的程序清单，并作注释。三、实验结果本原多项式，可纠正2位错误时，生成多项式为，写出生成矩阵，给出产生(15，7，2)BCH码的...
北斗卫星系统中几种常见编码的性能比较
qq_44905989的博客
04-24 1673
weil码的编码方法，相对于gold码的编码方式更为优秀，互相关特性和自相关特性都要好，但是设计方面更为复杂，程序设计更为复杂，gold码虽然并没有weil码的性能那么好，也就差一点点，但是也具有较好的自相关和互相关特性，编码方式更为简单。
北斗网格二维编码
Z_W_H_的博客
12-07 801
东经为 E —— east:东方 西经为 W —— west:西方 南纬为 S —— south:南方 北纬为 N —— north。其他级别网格计算见国标（有详细的计算过程和示例）将计算好的编码添加到列表中。
北斗三代卫星编号和发射数据
pirricy的博客
12-10 5510
出于工程需要，我需要北斗三代系统的所有卫星编号和发射日期这种数据。然而找来找去始终只有一些翻来覆去的新闻，北斗官网上也没有对卫星发射编号做详细的整理。明明是自己国家发射的卫星，却让我偏偏去北美国防部的公开数据里面查到北斗卫星的编号。发射时间从航天爱好者网整理的。 | 这是北斗三代系统的数据 | 用不同的颜色标注指的是同样的颜色是在相近的轨道上的。 颜色分类是按照以下STK做出的仿真图分类的。 另外需要全部北斗系统信息的可以去维基百科找。哎真是醉了，自己国家的东西全部要去国外找资料。 维基百科——北斗 ..
北斗卫星导航系统BCH解码的实现结构
王大锤的专栏
07-19 1805
如下：
北斗短报文通信实现源码
热门推荐
简与美
05-29 3万+
/****************************************************************************************************** 1、北斗协议和GPS协议不一样，不是以"\r\n"为一条协议的结束. 2、printf函数中的打印信息是const字符串常量，放在cpu内部flash，北斗模块printf打印
北斗B1I测距码特性及程序生成
红茶绿茶和奶茶的博客
09-28 7709
1测距码特性 测距码就是伪随机噪声( pseudo random noise，PRN）码，也称作伪码，常作为扩频系统的扩频码，工程上常用一系列0或1组成的二进制序列码，按照一定的编码规则来产生伪噪声码。伪码具有确定性、周期性以及良好的自相关性，大体被分为三类:m序列、组合码和非线性码，北斗伪卫星信号的Ⅰ支路测距码是一种由m序列构成的组合码。 北斗伪卫星系统中，Ⅰ支路使用的测距码是Gold 码，由两个11级的m序列G1、G2,通过模二和，并截短一个码片后产生。北斗伪卫星系统中Gold码的速率为2.046M.
基于北斗的气象自动观测系统数据编码设计
10-16
针对气象自动观测网数据传输需要，提出了基于北斗卫星导航系统的气象自动观测数据编码设计方案。该方案采用北斗短报文通信功能传输帧数据，从而实现自动观测数据的传输。采用该方案使得指挥终端能实时接收自动站发送...
基于北斗通信的工业过程数据压缩方法
10-16
该压缩方法分两步，第一步为有损压缩，以工业数据库压缩技术中的旋转门算法为基础，同时为了实现有损压缩的精度可调，用改进的BP神经网络PID控制器对旋转门算法的参数进行在线调整；第二步为无损压缩，提出了前置...
北斗网格位置码+图形分幅详解.zip
10-18
主要包括：GBT-13989-2012 国家基本比例尺地形图分幅和编号.pdf、GeoSOT编码.pdf、北斗网格位置码.pdf、国家大地测量基本技术规定.pdf、专题一全球等经纬度剖分网格GeoSOT1专题二.pdf 下一篇博文我将附上代码实现，...
北斗返回语句解析代码
05-07
北斗返回信息的解析，将北斗指挥机接收到的BDTXR语句中的信息提取出来，通过utf8编码来实现解码。
基于北斗Ⅱ_GPS车载定位终端的研究与实现.pdf
09-08
4. **数据传输协议**：研究分析了北斗二号和GPS的数据传输协议，理解它们的帧结构、数据编码方式和同步机制，这对于正确接收和解析卫星信号至关重要。 5. **Android串口通信技术**：在Android系统中，通过串口...
C语言经典指针运算笔试题图文解析
喜羊羊的博客
06-12 675
指针运算常常现在面试题中，画图解决是最好的办法。
CSAPP大作业 hello的一生
m0_73927795的博客
06-14 1060
计算机系统大作业 hello的一生
C++ Primer Chapter 10 Generic Algorithms
qq_22609521的博客
06-12 240
其中，capture list（捕获列表）是一个lambda所在函数中定义的局部变量的列表（通常为空）；单数，与普通函数不同，lambda必须使用尾置返回来指定返回类型。一个lambda表达式表示一个可调用的代码单元。与任何函数类似，一个lambda具有一个返回类型、一个参数列表和一个函数体。即，如果e是一个可调用的表达式，则我们可以编写代码e(args)，其中args是一个逗号分隔的一个或多个参数的列表。对于一个对象或一个表达式，如果对其使用调用运算符，则称它为可调用的。
二维数组与指针【C语言】
最新发布
2301_78085839的博客
06-15 714
​ 综上，一维数组的数组名是首元素的地址，对数组名解引用可以得到首元素（*arr = 首元素）；​ 二维数组的数组名是第一行所表示的一维数组的地址，对二维数组名第一次解引用得到的是第一行所表示的一维数组，即*arr = arr[0]，第一次解引用得到的是一维数组，那么只需在解引用一次就可得到首元素，即*(*arr) = 首元素**（以此类推，三维数组需对数组名解引用三次才能得到首元素）
16-C语言中的作用域——函数声明，函数头，局部与全局作用域以及作用域的掩盖
everything is possible，anything is OK.
06-11 576
16-C语言中的作用域——函数声明，函数头，局部与全局作用域以及作用域的掩盖
使用verilog语言实现北斗C/A码扩频与解扩电路设计编码
06-08
在Verilog中实现北斗C/A码扩频与解扩电路设计需要编写多个模块，下面是一个可能的编码示例： ``` // 扩频码模块 module prn_generator(clk, reset, prn); input clk, reset; output [10:0] prn; // 11位的伪随机码 reg [10:0] shift_reg; always @(posedge clk) begin if (reset) begin shift_reg <= 11'b11111111111; // 初始状态为全1 end else begin shift_reg <= {shift_reg[9:0], shift_reg[10] ^ shift_reg[2] ^ shift_reg[1] ^ shift_reg[0]}; end end assign prn = shift_reg; endmodule // 扩频器模块 module prn_multiplier(clk, reset, prn, signal_in, signal_out); input clk, reset; input [10:0] prn; input [7:0] signal_in; output [7:0] signal_out; reg [7:0] product; always @(posedge clk) begin if (reset) begin product <= 8'b0; end else begin product <= prn[10] ? signal_in : ~signal_in; // 根据扩频码的值进行相应的乘法运算 end end assign signal_out = product; endmodule // 解扩器模块 module prn_despreader(clk, reset, prn, signal_in, signal_out); input clk, reset; input [10:0] prn; input [7:0] signal_in; output [7:0] signal_out; reg [7:0] product; always @(posedge clk) begin if (reset) begin product <= 8'b0; end else begin product <= prn[10] ? signal_in : 8'b0; // 根据扩频码的值进行相应的乘法运算 end end assign signal_out = product; endmodule // 信号发生器模块 module signal_generator(clk, reset, signal_out); input clk, reset; output [7:0] signal_out; reg [31:0] counter; always @(posedge clk) begin if (reset) begin counter <= 32'b0; end else begin counter <= counter + 32'h1; // 以固定的频率递增计数器 end end assign signal_out = {8{counter[22]}}; // 将计数器的第22位复制8次作为输出信号 endmodule // 接收器模块 module signal_receiver(clk, reset, received_signal, prn, signal_out); input clk, reset; input [7:0] received_signal; input [10:0] prn; output [7:0] signal_out; wire [7:0] despread_signal; prn_despreader despreader(clk, reset, prn, received_signal, despread_signal); assign signal_out = despread_signal; endmodule ``` 以上是一个简单的Verilog实现示例，仅供参考。实际应用中需要根据具体的需求进行调整。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 3
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值