% plot the data from the Innolight Temperature sweep % Innolight temperature InnTemp = [0.60 .59 .56 .52 .65] + 39; FreqDiff = [412 463 564 737 306]; % fit the data p = polyfit(InnTemp, FreqDiff, 1); p1 = p(1)/1000; % plot the data close all plot(InnTemp, FreqDiff, 'o', InnTemp, p(1)*InnTemp + p(2), '-') xlabel('Innolight temperature (deg C)') ylabel('Frequency Difference (MHz)') legend('Data', 'Linear fit') title('Innolight - Lightwave frequency difference versus Innolight temperature. LW temp = 37.2602 C') text(39.62, 600, ['Slope = ', num2str(p1, '%5.3f'), 'GHz/K'],... 'HorizontalAlignment','center',... 'BackgroundColor',[.7 .9 .7]) grid on orient landscape print('-dpdf', 'Innolight_temp_sweep.pdf')