未名博客 - 未名空间
当前在线人数10057
首页 - 博客首页 - let comch pute be aos me you [博客首页] [博客论坛] [博客搜索]
欢迎访问heteroclinic的博客 - let comch pute be aos me you
[收藏该博客]
turb deter the ulence mine diaby EgramR
作者介绍
heteroclinic
上站次数:3226
经验值 :[14622](元始天尊)
表现值 :[26](还不错)
创建博客时间:2008-12-15
共发表日志: 522篇
栏目分类
每月档案
2020年09月 2020年8月
2020年7月 2020年6月
2020年5月 2020年4月
2020年3月 2020年2月
2020年1月 2019年12月
2019年11月 2019年10月
2019年9月 2019年8月
2019年7月 2019年6月
2019年5月 2019年4月
2019年3月 2019年2月
2019年1月 2018年12月
2018年11月 2018年10月
2018年9月 2018年8月
2018年7月 2018年6月
2018年5月 2018年4月
2018年3月 2018年2月
2018年1月 2017年12月
2017年11月 2017年10月
2017年9月 2017年8月
2017年7月 2017年6月
2017年5月 2017年4月
2017年3月 2017年2月
2017年1月 2016年12月
2016年11月 2016年10月
2016年9月 2016年8月
2016年7月 2016年6月
2016年5月 2016年4月
2016年3月 2016年2月
2016年1月 2015年12月
2015年11月 2015年10月
2015年9月 2015年8月
2015年7月 2015年6月
2015年5月 2015年4月
2015年3月 2015年2月
2015年1月 2014年12月
2014年11月 2014年10月
2014年9月 2014年8月
2014年7月 2014年6月
2014年5月 2014年4月
2014年3月 2014年2月
2014年1月 2013年12月
2013年11月 2013年10月
2013年9月 2013年8月
2013年7月 2013年6月
2013年5月 2013年4月
2013年3月 2013年2月
2013年1月 2012年12月
2012年11月 2012年10月
2012年9月 2012年8月
2012年7月 2012年6月
2012年5月 2012年4月
2012年3月 2012年2月
2012年1月 2011年12月
2011年11月 2011年10月
2011年9月 2011年8月
2011年7月 2011年6月
2011年5月 2011年4月
2011年3月 2011年2月
2011年1月 2010年12月
2010年11月 2010年10月
2010年9月 2010年8月
2010年7月 2010年6月
2010年5月 2010年4月
2010年3月 2010年2月
2010年1月 2009年12月
2009年11月 2009年10月
2009年9月 2009年8月
2009年7月 2009年6月
2009年5月 2009年4月
2009年3月 2009年2月
2009年1月 2008年12月
友情链接
XML
[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 ]
 
2010-10-24 [发表评论] [写信问候]
  继续我们计算non-prime number 的探险

先前提到过 y = sin(a * pi /x), 对于质数a 在(1,a)是没有解得。

写了段程序,用牛顿法算 x, 跑了跑发现,原来这个一维系统有很多folds,又想了想能不能用二维的系统来解这个问题呢?于是又了如下的构造
f = (f1,f2)^t
f1(x1,x2) = (x1 + exp(x2))* sin(pi * a /x1);
f2(x1,x2) = (x1*x1 -x2*x2)
求 f --〉0
这个系统看起来有个不错的jacobian,跑了跑程序,居然收敛了,但不知道收敛到哪里去了,您有时间帮忙看看吧,能不能有更好的构造。或者用什么论证着干脆就是不可行的,免得骑自行车去月球。

附程序c++


#include <iostream>
#include <math.h>
using namespace std;
#define PI (atan(1.0)*4.0)
#define __VERYSM


...阅读全文

heteroclinic 发布于2010-10-24 01:51:34  |  浏览[344]  |  评论[0]
 
2010-10-21 [发表评论] [写信问候]
  compute non-prime number

这两天胡乱想到一个问题,放到这里请教专家,我自己分叉太多还是老老实实写程序吧。

正整数的a的特性是sin(a PI/ x) = 0 if x is its divisor.

那么用taylor 把 sin(a PI/ x) 展开,有什么办法证明r(x)收敛还是不收敛。

如果收敛的话sin(a PI/ x) taylor 的 特征根如果有实数的话,则必定是整数,且未a 的因子。处吗?

由此可以判断a 是质数还是合数?

提示: 本博文来自于 Mathematics 版

heteroclinic 发布于2010-10-21 16:47:31  |  浏览[453]  |  评论[0]
 
2010-10-10 [发表评论] [写信问候]
  我现在在搞一个解码器,急需多项式eigenvalue 求法

1。不能用matlab
2。可以用fortran c/c++

比方 x^2 + m x + n = 0
矩阵
| m/2 sqrt(m^2/4 -n)|
|sqrt(m^2/4 -n) m/2 |
的eigenvalues 应该就是 x^2 + m x + n = 0 的 eigenvalues.
如果polynomial 的维数高了,怎么得到这样矩阵,因为矩阵式是numerical friendly 的。



提示: 本博文来自于 Mathematics 版

heteroclinic 发布于2010-10-10 23:22:26  |  浏览[702]  |  评论[1]
 
2010-08-27 [发表评论] [写信问候]
  有没有这样的算法

既定的条件:3d 空间里的一个solid,of course closed connected, bounded by 2d manifold (不是故意卖弄英文,这点儿都是用英语学的,的确不知道中文怎么说),this 2-manifold of course can be presented as a triangulated surface.

问题:我们用球体填充这个体,最大半径不限,最小半径要符合给定的精度要求,
那么
问题1. 是否存在最少球体的数量,从最大的球开始放,直到不能放进最小的球。
问题2. 如果是椭球呢?
问题3. 我们实在太需要一个方案,那么就限定球体的半径,那么如何把solid 充填满?

提示: 本博文来自于 CS 版

heteroclinic 发布于2010-08-27 20:43:51  |  浏览[662]  |  评论[0]
 
2010-08-14 [发表评论] [写信问候]
  如何处理这样的精度?

如何处理这样的精度?
double r1 = 1.0*2.0 + 1.0*2.0 + 1.0 *2.0;
double r10 = r1/sqrt(1.0+1.0+1.0)/sqrt(2.0*2.0+2.0*2.0+2.0*2.0);

VC++ 2010 express, Vista 64 结果
r10 = 1.0000000000000002;
求acos(r10)的时候溢出,用float也一样戴一个尾巴。

这样的运算怎样避免溢出?
非express 版也有这个问题吗?
没有试别的compler.


提示: 本博文来自于 CS 版

heteroclinic 发布于2010-08-14 20:52:07  |  浏览[464]  |  评论[0]
 
[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 ]

访问量 222354 更新时间: 2020-09-28 02:36:59
©All Rights Reserved http://blog.mitbbs.com/heteroclinic

 
Site Map - Contact Us - Terms and Conditions - Privacy Policy

版权所有BBS 未名空间站(mitbbs.com) since 1996