hello
this would be my suggestion - others methods I tried have not been successfull
flowrate = [ 0 93 100 120 140 160 172 180 200 220 232 ] ;
head = [ 9.51 8.76 8.68 8.41 8.01 7.45 7.03 6.70 5.73 4.49 3.61 ] ;
cv20 = [ 0.8 1.1 2.2 4.1 6.6 9.9 13.9 18.6 24.1 30.3 37.2] ;
cv30 = [ 0.8 1 1.5 2.4 3.7 5.3 7.3 9.6 12.3 15.3 18.7 ] ;
cv40 = [ 0.8 0.9 1.3 1.8 2.7 3.7 5 6.4 8.2 10.1 12.3 ] ;
cv50 = [ 0.8 0.9 1.1 1.6 2.2 2.9 3.9 5 6.3 7.7 9.3 ] ;
cv60 = [ 0.8 0.8 1.1 1.4 1.9 2.5 3.3 4.2 5.2 6.4 7.7 ] ;
cv70 = [ 0.8 0.8 1 1.3 1.8 2.3 2.9 3.7 4.6 5.6 6.7] ;
cv80 = [ 0.8 0.8 1 1.3 1.7 2.1 2.7 3.4 4.2 5.1 6.1 ] ;
% plot
figure(1)
plot(flowrate,head,flowrate,cv20,flowrate,cv30,...
flowrate,cv40,flowrate,cv50,flowrate,cv60,...
flowrate,cv70,flowrate,cv80);
title('Pump Performance Curves Estimates')
ylim([0.0 10.0])
% plot
flowrate2 = linspace(min(flowrate),max(flowrate),100);
[cv20new] = my_opti_fminsearch(flowrate,cv20,flowrate2);
[cv20new] = my_opti_fminsearch(flowrate,cv20,flowrate2);
[cv30new] = my_opti_fminsearch(flowrate,cv30,flowrate2);
[cv40new] = my_opti_fminsearch(flowrate,cv40,flowrate2);
[cv50new] = my_opti_fminsearch(flowrate,cv50,flowrate2);
[cv60new] = my_opti_fminsearch(flowrate,cv60,flowrate2);
[cv70new] = my_opti_fminsearch(flowrate,cv70,flowrate2);
[cv80new] = my_opti_fminsearch(flowrate,cv80,flowrate2);
figure(2)
plot(flowrate,head,flowrate2,cv20new,flowrate2,cv30new,...
flowrate2,cv40new,flowrate2,cv50new,flowrate2,cv60new,...
flowrate2,cv70new,flowrate2,cv80new);
title('Pump Performance Curves Estimates')
ylim([0.0 10.0])
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [cvnew] = my_opti_fminsearch(flowrate,cv,flowrate2)
cv = interp1(flowrate,cv,flowrate2);
% fminsearch method
x = flowrate2;
a = cv(1);
f = @(b,x) a + b.*(x.^2);
b_init = (cv(end) - a)./(flowrate(end).^2);
obj_fun = @(params) norm(f(params(1),flowrate2)-cv);
sol = fminsearch(obj_fun, [b_init]);
b_sol = sol(1);
cvnew = f(b_sol, flowrate2);
end
