日志
- 分享 圆周率的分数近似
- pi =1881244168 / 598818617 ~3.141592653589793117604424780267 pi =245850922 / 78256779~ 3.1415926535897931602832771842041
- 分享 勾股数的python代码
- 在CSDN看到一个研究够股数的帖子，感觉他的代码麻烦，就写了以下代码 y = 1 while True: a = 2 * y + 1 b = 2 * y * ( y + 1 ) c = b + 1 print a, b, c ...
- 分享 无限长导体方柱散射的计算
-
TM模式，所谓2维散射，导体截面边长3个波长 距离导体柱0.1个波长电场强度【绝对值】 圆圈是金建铭书上的数据，蓝线是算得的数据，电子书上的图形不是很规则，但是可以看出，数据复合的非常好 电场3D分布【绝对值】 实部3D 虚部3D 边界积分+某种算法 下面图形是今天晚上 用MOM算的 ...
- 分享 波前法
- 波前法是一种利用较小内存求解大型线性方程组的算法。其本质是分块高斯消去法的更灵活应用，它不形成体系总刚度矩阵，而只是形成一个波前内相关单元的“分块刚度阵”，分解后即记入硬盘；依次波前在整个结构中遍历即完成体系的分块总刚的形成和消元，回代时逆序进行即可。国内有文献根据这个特点，将波前法（front ...
- 分享 稀疏矩阵笔记
-
interl编译器，32400未知量，非零元数目1578197，Intel酷睿2，1.86GHz interl编译器，32400未知量，非零元数目1578197， AMD 620CPU4核 同上 Intel 酷睿2，1.86GHz，2G内存，VS2008 C/C++编译器 Intel 酷睿2，1.86GHz，2G内存，GNU GCC C/C++编译器 GNU GCC C/C++编译器, AMD 620CPU4核 Spar ...
- 分享 今天终于自己写了一个稀疏矩阵求解器
-
今天终于自己写了一个稀疏矩阵求解器，数据是链表【不是十字链表】存储，方法是高斯消元回带。性能如何有待于测试，回家测试一下。 2011-08-16:19:26 虽然可以运行，但是或许是由于链表操作过于频繁，性能并不好。或许还有改进的余地。 2011-08-17:20:40 稀疏矩阵求解器可以工作，7200未知数工作正常，但是求解时间比原 ...
- 分享 电波传播研究70年
- 电波传播研究70年 沙 踪 （中国电波传播研究所，北京 100040） 摘 要 回顾了我国电波传播研究从1936年到2007年，约70年的发展概况。全文分为我国学者早期的一些出色研究贡献、电离层结构探测、对流层结构探测、电离层传播特性研究、对流层传播特性研究、长波和超长波传播研究、地面地下和水面水下传播特性研究、特殊媒 ...
- 分享 翻硬币问题
- 翻硬币问题 翻硬币问题有好几种。 其中的一种是这样的： 桌子上有q = m + n枚硬币，m正面朝上，n枚反面朝上，每一轮翻p枚，在每一轮翻币的时候，被翻的同一枚硬币只能翻一次。问最少多少次能把所有的硬币翻成全部正面或者反面朝上？根据问题的描述，问题实际上隐含：m 、n、 p 0，m + n p。 这个问题往往是计 ...
- 分享 一个数值逼近问题
- 这是在www.CSDN.net的一个问题： 写一个函数，有三个输入参数，first, second, diff. 要求函数能计算出一个数值result，这个result应该是second的最小整数倍，且与first的整数倍的差的绝对值在diff以内的， 比如，first为3.6，second为35，diff为0.5的情况下，那么result就应该为140. 因为 140 = 35（second）*4; ...
- 分享 环形区域拉普拉斯/泊松方程求解
-
用Matlab很容易求解，不过下面的结果不是用Matlab做的，也不是用FD做的，方形区域FD做起来很方便，精度也较高（h^2）,下面的是用高阶方法做的，比FD麻烦很多，主要是图形方面 中心区域占1/3，挖空，泊松方程右端 = -1 中心区域占1/3，挖空，内边界边值 = 1，外边界边值 = 0，泊松方程右端 = -20 中心区域占 ...
