Matlab处理log文件

发表于 由

一般文本类文件都用java或python处理方便很多，此处记录一下matlab 处理文本文件的一些语法。

例如，程序记录下的日志文件为：dl4ss_output.log_2017_03_10_210643，其部分内容为：

  1. [Epoch: 1] - SDR:1.308344, SIR:2.063156, SAR:-1.509836, NSDR:2.908430  
  2. Epoch 1/1  
  3. 50/50 [==============================] - 178s - loss: 1.9203 - val_loss: 1.7383  
  4. Epoch 1/1  
  5. 50/50 [==============================] - 176s - loss: 1.6208 - val_loss: 1.5113  
  6. Epoch 1/1  
  7. 50/50 [==============================] - 176s - loss: 1.4212 - val_loss: 1.3334  
  8. Epoch 1/1  
  9. 50/50 [==============================] - 175s - loss: 1.2612 - val_loss: 1.1911  
  10. Epoch 1/1  
  11. 50/50 [==============================] - 176s - loss: 1.1320 - val_loss: 1.0740  
  12. [Epoch: 6] - SDR:2.128937, SIR:3.154744, SAR:0.774681, NSDR:3.729023  
  13. Epoch 1/1  
  14. 50/50 [==============================] - 175s - loss: 1.0312 - val_loss: 0.9913  

现在想对该日志文件进行处理，并将Epoch迭代过程中评估的结果值绘制成图，方法如下：

  1. clc   
  2. clear  
  3. colours = ’rbckgcmbkgcmbkgcmbgrcmykbgrcmykb’;  % ’bgrcmykw’ modified by jacoxu@msn.com   
  4. symbols = ’.x^os*.dvph><x+^os*.dvph><+’;  % ’.ox+*sdv^<>ph’ %参考 http://blog.sina.com.cn/s/blog_618af1950100kdi2.html   
  5. linetypes = {‘-’,'-’,'-’,'-’,'-’,'-’,'-’,'-’,'-’,'-’,'-’,'-’,'-’,'-’,'-’,'-’,'-’,'-’,'-’,'-’,'-’,'-’};  % {‘-’,':’,'-.’,'–’}   
  6. figureSize = [300,200,500,300];   
  7. lineWidth = 1.5;   
  8. para_list = {‘SDR’, ’SIR’, ’SAR’, ’NSDR’};   
  9.   
  10. log_file = ’dl4ss_output.log_2017_03_10_210643′;   
  11. log_file_fd = fopen(log_file);   
  12. epoch_idx = 0;   
  13. epoch_num = [];   
  14. sdr_value = [];   
  15. sir_value = [];   
  16. sar_value = [];   
  17. nsdr_value = [];   
  18.   
  19. % 开始循环读取log文件   
  20. while true   
  21.     % 获取当前行内容   
  22.     line = fgets(log_file_fd);   
  23.     % 如果当前行没有内容了   
  24.     if ischar(line) == false   
  25.         break  
  26.     end   
  27.     % 判断当前行中是否有评估值   
  28.     if isempty(strfind(line, ’NSDR:’))   
  29.         continue  
  30.     end   
  31.     epoch_idx = epoch_idx + 1;   
  32.     results = regexp(line, ’ - ’, ’split’);   
  33.     % 先获取当前的Epoch_num   
  34.     epoch_num(epoch_idx) = str2double(results{1}(8:(length(results{1})-1)));   
  35.     % 开始解析SDR, SIR, SAR 和 NSDR 值   
  36.     results = results{2};   
  37.     results = regexp(results, ’, ’, ’split’);   
  38.     for idx = 1:length(results)   
  39.         tmp_results = regexp(results{idx}, ’:', ’split’);   
  40.         if strcmp(tmp_results{1}, ’SDR’)   
  41.             sdr_value(epoch_idx) = str2double(tmp_results{2});   
  42.         elseif strcmp(tmp_results{1}, ’SIR’)   
  43.             sir_value(epoch_idx) = str2double(tmp_results{2});   
  44.         elseif strcmp(tmp_results{1}, ’SAR’)   
  45.             sar_value(epoch_idx) = str2double(tmp_results{2});   
  46.         elseif strcmp(tmp_results{1}, ’NSDR’)   
  47.             nsdr_value(epoch_idx) = str2double(tmp_results{2});   
  48.         else  
  49.             print (‘Error eval item’)   
  50.         end   
  51.     end   
  52. end   
  53.   
  54. % 绘图   
  55. figure(‘position’, figureSize)   
  56. plot(epoch_num, sdr_value, [colours(1), symbols(1), char(linetypes(1))], ’linewidth’, lineWidth);   
  57. hold on;   
  58. plot(epoch_num, sir_value, [colours(2), symbols(1), char(linetypes(2))], ’linewidth’, lineWidth);   
  59. hold on;   
  60. plot(epoch_num, sar_value, [colours(3), symbols(1), char(linetypes(3))], ’linewidth’, lineWidth);   
  61. hold on;   
  62. plot(epoch_num, nsdr_value, [colours(4), symbols(1), char(linetypes(4))], ’linewidth’, lineWidth);   
  63. title(‘BSS\_Eval results’)   
  64. grid on;   
  65. xlabel(‘Epoch Number’)   
  66. ylabel(‘Value’)   
  67. axis([0 epoch_num(end) -2 5]);   
  68. legend(para_list, ’Location’, ’SouthEast’)   

绘制的结果如下：

发表评论

电子邮件地址不会被公开。 必填项已用 * 标注

*

您可以使用这些 HTML 标签和属性： <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <del datetime=""> <em> <i> <q cite=""> <strike> <strong>