导出CCS3.3数据及使用matlab处理的方法

piday123

CCS3.3是一款DSP的集成开发环境（IDE）。在做DSP开发时，为验证算法，常常需要使用matlab进行算法验证，验证算法就需要数据。因此，一种交互的方法是：

使用DSP开发板连接CCS
- l9 {8 X- ?8 d/ f用CCS3.3菜单中的“File->Data->Save..”功能导出DSP内存中的数据到PC的.dat文件中
" S( F! `5 H+ }! V: i使用Matlab读取.dat文件，进行数据处理
; b1 o, b2 ?! L! T, W& H$ X下面给出使用Matlab读取CCS导出的dat文件的方法：
7 D7 U, @. z( L! M
; t2 B( K9 k: U7 l% =========================================================================
% Read *dat file of CCS3.3
! H. r8 R0 S) G) p% xiahouzuoxin
  y$ `& e* @& ^, f% 2014.04.21
/ U5 b. n8 Z1 J4 o! ~+ f& |% =========================================================================
2 i) R) f7 f: z# x
0 U5 Y8 u5 n( P, d* W& X% clc;clear all;close all;

$ w) g: y. R  x5 [% 对话框选择*.dat文件
) T( }5 e" x6 g' T: G[fname,pname]=uigetfile(...
4 J% J5 \& l) B% u# {    {'*.dat';'*.*'},'Input *.dat File');

  s$ |, C+ Y/ S. v; @fid = fopen(fullfile(pname,fname));
fseek(fid, 21,-1);  % 去文件头, 21字节，如 1651 1 80000006 0 100
% T8 g* C3 D. E1 K# U
% 按指定格式读取解析
fm = 4;
switch (fm)
% j) ~# B1 H: z9 ^) H. b: G    case 4  % 按4Byte格式读，如 0x 80000000
        x = textscan(fid, '%2s %8s');  
) V4 m+ n) H) E- e4 q5 ~% Y        z(:,1) = hex2dec(x{2});        
    case 2 % 按2Byte格式读，如 0x 8000 0000
2 b8 W2 ^; e, l/ u        x = textscan(fid, '%2s %4s %4s');  
        z(:,1) = hex2dec(x{3});
        z(:,2) = hex2dec(x{2});
    case 1 % 按1Byte格式读，如 0x 80 00 00 00
2 Q5 r  p8 J; V9 z: b        x = textscan(fid, '%2s %2s %2s %2s %2s');  
        z(:,1) = hex2dec(x{5});
        z(:,2) = hex2dec(x{4});
, r8 d1 ^- z6 ?  c        z(:,3) = hex2dec(x{3});
5 J3 E9 Z( s3 m" w4 w        z(:,4) = hex2dec(x{2});
    otherwise
        z = [];
end
0 |. T9 D9 q5 G. h  f
0 K$ m) A6 }4 o- ]if ~isempty(z)
9 u# y) S5 w7 j7 o! t$ s# I/ g) a    % 将数据处理代码放在这里
. k6 @; ^: E& G% ]9 l3 Y. a
end

' L& n; }, ^+ C# e6 R. g  {6 gfclose(fid);
程序中默认DSP内存中数据的存储格式为整数类型，而对于浮点型DSP如TMS320C6713B，常常在内存中是float或double格式存储的，
5 z3 F; Z1 \2 O( I8 f6 v
3 X! U+ [) O' y# C' c/ z/ |

8 Q1 Y7 z. p; z: }这时，我们还需要一个转化为浮点数的操作，因此我写了一个转化为float类型的函数，
5 R$ ]/ H2 y1 d4 L1 n: A' P
! N1 d0 E4 l: Y) w' `function y = integer2float(x, ishex)
% =========================================================================
% ^9 G! P$ h& y; k* `% 有些数据y原本是float类型，在内存中是按float格式存储
% 而现在按整数将其从内存中读出变成了x
% 这个函数就是用于将按整数格式读出的x转变成y
% ishex~=0表示输入为十六进制格式
% 注：要完成该函数，你必须了解IEEE浮点数的格式
%
% xiahouzuoxin
& K# _, p* x' W, O% 2014.04.21
% =========================================================================
7 n/ Z  _/ k( k) V( Bif nargin == 2
9 P' g9 R7 v( v" ]# E5 _# D( H    if (ishex)
        x = hex2dec(x);
; X3 R2 l/ g5 S, s    end
end
+ J) p* T5 F8 v" N2 Y% w
. r4 B2 h' q9 |( E[h w] = size(x);
, L+ F/ l- |, B5 A% t2 Ry = zeros(h,w);
& w( x/ @. l4 O6 V
for i = 1:h
    for j = 1:w
) L9 {7 C, D/ [, v8 L1 T        sign = bitget(x(i,j),32);
; g- G- l( d3 v# O# N' Q        exponent = bitget(x(i,j),24:31) * 2.^(0:7).';
        fraction = bitget(x(i,j),1:23) * 2.^(-23:-1).';
        y(i,j) = (-1)^sign * (1+fraction) * 2^(exponent-127);
    end
end
根据不同的输入情况（可以为十六进制的字符串表示或十进制数，但必须≤8Bytes），使用方法举个例子：

2 ]( T, Z; m  |' V" P4 x>> x='4565A012'
4 t1 c) Q  {* B1 @1 O>> y=integer2float(x,1);

& t+ f0 ?! R/ l. ~; I: w>> x = 23974881923;
>> y = integer2float(x)
: V/ `* [1 S" ?  B" [! h7 H后来发现，Matlab已经早有对应的方法了，
% D9 {1 R- t( @* e; k: W6 b
6 Y" b$ o1 r( `1 ^! jy = typecast(uint32(z),'single');
  Y9 g+ @7 d9 B0 n" m可以直接把z转化为解析为单精度浮点数，具体使用请参见

# O+ V& c6 w5 ~+ r) e) L>> help typecast

/ W; o3 i2 b& Y# [' ?" I* Z

$ ], K( `1 w. J  }! h

Heaven_1

typecast这个函数用处比较多