SOS: gnuplot fdtd的一个问题求助 perl vs python

我用perl和python写了相同功能的一段程序，计算一维fdtd,用gnuplot动态显示，可是python的数据没有显示出来，看横纵坐标的变化数据是正确收到了的，如最后的图片，求大神指点，谢谢。

 #!/usr/bin/perl

 # Author : Leon  Email: yangli0534@gmail.com
 # fdtd simulation , plotting with gnuplot, writting in perl
 # perl and gnuplot software packages should be installed before running this program
 # 1d fdtd with absorbing boundary and TFSF boundary between [49] and [50]
 # lossy dielectric material localted at > ez[150]
 #use Time::HiRes qw(sleep);
 #use autodie qw(:all);
 use strict;
 use warnings;

 #print "\@\n";
 #my $terminal = "";
 #    open GNUPLOT_TERM, "echo 'show terminal;' | gnuplot 2>&1 |";
 #    while (<GNUPLOT_TERM>) {
 #    if (m/terminal type is (\w+)/) {
 #        $terminal=$1;
 #    }
 #    }
 #    close GNUPLOT_TERM;

     # unfortunately, the wxt terminal type does not support positioning.
     # hardcode it...
 #    $terminal  = "x11";

 open my $PIPE ,"| gnuplot " || die "Can't initialize gnuplot number \n";

 print $PIPE "set size 0.85, 0.85\n";
 print $PIPE "set term png size 600, 400\n";

 my $title = "fdtd simulation by leon,yangli0534\\\\@"."gmail.com";
 print $PIPE "set terminal gif animate\n";# terminal type: png
 print $PIPE "set output \"fdtd_simulation_abc_tfsf_match_diel.gif\"\n";#output file name
 print $PIPE  "set title \"{/Times:Italic $title}\"\n";# title name and font
 #print $PIPE  "set title  \"fdtd simulation by leon,yangli0534\\\\@ gmail.com\"\n";# title name and font
 print $PIPE  "set title  font \",10\"  norotate tc rgb \"white\"\n";
 print $PIPE "unset key\n";
 print $PIPE "set tics textcolor rgb \"white\"\n";# text color
 print $PIPE "set border lc rgb \"orange\"\n";
 print $PIPE "set grid lc rgb\"orange\"\n";
 print $PIPE "set object 1 rectangle from screen 0,0 to screen 1,1 fc rgb \"gray10\" behind\n";#background color
 print $PIPE "set xlabel\" {/Times:Italic distance: wave length}\" tc rgb \"white\" \n";# xlabel
 print $PIPE "set ylabel\"{/Times:Italic amplitude: v}\" tc rgb \"white\"\n";#ylabel
 print $PIPE "set autoscale\n"; 

 my $size = ;#physical distance
 my @ez;#electric field
 my @hy;#magnetic field
 my @ceze;#
 my @cezh;#
 my @chye;#
 my @chyh;#
 my $LOSS_LAYER = ;

 #my @epsR;
 my $LOSS = 0.01;
 my $imp0 = 377.0;
 #initalization
 for (my $i = ; $i < $size; $i++){
     $ez[$i] = 0.0;
     $hy[$i] = 0.0;
     if ($i < )    {
         #$epsR[$i] = 1.0;
         $ceze[$i] = 1.0;
         $cezh[$i] = $imp0;
     }
     elsif($i < $LOSS_LAYER){
         #$epsR[$i] = 1.0;
         $ceze[$i] = 1.0;
         $cezh[$i] = $imp0/9.0;
     }
     else {
         #$epsR[$i] = 9.0;
         $ceze[$i] = (1.0-$LOSS)/(1.0+$LOSS);
         $cezh[$i] = $imp0 /  /(1.0+$LOSS);
     }
 } 

 for ( my $i = ; $i < $size; $i++){

     if($i < $LOSS_LAYER){
         #$epsR[$i] = 1.0;
         $chye[$i] = 1.0/$imp0;
         $chyh[$i] = 1.0;
     }
     else {
         #$epsR[$i] = 9.0;
         $chye[$i] = 1.0/ $imp0 /(1.0+$LOSS);
         $chyh[$i] = (1.0-$LOSS)/(1.0+$LOSS);
     }
 }
 my $qTime;
 my $MaxTime = ;
 my $pi = 3.141592563589793;
 print $PIPE "set xrange [0:$size-1]\n";
 my $mm = ;

 #do time stepping
 for($qTime = ; $qTime < $MaxTime; $qTime+=){

     # update magnetic field
     #$hy[$size - 1] = $hy[$size - 2 ];#abc
     for( $mm = ; $mm < $size - ; $mm++){
         $hy[$mm] = $hy[$mm]*$chyh[$mm] + ($ez[$mm+] - $ez[$mm])*$chye[$mm];
     }
     $hy[] -=exp(-($qTime - 30.0)*($qTime - 30.0)/100.0)/$imp0;
     # update electric field
     $ez[] = $ez[];#abc
     #$ez[$size-1] = $ez[$size-2];
     for( $mm = ; $mm < $size- ; $mm++){
         $ez[$mm] = $ez[$mm]*$ceze[$mm] + ($hy[$mm] - $hy[$mm-])*$cezh[$mm];
     }

     if($qTime %  == ){

         print $PIPE "plot \"-\" w l  lw 1.5 lc rgb \"green\"\n";
         my $cnt = ;
         for my $elem ( @ez) {
         #print " ".$elem;
             print $PIPE $cnt." ".$elem."\n";
             $cnt += ;
         }
         print $PIPE "e\n";
     }
     #hardwire a source
     $ez[] += exp(-($qTime +0.5-(-0.5)- 30.0)*($qTime +0.5-(-0.5)- 30.0)/100.0);
 }

 #print $PIPE "set terminal x11\n";

 print $PIPE "set output\n";

 close($PIPE);

正常结果

python的程序

 #!/usr/bin/env python

 import sys
 import math
 import os
 # Author : Leon  Email: yangli0534@gmail.com
 # fdtd simulation , plotting with gnuplot, writting in python
 # perl and gnuplot software packages should be installed before running this program
 # 1d fdtd with absorbing boundary and TFSF boundary between [49] and [50]
 # lossy dielectric material localted at > ez[150]
 gp = os.popen('gnuplot','w')
 #gp .write('plot sin(x)\n');
 #gp .flush();
 #p .close();
 gp.write('set size 0.85, 0.85\n')
 gp.write('set term png size 600, 400\n')
 title = 'fdtd simulation by leon,yangli0534\\\\@gmail.com'
 gp.write('set terminal gif animate\n')
 gp.write('set output "fdtd_simulation_abc_tfsf_match_diel_py.gif"\n')
 gp.write(''.join(['set title "{/Times:Italic ',title, '}"\n']))
 gp.write('set title  font ",10"  norotate tc rgb "white"\n')
 gp.write('unset key\n')
 gp.write('set tics textcolor rgb "white"\n')
 gp.write('set border lc rgb "orange"\n')
 gp.write('set grid lc rgb"orange"\n')
 gp.write('set object 1 rectangle from screen 0,0 to screen 1,1 fc rgb "gray10" behind\n')
 gp.write('set xlabel" {/Times:Italic distance: wave length}" tc rgb "white" \n')
 gp.write('set ylabel"{/Times:Italic amplitude: v}" tc rgb "white"\n')
 gp.write('set autoscale\n');

 size = 400#physical distance
 ez=size * [0.00]#electric field
 hy=size * [0.00]#magnetic field
 ceze=size * [0.00]#
 cezh=size * [0.00]#
 chye=size * [0.00]#
 chyh=size * [0.00]#
 imp0 = 377.00
 LOSS = 0.01
 LOSS_LAYER = 250
 MaxTime = 18000
 cnt = 0
 elem = 0.00000
 gp.write(''.join(['set xrange [0:',str(size),'-1]\n']));
 for i in range(0,size):
     ez[i] = 0.0
     hy[i] = 0.0
     if (i < 100):
         #$epsR[$i] = 1.0;
         ceze[i] = 1.0
         cezh[i] = imp0
     elif(i < LOSS_LAYER):
         #$epsR[$i] = 1.0;
         ceze[i] = 1.0
         cezh[i] = imp0/9.0
     else :
         #$epsR[$i] = 9.0;
         ceze[i] = (1.0-LOSS)/(1.0+LOSS)
         cezh[i] = imp0 / 9 /(1.0+LOSS)
     if( i < LOSS_LAYER):
         chye[i] = 1.0/imp0
         chyh[i] = 1.0
     else:
         chye[i] = 1.0/imp0/(1.0+LOSS)
         chyh[i] = (1.0-LOSS)/(1.0+LOSS)
 for qTime in range(0, MaxTime):
     # update magnetic field
     for mm in range(0, size-1):
         hy[mm] = hy[mm]*chyh[mm] + (ez[mm+1]-ez[mm])*chye[mm]
     hy[49] = hy[49]-math.exp(-(qTime - 30.0)*(qTime - 30.0)/100.0)/imp0
     # update electric field
     ez[0] = ez[1]#abc
     #$ez[$size-1] = $ez[$size-2];
     for mm in range(1, size-1):
         ez[mm] = ez[mm]*ceze[mm] + (hy[mm] - hy[mm-1])*cezh[mm]
     if(qTime % 20 == 0):
         gp.write('plot "-" w l  lw 1.5 lc rgb "white"\n')
         cnt = 0
         for elem in ez:
             #gp.write([cnt,' ',elem,'\n'])
             gp.write(''.join([str(cnt),' ',str(elem),'\n']))
             #gp.write(' ')
             #gp.write(str(elem))
             #gp.write('\n')
             cnt += 1
         gp.write('e\n')
     ez[50] = ez[50]+math.exp(-(qTime +0.5-(-0.5)- 30.0)*(qTime +0.5-(-0.5)- 30.0)/100.0);
 gp.write('set output\n')
 gp.close()

结果就是