matlab 解码 NMEA0183格式GGA数据

function [data] = decode_NMEA_GGA（file）
    format long g

    fid = fopen（file‘r‘）;
    % 只读到最后一个M
    rawdata = textscan（fid‘%s%f%f%s%f%s%f%f%f%f%s%f%s%*[^\n]‘‘Delimiter‘‘‘）;

    modeall = rawdata{7};
    idx = modeall ~=0;
    nsol = length（rawdata{4}）;
    
    % 默认为北半球
    latsig = ones（nsol1）;
    for i= 1:nsol
        if cell2mat（strfind（rawdata{14}（i）‘S‘））
            latsig（i） = -1;
        end 
    end 
    % 默认为东半球
    nsol = length（rawdata{6}）;
    lonsig = ones（nsol1）;
    for i= 1:nsol
        if cell2mat（strfind（rawdata{16}（i）‘W‘））
            lonsig（i） = -1;
        end 
    end     

    TT = cell2mat（rawdata（2））;
    hour = floor（TT/10000）;
    minute = floor（（TT-hour*10000）/100）;