% Исходные данные
R_c=145;
R_v=145;
U_c=220 ;
w=4000;
Potok_nom=0.007;
F=4000*[0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8];
Potok= Potok_nom*[0 0.3 0.52 0.67 0.78 0.86 0.92 0.96 1.01 1.02] ;
%Нормирование
F_nom=U_c*w/R_v;
F_norm= F/F_nom;
Potok_norm=Potok/Potok_nom;
% Поиск коэффициентов полинома 5 степени
G5=[Potok_norm'  Potok_norm'.^3 Potok_norm'.^5 ];
A5=G5'*G5;
B5=G5'*F_norm';
C5=inv(A5)*B5;
p_norm5= [C5(3), 0, C5(2), 0, C5(1), 0];
p5= F_nom*p_norm5% получаем полином p(Ф) 5 степени
% Поиск коэффициентов полинома 3 степени
G3=[Potok_norm'  Potok_norm'.^3 ];
A3=G3'*G3;
B3=G3'*F_norm';
C3=inv(A3)*B3;
p_norm3= [C3(2), 0, C3(1), 0];
p3= F_nom*p_norm3% получаем полином p(Ф) 3 степени
% Ищем значения  мдс из полинома
x=0:0.001:1.02;
F_raschet5=polyval(p5,x);
% Вычислим значение невязки для полинома 5 степени
j=1;
I5=0;
while  j<11
    I5=I5+(F(j)-polyval(p5,Potok_norm(j)))^2;
    j=j+1;
end;
F_raschet3=polyval(p3,x);
% Вычислим значение невязки для полинома 3 степени
j=1;
I3=0;
while  j<11
    I3=I3+(F(j)-polyval(p3,Potok_norm(j)))^2;
    j=j+1;
end;
% Строим графики
x1=0:0.001*Potok_nom:1.02*Potok_nom;
subplot(2,2,1)
plot(x1,F_raschet5,Potok,F, '*' );
xlabel('Ф ');
ylabel('F , p5 ');
legend ('p5(Ф)','F(Ф)')
grid on;
subplot(2,2,3);
plot(x1,F_raschet3,Potok,F, '*'  );
grid on;
xlabel('Ф ');
ylabel('F, p3 ')
legend ('p3(Ф)','F (Ф )')
F_raschet5=polyval(p_norm5,x);
subplot(2,2,2);
plot(x,F_raschet5,Potok_norm,F_norm, '*' );
xlabel('Ф норм');
ylabel('F норм, p5 норм')
legend ('p5(Ф норм)','F норм(Ф норм)')
grid on;
F_raschet3=polyval(p_norm3,x);
subplot(2,2,4);
plot(x,F_raschet3,Potok_norm,F_norm , '*' );
xlabel('Ф норм');
ylabel('F норм, p3 норм')
legend ('p3(Ф норм)','F норм(Ф норм)')
grid on;