版塊導航: 正在加載中...

登錄注冊

應《網絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

>論壇更新日志 (6051)
>導師招生 (2423)
>考研 (2286)
>蟲友互識 (491)
>文獻求助 (425)
>考博 (332)
>碩博家園 (230)
>招聘信息布告欄 (133)
>論文投稿 (113)
>休閑灌水 (107)
>博后之家 (85)
>基金申請 (79)
>公派出國 (49)
>找工作 (41)
>教師之家 (40)
>綠色求助(高懸賞) (39)

北京石油化工學院2026年研究生招生接收調劑公告

返回列表

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

ryan111a

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 9.5
帖子: 26
在線: 33.7小時
蟲號: 2475831
注冊: 2013-05-22
性別: GG
專業(yè): 電磁場與電路

[求助] 關于應用fortran編寫高斯消去法程序來求解方程組的問題

關于應用fortran編寫高斯消去法程序來求解方程組的問題。應用高斯消去法來求解方程組的fortran編程，我自己可以編寫，但如果方程的系數矩陣不便，而右邊的解矩陣發(fā)生了同樣的變化，怎樣才能通過循環(huán)，來求這個方程組的所有解。
   我當時在高斯消去法上加了個循環(huán)，但是在變換右邊解矩陣時，出現了如下錯誤，怎么解決？
   E:\machao reactor\assitance\restart\one.f90(399) : Error: A constant or named constant is required in this context. [IM]
      ImCmplx(Row,j)=(Im(j),0)
   Im(j)我想表達的意思是每次循環(huán)對應的的右邊的解矩陣。j是高斯消去法循環(huán)的次數。

回復此樓

» 猜你喜歡

11408總分309，一志愿東南大學求調劑，不挑專業(yè) 已經有6人回復
求調劑，一志愿南京航空航天大學，080500材料科學與工程學碩，總分289分已經有8人回復
環(huán)境工程 085701，267求調劑已經有7人回復
材料工程專碩求調劑已經有8人回復
297求調劑已經有3人回復
材料科學與工程求調劑已經有7人回復
281求調劑已經有7人回復
288資源與環(huán)境專碩求調劑，不限專業(yè)，有學上就行已經有8人回復
求收留已經有6人回復
085600 295分求調劑已經有9人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

【求助】用最小二乘法求解線性方程組的Fortran代碼已經有4人回復
【求助】用fortran求解大型線性方程組時出現的錯誤【已解決】已經有11人回復

live as if you will die today

1樓 2013-06-05 15:22:16

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

ryan111a

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 9.5
帖子: 26
在線: 33.7小時
蟲號: 2475831
注冊: 2013-05-22
性別: GG
專業(yè): 電磁場與電路

引用回帖:

2樓: Originally posted by pippi6 at 2013-06-05 15:54:26
給你個程序

! How to use
!  real*8 :: c(n_dim,n_dim),e(n_dim,n_dim)
!  integer :: n=4
!    n          size of the matrix to be inversed
!    n_dim    array size
!    e          unit ma ...

十分感謝您的回答，相信用您的程序肯定可以實現此功能，因為我只是注重fortran程序實現我的功能，不知您是否有時間，幫忙修改下我的這塊程序。

!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
SUBROUTINE ICGaussA(ICmplx,NumT,NumI,MInd,R,Im) !復數的高斯消去法(用來求包封電流或層電流)
IMPLICIT NONE

INTEGER NumT,NumI !總包封數或層數
COMPLEX ICmplx(NumT,NumI) !復數電流值，待輸出量
REAL MInd(NumT,NumT) !包封或層的互感矩陣
REAL R(NumT) !包封或層的電阻向量
REAL Im(NumI)

COMPLEX MZ[ALLOCATABLE](:,

!電壓阻抗矩陣
COMPLEX TZ[ALLOCATABLE](:,

!電流阻抗矩陣
COMPLEX ImCmplx[ALLOCATABLE](:,

!右端復數電流向量

COMPLEX Temp,Sum !中間變量
INTEGER Row,Col,Num !行號和列號
INTEGER ME !主元
INTEGER i
INTEGER j

REAL,PARAMETER:

mga=314.159
REAL,PARAMETER::pai=3.14159

      ALLOCATE(TZ(NumT,NumT))
ALLOCATE(MZ(NumT,NumT))
ALLOCATE(ImCmplx(NumT,NumT))
DO j=1,NumI
      DO Row=1,NumT
      DO Col=1,NumT
         IF(Row==Col) THEN
            MZ(Row,Col)=CMPLX(R(Row),omga*MInd(Row,Col))
      ELSE
      MZ(Row,Col)=CMPLX(0,omga*MInd(Row,Col))
      END IF
      END DO
      END DO

      DO Row=1,NumT
      IF(Row==NumT) THEN
   DO Col=1,NumT
      TZ(Row,Col)=CMPLX(1,0)
   END DO
      ELSE
            DO Col=1,NumT
               TZ(Row,Col)=MZ(Row,Col)-MZ(Row+1,Col)
            END DO
   END IF
      END DO
      DO Row=1,NumT
         IF(Row==NumT) THEN
         ImCmplx(Row,j)=(Im(j),0)
      ELSE
         ImCmplx(Row,j)=(0,0)
      END IF
      END DO
      DO Row=1,NumT
      ME=Row
      DO i=Row+1,NumT
      IF(ABS(TZ(i,Row))>ABS(TZ(ME,Row))) THEN
   ME=i
      END IF
      END DO
      DO i=1,NumT
      Temp=TZ(ME,i)
      TZ(ME,i)=TZ(Row,i)
      TZ(Row,i)=Temp
      END DO
      Temp=ImCmplx(Row,i)
      ImCmplx(Row,i)=ImCmplx(ME,i)
      ImCmplx(ME,i)=Temp
      DO Col=Row+1,NumT
      TZ(Col,Row)=TZ(Col,Row)/TZ(Row,Row)
      DO i=Row+1,NumT
         TZ(Col,i)=TZ(Col,i)-TZ(Col,Row)*TZ(Row,i)
      END DO
      ImCmplx(Col,j)=ImCmplx(Col,j)-TZ(Col,Row)*ImCmplx(Row,j)
      END DO
      END DO
      ICmplx(NumT,j)=ImCmplx(NumT,j)/TZ(NumT,NumT)
      DO Row=NumT-1,1,-1
      Sum=ImCmplx(Row,j)
      DO Num=Row+1,NumT
      Sum=Sum-TZ(Row,i)*ImCmplx(i,Num)
      END DO
      ICmplx(Row,j)=Sum/TZ(Row,Row)
      END DO
END DO
DEALLOCATE(MZ)
DEALLOCATE(TZ)
DEALLOCATE(ImCmplx)
RETURN
END SUBROUTINE
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@

贊一下

回復此樓

live as if you will die today

3樓2013-06-06 08:59:26

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 6 個回答

pippi6

鐵桿木蟲 (著名寫手)

工程和科學數值計算咨詢

應助: 413 (碩士)
貴賓: 0.002
金幣: 7116.5
散金: 15
紅花: 63
帖子: 1639
在線: 798.9小時
蟲號: 2469437
注冊: 2013-05-14
專業(yè): 計算數學與科學工程計算

【答案】應助回帖

★ ★ ★ ★ ★
感謝參與，應助指數 +1
nono2009: 編輯內容 2013-06-05 22:01
ryan111a: 金幣+5, ★★★很有幫助 2013-06-06 08:59:59

給你個程序

CODE:

! How to use

!  real*8 :: c(n_dim,n_dim),e(n_dim,n_dim)

!  integer :: n=4

!     n           size of the matrix to be inversed

!     n_dim       array size

!     e           unit matrix

!     c           matrix to be inversed

!     call gaussr_inverse(c,n,n_dim,e,n,n_dim)

!flag

subroutine gaussr_inverse(a,n,np,b,m,mp)

  !     Commented, reorganized and understood by xxxx, 1999-11-03

  !     1.   very good designed scheme. 

  !          seems No room for performance improve

  !     2.   row normalization (implicite pivoting) has not much  effect

  !           for full pivoting procedure

  !     991112  ASY/FJT, xxxx

  !     Test again. even important for full pivoting procedure

  !     3.   row normalization (implicite pivoting) has significant 

  !           accuracy enhancement on partial pivoting procedure

  !              7.e-5 -> 2.e-11

  !     000623  ASY/FJT, xxx

  !     Real version

  implicit none

  integer :: n,np,m,mp,i,j,k,irow,icol,l

  real*8 a(np,np),b(np,mp),pivinv,cc

  real*8, allocatable, dimension(:) :: dum

  real*8 big

  integer, allocatable, dimension(:) :: ipiv,indxr,indxc

  real*8, allocatable, dimension(:) :: c1

  allocate(ipiv(1:n),indxr(1:n),indxc(1:n),dum(1:max(n,m)),c1(1:n))

  ipiv=0

  !     row normalization (implicite pivoting)

  do i=1,n

     c1(i)=1/sqrt(sum(abs(a(i,1:n))**2))

     a(i,1:n)=a(i,1:n)*c1(i)

     b(i,1:m)=b(i,1:m)*c1(i)

  end do

  do i=1,n

     !     Full pivoting, however, not necessarily starting from 1st col

     big=0.d0

     do j=1,n

        if(ipiv(j)==1)cycle 

        do k=1,n

           if (ipiv(k)==1) cycle

           if (abs(a(j,k)) <= big) cycle

           big=abs(a(j,k))

           irow=j

           icol=k

        end do

     end do

     !     ipiv(icol)=1 means the column icol has been selected once

     ipiv(icol)=ipiv(icol)+1

     if (irow.ne.icol) then

        !     exchange a(irow,:) with a(icol,:)

        !     so a(irow,irow) is the selected pivot

        dum(1:n)= a(irow,1:n)  

        a(irow,1:n)= a(icol,1:n)  

        a(icol,1:n)=dum(1:n)

        dum(1:m)= b(irow,1:m)  

        b(irow,1:m)= b(icol,1:m)  

        b(icol,1:m)=dum(1:m)

     endif

     !     bookkeeping row- and column-indices

     indxr(i)=irow

     indxc(i)=icol

     if (abs(a(icol,icol))==0.d0) pause 'singular matrix.'

     pivinv=1.d0/a(icol,icol)

     a(icol,icol)=1.d0

     a(icol,1:n)=a(icol,1:n)*pivinv

     b(icol,1:m)=b(icol,1:m)*pivinv

     do l=1,n

        if(l==icol)cycle  ! only for row l not equal to icol

        cc=a(l,icol)

        a(l,icol)=0.d0    ! clever design

        a(l,1:n)=a(l,1:n)-a(icol,1:n)*cc

        b(l,1:m)=b(l,1:m)-b(icol,1:m)*cc

     end do

  end do

  do l=1,n

     if(indxr(l) == indxc(l)) cycle

     dum(1:n)=a(1:n,indxr(l))

     a(1:n,indxr(l))=a(1:n,indxc(l))

     a(1:n,indxc(l))=dum(1:n)

  end do

  do i=1,n

     a(1:n,i)=a(1:n,i)*c1(i)

  end do

  deallocate(ipiv,indxr,indxc,dum,c1)

end subroutine gaussr_inverse

[ Last edited by nono2009 on 2013-6-5 at 22:01 ]

贊一下(1人)

回復此樓

2樓2013-06-05 15:54:26

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

pippi6

鐵桿木蟲 (著名寫手)

工程和科學數值計算咨詢

應助: 413 (碩士)
貴賓: 0.002
金幣: 7116.5
散金: 15
紅花: 63
帖子: 1639
在線: 798.9小時
蟲號: 2469437
注冊: 2013-05-14
專業(yè): 計算數學與科學工程計算

你這兩句是想干什么？
                     ImCmplx(Row,j)=(Im(j),0)
                     ImCmplx(Row,j)=(0,0)
如果是想表示一個復數由實部和虛部構成，應這么寫
                  ImCmplx(Row,j)=cmplx(Im(j),0)
                     ImCmplx(Row,j)=cmplx(0,0)
anyway， (0,0) 這樣的寫法在fortran里沒有意義

贊一下(1人)

回復此樓

4樓2013-06-06 17:00:51

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

pippi6

鐵桿木蟲 (著名寫手)

工程和科學數值計算咨詢

應助: 413 (碩士)
貴賓: 0.002
金幣: 7116.5
散金: 15
紅花: 63
帖子: 1639
在線: 798.9小時
蟲號: 2469437
注冊: 2013-05-14
專業(yè): 計算數學與科學工程計算

引用回帖:

4樓: Originally posted by pippi6 at 2013-06-06 17:00:51
你這兩句是想干什么？
                     ImCmplx(Row,j)=(Im(j),0)
                     ImCmplx(Row,j)=(0,0)
如果是想表示一個復數由實部和虛部構成，應這么寫
                  ImCmplx(Row ...

sorry，（0,0） ok，但 (Im(j),0) 不行

贊一下(1人)

回復此樓

5樓2013-06-06 17:02:08

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 6 個回答

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 083000環(huán)境科學與工程調劑，總分281 +3	橙子（勝意） 2026-03-30	3/150	2026-03-30 20:07 by piklet
[考研] 309求調劑 +15	誰不是少年 2026-03-29	15/750	2026-03-30 19:28 by 2026材料調劑
[考研] 105500藥學求調劑，一志愿山東大學藥學，348分 +3	gr哈哈哈 2026-03-28	3/150	2026-03-30 18:56 by 源_2020
[考研] 一志愿廈門大學材料工程專碩354找調劑�。�！ +5	貝唄鋇鋇 2026-03-30	5/250	2026-03-30 18:16 by 無際的草原
[考研] 0703 化學求調劑，一志愿山東大學 342 分 +7	Shern—- 2026-03-28	7/350	2026-03-30 16:31 by nothing投稿中
[考研] 284求調劑 +14	junqihahaha 2026-03-26	15/750	2026-03-30 14:12 by 探123
[考研] 考研調劑 +5	Sanmu-124 2026-03-26	5/250	2026-03-30 13:36 by chemdavid
[考研] 283求調劑（080500） +14	A child 2026-03-27	14/700	2026-03-30 12:06 by 探123
[考博] 26申博自薦 +6	whh869393 2026-03-24	6/300	2026-03-29 21:05 by nxgogo
[考研] 085600，專業(yè)課化工原理，321分求調劑 +5	大饞小子 2026-03-28	5/250	2026-03-29 08:56 by qingfeng258
[考研] 0703化學調劑，求導師收 +9	天天好運來上岸?/a> 2026-03-24	10/500	2026-03-28 22:17 by chemzp
[考研] 一志愿太原理工安全工程300分，求調劑 +5	0857求調劑. 2026-03-24	6/300	2026-03-28 22:04 by zhq0425
[考研] 071000生物學求調劑，初試成績343 +7	小小甜面團 2026-03-25	7/350	2026-03-28 20:25 by 唐沐兒
[考研] 復試調劑 +3	raojunqi0129 2026-03-28	3/150	2026-03-28 15:27 by 落�？伤�
[考研] 283求調劑 +7	A child 2026-03-28	7/350	2026-03-28 12:05 by zllcz
[考研] 315分求調劑 +7	26考研上岸版26 2026-03-26	7/350	2026-03-28 04:05 by fmesaito
[考研] 295求調劑 +5	1428151015 2026-03-27	6/300	2026-03-28 04:04 by fmesaito
[考研] 308求調劑 +7	墨墨漠 2026-03-25	7/350	2026-03-27 14:47 by 狂炫麥當當
[考研] 材料求調劑 +5	.m.. 2026-03-25	5/250	2026-03-27 11:08 by 不吃魚的貓
[考研] 085602化學工程求調劑。 +4	平樂樂樂 2026-03-26	4/200	2026-03-26 17:57 by fmesaito

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

ryan111a

[求助] 關于應用fortran編寫高斯消去法程序來求解方程組的問題

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

ryan111a

pippi6

【答案】應助回帖

pippi6

pippi6

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

» 本主題相關價值貼推薦，對您同樣有幫助: