Auringon korkeuden yhtälö

Mä tein tämmösen aurinkopaneelin kääntelyä varten. En vaan osaa lähettää noita kulmia pc:ltä aaveeärrälle, kun oon niin tumpelo.

Turha tämmösen kanssa on kyllä mitään sekunteja ja asteen osia kyttäillä...

#include    "stdafx.h"
#include    "stdio.h"
#include    "stdlib.h"
#include    "windows.h"
#include    "conio.h"
#include    "math.h"
#include    "ctime"
using namespace std;

int main ()
{
        float pi = 3.14159265;
        float latitude = 0;
        float lognitude = 0;
        float x_angle=0;
        float y_angle=0;

        printf("Latitude: ");
        scanf("%f",&latitude);
        printf("Lognitude: ");
        scanf("%f",&lognitude);

    while(1)
    {       
        time_t now = time(0);
        tm *ltm = localtime(&now);

        float day;
        int this_day = ltm->tm_mday;                            // tämä päivä
        int this_month = ltm->tm_mon + 1;                    // tämä kuukausi
        int this_hour = ltm->tm_hour;                            // tämä tunti
        float this_min = ltm->tm_min;                            // tämä minuutti
        float this_sec = ltm->tm_sec;                            // tämä sekunti

        /* tunteja tänään kulunut */
        float hours = this_hour + this_min / 60 + this_sec / 3600;   
        /* korjaus normaaliaikaan tunteina */
        float fix = (30 - lognitude) * 0.06666667;               

        struct tm setdate;
        setdate = *localtime(&now);

        setdate.tm_hour = 13;                                    // asetettu aika kevätpäivän tasaus
        setdate.tm_min = 41;
        setdate.tm_sec = 0;
        setdate.tm_mon = 2;                                      // maaliskuu
        setdate.tm_mday = 20;                                    // 20 pv.

        double seconds = difftime(now,mktime(&setdate));                // sekunnit asetetusta ajasta

        printf("%.1f %.1f    ", latitude, lognitude);
        printf("%2d.%2d  ", this_day, this_month);               
        printf("%2d.%2.f    ", this_hour, this_min);           

        day = seconds / 86400;                                    // monesko päivä sitten asetetun ajan
       
        printf ("day %.f ", day);
        printf(" hours %.2f    ", hours);

        y_angle = (-cos(2*pi*(hours-fix)/24)) * ((180-2*latitude)/2) + (sin((2*pi/365)*day)*23.44);

        x_angle = ((hours-fix)*15-180);

        printf("y = %.1f  ", y_angle);
        printf("x = %.1f\n", x_angle);

        Sleep(10000);
    }
  return 0;
}

sieniukko

15.03.2007 klo 00:01

Päiväntasaajalla auringon kuuluisi tasauspäivänä nimittäin oikeasti edetä taivaan halki tasaisesti 1 aste / 4min.

Aurinko vaeltaa kyllä taivaalla tuota vauhtia mutta ei leveyspiireihin nähden vaan pituuspiireihin nähden Maan pyörimisen johdosta. Kyllä aurinko etenee vuorokauden aikana kevät- tai syyspäivän tasauksen aikaan hyvin pienen kulman päiväntasaajaan nähden, lienee vain jokunen kulmaminuutti.