% Исходные данные
clear;
clc;
Ra=0.3;
R_v=145;
U_c=220 ;
w=4000;
Ce=205;
Cm=200;
%Номинальные данные
F_nom=U_c*w/R_v;
Potok_nom=0.007;
w_nom=100;
i_nom=50;
M_nom=70;
x0=[1; 1;1];
u1=1.2:-0.2:0;
xx=[];
global a p dp;
a=[U_c*w/(F_nom*R_v) Ce*w_nom*Potok_nom i_nom*Ra U_c M_nom Cm*Potok_nom*i_nom];
dp=[3.7395 0 -0.5688 0 0.5022];
u=[u1; ones(size(u1)); ones(size(u1))];
%u=[ones(size(u1)); u1; ones(size(u1))];
%u=[ones(size(u1)); ones(size(u1)); 1.5    1.2    0.8000    0.6000    0.4000    0.2000         0.1];
% Найдем статические характеристики для различных входов
for i=1:length(u1),
    x=newton(x0,u(:,i));
    xx=[xx x];
    x0=x;
end
% Изобразим их на графике
hPlot=plot(u1,xx(1,:),'>:',u1,xx(2,:),'*-.',u1,xx(3,:), 'o-');
set( hPlot, 'LineWidth', 2 );
xlabel('Uв');
ylabel('Ф, Iг, w')
legend ('Ф','Iг','w' )
grid on;