±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 350  |  »Ø¸´: 1

wdphysics

½ð³æ (СÓÐÃûÆø)

[½»Á÷] ¡¾ÇóÖú¡¿ÈçºÎµ÷Óô˺¯Êý·Ö½â¾ßÌ徨Õó£¿ ÒÑÓÐ1È˲ÎÓë

¶Ô³Æ¾ØÕóµÄcholesky·Ö½â£¬aΪ¾ØÕó£¬nΪÆäάÊý¡£Ö÷ÒªÊDz»ÖªµÀÔõôŪµ÷Óýӿڣ¡
void choldc(float **a, int n, float p[])
{
        void nrerror(char error_text[]);
        int i,j,k;
        float sum;
        for (i=1;i<=n;i++) {
                for (j=i;j<=n;j++) {
                        for (sum=a[j],k=i-1;k>=1;k--) sum -= a[k]*a[j][k];
                        if (i == j) {
                                if (sum <= 0.0)
                                        nrerror("choldc failed";
                                p=sqrt(sum);
                        } else a[j]=sum/p;
                }
        }
}
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

yalefield

½ð³æ (ÎÄ̳¾«Ó¢)

ÀϺºÒ»Ã¶

¡ï ¡ï ¡ï ¡ï
Сľ³æ(½ð±Ò+0.5):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû½»Á÷
ÓàÔó³É(½ð±Ò+3):ÐÁ¿àÁË£¡ 2010-08-20 12:18:33
£¨1£©
   aÊÇÖ¸ÕëµÄÖ¸Õ룬ÔÚÕâÀËü¾ÍÊÇÖ¸ÏòÒ»¸ö2άµÄÊý×é
   C/C++ÀïÃæ£¬Êý×é¶¼ÊÇһάµÄ£¬Óöµ½¶þάÊý×飬¾Í°ÑËü¡°Õ¹¿ª¡±³ÉһάµÄ¡£
   Òò´Ë£¬ÄãÒªÖ¸¶¨¶þάÊý×éÖУ¬ÄǸöµÚһάÓжàÉÙ¸öÔªËØ£¬Ò²¾ÍÊDzÎÊýn
     µÚ0µ½µÚn-1¸ö£¬ÊǶþάÊý×éÖеÄǰn¸ö£¬Ò²¾ÍÊǵÚÒ»ÐУ¨»òÁУ©
   pÔòÊÇһάÊý×é

£¨2£©
   Òò´Ë£¬Òªµ÷ÓÃÕâ¸öº¯Êý£¬ÐèҪ׼±¸ºÃÒ»¸ö¶þάÊý×éºÍÒ»¸öһάÊý×é
   ÀýÈ磺
   int f ()
    {
        float X[100][200];
        float Y[100];
        return choldc(X, 100, Y);
    }

£¨3£©
   °´ÕÕĿǰµÄд·¨£¬Ó¦¸ÃÊÇ£ºint choldc(float **a, int n, float *p)
   È»ºó£¬·µ»Ø1±íʾ³É¹¦£¬·µ»Ø0±íʾ³ö´í£¨»òÕßÏà·´£©¡£
   ÔÚÕâÑùµÄº¯ÊýÀֱ½Óµ÷ÓôíÎóÊä³ö£¬ÊµÔÚËã²»ÉϵĵÂÒÕ˫ܰ¡£
CODE:
int choldc(float **a, int n, float p[])
{
        int i, j, k;
        float sum = 0.0F;
        for (i=1; i<=n; i++)
        {
                for (j=i; j<=n; j++)
                {
                        for (sum=a[i][j], k=i-1; k>=1; k--)
                        {
                                sum -= a[i][k]*a[j][k];
                        }
                        if (i == j)
                        {
                                if (sum <= 0.0)
                                {
                                        return 0;
                                }
                                p[i]=sqrt(sum);
                        }
                        else
                        {
                                a[j][i]=sum/p[i];
                        }
                }
        }
        return 1;
}

-------------------
CÓïÑÔ¾ÅÐÇ´óʦ£¬Å¼·¢ÈýË×´úÂë
2Â¥2010-08-20 11:20:52
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ wdphysics µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 326Çóµ÷¼Á +10 áÌáÌ×Ð 2026-04-02 10/500 2026-04-03 09:08 by ÅÁ¶ûÂíÀ­ÌØ
[¿¼ÑÐ] 306·Ö²ÄÁÏÓ뻯¹¤Çóµ÷¼Á +3 Àè°ÉÀ²À²ÄãºÜÓÐà 2026-04-03 3/150 2026-04-03 08:14 by laura9751
[¿¼ÑÐ] 372·Ö²ÄÁÏÓ뻯¹¤£¨085600£©Ò»Ö¾Ô¸ºþÄÏ´óѧÇóµ÷¼Á +5 À¶¼ãƬ 2026-04-02 6/300 2026-04-02 21:37 by dongzh2009
[¿¼ÑÐ] 070305¸ß·Ö×Ó»¯Ñ§ÓëÎïÀí 304·ÖÇóµ÷¼Á +16 c297914 2026-03-28 16/800 2026-04-02 20:50 by dongzh2009
[¿¼ÑÐ] 266·Ö£¬Ò»Ö¾Ô¸µçÆø¹¤³Ì£¬±¾¿Æ²ÄÁÏ£¬Çó²ÄÁÏרҵµ÷¼Á +7 ÍÛºôºßºôºß 2026-04-02 7/350 2026-04-02 20:13 by yulian1987
[¿¼ÑÐ] 349Çóµ÷¼Á +10 zwjjjjjj 2026-03-31 10/500 2026-04-02 20:13 by dongzh2009
[¿¼ÑÐ] 298Çóµ÷¼Á +4 zzz£¬£¬r 2026-04-02 7/350 2026-04-02 18:57 by 1939136013¹·×³
[¿¼ÑÐ] һ־Ըͬ¼Ã´óѧ323·Ö£¨080500£©Çóµ÷¼Á +6 yikeniu 2026-04-01 6/300 2026-04-02 14:19 by smileboy2006
[¿¼ÑÐ] Ò»Ö¾Ô¸±±¾©¿Æ¼¼´óѧ085601²ÄÁϹ¤³ÌÓ¢Ò»Êý¶þ³õÊÔ×Ü·Ö335Çóµ÷¼Á +9 Ë«ÂíβƦÀϰå2 2026-04-01 9/450 2026-04-02 12:14 by oooqiao
[¿¼ÑÐ] 337Çóµ÷¼Á +11 ¡¶Ê÷¡· 2026-03-29 11/550 2026-04-02 10:20 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] 070300»¯Ñ§Çóµ÷¼Á +14 С»ÆÑ¼±¦ 2026-03-30 14/700 2026-04-02 10:07 by oooqiao
[¿¼ÑÐ] 348»·¾³¹¤³Ìµ÷¼Á +3 ÎâÑå׿24k 2026-04-01 3/150 2026-04-02 09:14 by nanaliuyun
[¿¼ÑÐ] Ò»Ö¾Ô¸±±¾©¿Æ¼¼£¬085601×Ü·Ö305Çóµ÷¼Á +9 °ëÉú¹Ï£¡ 2026-04-01 11/550 2026-04-02 08:28 by Wang200018
[¿¼ÑÐ] ²ÄÁÏ¿ÆÑ§Ó빤³Ì339Çóµ÷¼Á +11 hyz0119 2026-03-31 12/600 2026-04-01 18:40 by ΰ´óºÓ±±
[¿¼ÑÐ] ʡ˫һÁ÷ÖØµãÒ»±¾´óѧÕÐÊÕµ÷¼Á +4 wwwwffffff 2026-03-31 7/350 2026-04-01 15:23 by wwwwffffff
[¿¼ÑÐ] ²ÄÁϵ÷¼Á +10 Eujd1 2026-03-31 11/550 2026-04-01 11:23 by ivanqyq
[¿¼ÑÐ] Ò»Ö¾Ô¸ ÄϾ©º½¿Õº½Ìì´óѧ £¬080500²ÄÁÏ¿ÆÑ§Ó빤³Ìѧ˶ +10 @taotao 2026-03-31 11/550 2026-04-01 09:43 by xiayizhi
[¿¼ÑÐ] 274Çóµ÷¼Á +6 xiao°®Í¬Ñ§ 2026-03-30 6/300 2026-03-31 10:04 by cal0306
[¿¼ÑÐ] 08¹¤¿Æ£¬295£¬½ÓÊÜ¿çרҵµ÷¼Á +6 lmnlzy 2026-03-30 6/300 2026-03-31 10:04 by cal0306
[¿¼ÑÐ] 22408 359·Öµ÷¼Á +4 Qshers 2026-03-27 8/400 2026-03-31 08:53 by Qshers
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û