Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- clear
- close all
- clc
- %%
- labels{1} = 'pristine';
- labels{2} = 'gold-doped';
- labels{3} = 'platinum-doped';
- data{1} = load('pris.txt');
- data{2} = load('Au.txt');
- data{3} = load('Pt.txt');
- faint = 0.25;
- clr_f{1} = [0.6350 0.0780 0.1840 faint];
- clr_f{2} = [0.4660 0.6740 0.1880 faint];
- clr_f{3} = [0 0.4470 0.7410 faint];
- clr_f{4} = [0.4940 0.1840 0.5560 faint];
- % data_smooth = cell(1,4);
- % for i = 1:length(data)
- % x = data{i}(:,1);
- % y = data{i}(:,2);
- % yy = smooth(y,'rlowess')
- % data_smooth{i} = yy;
- % end
- clr{1} = [0.6350 0.0780 0.1840];
- clr{2} = [0.4660 0.6740 0.1880];
- clr{3} = [0 0.4470 0.7410];
- clr{4} = [0.4940 0.1840 0.5560];
- figure
- val_range{1} = [100 240];
- val_range{2} = [1250 1400];
- val_range{3} = [1540 1640];
- val_range{4} = [2500 2750];
- span = [0.05 0.5 0.25 0.4];
- ylims{1} = [40 140];
- ylims{2} = [40 90];
- ylims{3} = [0 400];
- ylims{4} = [40 100];
- for i = 1:length(val_range)
- A = find(data{1}(:,1) >= val_range{i}(1), 1);
- B = find(data{1}(:,1) >= val_range{i}(2), 1);
- ind_range{i} = [A:B];
- end
- for i = 1:length(data)
- subplot(2,length(val_range),[1:length(val_range)])
- hold on
- grid on
- plot(data{i}(:,1),data{i}(:,2),'-','color',clr{i})
- xlabel('wavenumber, cm^{-1}')
- ylabel('intensity, a.u.')
- legend(labels)
- ftype()
- ylim([0 400])
- for j = 1:length(val_range)
- subplot(2,length(val_range),length(val_range)+j)
- [x1,y1,yy1] = smoother(data{i},ind_range{j},span(j));
- hold on
- grid on
- plot(x1,y1,'color',clr_f{i})
- plot(x1,yy1,'color',clr{i},'Linewidth',2)
- xlim(val_range{j})
- xlabel('wavenumber, cm^{-1}')
- ylabel('intensity, a.u.')
- ylim(ylims{j})
- ftype()
- end
- % for j = 1:length(val_range)
- % subplot(2,length(val_range),length(val_range)+j)
- % hold on
- % grid on
- % [x1,y1,yy1] = fitter(data{i},ind_range{j},fit_types(j));
- % xlim(val_range{j})
- % xlabel('wavenumber, cm^{-1}')
- % ylabel('intensity, a.u.')
- % ftype()
- % end
- % subplot(2,4,5)
- % [x1,y1,yy1] = smoother(data{i},ind_range{1});
- % hold on
- % grid on
- % plot(x1,y1,'color',clr_f{i})
- % plot(x1,yy1,'color',clr{i})
- % xlim(val_range{1})
- % xlabel('wavenumber, cm^{-1}')
- % ylabel('intensity, a.u.')
- %
- % subplot(2,4,6)
- % [x,y,yy] = smoother(data{i},ind_range{2});
- % hold on
- % grid on
- % plot(x,y,'color',clr_f{i})
- % plot(x,yy,'color',clr{i})
- % xlim(val_range{2})
- % xlabel('wavenumber, cm^{-1}')
- % ylabel('intensity, a.u.')
- %
- % subplot(2,4,7)
- % [x,y,yy] = smoother(data{i},ind_range{3});
- % hold on
- % grid on
- % plot(x,y,'color',clr_f{i})
- % plot(x,yy,'color',clr{i})
- % xlim(val_range{3})
- % xlabel('wavenumber, cm^{-1}')
- % ylabel('intensity, a.u.')
- %
- % subplot(2,4,8)
- % [x,y,yy] = smoother(data{i},ind_range{4});
- % hold on
- % grid on
- % plot(x,y,'color',clr_f{i})
- % plot(x,yy,'color',clr{i})
- % xlim(val_range{4})
- % xlabel('wavenumber, cm^{-1}')
- % ylabel('intensity, a.u.')
- end
- %%
- function [x,y,yy] = smoother(data,range,span);
- x = data(range,1);
- y = data(range,2);
- yy = smooth(y,span,'loess');
- end
- %{
- function [x,y,yy] = fitter(data,range,degree)
- faint = 0.25;
- clr_f{1} = [0.6350 0.0780 0.1840 faint];
- clr_f{2} = [0.4660 0.6740 0.1880 faint];
- clr_f{3} = [0 0.4470 0.7410 faint];
- clr_f{4} = [0.4940 0.1840 0.5560 faint];
- clr{1} = [0.6350 0.0780 0.1840];
- clr{2} = [0.4660 0.6740 0.1880];
- clr{3} = [0 0.4470 0.7410];
- clr{4} = [0.4940 0.1840 0.5560];
- x = data(range,1);
- y = data(range,2);
- F = fit(x,y,degree);
- plot(x,y,'color',clr_f{i})
- plot(F,x,y,'color',clr{i})
- end
- %}
- function ftype()
- set(gca,'FontSize',13,'fontWeight','bold')
- set(findall(gcf,'type','text'),'FontSize',13,'fontWeight','bold')
- end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement
