Math

格林威治时间转换为北京时间（单片机+GPS时使用）

杨老师的茅屋 — Sun, 27 Apr 2008 14:54:00 GMT

// 格林威治时间转换为北京时间（2000年以后）
// 输入参数年(2位)，月，日，时，分,秒
void Calc( int y, int m, int d, int hh, int mm, int ss )
{
    int w;        // 星期
    hh += 8;    // 格林威治时间 + 8 小时 = 北京时间

    if( hh < 24 )    goto l_out;        // 没有跨天，则计算完成

    ///////////// 下面是跨天后的计算 /////////////////////////////
    hh -= 24;
    d ++;        // 日期加一天

    switch( m )        // 按月判断
    {
    case 4:        case 6:                // 跨小月的判断
    case 9:        case 11:
        if( d > 30 )
        {
            d = 1;    m ++;
        }
        break;

    case 1:        case 3:                // 跨大月的判断
    case 5:        case 7:
    case 8:        case 10:
        if( d > 31 )
        {
            d = 1;    m ++;
        }
        break;

    case 12:                        // 12 月，要判断是否跨年
        if( d > 31 )
        {
            y ++;    d = 1;    m = 1;
        }
        break;

    case 2:                            // 2 月，要判断是否是闰年
        if( ( (y+2000)%400 == 0 ) ||                        // 能被400整除，一定是闰年
            ( (y+2000)%4 ==0 ) && ( (y+2000)%100 !=0 ) )    // 能被4整除，但不能被100整除，一定是闰年
        {
            if( d>29 )    // 闰年2月，可以有29号
            {
                m = 3;    d = 1;
            }
        }
        else if( d>28 )    // 非闰年2月，可以有28号
        {
            m = 3; d = 1;
        }
        break;
    }

l_out:    // 计算完成，开始输出
    printf( "%04d.%02d.%02d - %02d:%02d:%02d", y+2000, m, d, hh, mm, ss);
    if( 1 == m )
    {
        y --;
        m = 13;
    }
    else if( 2 == m )
    {
        y --;
        m = 14;
    }

    w = y + y/4 + 26 * ( m + 1 ) / 10 + d - 1;
    w %= 7;        // 0 表示星期日

    printf ( " week=%d ", w );
}

闰年的计算

杨老师的茅屋 — Sun, 15 Oct 2006 18:41:00 GMT

两天前我解答 [shutle] 网友计算闰年方法的时候，回答错误，在这里一并更正并致歉。后来查阅了资料，现把计算方法和原理粘贴一下：

我们知道地球自转一圈为1天，地球围绕太阳公转一圈为1年。那1年到底等于多少天哪？答案是：1年=365天5小时48分46秒，即 365.242199[074] 天(注：用[xxx]表示无限循环)

因此，4年 = 4*365.242199[074] = 4*365 + 0.9687[962] 天
显然，每隔4年，就要多出了将近1天，因此公元历法就定义㈠：每4年加一天，即闰年(366天)
那么，也就是说4年就多算了：1-0.9687[962] = 0.0312[037]天

        照这样算下去，100年中有25次闰年，那么就会多计算了 25*0.0312[037]=0.7800[925] 天，这又太接近1天了，因此公元历法就定义㈡：每100年的那年不算闰年。
于是：
    100年中就又会少算 1-0.7800[925]=0.2199[074] 天
    200年会少算 2*0.2199[074] = 0.439[814] 天
    300年会少算 3*0.2199[074] = 0.6597[2] 天
    400年会少算 4*0.2199[074] = 0.87[962] 天

停！400年后的误差又接近1天了，于是公元历法进行修正㈢：400年的那年(虽然是100年的倍数)，但又算闰年啦。于是，400年的时候，由于多算了一次闰年，那么，
每隔400年就会多算：1-0.87[962] = 0.12[037] 天

显然，等到4000年的时候（那时，看到这个帖子的人都死了2000年了），误差又到1天多了，也许那个时候公元历法有会进行修改了“每4000年少一个闰年”。哈哈，反正不管我们的事了。但也不尽然，因为地球的自转很不稳定，但总的趋势是每年要慢一丁点（大家是否还记得今年2006年元旦，全世界的人都要把表拨慢1秒钟？）由于阴差阳错，也许4000年后并不需要少算一个闰年也不一定呐。总结一下闰年的算法：
㈠每4年一闰㈡每100年的那年不闰㈢每400年的那年又闰
真罗嗦，用咱们IT人的话说，就是能被4整除但不能被100整除，可是却能被400整除的年是闰年

我所见过的最牛的，但我看不懂的程序

杨老师的茅屋 — Fri, 22 Sep 2006 19:13:00 GMT

void main()
{
    long a=10000,b=0,c=2800,d=0,e=0,f[2801],g=0;
    for( ; b-c; )
     f[ b++ ]=a/5;
    for( ;d=0, g=c*2; c-=14, printf( "%.4d",e+d/a ), e = d%a )
        for( b=c; d += f[b] * a, f[b] = d%--g, d /= g--, --b; d *=b );
}
运行一下吧，猜猜是什么结果。

算法征集

杨老师的茅屋 — Sun, 02 Apr 2006 03:46:00 GMT

本人在一个工程任务中,需要完成“在9个整数中，找到中值”。
比如：1,3,5,7,9,2,4,6,8 的中值是5。

看似非常简单的一个小问题，但由于每秒钟需要调用大约10万多次，因此我需要一个快速算法，不考虑空间浪费，只要求比较语句尽可能的少。

有好的想法吗？欢迎抛砖引玉......

YUV 与 RGB 的相互转换

杨老师的茅屋 — Wed, 28 Dec 2005 12:01:00 GMT

Y =   0.299R + 0.587G + 0.114B
U = -0.1684R - 0.3316G +    0.5B + 128
V =     0.5R - 0.4187G - 0.0813B + 128
--------------------------------------------------
保留小数后6位：
R = Y - 0.000037U + 1.401988V - 179.449700
G = Y - 0.344078U - 0.713885V + 135.419284
B = Y + 1.771795U - 0.001259V - 226.628590
--------------------------------------------------
保留小数后2位：
R = Y          + 1.4V - 179.45
G = Y - 0.34U - 0.71V + 135.42
B = Y + 1.77U         - 226.63

友情提示：彩色转灰度的算法

杨老师的茅屋 — Thu, 24 Jun 2004 07:52:00 GMT

一个最基本的分形图象。曼德布罗集

杨老师的茅屋 — Tue, 15 Jun 2004 08:22:00 GMT