| Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) | 
![]()  | 
	
| 		
			
			 | 
		#1 | 
| 
			
			 初级会员 
			
			
			
			注册日期: 2015-05-31 
				
				
				
					帖子: 1
				 
				
				
				声望力: 0 ![]()  | 
	
	
	
		
		
			
			 代码: 
	for npt=1:nt
 
        offset = ( ( npt - 1 ) * nx * ny * nz ) * 4 + hoffset ;
        fseek(fid,offset,'bof') ;
        for j=1:ny
            for i=1:nx
                val(j,i)=fread(fid,1,fbit);
            end
        end
        pcolor(x,y,val)
        colorbar('vert');
        caxis(colrange);
        colormap colmap
        titlestr=['Snapshot of ' vnm ' at ' ... 
                  '{\fontsize{16}{\bf ' ... 
                  num2str(double(t(npt)),'%06.3f') '}}s'];
        title(titlestr);
        colorbar('vert');
        caxis(colrange);
        colormap(colmap);
        drawnow
        F  = getframe(gcf);
        im = frame2im(F);
        [im, map] = rgb2ind(im, 256);
        if npt == 1;
            imwrite(A, map, fout, 'gif', 'LoopCount', Inf, 'DelayTime', 0.1);
        else
            imwrite(A, map, fout, 'gif', 'WriteMode', 'append', 'DelayTime', 0.1);
        end
    end
 | 
| 
		 | 
	
	
	
		
		
		
		
			 
		
		
		
		
		
		
		
			
		
		
		
	 |