(完整word版)高斯投影正反算 代码
(完整word版)Python题库(可编辑修改word版)
《Python 程序设计》填空题1-280 题
2017-03-23 董付国Python 小屋
1、Python 安装扩展库常用的是工具。(pip)
2、Python 标准库math 中用来计算平方根的函数是。(sqrt)
3、Python 程序文件扩展名主要有和两种,其中后者常用于GUI 程序。(py、pyw)
4、Python 源代码程序编译后的文件扩展名为。(pyc)
5、使用pip 工具升级科学计算扩展库numpy 的完整命令是
。(pip install --upgrade numpy)
6、使用pip 工具查看当前已安装的Python 扩展库的完整命令是
。(pip list)
7、在IDLE 交互模式中浏览上一条语句的快捷键是。(Alt+P)
8、使用pip 工具查看当前已安装Python 扩展库列表的完整命令是。(pip list)
9、在Python 中表示空类型。(None)
10、列表、元组、字符串是Python 的(有序?无序)序列。(有序)
11、查看变量类型的Python 内置函数是。(type())
12、查看变量内存地址的Python 内置函数是。(id())
13、以3 为实部4 为虚部,Python 复数的表达形式为或。(3+4j、3+4J)
14、Python 运算符中用来计算整商的是。(//)
15、Python 运算符中用来计算集合并集的是。(|)
16、使用运算符测试集合包含集合A 是否为集合B 的真子集的表达式可以写作
。(A<B)
17、表达式[1, 2, 3]*3 的执行结果为。([1, 2, 3, 1, 2, 3, 1, 2, 3])
(完整word版)ASCII码对照表.doc
(完整word版)ASCII码对照表.doc
十进制代
十六进制代码MCS 字符或缩写DEC 多国字符名码ASCII 控制字符 1
0 0 NUL 空字符
1 1 SOH 标题起始 (Ctrl/A)
2 2 STX 文本起始 (Ctrl/B)
3 3 ETX 文本结束 (Ctrl/C)
4 4 EOT 传输结束 (Ctrl/D)
5 5 ENQ 询问 (Ctrl/E)
6 6 ACK 认可 (Ctrl/F)
7 7 BEL 铃 (Ctrl/G)
8 8 BS 退格 (Ctrl/H)
9 9 HT 水平制表栏 (Ctrl/I)
10 0A LF 换行 (Ctrl/J)
11 0B VT 垂直制表栏 (Ctrl/K)
12 0C FF 换页 (Ctrl/L)
13 0D CR 回车 (Ctrl/M)
14 0E SO 移出 (Ctrl/N)
15 0F SI 移入 (Ctrl/O)
16 10 DLE 数据链接丢失 (Ctrl/P)
17 11 DC1 设备控制 1 (Ctrl/Q)
18 12 DC2 设备控制 2 (Ctrl/R)
19 13 DC3 设备控制 3 (Ctrl/S)
20 14 DC4 设备控制 4 (Ctrl/T)
21 15 NAK 否定接受 (Ctrl/U)
22 16 SYN 同步闲置符 (Ctrl/V)
23 17 ETB 传输块结束 (Ctrl/W)
24 18 CAN 取消 (Ctrl/X)
25 19 EM 媒体结束 (Ctrl/Y)
26 1A SUB 替换 (Ctrl/Z)
27 1B ESC 换码符
28 1C FS 文件分隔符
(完整word版)角谱法分析高斯光
矢量高斯光束传播分析和近轴球面近似有
效性
卡尔G.陈,保罗T.康科拉,胡安费雷拉,拉尔夫·K.海尔曼,和Mark
L.影子城堡
麻省理工大学,剑桥,马萨诸塞州02139
收稿2001年3月5日;五月接受24,2001;修改稿收到2001年6月20日
许多系统在光通信和测量利用高斯光束,如从单模光纤点衍射干涉自由空间传播和干涉光刻的分析,将受益于高斯光束传输的精确分析模型. 我们提出了高斯光束传播通过使用平面波的角谱的众所周
知的方法的完整矢量分析。高斯光束假定遍历一个自由,均匀,各向同性的线性非磁性电介质。角谱表示在其载体形式,被施加到一个问题高斯光强的边界条件。经过一些数学运算,每个非零传播电场分量被表示为一个幂级数展开项。先前导出的分析工作幂级数的横向场,其中第一项(零阶)中的膨胀对应于通常的标量傍轴近似。我们确认这个结果,并得出相应的纵向幂系列。我们证明了领先的纵向期限相当于其数值第一学期横标量傍轴术语以上,从而表明当超越了一个完整的矢量理论需要标量傍轴近似。尽管一个紧凑的分析形式主义的优点,从而实现快速和高斯光束系统的精确建模,这种方法有一个显着的缺点。高阶条件在分歧是从最初的边界足够远的位置,产生非物理结果。因此,任何有意义的使用扩展方法的要求进行了认
真研究适用性的范围。通过考虑到从傍轴高斯波球方法的过渡,我们能够得出一个简单的表达式在其中产生一系列数值令人满意的答案的范围。©2002 美国光学学会
OCIS代码:0260.2110,000.4430,350.5500
1.引言
由于其简洁的物理直观表示,角谱表示法已经被用来解决多种问题包括传播的问题和高斯波的反射。它的理论基础已经有大量的作者证明非常稳固。一般使用它的标量形式来表示在近轴结合以及传播领域的计算【1-3】。虽然研究人员如阿格拉瓦尔和Pattanayak【4】已延长到解决方案傍轴结果,他们解决问题的方法仍然是标量的性质。
(完整word版)LISP函数(分类)大全
AutoLisp函数
一、数学运算功能函数
1.l(十数值数值…)返回:累计实数或整数数值
1.2(一数值数值…)返回:差值
1.3(*数值数值…)返回:所有数值乘积
1.4(/ 数值数值…)返回:第一个数值除以第二个以后数值的商 1.5(l十数值)返回:数值十l
l. 6(1—数值)返回:数值一l
l.7(abs 数值)返回:数值的绝对值
1.8(atan 数值)返回:反正切值
1.9(cos 角度)返回:角度的余弦值,角度值为弧度
1.10(exp 数值)返回:数值的指数
1.11(expt 底数指数)返回:底数的指数值
1.12(fix 数值)返回:将数值转换为整数值
1.14(gcd 数值1 数值2)返回:两数值的最大公因数
1.15(log 数值)返回:数值的自然对数值
1.16(max 数值数值…)返回:数值中的最大值
1.17(min 数值数值…)返回:数值中的最小值
1.18 pi 常数∏,其值约为3.1415926
1.19(rem 数值 1数值 2)返回:M数值的相除的余数l.20(sin 角度)返回:角度的正旋值,角度值为弧度
1.21(sqrt 数值)返回:数值的平方根
二、检验与逻辑运算功能函数
2.l(= 表达式1 表达式2)比较表达式1是否等于式2,适用数值及字符串
2.2 (/= 表达式1 表达式2)比较表达式1是否大于等于表达式2
2.3(<表达式1 表达式2) 比较表达式1是否<小于表达式2
2.4(<= 表达式1 表达式2)比较表达式1是否<一小于等于表达式2
2.5(>表达式1 表达式2)比较表达式1是否>大于表达式2
(完整word版)Mathematica习题
复习题
1、写出求极限1
1261lim 22+--∞→x x x x 的mathematica 程序和结果。 2、写出求函数2sin()arctan y x x =的导数的mathematica 程序和结果。
3、写出已知2sin )2cos(x x x y -=,写出求)4(y 的mathematica 程序和计算结果。
4、写出求dx e x x
⎰3的mathematica 程序和结果。
5、计算定积分dx x x
⎰--1145,写出mathematica 的程序和结果。
6、写出绘制函数12/)(23x e x x f x +=- 在区间[-3,3]的图形的mathematica
程序。
7、设计一段mathematica 程序,用一个while 循环实现求2215
1...211+++。请写出该程序。 8、设有函数⎪⎩⎪⎨⎧≥+<+=0
,3sin 0,1)(2x x e x x x f x ,试在mathematica 中定义该函数,并写出绘制该函数在区间[-5,5]上的图形,写出定义函数的程序和绘制图形的程序。
9、编写程序在同一窗口绘制曲线1sin(2x 0.3)y =-;23cos(x 0.5)y =+在
[0,2]π区间上的图像。
10、某药店对顾客所购买营养品实行打折销售,标准如下(药品价格用price 来表示):
200price < 没有折扣
200500price ≤< 3%折扣
5001000price ≤< 5%折扣
10002500price ≤< 8%折扣
(完整word版)ASCII码对照表
[保留]2
174
AE
[保留]2
175
AF
[保留]2
176
B0
°
温度符
177
B1
±
加/减号
178
B2
²
上标 2
179
B3
³
上标 3
180
B4
[保留]2
181
B5
µ
微符
182
B6
¶
段落符,pilcrow
183
B7
·
中点
184
B8
[保留]2
185
B9
¹
上标 1
186
BA
º
阳性顺序指示符
187
BB
»
右角引号
左大括号
124
7C
|
垂直线
125
7D
}
右大括号 (ALTMODE)
126
7E
~
代字号 (ALTMODE)
127
7F
DEL
擦掉 (DELETE)
控制字符
128
80
[保留]
129
81
[保留]
130
82
[保留]
131
83
[保留]
132
84
IND
索引
(完整word版)ASCII码对照表
10
DLE
数据链接丢失 (Ctrl/P)
17
11
DC1
设备控制 1 (Ctrl/Q)
18
12
DC2
设备控制 2 (Ctrl/R)
19
13
DC3
设备控制 3 (Ctrl/S)
20
14
DC4
设备控制 4 (Ctrl/T)
21
15
NAK
否定接受 (Ctrl/U)
22
16
SYN
同步闲置符 (Ctrl/V)
32
2
2
51
33
3
3
52
34
4
4
53
35
5
5
54
36
6
6
55
37
7
7
56
38
8
8
57
39
9
9
58
3A
:
冒号
59
3B
;
分号
60
3C
<
小于
61
3D
=
等于
62
3E
>
大于
63
3F
?
问号
ASCII 字母字符
64
40
@
商业 at 符号
65
(完整word版)数控技术指令集合
数控铣床
常用代码类型及其含义
G代码—-准备功能
F代码—-加工功能
M代码—-辅助功能
S代码——主轴功能
T代码——刀具功能
N代码—-顺序号
G指令
(1)G0 快速移动(定位)指令:
格式:G0 X__Y__Z__B__
刀具以快速移动速度从当前点运动到坐标指定的位置(快速移动速度为6米/分).
G0指令执行时,主轴可以不转动.
(2)G1 直线移动(线性插补)指令:
格式:G1 X__Y__Z__B__(F__S__M__)
刀具以给定的进给速度、转速,从当前点移动到坐标所指定的点。运动时4个坐标同时移动,同时到达终点.
G1指令运动时,主轴必须转动,主轴转动状态由M指令指定.
G1指令的目标点坐标也可用极坐标来表示.在绝对坐标(G90)下,用B2、L2代表与水平轴的夹角和长度;在相对坐标(G91)下,用B1、L1代表与水平轴的夹角和长度。
也可用X坐标和B1表示。
若前面已定了F、S、M,则此处可省略.
(3) G2、G3 顺时针、逆时针圆弧移动(圆弧插补)指令:
格式:G2(G3) X__(Y__)Z___I__(J__)K__(R__) (F__S__M__)
刀具从当前点顺(逆时钟)运动,以给定的圆弧中心坐标(I、K)或圆弧半径R,
移动到坐标给定的目标点.
若前面已定了F、S、M,则此处可省略。
刀具坐标使用X-Y组合或是X—Z组合,还是Y-Z组合,取决于所选定的工作平面(G17、G18、G19)。圆心坐标也相应为I—J、I—K和J—K.
(4) G17、G18、G19工作平面定义指令:
格式:G17(G18、G19) 无参数
(完整word版)计算机视觉期末复习
一、
1.什么是计算机视觉?理解计算机视觉问题的产生原理。
研究用计算机来模拟生物视觉功能的技术学科。具体来说,就是让计算机具有对周围世界的空间物体进行传感、抽象、分析判断、决策的能力,从而达到识别、理解的目的。
2.直方图的均衡化
处理的“中心思想”是把原始图像的灰度直方图从比较集中的某个灰度区间变成在全部灰度范围内的均匀分布。直方图均衡化就是对图像进行非线性拉伸,重新分配图像像素值,使一定灰度范围内的像素数量大致相同。直方图均衡化就是把给定图像的直方图分布改变成“均匀”分布直方图分布。
是将原图像通过某种变换,得到一幅灰度直方图更为均匀分布的新图像的方法。设图像均衡化处理后,图像的直方图是平直的,即各灰度级具有相同的出现频数,那么由于灰度级具有均匀的概率分布,图像看起来就更清晰了。
二、 1.常见的几何变换:平移T x 为点(x ,y )在x 方向要平移的量。 旋转 变尺度:x 轴变大a 倍,y 轴变大b 倍。 2.卷积掩膜技术:(,)(,)(,)(,)m n f i j h i m j n g m n =--∑∑ 对应相乘再相加 掩膜的有效应用——去噪问题
3. 均值滤波器(低通):抑制噪声
主要用于抑制噪声,对每一个目标像素值用其局部邻域内所有像素值的加权均值置换。con 命令
高斯滤波器:一个朴素的道理,距离目标像素越近的点,往往相关性越大,越远则越不相干。所以,高斯
滤波器根据高斯函数选择邻域内各像素的权值 medfilt1 。
区别方法是:高通滤波器模板的和为0,低通滤波器模板的和为1
常用的非线性滤波器:中值滤波;双边滤波;非局部滤波
(完整word版)FLAC3D5.0模型及输入参数说明
1.1模型参数代码
可参考manual中各个章节的command命令及说明,注意单位.用prop 赋值。
1.1.12
经典粘弹性模型的材料参数(Classical Viscoelastic (Maxwell Substance) –MODEL mechanical viscous)
1bulk弹性体积模量,K
2shear弹性剪切模量,G
3viscosity动力粘度,η
1.1.13粘弹性模型
粘弹性模型的材料参数(Burgers Model –MODEL mechanical burgers)
1bulk弹性体积模量,K
2kshear Kelvin弹性剪切模量,G K
3kviscosity Kelvin动力粘度,ηK
4mkshear Maxwell切边模量,G M
5mviscosity Maxwell动力粘度,ηM
1.1.14二分幂律模型
二分幂律模型的材料参数(Power Law –MODEL mechanical power)
1a_1常数,A1
2a_2常数,A2
3bulk弹性体积模量,K
4n_1指数,n1
5n_2指数,n2
6rs_1参考应力,σ1ref
7rs_2参考应力,σ2ref
8shear弹性剪切模量,G
1.1.15蠕变模型
蠕变模型材料参数(WIPP Model –MODEL mechanical wipp)
1act_energy活化能,Q
2a_wipp常数,A
3b_wipp常数,B
4bulk弹性体积模量,K
5d_wipp常数,D
6e_dot_star
临界稳定状态蠕变率,
7gas_c气体常数,R
(完整word版)Word通配符用法详解
Word通配符用法详解
1、任意单个字符:
“?”可以代表任意单个字符,输入几个“?”就代表几个未知字符。如:
输入“? 国”就可以找到诸如“中国”、“美国”、“英国”等字符;
输入“国”可以找到“孟加拉国”等字符。
2、任意多个字符:
“*”可以代表任意多个字符。如:
输入“*国”就可以找到“中国”、“美国”、“孟加拉国”等字符。
3、指定字符之一:
“[]”框内的字符可以是指定要查找的字符之一,如:
输入“[中美]国”就可以找到“中国”、“美国”。又如:
输入“th[iu]g”,就可查找到“thigh”和“thug”。
输入“[学硕博]士”,查找到的将会是学士、士、硕士、博士。
输入“[大中小]学”可以查找到“大学”、“中学”或“小学”,但不查找“求学”、“开学”等。
输入“[高矮]个”的话,Word查找工具就可以找到“高个”、“矮个”等内容。
4、指定范围内的任意单个字符:
“[x-x]”可以指定某一范围内的任意单个字符,如:
输入“[a-e]ay”就可以找到“bay”、“day”等字符,要注意的是指定范围内的字符必须用升序。用升序。如:
输入“[a-c]mend”的话,Word查找工具就可以找到“amend”、“bmend”、“cmend”等字符内容。
5、排除指定范、排除指定范围内的任意单个字符:
“[!x-x]”可以用来排除指定范围内的任意单个字符,如:
输入“[!c-f]”就可以找到“bay”、“gay”、“lay”等字符,但是不等字符,但是不会找到“cay”、“day”等字符。要注意范围必须用升序。
又如:
输入“[!a-c]”的话,word程序就可以找到“good”、“see”、“these”等目标字符,而所有包含字符a、b、c之类的内容都不会在查找结果中出现。
(完整word版)湘源控规快捷键以及相关操作
(完整word版)湘源控规快捷键以及相关操作
湘源控规的快捷命令
AAA:计算面积AD:选数求和BB:测量边界
BD:交点打断BG:标高标注CB:修改标高
CD:设当前层CG:图层相同CL:改实体色
CS:改线型比CW:改曲线宽CY:改线型
CT:修改内容DB:显示后置DF:显示前置
DD:关用地层DDD:按类显示EE:删同实体
FF:关指定层GG:图层相同HB:填充关联
FT:构选择集HHH:选择填充PPP:选多义线
TTT:选择文字II:对象查询JJ:线变复线
LL:查询LK:按次选线OL:图层孤立
PU:清理SS:设当前层SSS:按层显示
TF:边界打断TT:快速修剪VV:弯道圆角
WL:融合直线WWW:拷贝III:粘贴
PP:返回前图ZB:注坐标ZD:放大一倍
ZX:缩小一倍ZE:充满全屏ZZ:缩小0.8倍
湘源控规参数设置命令集合
一、高程分析或者坡度分析的面积统计问题
1、ft命令回车命令行:选择参照块(即在图例上选择一个色块)回车
2、命令行:构造实体集类型[0-同层实体1-同类实体2-同层及同类实体3-同色实体]<3>选择3回车
3、输入命令aaa 回车命令行:单位P/<输入字体高度><3.50>输入p回车选择面积的单位公顷、平方米或者亩回车输入字体的高度回车
4、命令行:选择计算方法[0-点选1-选实体2-描边界3-按次选线]<1>:选择1回车
5、命令行:选择数字、闭合曲线或填充图案: 输入p回车
6、命令行:请输入位置点: 在图上选择要输入的位置即可。
7、重复上面的工作可以把所有的不同颜色的高程区域或坡度区域的面积统计出来。(地形分析的高程分析图和坡度分析图的操作步骤:1地形-字转高程。把地形图的高程点转为湘源可识别的离散点;2地形-地形分析-高程分析或坡度分析选择对话框的等间距,确定即刻。即可绘出高程分析图或者是坡度分析图,图例也自动生成)
(完整word版)高斯平滑滤波器(含matlab代码)
Gaussian Smoothing Filter
高斯平滑滤波器
一、图像滤波的基本概念
图像常常被强度随机信号(也称为噪声)所污染.一些常见的噪声有椒盐(Salt & Pepper)噪声、脉冲噪声、高斯噪声等.椒盐噪声含有随机出现的黑白强度值.而脉冲噪声则只含有随机的白强度值(正脉冲噪声)或黑强度值(负脉冲噪声).与前两者不同,高斯噪声含有强度服从高斯或正态分布的噪声.研究滤波就是为了消除噪声干扰。
图像滤波总体上讲包括空域滤波和频域滤波。频率滤波需要先进行傅立叶变换至频域处理然后再反变换回空间域还原图像,空域滤波是直接对图像的数据做空间变换达到滤波的目的。它是一种邻域运算,即输出图像中任何像素的值都是通过采用一定的算法,根据输入图像中对用像素周围一定邻域内像素的值得来的。如果输出像素是输入像素邻域像素的线性组合则称为线性滤波(例如最常见的均值滤波和高斯滤波),否则为非线性滤波(中值滤波、边缘保持滤波等)。
线性平滑滤波器去除高斯噪声的效果很好,且在大多数情况下,对其它类型的噪声也有很好的效果。线性滤波器使用连续窗函数内像素加权和来实现滤波。特别典型的是,同一模式的权重因子可以作用在每一个窗口内,也就意味着线性滤波器是空间不变的,这样就可以使用卷积模板来实现滤波。如果图像的不同部分使用不同的滤波权重因子,且仍然可以用滤波器完成加权运算,那么线性滤波器就是空间可变的。任何不是像素加权运算的滤波器都属于非线性滤波器.非线性滤波器也可以是空间不变的,也就是说,在图像的任何位置上可以进行相同的运算而不考虑图像位置或空间的变化。
(完整word版)keil C语言与汇编语言混合编程
keil C语言与汇编语言混合编程
1. C语言中嵌入汇编
1、在C 文件中要嵌入汇编代码片以如下方式加入汇编代码:
#pragma ASM;
Assembler Code Here
#pragma ENDASM
2、在Project 窗口中包含汇编代码的C 文件上右键,选择“Options for ...”,点击右边的“Generate Assem bler SRC File”和“Assemble SRC File”,使检查框由灰色变成黑色(有效)状态;
3、根据选择的编译模式,把相应的库文件(如Small 模式时,是Keil\C51\Lib\C51S.Lib)加入工程中, 该文件必须作为工程的最后文件;
4、编译,即可生成目标代码。
来个实例吧:
#i nclude
void main(void)
{
P2=1;
#pragma asm
MOV R7,#10
DEL:MOV R6,#20
DJNZ R6,$
DJNZ R7,DEL
#pragma endasm
P2=0;
}
2 . 无参数传递的函数调用
C51调用汇编函数
1.无参数传递的函数调用
先来个例子:其中example.c和example.a51为项目中的两个文件
***********************example.c***********************************************
extern void delay100();
main()
{delay100;}
***********************example.a51***********************************************
测井曲线代码大全word版本
测井曲线代码大全
测井曲线代码
RD、RS—深、浅侧向电阻率
RDC、RSC—环境校正后的深、浅侧向电阻率VRD、VRS—垂直校正后的深、浅侧向电阻率DEN—密度
DENC—环境校正后的密度
VDEN—垂直校正后的密度
CNL—补偿中子
CNC—环境校正后的补偿中子
VCNL—垂直校正后的补偿中子
GR—自然伽马
GRC—环境校正后的自然伽马
VGR—垂直校正后的自然伽马
AC—声波
VAC—垂直校正后声波
PE—有效光电吸收截面指数
VPE—垂直校正后的有效光电吸收截面指数SP—自然电位
VSP—垂直校正后的自然电位
CAL—井径
VCAL—垂直校正后井径
KTh—无铀伽马
GRSL—能谱自然伽马
U—铀
Th—钍
K—钾
WCCL—磁性定位
TGCN—套管中子
TGGR—套管伽马
R25—2.5米底部梯度电阻率
VR25—环境校正后的2.5米底部梯度电阻率DEV—井斜角
AZIM—井斜方位角
TEM—井温
RM—井筒钻井液电阻率
POR2—次生孔隙度
POR—孔隙度
PORW—含水孔隙度
PORF—冲洗带含水孔隙度
PORT—总孔隙度
PERM—渗透率
SW-含水饱和度
SXO—冲洗带含水饱和度
SH—泥质含量
CAL0—井径差值
HF—累计烃米数
PF—累计孔隙米数
DGA—视颗粒密度
SAND,LIME,DOLM,OTHR—分别为砂岩,石灰岩,白云岩,硬石膏含量VPO2—垂直校正次生孔隙度
VPOR—垂直校正孔隙度
VPOW—垂直校正含水孔隙度
VPOF—垂直校正冲洗带含水孔隙度
VPOT—垂直校正总孔隙度
VPEM—垂直校正渗透率
VSW-垂直校正含水饱和度
VSXO—垂直校正冲洗带含水饱和度
(完整word版)模式作业设计
(完整word 版)模式作业设计
1.1题目
给出K —均值算法的程序框图,编写程序,自选一组分别属于三类的三维模式样本(图2.13),并给它们进行聚类分析。
1。2算法原理
K-均值算法也称C —均值算法,是根据函数准则进行分类的聚类算法,基于使聚类准则最小化.此处所用的聚类准则函数是聚类集中每一个样本点到该类聚类中心的距离平方和,对于第j 个聚类集,准则函数定义为
式中,S j 表示第j 个聚类集,也称聚类域,其聚类中心为j Z ;j N 为第j 个聚类集j S 中所包含的样本个数。
(1)任选K 个初始聚类中心Z 1(1),Z 2(1),…,Z K (1),K N <.括号内的序号代表寻找聚类中心的迭代运算的次序号。一般可选择样本集中前K 个样本作为初始聚类中心.
(2)按最小距离原则将其余样本分配到K 个聚类中心中的某一个,即:
若min{(),1,2,K}()()i j X Z k i X Z k D k -=⋅⋅⋅=-=,则i j X S ∈.式中,k 代表迭代运算次序号;K 代表聚类中心的个数。
(3)计算各个聚类中心的新向量值(1)j Z k +,1,2,,j K =⋅⋅⋅
即以均值向量作为新的聚类中心。这一步要分别计算K 个聚类中心的样本均值向量,故该算法被称为K-均值算法.
(4)如果(1)()j j Z k Z k +≠,1,2,,j K =⋅⋅⋅,则回到(2),将模式样本逐个重新分类,并重复迭代计算;如果(1)()j j Z k Z k +=,1,2,,j K =⋅⋅⋅,算法收敛,计算完毕。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
a4=(0.25+0.00252*cB2)*cB2-0.04166;
a6=(0.166*cB2-0.084)*cB2;
a3=(0.3333333+0.001123*cB2)*cB2-0.1666667;
a5=0.0083-(0.1667-(0.1968+0.0040*cB2)*cB2)*cB2;
// x=X+N*sin(B)*cos(B)*l*l/2+N*sin(B)*pow(cos(B),3)*(5-t*t+9*n2+4*n2*n2)*pow(l,
4)/24+N*sin(B)*pow(cos(B),5)*(61-58*t*t+pow(t,4))*pow(l,6)/720;
y=(1+(a3+a5*l*l)*l*l)*l*N*cos(B);
printf("输出高斯邻带坐标换算\n");
printf("x=%f\ny=%f\n",x,y);
}
void main()
{
int i;
printf("输入1 2 3选择投影的计算方法\n");
scanf("%d",&i);
switch (i)
// y=N*cos(B)*l+N*pow(cos(B),3)*(1-t*t+n2)*pow(l,5)/6+N*pow(cos(B),5)*(5-18*t*t
+pow(t,4)+14*n2-58*n2*t*t)*pow(l,5)/120;
x=6367558.4969*B-(a0-(0.5+(a4+a6*l*l)*l*l)*l*l*N)*sin(B)*cos(B);
l=L-111*3600/P // l=((m%6)*3600+n*60+h)/P;
t=tan(B);
n2=e2*cos(B)*cos(B);
V=sqrt(1+n2);
cB2=pow(cos(B),2);
N=6399698.902-(21562.267-(108.973-0.612*cB2)*cB2)*cB2; // N=c/V;
double a0,a4,a6,a3,a5,cB2;
e2=0.006738525414683;
c=6399698.901782711;
// printf("输入x y的值\n");
// cout<<"x=";
// cin>>x;
// cout<<"y=";
// cin>>y;
// scanf("x=%f",&x); ///////////// ?????
{
case 1: GaosZ_fun(); break;
case 2: GaosF_fun(); break;
case 3: GaosLT_fun(); break;
default : break;
}
}
b3=0.333333-(0.166667-0.001123*cos(Bf)*cos(Bf))*cos(Bf)*cos(Bf);
b4=0.25+(0.16161+0.00562*cos(Bf)*cos(Bf))*cos(Bf)*cos(Bf);
b5=0.2-(0.167-0.0088*cos(Bf)*cos(Bf))*cos(Bf)*cos(Bf);
// y=N*cos(B)*l+N*pow(cos(B),3)*(1-t*t+n2)*pow(l,5)/6+N*pow(cos(B),5)*(5-18*t*t
+pow(t,4)+14*n2-58*n2*t*t)*pow(l,5)/120;
x=6367558.4969*B-(a0-(0.5+(a4+a6*l*l)*l*l)*l*l*N)*sin(B)*cos(B);
a0=32140.404-(135.3302-(0.7092-0.004*cB2)*cB2)*cB2;
a4=(0.25+0.00252*cB2)*cB2-0.04166;
a6=(0.166*cB2-0.084)*cB2;
a3=(0.3333333+0.001123*cB2)*cB2-0.1666667;
// scanf("x=%f",&x);
// scanf("y=%f",&y);
x=3380330.875;
y=320089.976;
L0=111;
b=x/6367558.4969;
Bf=b+(50221746+ ( 293622+ (2350+22*cos(b)*cos(b))*cos(b)*cos(b) )
double l,L,B,n2,x,y,N,t,V,c,e2;
double i,j,k,n,h,a0,a4,a6,a3,a5,cB2;
int m;
e2=0.006738525414683;
c=6399698.901782711;
B=17.33557339*3600/P;
L=119.15521159*3600/P;
y=(1+(a3+a5*l*l)*l*l)*l*N*cos(B);
printf("x=%f\ny=%f\n",x,y);
}
void GaosF_fun()
{
printf("高斯投影的反算\n");
double B,Bf,Nf,b,b2,b3,b4,b5,Z,x,y,L0,l;
// printf("输入x y的值\n");
b3=0.333333-(0.166667-0.001123*cos(Bf)*cos(Bf))*cos(Bf)*cos(Bf);
b4=0.25+(0.16161+0.00562*cos(Bf)*cos(Bf))*cos(Bf)*cos(Bf);
b5=0.2-(0.167-0.0088*cos(Bf)*cos(Bf))*cos(Bf)*cos(Bf);
l=(1-(b3-b5*Z*Z)*Z*Z)*Z;
L=l+L01; ///
反算就出
L B
l=L-L02;
B=Bf-(1-(b4-0.12*Z*Z)*Z*Z)*Z*Z*b2;
t=tan(B);
n2=e2*cos(B)*cos(B);
V=sqrt(1+n2);
cB2=pow(cos(B),2);
N=6399698.902-(21562.267-(108.973-0.612*cB2)*cB2)*cB2; //N=c/V;
*cos(b)*cos(b)) *sin(b)*cos(b)*1e-10;
Nf=6399698.902-(21562.267-(108.973-0.612*cos(Bf)*cos(Bf))*cos(Bf)*cos(Bf))*cos(
Bf)*cos(Bf);
Z=y/Nf/cos(Bf);
b2=(0.5+0.003369*cos(Bf)*cos(Bf))*sin(Bf)*cos(Bf);
a5=0.0083-(0.1667-(0.1968+0.0040*cB2)*cB2)*cB2;
// x=X+N*sin(B)*cos(B)*l*l/2+N*sin(B)*pow(cos(B),3)*(5-t*t+9*n2+4*n2*n2)*pow(l,
4)/24+N*sin(B)*pow(cos(B),5)*(61-58*t*t+pow(t,4))*pow(l,6)/720;
// scanf("y=%f",&y);
// x=3380330.875;
// y=320089.976;
x=1Leabharlann Baidu44359.607;
y=240455.4563;
L01=117*3600/P;
L02=120*3600/P;
b=x/6367558.4969;
Bf=b+(50221746+ ( 293622+ (2350+22*cos(b)*cos(b))*cos(b)*cos(b) )
B=Bf-(1-(b4-0.12*Z*Z)*Z*Z)*Z*Z*b2;
l=(1-(b3-b5*Z*Z)*Z*Z)*Z;
printf("B=%f\nL=%f\n",B*180/PI,(L0+l*180/PI));
}
void GaosLT_fun()
{
double B,Bf,Nf,b,b2,b3,b4,b5,Z,x,y,L01,L02,n2,l,V,L,N,c,t,e2;
*cos(b)*cos(b)) *sin(b)*cos(b)*1e-10;
Nf=6399698.902-(21562.267-(108.973-0.612*cos(Bf)*cos(Bf))*cos(Bf)*cos(Bf))*cos(
Bf)*cos(Bf);
Z=y/Nf/cos(Bf);
b2=(0.5+0.003369*cos(Bf)*cos(Bf))*sin(Bf)*cos(Bf);
#include "stdafx.h"
#include "iostream.h"
#include "math.h"
#include "stdio.h"
#define P 206264.806247096355
#define PI 3.141592653589793
void GaosZ_fun()
{
printf("高斯投影的正算\n");