fortran常用演算法程序集
1. 如何用FORTRAN 2011內部數學函數求解線性方程組
MKL的建議看下幫助文檔,國內用Fortran的用到數學庫的極少,寫書的都盯著錢,就算出了也沒幾個人買的
2. fortran中階乘演算法
http://www.fcode.cn/algorithm-50-1.html
我也不太理解,我是問的大神,論壇那有專門的講解,你看看吧
3. 二元三點插值法
不知道你還需要知道什麼。
4. 求fortran實現FFT(快速傅里葉變換)程序。
徐士良的《fortran常用演算法程序集》中有fft代碼,直接可以用的那種
5. 求一個fortran主程序,用來計算矩陣方程AX+XB=C的。方程的演算法已經有了,是個子例行程序。
有一個叫做徐士良 Fortran常用演算法程序集的東西,裡面有很多編好的子程序,用不同的方式解矩陣方程的子程序也有,很好用,你可以參考。
6. 怎麼用fortran編寫解非線性方程
也可參看《Fortran常用演算法程序集》
7. 求FORTRAN程序設計的電子書,最好是權威出版社的,清晰版本
彭國倫《Fortran95程序設計》用於入門,學習基本的語法
徐士良_Fortran常用演算法程序集-第二版,用於具體的數值計算,可查詢常用的演算法
這兩本書的電子書和代碼都給你發網路消息了,請查閱。
8. 用直接消去法解方程組的程序如何編寫(Fortran程序)
!高斯消去法
subroutine agaus(a,b,n,x,l,js)
dimension a(n,n),x(n),b(n),js(n)
double precision a,b,x,t
l=1 !邏輯變數
do k=1,n-1
d=0.0
do i=k,n
do j=k,n
if (abs(a(i,j))>d) then
d=abs(a(i,j))
js(k)=j
is=i
end if
end do
end do !把行絕對值最大的元素換到主元位置
if (d+1.0==1.0) then
l=0
else !最大元素為0無解
if(js(k)/=k) then
do i=1,n
t=a(i,k)
a(i,k)=a(i,js(k))
a(i,js(k))=t
end do !最大元素不在K行,K行
end if
if(is/=k) then
do j=k,n
t=a(k,j)
a(k,j)=a(is,j)
a(is,j)=t !交換到K列
end do
t=b(k)
b(k)=b(is)
b(is)=t
end if !最大元素在主對角線上
end if !消去
if (l==0) then
write(*,100)
return
end if
do j=k+1,n
a(k,j)=a(k,j)/a(k,k)
end do
b(k)=b(k)/a(k,k) !求三角矩陣
do i=k+1,n
do j=k+1,n
a(i,j)=a(i,j)-a(i,k)*a(k,j)
end do
b(i)=b(i)-a(i,k)*b(k)
end do
end do
if (abs(a(n,n))+1.0==1.0) then
l=0
write(*,100)
return
end if
x(n)=b(n)/a(n,n)
do i=n-1,1,-1
t=0.0
do j=i+1,n
t=t+a(i,j)*x(j)
end do
x(i)=b(i)-t
end do
100 format(1x,'fail')
js(n)=n
do k=n,1,-1
if (js(k)/=k) then
t=x(k)
x(k)=x(js(k))
x(js(k))=t
end if
end do
return
end
program main
dimension a(4,4),b(4),x(4),js(4)
double precision a,b,x
real m1,m2,j
open(1,file="laiyi.txt")
read(1,*)m1,m2,j
close(1)
n=4
print*,m1,m2,j
a(1,1)=m1*cos(3.14159*j/180)
a(1,2)=-m1
a(1,3)=-sin(3.14159*j/180)
a(1,4)=0
a(2,1)=m1*sin(3.14159*j/180)
a(2,2)=0
a(2,3)=cos(3.14159*j/180)
a(2,4)=0
a(3,1)=0
a(3,2)=m2
a(3,3)=-sin(3.14159*j/180)
a(3,4)=0
a(4,1)=0
a(4,2)=0
a(4,3)=-cos(3.14159*j/180)
a(4,4)=1
b(1)=0
b(2)=m1*9.8
b(3)=0
b(4)=m2*9.8
call agaus(a,b,n,x,l,js)
if (l/=0) then
write(*,*)"a1=",x(1),"a2=",x(2) ,"n1=",x(3),"n2=",x(4)
end if
end
!逆矩陣求解
SUBROUTINE qiuni(A,N,L,IS,JS)
DIMENSION A(N,N),IS(N),JS(N)
DOUBLE PRECISION A,T,D
L=1
DO K=1,N
D=0.0
DO I=K,N
DO J=K,N
IF(ABS(A(I,J)).GT.D) THEN !把最大的元素給D
D=ABS(A(I,J))
IS(K)=I
JS(K)=J
END IF
END DO
END DO
IF (D+1.0.EQ.1.0)THEN
L=0
WRITE(*,200)
RETURN
END IF
200 FORMAT(1X,'ERR**NOT INV')
DO J=1,N
T=A(K,J)
A(K,J)=A(IS(K),J)
A(IS(K),J)=T
END DO
DO I=1,N
T=A(I,K)
A(I,K)=A(I,JS(K))
A(I,JS(K))=T
END DO
A(K,K)=1/A(K,K)
DO J=1,N
IF(J.NE.K)THEN
A(K,J)=A(K,J)*A(K,K)
END IF
END DO
DO I=1,N
IF(I.NE.K)THEN
DO J=1,N
IF(J.NE.K)THEN
A(I,J)=A(I,J)-A(I,K)*A(K,J)
END IF
END DO
END IF
END DO
DO I=1,N
IF(I.NE.K)THEN
A(I,K)=-A(I,K)*A(K,K)
END IF
END DO
END DO
DO K=N,1,-1
DO J=1,N
T=A(K,J)
A(K,J)=A(JS(K),J)
A(JS(K),J)=T
END DO
DO I=1,N
T=A(I,K)
A(I,K)=A(I,IS(K))
A(I,IS(K))=T
END DO
END DO
RETURN
END
SUBROUTINE BRMUL(A,B,N,C)
DIMENSION A(N,N),B(N),C(N)
DOUBLE PRECISION A,B,C
DO I=1,N
DO J=1,N
C(I)=0.0
DO L=1,N
C(I)=C(I)+A(I,L)*B(L)
END DO
END DO
END DO
RETURN
END
program main
DIMENSION A(4,4),B(4,1),C(4,1),IS(4),JS(4)
DOUBLE PRECISION A,B,C
REAL M1,M2,JD
OPEN(1,FILE='LAIYI.TXT')
READ(1,*) M1,M2,JD
PRINT*,M1,M2,JD
CLOSE(1)
A(1,1)=M1*COS(3.14*JD/180)
A(1,2)=-M1
A(1,3)=-SIN(3.14*JD/180)
A(1,4)=0
A(2,1)=M1*SIN(3.14*JD/180)
A(2,2)=0
A(2,3)=COS(3.14*JD/180)
A(2,4)=0
A(3,1)=0
A(3,2)=M2
A(3,3)=-SIN(3.14*JD/180)
A(3,4)=0
A(4,1)=0
A(4,2)=0
A(4,3)=-COS(3.14*JD/180)
A(4,4)=1
B(1,1)=0
B(2,1)=M1*9.8
B(3,1)=0
B(4,1)=M2*9.8
CALL QIUNI(A,4,L,IS,JS)
CALL BRMUL(A,B,4,C)
WRITE(*,*) (C(I,1),I=1,4)
END
畫圖
USE MSFLIB
INTEGER status
TYPE(xycoord) xy
status=SETCOLORRGB(#FFFFFF)
status1=SETCOLORRGB(#0000FF)
OPEN(1,FILE="G.TXT")
READ(1,*) G1,G2,G3,G4
OPEN(2,FILE="N.TXT")
READ(2,*) N1,N2,N3,N4
CALL MOVETO(INT(20),INT(20),xy)
status=LINETO(INT(40),INT(G1))
status=LINETO(INT(80),INT(G2))
status=LINETO(INT(120),INT(G3))
status=LINETO(INT(160),INT(G4))
CALL SETLINESTYLE(#FF00)
CALL MOVETO(INT(20),INT(20),xy)
status1=LINETO(INT(40),INT(N1))
status1=LINETO(INT(80),INT(N2))
status1=LINETO(INT(120),INT(N3))
status1=LINETO(INT(160),INT(N4))
READ(*,*)
END
如果對您有幫助,請記得採納為滿意答案,謝謝!祝您生活愉快!
9. 關於清華大學電子工程系
培養目標
電子工程各專業的本科學生應掌握扎實的基礎理論、專業基礎理論和專業知識及基本技能,並掌握一定的人文社會科學、經濟管理、環境工程等方面的基本理論,具有在專業領域跟蹤新理論、新知識、新技術的能力,能從事該領域的科學研究、技術開發、教學及管理等工作。
學位授予
電子信息工程
電子科學與技術
上述兩個專業均可授予學士學位。
電子信息工程
電子信息工程是信息產業的重要基礎和支柱之一。它以電路與系統、信號與信息處理、電磁場與微波等理論為基礎,研究各種信息如語音、文字、圖象、遙感信息等的處理、交換、及無線、電纜、光纜等的傳輸,在此基礎上研究和發展各種電子與信息系統。本專業主要內容為:
各種信息的處理、傳輸、交換、檢測的理論與技術、各種電信系統和網路;
電路理論、電子系統設計及應用、系統模擬和集成電路設計自動化的理論和技術;
微波理論與技術、天線與電磁兼容理論與技術;
電子電路應用技術、計算機應用技術、信息處理應用技術、電磁波應用技術。
在本專業設有微波與數字通信國家重點實驗室,智能技術與系統國家重點實驗室智能圖形圖象分室,電子設計自動化(EDA)實驗室等,擁有多種先進儀器設備。
本專業畢業生適應的工作范圍為:信號的處理、傳輸、交換及檢測技術的研究與教學工作,電子設備與系統的研製、生產與應用,電子技術及計算機技術應用與開發,微波技術的研究、應用及開發工作。
電子科學與技術
電子科學與技術是信息科學技術的前沿學科,它以近代物理學與數學為基礎,研究電磁波的產生、運動及在不同介質中的相互作用規律,在此基礎上發明和發展各種信息電子材料、元器件、集成電路乃至集成電子系統。電子科學與技術專業有兩個方向:物理電子學與光電子學、微電子學。
物理電子學與光電子學方向主要內容為:
信息光電子學與光子器件;
光纖通信系統與網路應用技術;
新型顯示器件和新型電光薄膜材料與器件;
微細技術和信息光電子材料評價與檢測技術;
電子電路應用技術、計算機應用技術。
微電子學方向主要內容為:
大規模、超大規模集成電路和集成系統的研究、設計、應用;
各種半導體器件、微電子機械繫統的研究、應用;
亞微米和深亞微米集成電路工藝;
集成電路的計算機輔助設計、製造和測試方法學。
本專業設有集成光電子學國家重點實驗室(清華大學實驗區),超凈工藝線和電子系統集成與專用集成電路技術研究中心。
本專業畢業生適應的工作范圍為:光子學、光電子學、物理電子學與微電子學的研究、教學工作和相關技術的開發、應用。光通信、光感測、信息顯示系統、集成電路計算機輔助設計、集成電子系統等方面的研究、設計、製造與開發工作。
要考研的話,看這些沒有,有個群85501851,你可以去多了解
教材名稱 主編姓名 出版社 版次 出版年月 適用專業 適用層次 統一書號(ISBN)
語音信號數字處理 楊行峻等 電子 1 1995.8 電子信息工程 本、碩 7-5053-3147-7
現代分析技術 陸家和 陳長彥 清華 1 1995.9 物理電子 光電子 真空技術 本、碩 7-302-01830-8
計算機常用演算法 徐士良 清華 2 1995.11 理工 本、碩 7-302-01958-4
FORTRAN常用演算法程序集 徐士良 清華 2 1995.12 理工 本、碩 7-302-01947-9
UNIX應用教程 楊華中 郵電 1 1996.9 電子科學與技術 本、碩 7-115-06214-5
C常用演算法程序集 徐士良 清華 2 1996.11 所有專業 本、碩 7-302-02290-9
微波工程基礎 李宗謙 東南大學 1 1996.12 電子信息工程 本科 7-81050-144-5
電子線路基礎 高文煥 高教 1 1997 所有專業 本科 7-04-005996-7
7-118-01730-2
程序設計實踐指導 朱明方 趙純善 雷田玉 清華 1 1997.3 電子信息工程
電子科學與技術 本科 7-302-02344-1
軟體技術基礎 徐士良 高教 1 1997.4 理工 本科 7-04-005950-9
C語言設計實驗指導 徐士良 清華 1 1997.9 所有專業 本科 7-302-02621-1
BASIC常用演算法程序集 徐士良 清華 1 1997.9 所有專業 本、碩 7-302-02662-9
模擬電路的計算機分析與設計 高文煥 汪 蕙 清華 1 1999 所有專業 本科 7-302-03158-4
模擬集成電路的自動綜合方法 楊華中 汪 蕙 科學 1 1999.2 電子科學與技術 碩士 7-03-006975-7
圖象工程 上冊 圖象處理和分析 章毓晉 清華 1 1999.3 電子信息工程 本科 7-302-03343-9
超大規模集成電路方法學導論 楊之廉 清華 1 1999.3 電子科學與技術 本科 7-302-03275-0
程式控制交換與綜合業務通信網 樂正友 楊為理 清華 1 1999.5 電子信息工程 本科 7-302-03355-2
軟體工程 徐士良 清華 1999.5 本科 7-302-03453-2
信號與系統 (第二版)上冊 鄭君里 應啟珩 楊為理 高教 2 2000 工科、理科 7-04-007981-X
信號與系統 (第二版)下冊 鄭君里 應啟珩 楊為理 高教 2 2000 工科、理科 7-04-007983-6
半導體激光器 江劍平 電子 1 2000.2 本科 7-5053-5535-X
微機圖象處理系統 蘇光大 清華 1 2000.7 本科 7-302-03963-1
數字集成系統的結構化設計與高層次綜合 王志華 清華 1 2000.7 本科 7-302-03837-6
圖象工程 下冊 圖象處理與計算機視覺 章毓晉 清華 1 2000.8 電子信息工程 碩士 7-302-03980-1
FORTRAN程序設計上機指導 徐士良 清華 2 2001.1 7-302-03674-8
電磁場理論基礎 王薔 李國定 清華 1 2001.1 7-302-04251-9
半導體集成電路 朱正涌 清華 1 2001.1 7-302-040851-03333
應用資訊理論基礎 朱雪龍 清華 1 2001.3 7-302-04154-7
離散時間信號分析和處理 應啟珩 馮一雲 竇維蓓 清華 1 2001.9 本科 7-302-04679-4
可編程ASIC設計及應用 李廣軍 孟憲元 電子科大(成都) 2002
圖像並行處理技術 蘇光大 清華 1 2002.1 碩士 7-302-05282-4
通信電路原理 董在望主編 陳雅琴 雷有華 肖華庭 高教 2 2002.8 本科 7-04-010643-4
10. 想好好學FORTRAN
彭國倫<<Fortran 95程序設計>>是非常流行的一本教材,很多人都在看,市面上沒有更好的入門書了。
網上也有電子版,你可以網路一下。
如果找常用演算法,可以查徐士良《fortran常用演算法程序集》第二版。