徐士良fortran常用算法程序集

㈣ 二元三点插值法

不知道你还需要知道什么。

㈤ 求一个fortran主程序，用来计算矩阵方程AX+XB=C的。方程的算法已经有了，是个子例行程序。

有一个叫做徐士良 Fortran常用算法程序集的东西，里面有很多编好的子程序，用不同的方式解矩阵方程的子程序也有，很好用，你可以参考。

㈥ fortran语言矩阵求逆

! aa为原矩阵，b为存放aa的逆矩阵，n为矩阵aa的维数
subroutine nizhen(aa,b,n)
integer n,i,j,k
real:: aa(n,n),b(n,n),a(n,n)
a=aa
do i=1,n
b(i,i)=1
enddo
do i=1,n
b(i,:)=b(i,:)/a(i,i)
a(i,i:n)=a(i,i:n)/a(i,i)
do j=i+1,n
do k=1,n
b(j,k)=b(j,k)-b(i,k)*a(j,i)
enddo
a(j,i:n)=a(j,i:n)-a(i,i:n)*a(j,i)
enddo
enddo
do i=n,1,-1
do j=i-1,1,-1
do k=1,n
b(j,k)=b(j,k)-b(i,k)*a(j,i)
enddo
enddo
enddo
end

㈦ 关于清华大学电子工程系

培养目标
电子工程各专业的本科学生应掌握扎实的基础理论、专业基础理论和专业知识及基本技能，并掌握一定的人文社会科学、经济管理、环境工程等方面的基本理论，具有在专业领域跟踪新理论、新知识、新技术的能力，能从事该领域的科学研究、技术开发、教学及管理等工作。

学位授予
电子信息工程
电子科学与技术

上述两个专业均可授予学士学位。

电子信息工程
电子信息工程是信息产业的重要基础和支柱之一。它以电路与系统、信号与信息处理、电磁场与微波等理论为基础，研究各种信息如语音、文字、图象、遥感信息等的处理、交换、及无线、电缆、光缆等的传输，在此基础上研究和发展各种电子与信息系统。本专业主要内容为:

各种信息的处理、传输、交换、检测的理论与技术、各种电信系统和网络;
电路理论、电子系统设计及应用、系统仿真和集成电路设计自动化的理论和技术;

微波理论与技术、天线与电磁兼容理论与技术；

电子电路应用技术、计算机应用技术、信息处理应用技术、电磁波应用技术。
在本专业设有微波与数字通信国家重点实验室,智能技术与系统国家重点实验室智能图形图象分室,电子设计自动化(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

㈧ 求fortran实现FFT（快速傅里叶变换）程序。

徐士良的《fortran常用算法程序集》中有fft代码，直接可以用的那种

㈨ 徐士良fortran的 快速傅里叶变换 到底怎么用啊 求助 给分

我想，这可能不是一个Fortran的问题，也不是徐世良的算法问题。

而是你对 FFT 问题的数学物理意义的理解问题。

一个实序列的信号，N个时间域采样，经过FFT以后，第一个数值是直流分量的 N 倍，第 i 个点表示 (i-1)*Fs/N 频率振幅的 N/2 倍。

第一个频率就是 1/dx，其中dx=Fs/N，其他的频率按照 (i-1)/dx 计算

㈩ 用直接消去法解方程组的程序如何编写(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

如果对您有帮助，请记得采纳为满意答案，谢谢！祝您生活愉快！