求控制工程中的传递函数的动态误差系数

wd840611

<>各位仁兄，大家好。最近和本组同学为准备国际大赛作了几个程序。希望多多指教。</P>

<>求控制工程中的传递函数的动态误差系数：
function [G,d]=as(G)
    [num,den]=numden(G);</P>
<>disp('请确定是否为开环传递系统，是请按1，否请按0')
c=input('');
disp ('请确定所要的展开级数：')
N=input('');
disp('请输入所要选取的级数系数：')
M=input('')</P>
<P>if M&gt;N
     disp('所选取的级数超出最高的级数。')
     disp('级数的系数为：0');</P>
<P>     else
         if c==1
           G=1/(1+G);
        else
           G=G;
         end
         syms s
         sum=1;
         [num,den]=numden(G);
            if den==0
                disp('输入的方程式的分母为零！')
                 disp ('请重新输入：')
                 else
                       for i=1:1:N
                       G0=diff(G);
      
                       sum=sum*i;
                       [num,den]=numden(G0);
                         if M&lt;N
                            for j=M:1:N
                             disp('方程的输出的状态的系数为：')
                                    numde=limit(G0,s,0)/sum;
                                      
                                  end
                                   G=G0;
                                  if M==N
                                      for i=1:1:N
                                 
                                   disp('方程的输出的状态的系数为：')
                                    numde=limit(G0,s,0)/sum
                                    
                                   end  
                                     G=G0;
                                   end
end
end
end
end</P>