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

調(diào)劑小程序

登錄注冊

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

24小時熱門版塊排行榜

>論壇更新日志 (5528)
>導師招生 (2062)
>考研 (1856)
>蟲友互識 (416)
>文獻求助 (396)
>考博 (293)
>碩博家園 (186)
>招聘信息布告欄 (121)
>論文投稿 (99)
>休閑灌水 (82)
>博后之家 (80)
>基金申請 (72)
>公派出國 (43)
>綠色求助(高懸賞) (39)
>找工作 (38)
>教師之家 (35)

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

返回列表

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

ryan111a

新蟲 (初入文壇)

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

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

關(guān)于應用fortran編寫高斯消去法程序來求解方程組的問題。應用高斯消去法來求解方程組的fortran編程，我自己可以編寫，但如果方程的系數(shù)矩陣不便，而右邊的解矩陣發(fā)生了同樣的變化，怎樣才能通過循環(huán)，來求這個方程組的所有解。
   我當時在高斯消去法上加了個循環(huán)，但是在變換右邊解矩陣時，出現(xià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)的次數(shù)。

回復此樓

» 猜你喜歡

0703化學321分求調(diào)劑已經(jīng)有10人回復
0703化學求調(diào)劑已經(jīng)有7人回復
一志愿廈門大學材料工程專碩354找調(diào)劑！��！已經(jīng)有5人回復
266分，求材料相關(guān)專業(yè)調(diào)劑已經(jīng)有5人回復
292求調(diào)劑已經(jīng)有13人回復
材料化工340求調(diào)劑已經(jīng)有3人回復
295材料工程專碩求調(diào)劑已經(jīng)有9人回復
085600 材料與化工 329分求調(diào)劑已經(jīng)有18人回復
一志愿北京化工大學材料與化工（085600）296求調(diào)劑已經(jīng)有23人回復
08工科，295，接受跨專業(yè)調(diào)劑已經(jīng)有3人回復

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

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

live as if you will die today

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

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

pippi6

鐵桿木蟲 (著名寫手)

工程和科學數(shù)值計算咨詢

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

你這兩句是想干什么？
                     ImCmplx(Row,j)=(Im(j),0)
                     ImCmplx(Row,j)=(0,0)
如果是想表示一個復數(shù)由實部和虛部構(gòu)成，應這么寫
                  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

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

查看全部 6 個回答

pippi6

鐵桿木蟲 (著名寫手)

工程和科學數(shù)值計算咨詢

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

【答案】應助回帖

★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
nono2009: 編輯內(nèi)容 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

已閱回復此樓關(guān)注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 ...

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

!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
SUBROUTINE ICGaussA(ICmplx,NumT,NumI,MInd,R,Im) !復數(shù)的高斯消去法(用來求包封電流或?qū)与娏?
IMPLICIT NONE

INTEGER NumT,NumI !總包封數(shù)或?qū)訑?shù)
COMPLEX ICmplx(NumT,NumI) !復數(shù)電流值，待輸出量
REAL MInd(NumT,NumT) !包封或?qū)拥幕ジ芯仃?br /> REAL R(NumT) !包封或?qū)拥碾娮柘蛄?br /> REAL Im(NumI)

COMPLEX MZ[ALLOCATABLE](:,

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

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

!右端復數(shù)電流向量

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

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

pippi6

鐵桿木蟲 (著名寫手)

工程和科學數(shù)值計算咨詢

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

引用回帖:

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

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

贊一下(1人)

回復此樓

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

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

查看全部 6 個回答

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

[考研] 材料專碩 085600求調(diào)劑 +7	BBQ233 2026-03-30	7/350	2026-03-30 17:44 by oooqiao
[考研] 322求調(diào)劑：一志愿湖南大學材料與化工（085600），已過六級。 +9	XX小鄧 2026-03-29	9/450	2026-03-30 17:18 by limeifeng
[考研] 342求調(diào)劑 +4	加油a李zs 2026-03-26	4/200	2026-03-30 16:39 by 晶體之美
[考研] 332求調(diào)劑 +13	032500 2026-03-25	13/650	2026-03-30 16:25 by wangjy2002
[考研] 322求調(diào)劑 +9	宋明欣 2026-03-27	9/450	2026-03-30 16:15 by wangjy2002
[考研] 329求調(diào)劑，一志愿西北工業(yè)大學，材料工程（085601） +5	小小機靈蟲 2026-03-29	11/550	2026-03-30 15:02 by Wang200018
[考研] 考研調(diào)劑 +5	Sanmu-124 2026-03-26	5/250	2026-03-30 13:36 by chemdavid
[考研] 一志愿南開大學0710生物學359求調(diào)劑 +5	兔兔兔111223314 2026-03-29	5/250	2026-03-30 10:41 by malone0211
[考研] 295求調(diào)劑 +5	wei-5 2026-03-26	5/250	2026-03-30 08:34 by 探123
[考研] 327求調(diào)劑 +4	小卡不卡. 2026-03-29	4/200	2026-03-30 06:15 by wxiongid
[考研] 085600 286分材料求調(diào)劑 +10	麻辣魷魚 2026-03-27	11/550	2026-03-29 21:52 by nanaliuyun
[考博] 26申博自薦 +6	whh869393 2026-03-24	6/300	2026-03-29 21:05 by nxgogo
[考研] 327求調(diào)劑 +6	汲亦昊 2026-03-29	6/300	2026-03-29 13:40 by peike
[考研] 308求調(diào)劑 +7	墨墨漠 2026-03-27	7/350	2026-03-28 07:43 by 熱情沙漠
[考研] 315分求調(diào)劑 +7	26考研上岸版26 2026-03-26	7/350	2026-03-28 04:05 by fmesaito
[考研] 298調(diào)劑 +3	jiyingjie123 2026-03-27	3/150	2026-03-27 11:57 by wxiongid
[考研] 0703化學338求調(diào)劑！ +6	Zuhui0306 2026-03-26	7/350	2026-03-27 10:35 by shangxh
[考研] 考研一志愿蘇州大學初始315（英一）求調(diào)劑 +3	sbdksD 2026-03-24	4/200	2026-03-25 18:16 by xcjcqu
[考研] 296求調(diào)劑 +4	汪�。�！ 2026-03-25	7/350	2026-03-25 16:41 by 汪！？！
[考研] 277分求調(diào)劑，跨調(diào)材料 +3	考研調(diào)劑lxh 2026-03-24	3/150	2026-03-24 13:52 by JourneyLucky

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

24小時熱門版塊排行榜

ryan111a

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

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

pippi6

pippi6

【答案】應助回帖

ryan111a

pippi6

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

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