matlab生成HEX文件-任意信號 大于64K長度

HEX文件格式不贅述，寫里直接放上代碼。請批評改正。

%%convert a signal data into hex file format
% data format:16bit 
% signal length: less than 2^24-1
% author: Yang Li yangli0534@gmail.com
% data:2015.01.27

clear all;
close all;
clc;

%% fixed point or binary16 float point
fixed=1;
%% generate a signal waveform
%
T=10;%time length
fs=1e3;%sample rate
N=10000;
Period=N/fs/10;
f=1/Period;
t=linspace(0,(N-1)/fs,N);

%data=1:2^5-1; 
data= sin(2*pi*f*t) +(rand(1,N)-0.5)/10;
% data=mod(50000*t,1000);
% plot(t,data)
dataBits=16;
byteBits=8;

%% convert into 16 bits fixed point number
if fixed == 1
    maximum=max(abs(data));%maximum of absolute signal 
    scale=(2^15-1)/maximum;%scale 
    data=round(data*scale);
    for i=1:1:length(data)
        if data(i) <0 
            data(i) =2^16-abs(data(i)) ;
        end
    end
end
%% write
fh=fopen('data.hex','w');%open a new file (file path, file name)
bytesEveryRow=16;%write 16 bytes every row
dataEveryRow=bytesEveryRow/(dataBits/byteBits);%number of data every row
lengthData=length(data);%length of signal data
rowLengths=ceil(lengthData/dataEveryRow);%number of rows of the file to be written

for i=1:rowLengths % write a row each time
    if(mod(i,hex2dec('1000'))==1)
        baseaddr=dec2hex(floor(i/hex2dec('1000')),4);
        checkSum=6+hex2dec(baseaddr(1:2))+hex2dec(baseaddr(3:4));
         checkSum=(dec2hex(bitand(256-mod(checkSum,256),255),2));
        fprintf(fh,[':02000004',baseaddr,checkSum,13,10]);%產生HEX碼(end-1:end)
    
    end
    checkSum=0;
    addr=dec2hex(mod((i-1)*16,hex2dec('10000')),4);%the start address of this row
    checkSum=hex2dec(addr(1:2))+hex2dec(addr(3:4));
    if i*dataEveryRow <lengthData %  numbers of data to be written greater than dataEveryRow
        dataRow=dataEveryRow;
        dataHex='';
        for j=8*(i-1)+1:1:8*(i)            
            if fixed == 1%                
                dataHexCurr=dec2hex(data(j),4);%16 bits fixed point
             else
                dataHexCurr=num2binary16hex(data(j));%number format: 16bits floatIEEE 754 half precison float point standard,result is a hex format
             end
            checkSum=checkSum+hex2dec(dataHexCurr(1:2))+hex2dec(dataHexCurr(3:4));
            dataHex=strcat(dataHex,dataHexCurr);
        end    
            
    else
         dataHex='';
         dataRow =lengthData-(rowLengths-1)*dataEveryRow;
        for j=(rowLengths-1)*dataEveryRow+1:1:lengthData            
             if fixed == 1% 
                dataHexCurr=dec2hex(data(j),4);%16 bits fixed point
             else
                dataHexCurr=num2binary16hex(data(j));%number format: 16bits floatIEEE 754 half precison float point standard,result is a hex format
             end  
            checkSum=checkSum+hex2dec(dataHexCurr(1:2))+hex2dec(dataHexCurr(3:4));
            dataHex=strcat(dataHex,dataHexCurr);
        end
    end
   checkSum=checkSum+dataRow*dataBits/byteBits;
    checkSum=(dec2hex(bitand(256-mod(checkSum,256),255),2)); 
%     disp([num2str(i),'--',dataHex,'--',checkSum])
    fprintf(fh,[':',dec2hex(dataRow*dataBits/byteBits,2),addr,'00',dataHex,checkSum,13,10]);%產生HEX碼(end-1:end)
end

fprintf(fh,':00000001FF');    
fclose(fh);

代碼中有定點浮點兩種，其中浮點格式為IEEE754標準中的half-precision-floating point, 16bit, 1bit 浮號位，5bit指數位，10小數位，代碼如下：

function [ hex ] = num2binary16he( number )
%The IEEE 754 standard specifies a binary16 as having the following format:

%Sign bit: 1 bit
%Exponent width: 5 bits
%Significand precision: 11 bits (10 explicitly stored)
%S EEEEE MMMMMMMMMM
%INPUT: number is a random precision number

%OUTPUT : var binary is a 2's string of the binary16
%         var decimail is the value of the number in decimal

%Author :Yang Li .
%Email: yangli0534@gmail.com
if real(number) >= 2^14  
    error('beyond the maximum :-2^14--+2^24');
end

S=number<0 ;
E=0; 
M=abs(number);
while((E>-15 && E<15)&&( M>=2 || M<1 ))
    if(M>=2)
     E=E+1;
     M=M/2;
    else
        E=E-1;
        M=M*2;
    end    
end
if ((E==-15 || E==15)&& M>=2 || M<1)
    M=0;
else
    M=round((M-1)*2^10) ; 
end


  number1=(-1)^S*(2^ E   )*(1+M/2^10);
% binary=strcat(num2str(S),num2str(E),num2str(M)) 
E =E+15;
binary=strcat(num2str(S),dec2bin(E,5), dec2bin(M,10)) ;


sReal=(str2double(binary(1)));
eReal =-15; 
 for i=2:6
    eReal=eReal+str2double(binary(i))*2^(6-i);
 end    
%  eReal=eReal*(-1)^str2double(binary(2));
 mReal = 0; 
for i=7:16
    mReal=mReal+str2double(binary(i))*2^(16-i);
end 
 
numberReal=(-1)^sReal*2^eReal*(1+mReal/2^10);
err=num2str(abs(number-numberReal));
decimal=numberReal;
% disp(['the origin data is : ',num2str(number)]);
% disp(['the float format is : ',binary]); 
% disp(['so ,there is a error :',num2str(abs(number-numberReal))]);

num=0;
for i=1:16
    num=num+str2double(binary(i))*2^(16-i);
end
hex=dec2hex(num,4);
end

?