[gcf5,gcf6,gcf7,gcf8,gcf9]=plot_ampli(ntotal_colum,mm_index,fop_out,ampli_out,phase_out);
%figure(10),figure(11),figure(12),figure(13)是拟合数据与原始数据的比较图
t=150:dt:(n_row_origin-n_cut_e)*dt;
[gcf10,gcf11,gcf12,gcf13,y]=plot_compare(x,dt,t,ntotal_row,mm_index,ntotal_colum,fop_out,ampli_out,phase_out);
y;
save fitdata_im2053_2.dat y -Ascii%将拟合后数据存起来
%谱图
saveas(gcf1,'pu1','emf')
saveas(gcf2,'pu2','emf')
saveas(gcf3,'pu3','emf')
saveas(gcf4,'pu4','emf')
%恢复幅值
saveas(gcf5,'am1','emf')
saveas(gcf6,'am2','emf')
saveas(gcf7,'am3','emf')
saveas(gcf8,'am4','emf')
saveas(gcf9,'am_all','emf')
%原始数据与拟合数据比较
saveas(gcf10,'compare_1','emf')
saveas(gcf11,'compare_2','emf')
saveas(gcf12,'compare_3','emf')
saveas(gcf13,'compare_4','emf')
fclose('all') |