You cannot select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
158 lines
5.8 KiB
Matlab
158 lines
5.8 KiB
Matlab
function collision_physics_data()
|
|
figoff=0;
|
|
|
|
% Expt 1 of Todd & Warren
|
|
clear est
|
|
est(1,:)=[0.7408 0.8911 0.9917 0.9964];
|
|
est(2,:)=[0.7020 0.8450 0.9666 0.9938];
|
|
est(3,:)=[0.6193 0.7802 0.8946 0.9590];
|
|
Ua=[1.19:0.13:4.45];
|
|
mb=1;
|
|
Ma=[1.25,1.5,2,3];
|
|
E=[0.9,0.5,0.1];
|
|
figure(figoff+1),clf
|
|
semilogx(Ma,est(1,:),'m-d','LineWidth',2,'MarkerFaceColor','m','MarkerSize',9),hold on
|
|
semilogx(Ma,est(2,:),'b-s','LineWidth',2,'MarkerFaceColor','b','MarkerSize',9)
|
|
semilogx(Ma,est(3,:),'g-o','LineWidth',2,'MarkerFaceColor','g','MarkerSize',9)
|
|
semilogx([1,3],[0.5,0.5],'k--')
|
|
set(gca,'XTick',[1,1.5,3],'XTickLabel',['1/1';'3/2';'3/1'],'FontSize',12)
|
|
xlabel('m_A/m_B');ylabel('proportion choosing object A')
|
|
set(gcf,'PaperSize',[6 9],'PaperPosition',[0 0.25 6 9],'PaperOrientation','Portrait');
|
|
print(gcf, '-dpdf','collision_physics1_data.pdf');
|
|
|
|
|
|
% Expt 2 of Todd & Warren - moving condition
|
|
clear est
|
|
est(1,:)=[1-0.99,1-0.95,1-0.95,1-0.89,0.90,0.98,0.96,0.98];
|
|
est(2,:)=[1-0.99,1-0.95,1-0.95,1-0.91,0.85,0.94,0.95,0.98];
|
|
est(3,:)=[1-0.94,1-0.99,1-0.98,1-0.95,0.91,0.96,1,0.99];
|
|
ua=3.18;
|
|
ub=-3.18;
|
|
mb=1;
|
|
Ma=[1/3,0.5,2/3,1/1.25,1.25,1.5,2,3];
|
|
E=[0.9,0.5,0.1];
|
|
figure(figoff+2),clf
|
|
semilogx(Ma,est(1,:),'m-d','LineWidth',2,'MarkerFaceColor','m','MarkerSize',9),hold on
|
|
semilogx(Ma,est(2,:),'b-s','LineWidth',2,'MarkerFaceColor','b','MarkerSize',9)
|
|
semilogx(Ma,est(3,:),'g-o','LineWidth',2,'MarkerFaceColor','g','MarkerSize',9)
|
|
semilogx([1,1],[0,1],'k--')
|
|
semilogx([1/3,3],[0.5,0.5],'k--')
|
|
set(gca,'XTick',[1/3,2/3,1,1.5,3],'XTickLabel',['1/3';'2/3';'1/1';'3/2';'3/1'],'FontSize',12)
|
|
xlabel('m_A/m_B');ylabel('proportion choosing object A')
|
|
set(gcf,'PaperSize',[6 9],'PaperPosition',[0 0.25 6 9],'PaperOrientation','Portrait');
|
|
print(gcf, '-dpdf','collision_physics2_data.pdf');
|
|
|
|
|
|
|
|
% Expt 2 of Todd & Warren - stationary condition
|
|
clear est
|
|
est(1,:)=[1-0.9,1-0.75,1-0.48,1-0.38,0.93,0.96,0.95,0.96];
|
|
est(2,:)=[1-0.74,1-0.43,1-0.23,1-0.24,0.95,0.94,0.96,0.96];
|
|
est(3,:)=[1-0.56,1-0.45,1-0.36,1-0.39,0.65,0.76,0.68,0.78];
|
|
ua=3.2;
|
|
ub=0;
|
|
mb=1;
|
|
Ma=[1/3,0.5,2/3,1/1.25,1.25,1.5,2,3];
|
|
E=[0.9,0.5,0.1];
|
|
figure(figoff+3),clf
|
|
semilogx(Ma,est(1,:),'m-d','LineWidth',2,'MarkerFaceColor','m','MarkerSize',9),hold on
|
|
semilogx(Ma,est(2,:),'b-s','LineWidth',2,'MarkerFaceColor','b','MarkerSize',9)
|
|
semilogx(Ma,est(3,:),'g-o','LineWidth',2,'MarkerFaceColor','g','MarkerSize',9)
|
|
semilogx([1,1],[0,1],'k--')
|
|
semilogx([1/3,3],[0.5,0.5],'k--')
|
|
set(gca,'XTick',[1/3,2/3,1,1.5,3],'XTickLabel',['1/3';'2/3';'1/1';'3/2';'3/1'],'FontSize',12)
|
|
xlabel('m_A/m_B');ylabel('proportion choosing object A')
|
|
set(gcf,'PaperSize',[6 9],'PaperPosition',[0 0.25 6 9],'PaperOrientation','Portrait');
|
|
print(gcf, '-dpdf','collision_physics3_data.pdf');
|
|
|
|
|
|
|
|
% experiment shown in Fig.4 of Runesun and Vedeler
|
|
clear est
|
|
est(1,:)=1-[0.0016 0.0109 0 0.0328 0.0733 0.5460 0.7488 0.8440];
|
|
est(2,:)=1-[0.0016 0 0.0016 0.0031 0.0187 0.2543 0.5086 0.5211];
|
|
Mb=[0.25,0.37,0.55,0.82,1.22,1.81,2.69,4];
|
|
Ma=1;
|
|
ma=1;
|
|
e=0.9;
|
|
figure(figoff+4),clf
|
|
semilogx(1./Mb,est(1,:),'m-d','LineWidth',2,'MarkerFaceColor','m','MarkerSize',9),hold on
|
|
semilogx(1./Mb,est(2,:),'m--d','LineWidth',2,'MarkerFaceColor','m','MarkerSize',9)
|
|
semilogx([1,1],[0,1],'k--')
|
|
semilogx([1/4,4],[0.5,0.5],'k--')
|
|
set(gca,'XTick',[1/4,1/2,1,2,4],'XTickLabel',['1/4';'1/2';'1/1';'2/1';'4/1'],'FontSize',12)
|
|
xlabel('m_A/m_B');ylabel('proportion choosing object A')
|
|
set(gcf,'PaperSize',[6 9],'PaperPosition',[0 0.25 6 9],'PaperOrientation','Portrait');
|
|
print(gcf, '-dpdf','collision_physics4_data.pdf');
|
|
|
|
|
|
|
|
% experiment shown in Fig.3a of Sanborn
|
|
clear est
|
|
est(1,:)=[0.1585 0.5916 0.8265 0.8720];
|
|
est(2,:)=[0.3012 0.5492 0.6857 0.8122];
|
|
Ua=[1.19:0.13:4.45];
|
|
ub=0;
|
|
mb=1;
|
|
Ma=[1/3,1/1.25,1.25,3];
|
|
E=[1,0.1];
|
|
figure(figoff+5),clf
|
|
semilogx(Ma,est(1,:),'r-+','LineWidth',2,'MarkerFaceColor','r','MarkerSize',9),hold on
|
|
semilogx(Ma,est(2,:),'g-o','LineWidth',2,'MarkerFaceColor','g','MarkerSize',9)
|
|
semilogx([1,1],[0,1],'k--')
|
|
semilogx([1/3,3],[0.5,0.5],'k--')
|
|
set(gca,'XTick',[1/3,2/3,1,1.5,3],'XTickLabel',['1/3';'2/3';'1/1';'3/2';'3/1'],'FontSize',12)
|
|
xlabel('m_A/m_B');ylabel('proportion choosing object A')
|
|
set(gcf,'PaperSize',[6 9],'PaperPosition',[0 0.25 6 9],'PaperOrientation','Portrait');
|
|
print(gcf, '-dpdf','collision_physics5_data.pdf');
|
|
|
|
|
|
|
|
|
|
% experiment shown in Fig.3b of Sanborn (as above but with object B occluded after collision)
|
|
clear est
|
|
est(1,:)=[0.1769 0.1692 0.2788 0.7000];
|
|
est(2,:)=[0.3519 0.5442 0.6212 0.8481];
|
|
Ua=[1.19:0.13:4.45];
|
|
ub=0;
|
|
mb=1;
|
|
Ma=[1/3,1/1.25,1.25,3];
|
|
E=[1,0.1];
|
|
figure(figoff+6),clf
|
|
semilogx(Ma,est(1,:),'r-+','LineWidth',2,'MarkerFaceColor','r','MarkerSize',9),hold on
|
|
semilogx(Ma,est(2,:),'g-o','LineWidth',2,'MarkerFaceColor','g','MarkerSize',9)
|
|
semilogx([1,1],[0,1],'k--')
|
|
semilogx([1/3,3],[0.5,0.5],'k--')
|
|
set(gca,'XTick',[1/3,2/3,1,1.5,3],'XTickLabel',['1/3';'2/3';'1/1';'3/2';'3/1'],'FontSize',12)
|
|
xlabel('m_A/m_B');ylabel('proportion choosing object A')
|
|
set(gcf,'PaperSize',[6 9],'PaperPosition',[0 0.25 6 9],'PaperOrientation','Portrait');
|
|
print(gcf, '-dpdf','collision_physics6_data.pdf');
|
|
|
|
|
|
|
|
|
|
% experiment shown in Fig.3c of Sanborn (as above but with object A occluded after collision)
|
|
clear est
|
|
est(1,:)=[0.6442 0.8615 0.9192 0.9212];
|
|
est(2,:)=[0.3404 0.5904 0.7288 0.8462];
|
|
Ua=[1.19:0.13:4.45];
|
|
ub=0;
|
|
mb=1;
|
|
Ma=[1/3,1/1.25,1.25,3];
|
|
E=[1,0.1];
|
|
figure(figoff+7),clf
|
|
semilogx(Ma,est(1,:),'r-+','LineWidth',2,'MarkerFaceColor','r','MarkerSize',9),hold on
|
|
semilogx(Ma,est(2,:),'g-o','LineWidth',2,'MarkerFaceColor','g','MarkerSize',9)
|
|
semilogx([1,1],[0,1],'k--')
|
|
semilogx([1/3,3],[0.5,0.5],'k--')
|
|
set(gca,'XTick',[1/3,2/3,1,1.5,3],'XTickLabel',['1/3';'2/3';'1/1';'3/2';'3/1'],'FontSize',12)
|
|
xlabel('m_A/m_B');ylabel('proportion choosing object A')
|
|
set(gcf,'PaperSize',[6 9],'PaperPosition',[0 0.25 6 9],'PaperOrientation','Portrait');
|
|
print(gcf, '-dpdf','collision_physics7_data.pdf');
|
|
|
|
|
|
|
|
|
|
|
|
|