麻豆小视频在线观看_中文黄色一级片_久久久成人精品_成片免费观看视频大全_午夜精品久久久久久久99热浪潮_成人一区二区三区四区

首頁 > 學院 > 開發(fā)設計 > 正文

學數(shù)學也要學點偷懶的技術(二）

2019-11-17 05:31:10
字體：大 中 小
來源：轉載
供稿：網(wǎng)友
    老師剛教逆矩陣，我就覺得求逆矩陣很難算，假如求四階以上的逆矩陣，至少也要四五分鐘才能解決。于是我想到用C來解決了。順便為自己上次寫的文章升升級也好。說寫就寫，可我很快發(fā)現(xiàn)找求逆矩陣的算法并不輕易。想了很久，沒辦法，就去圖書館查資料。找是找到了“高斯全選主元法”在一本久版的《c常用算法程序集》上，可是看了半個小時，還是看不明白，想必大師之作不讓我這小人物看懂。最后我下決心要自己找一個好的算法。   工夫不負有心人，我終于找到求逆矩陣的算法，而且有兩種。這兩種的算法都要調用到第一次寫的《學數(shù)學也要學點偷懶的技術》里的函數(shù)。下面我只用了其中的一種算法。   數(shù)學依據(jù)是：逆矩陣=原矩陣的N階行列式值的倒數(shù)*原矩陣的伴隨陣。   條件：原矩陣的N階行列式的值不能為零     請大家多多指教。   QQ:258220980
#include <stdio.h>
#include <math.h>
#include <conio.h>
#include <stdlib.h>#define N  10
#define M 10
#define S 10void njie_valu();
void njief_value();
float Dvalue(float a[N][N],int n);
void juzhench();
void nijuzhen();void nijuzhen()
{
    int i,j,n,m1,m2,n1,n2 ;
    int k ;
    float x[(N-1)*(N-1)];
    float y,s ;
    float a1[N][N],b[N][N],a[N][N],A[N][N],p[N][N],c[N][N];
    clrscr();
        scanf("%d",&n);
    printf("請以正確的順序輸入矩陣:/n");
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
        {
            scanf("%f",&a[i][j]);
        }
    }
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
        {
            printf("%f  ",a[i][j]);
        }
        printf("/n");
    }
    getch();
 if(Dvalue(a,n)==0)
 {
  printf("D=0,沒有逆矩陣！");
  exit(1);
 }
 s=1.0/Dvalue(a,n);
    for(m1=0;m1<n;m1++)
    {
        for(n1=0;n1<n;n1++)
        a1[m1][n1]=a[m1][n1];
    }
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
        {
            /*位于i行j列的元素變?yōu)?*/
            for(m1=0;m1<n;m1++)
            {
                for(n1=0;n1<n;n1++)
                if((m1==i)(n1==j))
                a1[m1][n1]=0 ;
            }
            k=0 ; 
            for(m1=0;m1<n;m1++)
            {
                for(n1=0;n1<n;n1++)
                {
                    if(a1[m1][n1]!=0)
                    {
                        x[k]=a1[m1][n1];
                        k++;
                    }
                }
            }
            k=0 ;
            for(m2=0;m2<n-1;m2++)
            {
                for(n2=0;n2<n-1;n2++)
                {
                    b[m2][n2]=x[k];
                    k++;
                }
            }
            /*對數(shù)組a1從新賦值*/
            for(m1=0;m1<n;m1++)
            {
                for(n1=0;n1<n;n1++)
                a1[m1][n1]=a[m1][n1];
            }
            y=pow(-1,i+j);
            A[i][j]=y*Dvalue(b,n-1);
        }
    }
    printf("A:/n");
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
        printf("%f    ",A[j][i]);
        printf("/n");
    }
    getch();
    printf("1/A=%f/n  ",s);
 getch();
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
        c[i][j]=0 ;
    }
 for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
        c[i][j]=c[i][j]+A[i][j]*s ;
    }
    printf("要求的逆矩陣是：/n");
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
     &n


發(fā)表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發(fā)表
主站蜘蛛池模板： 国产手机国产手机在线 | 国产高潮好爽好大受不了了 | 原来神马影院手机版免费 | 欧美视屏一区二区 | 成人午夜视频在线观看 | 国产精品久久久久久久久久妇女 | 91九色网| 毛片免费观看完整版 | 逼特逼视频在线观看 | 全黄性性激高免费视频 | 99精品视频免费看 | 人禽l交免费视频观看 视频 | 免费一级毛片在线播放视频老 | 在线播放中文 | 国产羞羞视频在线观看免费应用 | 免费男女视频 | 成人做爰s片免费看网站 | 亚洲成人国产综合 | 免费黄色大片网站 | 中文字幕在线亚洲 | 精品国产精品久久 | 国产亚洲精品久久午夜玫瑰园 | 九九午夜| 久久我不卡 | 91精品国 | 国产精品视频久久久 | 成人精品 | 亚洲少妇诱惑 | 久夜草 | 激情视频导航 | 国产精品成人一区二区三区电影毛片 | 国产九色在线播放九色 | 一区国产精品 | 美女视频在线观看黄 | 九九热在线精品视频 | 黄色免费小视频网站 | av久草 | 日韩精品一区二区三区中文 | 成人在线视频免费播放 | 久久96国产精品久久秘臀 | 久久精品9|