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

piday123

CCS3.3是一款DSP的集成开发环境（IDE）。在做DSP开发时，为验证算法，常常需要使用matlab进行算法验证，验证算法就需要数据。因此，一种交互的方法是：
. f, S. I2 X+ d: f& q
使用DSP开发板连接CCS
# t1 l* ]' D9 {- G8 \: K用CCS3.3菜单中的“File->Data->Save..”功能导出DSP内存中的数据到PC的.dat文件中
% P" `# w- o4 W1 v使用Matlab读取.dat文件，进行数据处理
下面给出使用Matlab读取CCS导出的dat文件的方法：
& W8 n  C5 P1 k1 \4 o$ W0 H- q7 B
* t" H' J! b# k$ h# `! z( I) u% =========================================================================
+ M+ Z: z$ o' u- |9 P% ^% Read *dat file of CCS3.3
% xiahouzuoxin
/ X. ?* n0 Z- W2 G% 2014.04.21
5 `3 Z0 g: H3 J: A0 i% =========================================================================
' P9 ]; N# S1 }( m  a
% clc;clear all;close all;
" r) @0 U" a/ w6 k  B
5 T( Y- E' W' J2 t% 对话框选择*.dat文件
[fname,pname]=uigetfile(...
    {'*.dat';'*.*'},'Input *.dat File');

fid = fopen(fullfile(pname,fname));
6 S  z: W1 f8 j( ~fseek(fid, 21,-1);  % 去文件头, 21字节，如 1651 1 80000006 0 100

# l) R* N7 J( L8 E8 J% 按指定格式读取解析
fm = 4;
: J. B: N6 o! C7 U( w7 lswitch (fm)
    case 4  % 按4Byte格式读，如 0x 80000000
        x = textscan(fid, '%2s %8s');  
7 q$ @/ Z7 b- ?, v- o        z(:,1) = hex2dec(x{2});        
    case 2 % 按2Byte格式读，如 0x 8000 0000
4 c" N% c" `! Y3 V        x = textscan(fid, '%2s %4s %4s');  
        z(:,1) = hex2dec(x{3});
: Y; W% P  }$ b, L5 x8 I: V3 X+ r0 ^        z(:,2) = hex2dec(x{2});
. L* \$ u9 `0 J. B9 S- [0 B    case 1 % 按1Byte格式读，如 0x 80 00 00 00
        x = textscan(fid, '%2s %2s %2s %2s %2s');  
# Y0 g4 Y2 n; C: g, V- R" U        z(:,1) = hex2dec(x{5});
        z(:,2) = hex2dec(x{4});
% p7 u" Q1 p" X: u        z(:,3) = hex2dec(x{3});
        z(:,4) = hex2dec(x{2});
    otherwise
  `# \) }) W5 \( y8 R        z = [];
end

8 r) b7 a% s9 j& G  R; Nif ~isempty(z)
0 U. t4 w! \0 L- ^& t8 a    % 将数据处理代码放在这里
7 y; }7 y, V* `  _# y2 f* w
end

fclose(fid);
程序中默认DSP内存中数据的存储格式为整数类型，而对于浮点型DSP如TMS320C6713B，常常在内存中是float或double格式存储的，



这时，我们还需要一个转化为浮点数的操作，因此我写了一个转化为float类型的函数，
- L7 v; h3 B4 @' G* g' O( K
function y = integer2float(x, ishex)
$ J$ i" I/ b) g: g' C% =========================================================================
% 有些数据y原本是float类型，在内存中是按float格式存储
% 而现在按整数将其从内存中读出变成了x
; Y5 s8 P4 c  ]: X% B2 f: R8 Y, P% 这个函数就是用于将按整数格式读出的x转变成y
+ ^( g. o0 U5 m: B% ishex~=0表示输入为十六进制格式
% 注：要完成该函数，你必须了解IEEE浮点数的格式
%
& |7 D& D1 l$ n/ d/ A% xiahouzuoxin
3 a+ L0 a7 R) r5 w% 2014.04.21
% =========================================================================
if nargin == 2
4 e! d: v* W) U2 k$ |$ C0 D    if (ishex)
$ v" B0 j4 E, o        x = hex2dec(x);
    end
. `8 b2 W% b# x1 ^7 q0 D% k, _- nend

1 h- Q. \7 p9 J3 k4 E1 B0 D7 l[h w] = size(x);
y = zeros(h,w);
3 ~' S/ y5 f- L: n8 {
for i = 1:h
2 Q( U- z8 b8 i' w* n* `    for j = 1:w
        sign = bitget(x(i,j),32);
        exponent = bitget(x(i,j),24:31) * 2.^(0:7).';
        fraction = bitget(x(i,j),1:23) * 2.^(-23:-1).';
: J% B& l+ P& `# ]. M        y(i,j) = (-1)^sign * (1+fraction) * 2^(exponent-127);
* E! P# Y  a# ]/ G2 l- h: _    end
end
根据不同的输入情况（可以为十六进制的字符串表示或十进制数，但必须≤8Bytes），使用方法举个例子：
# I2 Y. Y1 f* u& c& x
. z  i" I, h: j$ G$ e>> x='4565A012'
0 j* B' F* ~- I% g>> y=integer2float(x,1);
8 N' a+ ~8 z2 P& ^, j# w
>> x = 23974881923;
0 w# k7 s% {' R& A>> y = integer2float(x)
后来发现，Matlab已经早有对应的方法了，

y = typecast(uint32(z),'single');
- _' g! l; [7 D0 j- i6 \8 F可以直接把z转化为解析为单精度浮点数，具体使用请参见
( f! O5 f2 a* ?& ?7 k8 \
>> help typecast



" N  l" J- u! B6 g8 U7 t

Heaven_1

typecast这个函数用处比较多