isingfidatat.m

matlat程序数据拟合例子 要自己再输入数据文件的

clc
clear
%%效果不好，没有用了，直接用原始数据也可以了。
Fid1=fopen('heatcapacity8.dat', 'r');
h8=fscanf(Fid1, '%e  %f',[2,inf]);
fclose(Fid1);

Fid1=fopen('heatcapacity16.dat', 'r');
h16=fscanf(Fid1, '%e  %f',[2,inf]);
fclose(Fid1);

Fid1=fopen('heatcapacity24.dat', 'r');
h24=fscanf(Fid1, '%e  %f',[2,inf]);
fclose(Fid1);

Fid1=fopen('heatcapacity48.dat', 'r');
h48=fscanf(Fid1, '%e  %f',[2,inf]);
fclose(Fid1);

h8=h8';
h16=h16';
h24=h24';
h48=h48';

p8=polyfit(h8(:,1),h8(:,2),12);
p16=polyfit(h16(:,1),h16(:,2),12);
p24=polyfit(h24(:,1),h24(:,2),12);
p48=polyfit(h48(:,1),h48(:,2),10);


%get the polyval
y8=polyval(p8,h8(:,1));
y16=polyval(p16,h16(:,1));
y24=polyval(p24,h24(:,1));
y48=polyval(p48,h48(:,1));

yy=[y8,y16,y24,y48];
B=[h8(:,1),yy]

fid8=fopen('3dheatfit.dat', 'w');
fprintf(fid8,'%5.3f %8.5f %8.5f %8.5f %8.5f \n',B');
fclose(fid8);

figure
plot(h8(:,1),yy(:,1),'b',h8(:,1),yy(:,2),'g',h8(:,1),yy(:,3),'r',...
    h8(:,1),yy(:,4),'c')  
legend('8x8x8','16x16x16','24x24x24','48x48x48')
axis([2,6.4, 0,3.5])