AblationRateCurveFit.m

Data = [40	4.4	200	6.42;
        42	1.8	200	3.73;
        45	4.9	240	9.45;
        65	1.6	200	6.955;
        73	8.9	120	7.38;
        73	13.3 120 8.61;
        107	2.3	55 2.3;
        114	20.9 80 8.685;
        133	31.6 80 10.07;
        143	3.8	200	12.66;
        143	3.8	400	24.72;
        154	13.3 33	2.925;
        154	13.3 66	5.61;
        165	6.8	50 3.96;
        169	5 33 2.235;
        169	5 66 4.28;
        175	13.9 50 5.22;
        183	26.7 50 6.765];
% Heat Flux(X);  Stag. Pressure(Y);  Time;  Ablation
% W/cm^2;             kPa;            s;       mm

% data from NTRS 20120017464

abRate = Data(:,4)./Data(:,3);

figure(1)
options = fitoptions('poly11');
options.Lower = [0 0.0002 0.002];
options.Upper = [0 0.0004 0.004];
f = fit([Data(:,1), Data(:,2)], abRate,'poly11',options);

plot(f,[Data(:,1), Data(:,2)], abRate);
xlabel('Heat Flux (W/cm^2)');
ylabel('Stagnation Pressure (kPa)');
zlabel('Ablation Rate (mm/s)');
ylim([0 35]);
xlim([0 200]);