sin

经纬度与平面坐标互转,经纬度与空间直角坐标互转(C++代码)

早过忘川 提交于 2019-11-27 21:28:34
在三维激光点云处理中,需经常用到经纬度与平面坐标、空间直角坐标互转的功能,有时只是临时写一个测试demo,不想调用gdal,太麻烦,希望有更简单的调用方式。 网上一通搜索,并没有找到很完整的代码,一些代码杂乱无章,正确性还需确认,于是自己动手写了这四个转换函数,在此与大家分享使用: 头文件: /******************************************************************* * * 作者: Sun Zhenxing * 创建日期: 20190819 * * 说明:实现经纬度与平面坐标互转,实现经纬度与空间直角坐标互转 * ******************************************************************/ #ifndef ZTGEOGRAPHYCOORDINATETRANSFORM_H #define ZTGEOGRAPHYCOORDINATETRANSFORM_H #include <math.h> struct EllipsoidParameter { double a, b, f; double e2, ep2; // 高斯投影参数 double c; double a0, a2, a4, a6; EllipsoidParameter() { // Default:

MATLAB 滤波

只愿长相守 提交于 2019-11-27 20:50:24
clear fs=500; t=0:pi/200:pi/1.5; x=sin(t)+sin(3*t)+1; x_=x'; wp=1.5/500;%截止频率,以下的频率可以通过 ws=5/500;%被滤除 Rp=1; As=10; subplot(311); plot(t,x); title('sin(t)+sin(5*t)+4'); %N滤波阶数,Wn代表滤波器的截止频率 [N,wc]=buttord(wp,ws,Rp,As); [B,A]=butter(N,wc,'low'); [H,W]=freqz(B,A); y=filter(B,A,x); subplot(312); test=sin(t); plot(t,test); %plot(W,abs(H)); title('y=sin(x)'); subplot(313); plot(t,y) title('低通滤波'); 来源: https://www.cnblogs.com/hsy1941/p/11376077.html

python DataScience数据分析笔记day02

房东的猫 提交于 2019-11-27 20:34:38
day02-数据可视化-坐标-图例-标注-子图-刻度-半对数坐标-散点图-条形图-饼状图-等高线 3、设置坐标范围 mp.xlim(水平坐标最小值,水平坐标最大值) mp.ylim(垂直坐标最小值,垂直坐标最大值) 示例: import numpy as np import matplotlib.pyplot as mp x = np.linspace(-np.pi,np.pi,1000) cos_y = np.cos(x) / 2 sin_y = np.sin(x) mp.xlim(x.min() * 1.2,x.max() * 1.2 ) mp.ylim(sin_y.min() * 1.2,sin_y.max() * 1.2 ) mp.plot(x,cos_y,linestyle='--',linewidth= 1 ,color = 'green') # -- 为虚线 mp.plot(x,sin_y,linestyle=':', linewidth = 4,color='blue') # : 为点线 mp.show() 4、设置坐标刻度 mp.xticks([位置序列],[标签序列]),标签序列可以不需要 mp.yticks([位置序列],[标签序列]) 示例: import numpy as np import matplotlib.pyplot as mp x = np

matlab基本语法

こ雲淡風輕ζ 提交于 2019-11-27 09:55:24
2019-08-15 13:12:07 %基本绘图 x = 0:0.5:2*pi; y = sin(x); plot(x,y); %第二个参数为矩阵 y1 = sin(x); y2 = cos(x); y3 = 0.02 * exp(x); y4 = x; y5 = 0.02 * tan(x); z = [y1;y2;y3;y4;y5]; plot(x,z); %两个参数都是矩阵 x1 = 0:0.01:2 * pi; x2 = -pi:0.01:pi; x = [x1;x2]; y1 = cos(x1); y2 = sin(x2); y = [y1;y2]; plot(x,y); %plot只有一个参数 x = linspace(0,2*pi,200); y = sin(x); plot(y); y2 = cos(x); plot(y2); y3 = y + i * y2; plot(y3); %plot有多个参数 x1 = linspace(0,2*pi,200); x2 = linspace(0,2*pi,100); y1 = cos(x1); y2 = sin(x2); plot(x1,y1,x2,y2); %plot含有的曲线选项 x = linspace(0,2*pi,100); y = sin(x); plot(x,y,'k');% r(红),g(绿),y(黄),m

深度学习基础 - 余弦定理

依然范特西╮ 提交于 2019-11-27 05:44:58
深度学习基础 - 余弦定理 flyfish A D = b cos ⁡ A , C D = b sin ⁡ A , A D=b \cos A, \\C D=b \sin A, A D = b cos A , C D = b sin A , B D = A B − A D B D = c − b cos ⁡ A B D=A B-A D \\ B D=c-b \cos A B D = A B − A D B D = c − b cos A 根据是 勾股定理 B C 2 = B D 2 + C D 2 = ( c − b cos ⁡ A ) 2 + ( b sin ⁡ A ) 2 = c 2 − 2 c b cos ⁡ A + b 2 整 理 得 &ThinSpace; a 2 = b 2 + c 2 − 2 b c cos ⁡ A \begin{aligned} B C^{2} &=B D^{2}+C D^{2} \\ &=(c-b \cos A)^{2}+(b \sin A)^{2} \\ &=c^{2}-2 c b \cos A+b^{2} \\ \mathbb{整理得} \, a^{2}=& b^{2}+c^{2}-2 b c \cos A \end{aligned} B C 2 整 理 得 a 2 = ​ = B D 2 + C D 2 = ( c − b cos A ) 2

2018ICPC 沈阳 训练赛

扶醉桌前 提交于 2019-11-27 03:27:29
C.1:33:41(-1) solved by zcz 推式子 找规律 #include<iostream> using namespace std; long long F(long long n,long long mod) { long long rec=1; for(long long i=1;i<=n;i++) { rec*=i; rec%=mod; } return rec; } long long F2(long long t) { return t*t-2*t+2; } int main() { int T; cin>>T; int cas=1; while(T--) { long long n,k,q; scanf("%lld %lld %lld",&n,&k,&q); if(k>n) k=n; long long ans=F(k,q)*(k*(n-k)+F2(n-k)+n-k-1)%q; if(k==1) ans=F2(n)%q; if(k==n) ans=F(n,q); if(k==2&&n==3) ans=6%q; printf("Case #%d: %lld\n",cas++,ans); } return 0; } C G.2:18:23(-4) solved by hl 预处理出1 - 1e7范围内的k可以被哪些i² + j²得出,显然对于每个k

利用subplot绘制多个图像

被刻印的时光 ゝ 提交于 2019-11-26 10:01:48
利用subplot绘制多个图像 subplot(m,n,p) subplot是将多个图画到一个平面上的函数,m是行,n是列,p是所要绘制图所在的位置 x = 0:0.1:100; sinY = sin(x); cosY = cos(Y); subplot(2, 1, 1) plot(x, sinY) %绘制第一个图像 y = sin(x) subplot(2, 1, 2) plot(x, cosY) %绘制第二个图像 y = cos(x) 来源: https://www.cnblogs.com/woxiaosade/p/11317437.html

如何理解傅立叶级数公式?

筅森魡賤 提交于 2019-11-26 07:34:21
文章目录 1. 对周期函数进行分解的猜想 2. 分解的思路 2.1 常数项 2.2 通过sin(x),cos(x)进行分解 2.3 保证组合出来周期为T 2.4 调整振幅 3. sin(x)的另外一种表示方法 3.1 $e^{i\omega t}$ 3.2 通过$e^{i\omega t}$表示sin(t) 4. 通过频域来求系数 4.1 函数是线性组合 4.2 如何求正交基的坐标 如何求sin(nt)基下的坐标 4.4 更一般的 5. 傅立叶级数的另外一种表现形式 1. 对周期函数进行分解的猜想 拉格朗日等数学家发现某些周期函数可以由三角函数的和来表示,比如下图中,黑色的斜线就是周期为2\pi的函数,而红色的曲线是三角函数之和,可以看出两者确实近似: 而另外一位数学家,让·巴普蒂斯·约瑟夫·傅里叶男爵(1768-1830)猜测任意周期函数都可以写成三角函数之和。 2. 分解的思路 假设f(x)是周期为T的函数,傅里叶男爵会怎么构造三角函数的和,使之等于f(x)? 2.1 常数项 对于 y = C , C ∈ R y=C,C\in\mathbb{R} y = C , C ∈ R 这样的常数函数: 根据周期函数的定义,常数函数是周期函数,周期为任意实数。 所以,分解里面得有一个常数项。 2.2 通过sin(x),cos(x)进行分解 首先,sin(x),cos(x)是周期函数