求解嵌套非线性方程组，改程序，求助

ttt

我要求A1，A3，可是只有含有A1,A3的表达式A，B,C,D的关系式，编程如下，求大师指导

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	program tu2

02

use IMSL

03	implicit none

04	external FCN

05	real*8 FCN

06	real,parameter::ERRREL=0.0001

07	integer,parameter::N=2

08	integer,parameter::ITMAX=100

09	real*8::XGUESS(N)=(/0.0467,0.0732/)

10	real X(N),FNORM

11	real*8 A1,A2,A3,A,B,C,D

12	real*8 Z0,Zi

13	real*8 E,F

14	real,parameter::pi=3.14159

15	integer,parameter::Zrs=5

16	integer,parameter::A4s=9

17	integer k,i

18	real*8::Zr(Zrs)=(/1.0,2.0,3.0,4.0,5.0/)

19	real*8::A4(A4s)=(/0.075,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.425/)

20	common /group2/ zi

21	common /group3/ A2

22	!!**************************************************

23	Zi=-3.0                  !!天线阻抗虚部的值

24	Z0=50.0                  !!传输线的特性阻抗

25	A2=0.25                  !!第二个支节的长度

26	!!***************************************************

27	do i=1,Zrs

28	do k=1,A4s

29	call  NEQNF(Zr(i),A4(k),FCN(i,k),ERRREL,N,ITMAX,XGUESS,X,FNORM)

30	write(*,*) X

31

stop

32

end

33	!!!**************************

34	subroutine A4_get_A(A4,A)

35	implicit none

36	real*8 A2,A3,A4,A

37	real,parameter::pi=3.14159

38	common /group3/ A2

39	A=cos(2.0*pi*A2)*cos(2.0*pi*A4)-sin(2.0*pi*A2)*(-cos(2.0*pi*A4)/tan(2.0*pi*A3)+sin(2.0*pi*A4))              !!A的函数式

40

41

return

42

end

43	!!!***************************

44	subroutine A4_get_B(A4,B)

45	implicit none

46	real*8 A2,A3,A4,z0,B

47	real,parameter::pi=3.14159

48	common /group3/ A2

49	common /group1/ Z0

50	B=Z0*(cos(2.0*pi*A2)*sin(2.0*pi*A4)+sin(2.0*pi*A2)*(sin(2.0*pi*A4)/tan(2.0*pi*A3)+cos(2.0*pi*A4)))            !!B的函数式

51

return

52

end

53	!!*****************************

54	subroutine A4_get_C(A4,C)

55	implicit none

56	real*8 A1,A2,A3,A4,Z0,C

57	real,parameter::pi=3.14159

58	common /group3/ A2

59	common /group1/ Z0

60	C=1.0/Z0*(cos(2.0*pi*A4)*((-cos(2.0*pi*A2))/tan(2.0*pi*A1)+sin(2.0*pi*A2))+(sin(2.0*pi*A2)/tan(2.0*pi*A1)+cos(2.0*pi*A2))*((-cos(2.0*pi*A4))/tan(2.0*pi*A3)+sin(2.0*pi*A4)))

61

return

62

end

63	!!*******************************

64	subroutine A4_get_D(A4,D)

65	implicit none

66	real*8 A1,A2,A3,A4,D

67	real,parameter::pi=3.14159

68

69	common /group3/ A2

70	D=-sin(2.0*pi*A4)*(-cos(2.0*pi*A2)/tan(2.0*pi*A1)+sin(2.0*pi*A2))+(sin(2.0*pi*A2)/tan(2.0*pi*A1)+cos(2.0*pi*A2))*(sin(2.0*pi*A4)/tan(2.0*pi*A3)+cos(2.0*pi*A4))

71

return

72

end

73	!!****************************************************

74

75	subroutine FCN(A,B,C,D,A4,Zr,XA,F,N)

76	implicit none

77

integer N

78	real,target::XA(N)

79

real F(N)

80	real,pointer::A1,A3

81	real*8 A,B,C,D,A4,Z0,Zr,Zi

82	common /group1/Z0

83	common /group2/Zi

84

85

A1=>XA(1)

86

A3=>XA(2)

87	call A4_get_A(A4,A)

88	call A4_get_B(A4,B)

89	call A4_get_C(A4,C)

90	call A4_get_D(A4,D)

91

92	F(1)=A*D+B*C-Z0*((D-C*Zi)*(D-C*Zi)+(C*Zr)*(C*Zr))

93	F(2)=(D*A-B*C)*Zi-A*C*(Zr*Zr+Zi*Zi)+D*B

94

return

95

end

以下是错误指示
Compiling Fortran...
C:\Program Files\Microsoft Visual Studio\MyProjects\tu2\tu2.f90
C:\Program Files\Microsoft Visual Studio\MyProjects\tu2\tu2.f90(29) : Error: The number of actual arguments cannot be greater than the number of dummy arguments.   [NEQNF]
call  NEQNF(Zr(i),A4(k),FCN(i,k),ERRREL,N,ITMAX,XGUESS,X,FNORM)
------^
C:\Program Files\Microsoft Visual Studio\MyProjects\tu2\tu2.f90(29) : Error: The type of the actual argument differs from the type of the dummy argument.   [A4]
call  NEQNF(Zr(i),A4(k),FCN(i,k),ERRREL,N,ITMAX,XGUESS,X,FNORM)
------------------^
C:\Program Files\Microsoft Visual Studio\MyProjects\tu2\tu2.f90(29) : Error: The type of the actual argument differs from the type of the dummy argument.   [FCN]
call  NEQNF(Zr(i),A4(k),FCN(i,k),ERRREL,N,ITMAX,XGUESS,X,FNORM)
------------------------^
C:\Program Files\Microsoft Visual Studio\MyProjects\tu2\tu2.f90(29) : Error: The type of the actual argument differs from the type of the dummy argument.   [0.0001]
call  NEQNF(Zr(i),A4(k),FCN(i,k),ERRREL,N,ITMAX,XGUESS,X,FNORM)
---------------------------------^
C:\Program Files\Microsoft Visual Studio\MyProjects\tu2\tu2.f90(29) : Error: The type of the actual argument differs from the type of the dummy argument.   [2]
call  NEQNF(Zr(i),A4(k),FCN(i,k),ERRREL,N,ITMAX,XGUESS,X,FNORM)
----------------------------------------^
C:\Program Files\Microsoft Visual Studio\MyProjects\tu2\tu2.f90(29) : Error: The type of the actual argument differs from the type of the dummy argument.   [100]
call  NEQNF(Zr(i),A4(k),FCN(i,k),ERRREL,N,ITMAX,XGUESS,X,FNORM)
------------------------------------------^
C:\Program Files\Microsoft Visual Studio\MyProjects\tu2\tu2.f90(29) : Error: The type of the actual argument differs from the type of the dummy argument.   [XGUESS]
call  NEQNF(Zr(i),A4(k),FCN(i,k),ERRREL,N,ITMAX,XGUESS,X,FNORM)
------------------------------------------------^
C:\Program Files\Microsoft Visual Studio\MyProjects\tu2\tu2.f90(29) : Error: This actual argument must be the name of a user subroutine or the name of intrinsic subroutine.   [ZR]
call  NEQNF(Zr(i),A4(k),FCN(i,k),ERRREL,N,ITMAX,XGUESS,X,FNORM)
------------^
C:\Program Files\Microsoft Visual Studio\MyProjects\tu2\tu2.f90(29) : Error: The shape matching rules of actual arguments and dummy arguments have been violated.   [XGUESS]
call  NEQNF(Zr(i),A4(k),FCN(i,k),ERRREL,N,ITMAX,XGUESS,X,FNORM)
------------------------------------------------^
C:\Program Files\Microsoft Visual Studio\MyProjects\tu2\tu2.f90(28) : Error: An unterminated block exists.
do k=1,A4s
^
C:\Program Files\Microsoft Visual Studio\MyProjects\tu2\tu2.f90(27) : Error: An unterminated block exists.
do i=1,Zrs
^
Error executing df.exe.

tu2.obj - 11 error(s), 0 warning(s)

在此谢过！！

fcode

你的代码问题比较多。

1.主程序中，少两个 End Do
2.我在 IMSL 里找到的 NEQNF 是 7 个参数，而你给了 9 个。且第一个参数应该为一个外部函数，你给的是 Zr(i)

你先改这两个问题吧。


另外，我很好奇，这代码是你自己书写的吗？
为什么要用 common 呢？而且还用了 pointer
一个很新潮的指针，和一个老掉牙的 common 语句，放一起真是...很另类

ttt

fcode 发表于 2014-4-26 12:45
你的代码问题比较多。

1.主程序中，少两个 End Do

见笑了，我改改，那个啥我是模仿来的，嗯，谢谢你

我是处手

fcode 发表于 2014-4-26 12:45
你的代码问题比较多。

1.主程序中，少两个 End Do

这是新手,我也是新学的,其中关于传递参数,指针有这个功能,common也有这个功能,但是common又有面临参数很多的问题,所以我们会几种方法混合用

111

您好，请问您求解嵌套非线性方程组的问题解决了吗？