回 帖 发 新 帖 刷新版面

主题:求教:这是否FORTRAN?求助转换为MATLAB

求教:这是BASIC还是FORTRAN? 能否转换为MATLAB或者标准的FORTRAN执行?

感谢指教。这是片段,有整个程序,在引用别人的发表资料。为博士论文的数学模型拟合,有酬谢

412 NEXT X:LPRINT:LPRINT

414 LPRINT TAB(6) "AIR TEMP. ";:FOR X=0 TO XT STEP XPR

416 LPRINT USING F1$;(TA(X)-32)/l.8;:NEXT X:LPRINT

418 LPRINT TAB(6) "PELLET TEMP.";:FOR X=0 TO XT STEP XPR

420 LPRINT USING F1$;(TC(X,0)-32)/1.8;:NEXT X:LPRINT

422 LPRINT TAB(5) " MOIS.CONT. ";

424 FOR X=0 TO XT STEP XPR:LPRINT USING F1$;MC(X,O)*100;:NEXT X

426 LPRINT:LPRINT TAB(5) " ABS. HUMID. ";

428 FOR X=0 TO XT STEP XPR:LPRINT USING Fl$;H(X);:NEXT X:LPRINT

430 LPRINT:LPRINT TAB(10) STRING;(53,"*"): LPRINT

432 REM

434 CPRI=0:PTI=1/(DT*60)

436 REM TIME LOOP

438 FOR L = DT TO TI STEP DT

440 CPRI=CPRI+L:SMC=0: STA=0

442 LOCATE 22,5

444 PRINT " ELAPSED TIME = ";:PRINT USING "##.##"; L*60;

446 PRINT " min.:";:PRINT " COMPUTING CONDITIONS AT "

448 REM DEPTH LOOP

450 FOR X = 0 TO XT

452 LOCATE 22,56

454 PRINT USING " ##.##"; ((X*30.48*BED)/XT);: PRINT " cm."

456 REM EQUILIBRIUM MOISTURE CONTENT EQUATIONS

458 TPC=(TAV(X)-32)*(5/9)

460 REM NELLIST EQUATION

462 REM EMCD(X)=.191-.055*LOG(1-RH(X))-.028*LOG(TPC)

464 REM HENDERSON EQUATION

466 REM EMCD(X)=(-LOG(1-RH(X))/(6.66*(TPC+22.12)))^(1/3.11)

468 REM CHUNG-PFOST EQUATION

470 EMCD(X)=.277-.042*LOG(-(TPC+13.3)*LOG(RH(X)))

472 REM DIFFUSION COEFFICIENT EQUATION

474 DIC=1.015E-05*EXP(-547/(TPC+273.15))

476 REM MOISTURE CONTENT DIFFERENTIAL EON.

478 SUMM=O:FDR I = 1 TO KL

480 SUMO=B(I)*EXP(AR(I)*DIC*L)

482 SUMM=SUKM+SUMC:NEXT I
484 MCEQ(X)=EMCD(X)+(MCID-EMCD(X))*SUMM

486 REM MOISTURE CONTENT - FINITE-DIFF. FORMULAS

488 MCT=MC(X,O):EMC=EMCD(X)

490 AM=(DIC*DT)/(DELR*DELR)

492 GOSUB 608

494 MC(X,1)=MT(OgNI-1):MS(X)=MTIM,NII:MAV(X,1)=MAVE

496 SMC=SMC+MAV(X,1)

498 REM ABSOLUTE HUMIDITY EQUATION

500 H(X+1)=H(X)-C1*(MC(X,1)-MC(X,0))

502 REM PELLET TEMPERATURE EQUATIONS

504 P=MAV(X,O)

506 K=.1133-2.936*(P*P)+25.44*(P*P*P)-38.71*(P*P*P*P)

508 CP=4180*(.343+P)/(1+P):CPE=CP/4186.69

510 A=(K*DELT)/(POEM*CP*(DELR*DELR))

512 REM

514 TPK=TPT(X,O)+459.69

516 HFG=(1056.5-.55*(TPK-520))*(1+23*EXP(-4*TPK*P))

518 DTEMP=TA(X)-TS(X):AUX=CPE+CW*P:AUXI=HFG+CV*DTEMP

520 DHT=H(X+1)-H(X)


回复列表 (共2个回复)

沙发

basic

很古老了

板凳

不是 Fortran 语言编写的,按您的说法,二选一的话,应该是 Basic。

我来回复

您尚未登录,请登录后再回复。点此登录或注册