% 編寫直流發電機外特性分析腳本函數
% 將該函數定義為dc_ge_load(dc_generator_load)
% 下面輸入電機基本數據：
nN=1000;IfN=2.5;Rf=291/12;Ra=1.4;
% 下面輸入750轉／分鐘時的空載特性試驗數據(Ifdata是勵磁電流,Eadata是感應電勢)：
 Ifdata=[ 0, 0.4, 1.0, 1.6,	2.0,  3.0, 8.5  12]; 

 Eadata=[15, 43, 88, 130,  150,  193.5,   290, 291 ]; 
% 下面進行空載特性曲線擬合：
p=polyfit(Ifdata,Eadata,3);
If=0:.01:12;% 進行勵磁電流參數重新設置：
 Ea=polyval(p,If);% 計算對應于擬合曲線的感應電勢
%plot(If,Ea)
%hold on
% 繪制擬合空載特性曲線：
% plot(If,Ea)
%hold on
%hold on
%xlabel('If[A]')
%ylabel('Ea[V]')
% 計算端電壓：
U=If*Rf;
Ia=(Ea-U)/Ra;
I=Ia-If;

% 計算負載電流：

 plot(I,U)
 hold on
%for j=1:9
% 對K和KK進行比較：
% for i=1:900
 % mmm(i)=Ea(i)/If(i);  
   %m(i)=abs(mmm(i)-K(j));
    %if m(i)<1
        %j
        %RL(j)
        %If(i)
       % If(i)*Rf
        %If(i)*Rf/RL(j)
        %Ea(i)
        % aa=If(i)
        %plot(If(i)*Rf/RL(j),If(i)*Rf,'*')
       % axis([0,50,0,350])
%hold on
  %      break;
  % end;
 
  %end
%end

% k=1:90;


% 計算他勵發電機外特性