https://github.com/perrysou/GNSS_SDR/blob/master/include/ephemeris.m line 95~96: eph.T_GD = twosComp2dec(subframe(197:204)) * 2^(-31); eph.IODC = bin2dec([subframe(83:84) subframe(197:204)]); IODC should be at 211:218 as IS-GPS-200N page 99: 
https://github.com/perrysou/GNSS_SDR/blob/master/include/ephemeris.m

line 95~96:
eph.T_GD = twosComp2dec(subframe(197:204)) * 2^(-31);
eph.IODC = bin2dec([subframe(83:84) subframe(197:204)]);
IODC should be at 211:218 as IS-GPS-200N page 99: