NASIM.1111
عضو جدید
سلام خواهش میکنم کمکم کنید معادله X[SUP]3[/SUP]/3+X[SUP]2[/SUP]/2+2X+1/3 هم حل بشه هم ریشه هاش به دست بیاد تو نرم افزار متلب
سلام خواهش میکنم کمکم کنید معادله X[SUP]3[/SUP]/3+X[SUP]2[/SUP]/2+2X+1/3 هم حل بشه هم ریشه هاش به دست بیاد تو نرم افزار متلب
چگونه می توان در متلب برنامه ای نوشت که با داشتن مشتق یک تابع در هر نقطه مقدار تابع را در هر نقطه حساب کند؟(در واقع انتگرال)
ممنون میشم اگه راهنماییم کنید.
خیلی ساده با sort کردن میشه این کار رو انجام داد.کسی هست که جواب بده؟؟
1 سوال برنامه نویسی توو محیط متلب دارم..
ﺳﻮال.ﺑﺪون اﺳﺘﻔﺎده از تابع max، تابعی ﺑﻨﻮﯾﺴﯿﺪ ﮐﻪ ﻧﺎم یک ﻣﺎﺗﺮﯾﺲ را ﺑﻪ ﻋﻨﻮان ﭘﺎراﻣﺘﺮ ورودی درﯾﺎﻓﺖ ﮐﻨﺪ و ﻋﻨﺼﺮ ﻣﺎﮐﺰﯾﻤﻢ و اﻧﺪﯾﺴﻬﺎیﺳﻄﺮی وﺳﺘﻮنی آنرا ﺑﺮﮔﺮداﻧﺪ.
ممنون
حلقه بذارید که درایههای هر سطر یا ستون رو دو به دو با هم مقایسه کرده و اونا رو مرتب بکنه؛ تا زمانی که تمام درایههای اون سطر یا ستون مرتب شده باشن. حالا بسته به مقایسهتون (از کمتر به بیشتر یا بالعکس)، اولین درایه یا آخرین درایه هر سطر یا ستون، میشه مقدار ماکزیمم یا مینیمم اون سطر یا ستونبا حلقه ی for چطور میشه انجام داد؟
sort که ن..
با دستور if...مقایسه رو چطور انجام میده؟ با چه دستوری؟باید شمارنده تعریف کنم؟
خواهش میکنم دوست عزیز!ممنون از راهنماییتون
SOLUTION OF A SET OF ORDINARY DIFFERENTIAL EQUATIONS
USING THE FOURTH ORDER ADAMS-BASHFORTH, ADAMS-MOULTON
PREDICTOR-CORRECTOR METHOD WITH AN ITERATIVE ALGORITHM.
THE INTEGRATION IS STARTED USING CLASSIC FOURTH ORDER
RUNGE-KUTTA METHOD AND THE ITERATION PROCEDURE IS THE
THE GAUSS-SEIDEL WITH RELAXATION.
PROGRAM NAME: ABAMPC4.FOR
EXAMPLE CASE: dY/dX = -.3 Y + .1 Z + .1 U
dZ/dX = -.2 Z + .1 U
dU/dX = -.1 U
WITH INITIAL CONDITIONS: AT X=0, Y=3, Z=2, U=1
ANALYTICAL SOLUTION IS : Y = EXP(-.1X) + EXP(-.2X) + EXP(-.3X)
Z = EXP(-.1X) + EXP(-.2X)
U = EXP(-.1X)
REAL X(51),Y(51),Z(51),U(51),DYDX(51),DZDX(51),DUDX(51)
REAL YS(51),ZS(51),US(51),DYDXS(51),DZDXS(51),DUDXS(51)
REAL YTRUE(51),ZTRUE(51),UTRUE(51),H,K1,K2,K3,K4
REAL L1,L2,L3,L4,M1,M2,M3,M4
OPEN(10,FILE='ABAMPC4.RES',STATUS='NEW')
SET TOLERANCE FOR ITERATION PROCEDURE
AND MAXIMUM NUMBER OF ITERATIONS
TOL=0.00000001
IMAX=50
SPECIFY THE RELAXATION FACTOR
WRITE(*,7)
FORMAT(' ENTER RELAXATION FACTOR')
READ(*,8) W
FORMAT(F4.2)
WRITE(10,9) W
FORMAT(' RELAXATION FACTOR: ',F4.2)
SPECIFY THE INITIAL CONDITION AND THE STEP SIZE.
X(1)=0.0
Y(1)=3.0
Z(1)=2.0
U(1)=1.0
H=0.01
USE RK FOR THE FIRST 4 INTERVALS TO START
THE PREDICTOR-CORRECTOR METHOD.
DO 100 J=2,5
I=J-1
XT=X(I)
YT=Y(I)
ZT=Z(I)
UT=U(I)
CALL DERIV (XT,YT,ZT,UT,F1,F2,F3)
DYDX(I) =F1
DZDX(I) =F2
DUDX(I) =F3
K1=H*F1
L1=H*F2
M1=H*F3
XT=X(I)+(H/2.)
YT=Y(I)+(K1/2.)
ZT=Z(I)+(L1/2.)
UT=U(I)+(M1/2.)
CALL DERIV(XT,YT,ZT,UT,F1,F2,F3)
K2=H*F1
L2=H*F2
M2=H*F3
XT=X(I)+(H/2.)
YT=Y(I)+(K2/2.)
ZT=Z(I)+(L2/2.)
UT=U(I)+(M2/2.)
CALL DERIV(XT,YT,ZT,UT,F1,F2,F3)
K3=H*F1
L3=H*F2
M3=H*F3
XT=X(I)+H
YT=Y(I)+K3
ZT=Z(I)+L3
UT=U(I)+M3
CALL DERIV(XT,YT,ZT,UT,F1,F2,F3)
K4=H*F1
L4=H*F2
M4=H*F3
X(J)=X(I)+H
Y(J)=Y(I)+((1./6.)*(K1+(2.*K2)+(2.*K3)+K4))
Z(J)=Z(I)+((1./6.)*(L1+(2.*L2)+(2.*L3)+L4))
U(J)=U(I)+((1./6.)*(M1+(2.*M2)+(2.*M3)+M4))
CONTINUE
J=5
XT=X(J)
YT=Y(J)
ZT=Z(J)
UT=U(J)
CALL DERIV(XT,YT,ZT,UT,F1,F2,F3)
DYDX(J) =F1
DZDX(J) =F2
DUDX(J) =F3
END OF THE STARTER SEGMENT
DO 800 N=5,50
ITER=1
N1=N-1
N2=N-2
N3=N-3
NN=N+1
X(NN)=X(N)+H
PREDICTOR
TEMP=(55.*DYDX(N))-(59.*DYDX(N1))+(37.*DYDX(N2))-(9.*DYDX(N3))
YS(NN)=Y(N)+((H/24.)*TEMP)
TEMP=(55.*DZDX(N))-(59.*DZDX(N1))+(37.*DZDX(N2))-(9.*DZDX(N3))
ZS(NN)=Z(N)+((H/24.)*TEMP)
TEMP=(55.*DUDX(N))-(59.*DUDX(N1))+(37.*DUDX(N2))-(9.*DUDX(N3))
US(NN)=U(N)+((H/24.)*TEMP)
USE GAUSS-SEIDEL ITERATION ON THE CORRECTOR
XT=X(NN)
YT=YS(NN)
ZT=ZS(NN)
UT=US(NN)
CALL DERIV(XT,YT,ZT,UT,F1,F2,F3)
DYDXS(NN)=F1
TEMP=(9.*DYDXS(NN))+(19.*DYDX(N))-(5.*DYDX(N1))+DYDX(N2)
Y(NN)=Y(N)+((H/24.)*TEMP)
Y(NN)=YS(NN)+(W*(Y(NN)-YS(NN)))
YT=Y(NN)
CALL DERIV(XT,YT,ZT,UT,F1,F2,F3)
DZDXS(NN)=F2
TEMP=(9.*DZDXS(NN))+(19.*DZDX(N))-(5.*DZDX(N1))+DZDX(N2)
Z(NN)=Z(N)+((H/24.)*TEMP)
Z(NN)=ZS(NN)+(W*(Z(NN)-ZS(NN)))
ZT=Z(NN)
CALL DERIV(XT,YT,ZT,UT,F1,F2,F3)
DUDXS(NN)=F3
TEMP=(9.*DUDXS(NN))+(19.*DUDX(N))-(5.*DUDX(N1))+DUDX(N2)
U(NN)=U(N)+((H/24.)*TEMP)
U(NN)=US(NN)+(W*(U(NN)-US(NN)))
ITER=ITER+1
CHECK FOR CONVERGENCE
EY=((Y(NN)-YS(NN))/Y(NN))*100.
EZ=((Z(NN)-ZS(NN))/Z(NN))*100.
EU=((U(NN)-US(NN))/U(NN))*100.
IF(ITER .GT. IMAX) GOTO 500
IF(EY .GT. TOL) GOTO 600
IF(EZ .GT. TOL) GOTO 600
IF(EU .GT. TOL) GOTO 600
GOTO 700
WRITE(10,550)
FORMAT(' CONVERGENCE NOT ACHIEVED IN SPECIFIED ITERATIONS')
STOP
YS(NN)=Y(NN)
ZS(NN)=Z(NN)
US(NN)=U(NN)
GOTO 400
XT=X(NN)
YT=Y(NN)
ZT=Z(NN)
UT=U(NN)
CALL DERIV(XT,YT,ZT,UT,F1,F2,F3)
DYDX(NN) =F1
DZDX(NN) =F2
DUDX(NN) =F3
CONTINUE
WRITE RESULTS
DO 900 I=1,51
YTRUE(I)=EXP(-0.1*X(I))+EXP(-0.2*X(I))+EXP(-0.3*X(I))
ZTRUE(I)=EXP(-0.1*X(I))+EXP(-0.2*X(I))
UTRUE(I)=EXP(-0.1*X(I))
WRITE(10,850) X(I),Y(I),YTRUE(I),Z(I),ZTRUE(I),U(I),UTRUE(I)
FORMAT(F4.2,' ',2F11.6,' ',2F11.6,' ',2F11.6)
CONTINUE
STOP
END
SUBROUTINE DERIV(XT,YT,ZT,UT,F1,F2,F3)
F1=(-0.3*YT)+(0.1*ZT)+(0.1*UT)
F2=(-0.2*ZT)+(0.1*UT)
F3=(-0.1*UT)
RETURN
END
به نام خدا
سلام
این برنامه تو جزوه خواهرم بوده که تایپش کرده. انگار روش حل معادله مرتبه چهارم هست. اما خیلی سینتکس ارور داره
ممکنه راهنمایی کنید و بگین اشکالاتش کجاست
با تشکر
کد:SOLUTION OF A SET OF ORDINARY DIFFERENTIAL EQUATIONS USING THE FOURTH ORDER ADAMS-BASHFORTH, ADAMS-MOULTON PREDICTOR-CORRECTOR METHOD WITH AN ITERATIVE ALGORITHM. THE INTEGRATION IS STARTED USING CLASSIC FOURTH ORDER RUNGE-KUTTA METHOD AND THE ITERATION PROCEDURE IS THE THE GAUSS-SEIDEL WITH RELAXATION. PROGRAM NAME: ABAMPC4.FOR EXAMPLE CASE: dY/dX = -.3 Y + .1 Z + .1 U dZ/dX = -.2 Z + .1 U dU/dX = -.1 U WITH INITIAL CONDITIONS: AT X=0, Y=3, Z=2, U=1 ANALYTICAL SOLUTION IS : Y = EXP(-.1X) + EXP(-.2X) + EXP(-.3X) Z = EXP(-.1X) + EXP(-.2X) U = EXP(-.1X) REAL X(51),Y(51),Z(51),U(51),DYDX(51),DZDX(51),DUDX(51) REAL YS(51),ZS(51),US(51),DYDXS(51),DZDXS(51),DUDXS(51) REAL YTRUE(51),ZTRUE(51),UTRUE(51),H,K1,K2,K3,K4 REAL L1,L2,L3,L4,M1,M2,M3,M4 OPEN(10,FILE='ABAMPC4.RES',STATUS='NEW') SET TOLERANCE FOR ITERATION PROCEDURE AND MAXIMUM NUMBER OF ITERATIONS TOL=0.00000001 IMAX=50 SPECIFY THE RELAXATION FACTOR WRITE(*,7) FORMAT(' ENTER RELAXATION FACTOR') READ(*,8) W FORMAT(F4.2) WRITE(10,9) W FORMAT(' RELAXATION FACTOR: ',F4.2) SPECIFY THE INITIAL CONDITION AND THE STEP SIZE. X(1)=0.0 Y(1)=3.0 Z(1)=2.0 U(1)=1.0 H=0.01 USE RK FOR THE FIRST 4 INTERVALS TO START THE PREDICTOR-CORRECTOR METHOD. DO 100 J=2,5 I=J-1 XT=X(I) YT=Y(I) ZT=Z(I) UT=U(I) CALL DERIV (XT,YT,ZT,UT,F1,F2,F3) DYDX(I) =F1 DZDX(I) =F2 DUDX(I) =F3 K1=H*F1 L1=H*F2 M1=H*F3 XT=X(I)+(H/2.) YT=Y(I)+(K1/2.) ZT=Z(I)+(L1/2.) UT=U(I)+(M1/2.) CALL DERIV(XT,YT,ZT,UT,F1,F2,F3) K2=H*F1 L2=H*F2 M2=H*F3 XT=X(I)+(H/2.) YT=Y(I)+(K2/2.) ZT=Z(I)+(L2/2.) UT=U(I)+(M2/2.) CALL DERIV(XT,YT,ZT,UT,F1,F2,F3) K3=H*F1 L3=H*F2 M3=H*F3 XT=X(I)+H YT=Y(I)+K3 ZT=Z(I)+L3 UT=U(I)+M3 CALL DERIV(XT,YT,ZT,UT,F1,F2,F3) K4=H*F1 L4=H*F2 M4=H*F3 X(J)=X(I)+H Y(J)=Y(I)+((1./6.)*(K1+(2.*K2)+(2.*K3)+K4)) Z(J)=Z(I)+((1./6.)*(L1+(2.*L2)+(2.*L3)+L4)) U(J)=U(I)+((1./6.)*(M1+(2.*M2)+(2.*M3)+M4)) CONTINUE J=5 XT=X(J) YT=Y(J) ZT=Z(J) UT=U(J) CALL DERIV(XT,YT,ZT,UT,F1,F2,F3) DYDX(J) =F1 DZDX(J) =F2 DUDX(J) =F3 END OF THE STARTER SEGMENT DO 800 N=5,50 ITER=1 N1=N-1 N2=N-2 N3=N-3 NN=N+1 X(NN)=X(N)+H PREDICTOR TEMP=(55.*DYDX(N))-(59.*DYDX(N1))+(37.*DYDX(N2))-(9.*DYDX(N3)) YS(NN)=Y(N)+((H/24.)*TEMP) TEMP=(55.*DZDX(N))-(59.*DZDX(N1))+(37.*DZDX(N2))-(9.*DZDX(N3)) ZS(NN)=Z(N)+((H/24.)*TEMP) TEMP=(55.*DUDX(N))-(59.*DUDX(N1))+(37.*DUDX(N2))-(9.*DUDX(N3)) US(NN)=U(N)+((H/24.)*TEMP) USE GAUSS-SEIDEL ITERATION ON THE CORRECTOR XT=X(NN) YT=YS(NN) ZT=ZS(NN) UT=US(NN) CALL DERIV(XT,YT,ZT,UT,F1,F2,F3) DYDXS(NN)=F1 TEMP=(9.*DYDXS(NN))+(19.*DYDX(N))-(5.*DYDX(N1))+DYDX(N2) Y(NN)=Y(N)+((H/24.)*TEMP) Y(NN)=YS(NN)+(W*(Y(NN)-YS(NN))) YT=Y(NN) CALL DERIV(XT,YT,ZT,UT,F1,F2,F3) DZDXS(NN)=F2 TEMP=(9.*DZDXS(NN))+(19.*DZDX(N))-(5.*DZDX(N1))+DZDX(N2) Z(NN)=Z(N)+((H/24.)*TEMP) Z(NN)=ZS(NN)+(W*(Z(NN)-ZS(NN))) ZT=Z(NN) CALL DERIV(XT,YT,ZT,UT,F1,F2,F3) DUDXS(NN)=F3 TEMP=(9.*DUDXS(NN))+(19.*DUDX(N))-(5.*DUDX(N1))+DUDX(N2) U(NN)=U(N)+((H/24.)*TEMP) U(NN)=US(NN)+(W*(U(NN)-US(NN))) ITER=ITER+1 CHECK FOR CONVERGENCE EY=((Y(NN)-YS(NN))/Y(NN))*100. EZ=((Z(NN)-ZS(NN))/Z(NN))*100. EU=((U(NN)-US(NN))/U(NN))*100. IF(ITER .GT. IMAX) GOTO 500 IF(EY .GT. TOL) GOTO 600 IF(EZ .GT. TOL) GOTO 600 IF(EU .GT. TOL) GOTO 600 GOTO 700 WRITE(10,550) FORMAT(' CONVERGENCE NOT ACHIEVED IN SPECIFIED ITERATIONS') STOP YS(NN)=Y(NN) ZS(NN)=Z(NN) US(NN)=U(NN) GOTO 400 XT=X(NN) YT=Y(NN) ZT=Z(NN) UT=U(NN) CALL DERIV(XT,YT,ZT,UT,F1,F2,F3) DYDX(NN) =F1 DZDX(NN) =F2 DUDX(NN) =F3 CONTINUE WRITE RESULTS DO 900 I=1,51 YTRUE(I)=EXP(-0.1*X(I))+EXP(-0.2*X(I))+EXP(-0.3*X(I)) ZTRUE(I)=EXP(-0.1*X(I))+EXP(-0.2*X(I)) UTRUE(I)=EXP(-0.1*X(I)) WRITE(10,850) X(I),Y(I),YTRUE(I),Z(I),ZTRUE(I),U(I),UTRUE(I) FORMAT(F4.2,' ',2F11.6,' ',2F11.6,' ',2F11.6) CONTINUE STOP END SUBROUTINE DERIV(XT,YT,ZT,UT,F1,F2,F3) F1=(-0.3*YT)+(0.1*ZT)+(0.1*UT) F2=(-0.2*ZT)+(0.1*UT) F3=(-0.1*UT) RETURN END
به نام خداسلام
دوست عزيز اين كد با زبان متلب نيست. vb هست احتمالا يا فورترن
به نام خدا
سلام
این برنامه تو جزوه خواهرم بوده که تایپش کرده. انگار روش حل معادله مرتبه چهارم هست. اما خیلی سینتکس ارور داره
ممکنه راهنمایی کنید و بگین اشکالاتش کجاست
با تشکر
کد:SOLUTION OF A SET OF ORDINARY DIFFERENTIAL EQUATIONS USING THE FOURTH ORDER ADAMS-BASHFORTH, ADAMS-MOULTON PREDICTOR-CORRECTOR METHOD WITH AN ITERATIVE ALGORITHM. THE INTEGRATION IS STARTED USING CLASSIC FOURTH ORDER RUNGE-KUTTA METHOD AND THE ITERATION PROCEDURE IS THE THE GAUSS-SEIDEL WITH RELAXATION. PROGRAM NAME: ABAMPC4.FOR EXAMPLE CASE: dY/dX = -.3 Y + .1 Z + .1 U dZ/dX = -.2 Z + .1 U dU/dX = -.1 U WITH INITIAL CONDITIONS: AT X=0, Y=3, Z=2, U=1 ANALYTICAL SOLUTION IS : Y = EXP(-.1X) + EXP(-.2X) + EXP(-.3X) Z = EXP(-.1X) + EXP(-.2X) U = EXP(-.1X) REAL X(51),Y(51),Z(51),U(51),DYDX(51),DZDX(51),DUDX(51) REAL YS(51),ZS(51),US(51),DYDXS(51),DZDXS(51),DUDXS(51) REAL YTRUE(51),ZTRUE(51),UTRUE(51),H,K1,K2,K3,K4 REAL L1,L2,L3,L4,M1,M2,M3,M4 OPEN(10,FILE='ABAMPC4.RES',STATUS='NEW') SET TOLERANCE FOR ITERATION PROCEDURE AND MAXIMUM NUMBER OF ITERATIONS TOL=0.00000001 IMAX=50 SPECIFY THE RELAXATION FACTOR WRITE(*,7) FORMAT(' ENTER RELAXATION FACTOR') READ(*,8) W FORMAT(F4.2) WRITE(10,9) W FORMAT(' RELAXATION FACTOR: ',F4.2) SPECIFY THE INITIAL CONDITION AND THE STEP SIZE. X(1)=0.0 Y(1)=3.0 Z(1)=2.0 U(1)=1.0 H=0.01 USE RK FOR THE FIRST 4 INTERVALS TO START THE PREDICTOR-CORRECTOR METHOD. DO 100 J=2,5 I=J-1 XT=X(I) YT=Y(I) ZT=Z(I) UT=U(I) CALL DERIV (XT,YT,ZT,UT,F1,F2,F3) DYDX(I) =F1 DZDX(I) =F2 DUDX(I) =F3 K1=H*F1 L1=H*F2 M1=H*F3 XT=X(I)+(H/2.) YT=Y(I)+(K1/2.) ZT=Z(I)+(L1/2.) UT=U(I)+(M1/2.) CALL DERIV(XT,YT,ZT,UT,F1,F2,F3) K2=H*F1 L2=H*F2 M2=H*F3 XT=X(I)+(H/2.) YT=Y(I)+(K2/2.) ZT=Z(I)+(L2/2.) UT=U(I)+(M2/2.) CALL DERIV(XT,YT,ZT,UT,F1,F2,F3) K3=H*F1 L3=H*F2 M3=H*F3 XT=X(I)+H YT=Y(I)+K3 ZT=Z(I)+L3 UT=U(I)+M3 CALL DERIV(XT,YT,ZT,UT,F1,F2,F3) K4=H*F1 L4=H*F2 M4=H*F3 X(J)=X(I)+H Y(J)=Y(I)+((1./6.)*(K1+(2.*K2)+(2.*K3)+K4)) Z(J)=Z(I)+((1./6.)*(L1+(2.*L2)+(2.*L3)+L4)) U(J)=U(I)+((1./6.)*(M1+(2.*M2)+(2.*M3)+M4)) CONTINUE J=5 XT=X(J) YT=Y(J) ZT=Z(J) UT=U(J) CALL DERIV(XT,YT,ZT,UT,F1,F2,F3) DYDX(J) =F1 DZDX(J) =F2 DUDX(J) =F3 END OF THE STARTER SEGMENT DO 800 N=5,50 ITER=1 N1=N-1 N2=N-2 N3=N-3 NN=N+1 X(NN)=X(N)+H PREDICTOR TEMP=(55.*DYDX(N))-(59.*DYDX(N1))+(37.*DYDX(N2))-(9.*DYDX(N3)) YS(NN)=Y(N)+((H/24.)*TEMP) TEMP=(55.*DZDX(N))-(59.*DZDX(N1))+(37.*DZDX(N2))-(9.*DZDX(N3)) ZS(NN)=Z(N)+((H/24.)*TEMP) TEMP=(55.*DUDX(N))-(59.*DUDX(N1))+(37.*DUDX(N2))-(9.*DUDX(N3)) US(NN)=U(N)+((H/24.)*TEMP) USE GAUSS-SEIDEL ITERATION ON THE CORRECTOR XT=X(NN) YT=YS(NN) ZT=ZS(NN) UT=US(NN) CALL DERIV(XT,YT,ZT,UT,F1,F2,F3) DYDXS(NN)=F1 TEMP=(9.*DYDXS(NN))+(19.*DYDX(N))-(5.*DYDX(N1))+DYDX(N2) Y(NN)=Y(N)+((H/24.)*TEMP) Y(NN)=YS(NN)+(W*(Y(NN)-YS(NN))) YT=Y(NN) CALL DERIV(XT,YT,ZT,UT,F1,F2,F3) DZDXS(NN)=F2 TEMP=(9.*DZDXS(NN))+(19.*DZDX(N))-(5.*DZDX(N1))+DZDX(N2) Z(NN)=Z(N)+((H/24.)*TEMP) Z(NN)=ZS(NN)+(W*(Z(NN)-ZS(NN))) ZT=Z(NN) CALL DERIV(XT,YT,ZT,UT,F1,F2,F3) DUDXS(NN)=F3 TEMP=(9.*DUDXS(NN))+(19.*DUDX(N))-(5.*DUDX(N1))+DUDX(N2) U(NN)=U(N)+((H/24.)*TEMP) U(NN)=US(NN)+(W*(U(NN)-US(NN))) ITER=ITER+1 CHECK FOR CONVERGENCE EY=((Y(NN)-YS(NN))/Y(NN))*100. EZ=((Z(NN)-ZS(NN))/Z(NN))*100. EU=((U(NN)-US(NN))/U(NN))*100. IF(ITER .GT. IMAX) GOTO 500 IF(EY .GT. TOL) GOTO 600 IF(EZ .GT. TOL) GOTO 600 IF(EU .GT. TOL) GOTO 600 GOTO 700 WRITE(10,550) FORMAT(' CONVERGENCE NOT ACHIEVED IN SPECIFIED ITERATIONS') STOP YS(NN)=Y(NN) ZS(NN)=Z(NN) US(NN)=U(NN) GOTO 400 XT=X(NN) YT=Y(NN) ZT=Z(NN) UT=U(NN) CALL DERIV(XT,YT,ZT,UT,F1,F2,F3) DYDX(NN) =F1 DZDX(NN) =F2 DUDX(NN) =F3 CONTINUE WRITE RESULTS DO 900 I=1,51 YTRUE(I)=EXP(-0.1*X(I))+EXP(-0.2*X(I))+EXP(-0.3*X(I)) ZTRUE(I)=EXP(-0.1*X(I))+EXP(-0.2*X(I)) UTRUE(I)=EXP(-0.1*X(I)) WRITE(10,850) X(I),Y(I),YTRUE(I),Z(I),ZTRUE(I),U(I),UTRUE(I) FORMAT(F4.2,' ',2F11.6,' ',2F11.6,' ',2F11.6) CONTINUE STOP END SUBROUTINE DERIV(XT,YT,ZT,UT,F1,F2,F3) F1=(-0.3*YT)+(0.1*ZT)+(0.1*UT) F2=(-0.2*ZT)+(0.1*UT) F3=(-0.1*UT) RETURN END
سلامسلام دوستان
من یک برنامه نوشتم که طی اون به 3معادله خطی منتهی میشه. حالا مخوام این 3معادله رو که بر حسب x y z هستن با شرط x^2+y^2+z^2=1 حل کنه. اما وقتی با تابع solve مخوام حل کنه خطا میده و میگه 4معادله و 3 مجهول رو حل نمی کنم.
حالا باید چه کار کنم؟ لطفا راهنماییم کنید ممنون.
سلام
من یه فایل اکسل دارم
میخوام در متلب ایمپورت کنم
در نرم افزار متلب file/open رو می زنم
بعد اون فایل اکسل رو انتخاب می کنم
ایمپورت میشه اما فایل اکسل داخل نرم افزار متلب بد نشون داده میشه
عدد ها رو نشون میده اما حروف رو نه
مشکل از کجاست
چکارش باید بکنم
کسی می دونه؟؟؟
Thread starter | عنوان | تالار | پاسخ ها | تاریخ |
---|---|---|---|---|
F | سوالات و مشکلات در متلب | MATLAB | 1 | |
N | راهنمایی پیاده سازی کلاسبند جنگل تصادفی و خوشه بند DBSCAN در متلب | MATLAB | 0 | |
ب | کمک در برنامه نویسی متلب...فوری | MATLAB | 0 | |
A | شبیه سازی متغییر گسسته وپیوسته باهم برای یک مسئله در متلب 2014 | MATLAB | 0 | |
A | درخواست کمک در پروژه متلب.فوری | MATLAB | 4 |