initial_cheenu.sci

//this is the actual file initial4
//work on this only

function[varargout]=initial(varargin)
    n=length(varargin);
    m=length(varargout);
    [lhs,rhs]=argn(0)
    if(typeof(varargin(1))<>'state-space') then
     error('wrong input argument .Enter state space model only' )
    end
 k=0;
 //////////determining no. of constants///////////////
 for j=1:rhs
     if(typeof(varargin(j))=='constant') then
         k=k+1;
     end
 end
 
 /////////////getting subplot credentiALA////////////

 /////////////////////////////////////////////////////////
 ///////////////////////////////////////////////////
 //disp(k);
 ///////if two constants///////////////////////////////////////time & x0 specification part//////////
 if (k==2) then
     x0=varargin(n-1);///////second last as initial condition///////////
      //disp(x0);
      if typeof(varargin($)) == 'constant' then /////time//////////
            if size(varargin($)) == [1 1] then
                if varargin($) <= 0 then
                    error(msprintf(gettext("%s: The final time value must be a positive real number.\n"),"stepplot"))
                end
                tFinal = varargin($)
                totalTime =0:0.01:tFinal; 
            elseif isequal(size(varargin($)),[1 1]) == %f then
                // finding that the time vector has positive time value
                tempTimeIndex = find(varargin($) >= 0)
                if isequal(size(varargin($)),size(tempTimeIndex)) == %t then
                    totalTime = varargin($)
                    for l=1:(n-2)
                    if(typeof(varargin(l))==['state-space'] & varargin(l).dt<>'c') then
                        if(varargin(l).dt=='d') then
                              t=1; /////////if 'd' then the sampling time is 1
                         else
                            t=varargin(l).dt
                        end
                        if(t<>(totalTime(2)-totalTime(1))) then
                            error('time interval must be equal to sampling time')
                        else
                            totalTime=0:t:totalTime($)
                        end
                    end
                 end           
                            
                else
                    ////////////here also we should check dt but lets do it later---""CHEENU""/////////////////
                    
                    tempTime = varargin($)
                    tempTime = tempTime(tempTimeIndex(1):tempTimeIndex($))
                    totalTime = tempTime
                end
            end
        end
    
elseif(k==1) then        //////////////////if time is not specified/////////////////////////
        x0=varargin($);
end
///////////////////////////////////////////////////////////////////////////////
ll=0;
 for j=1:rhs-k
     xx1(j)=size(varargin(j),'r');
   
////    if(typeof(varargin(j))=='state-space') then
////     bb=varargin(j).a
////    if(size(x0,'r')==size(bb,'r')) then
////        ll=ll+1;
////    else
////        error('wrong size of initial condition');
//// end
////end
////end
////tl=find(typeof(varargin(j))=='state-space');
////if (size(tl,'c')==ll) then
////    disp('paro')
////else
////    disp('cheenu')
end

index=1;
 for i=1:(n-k)
     index=index+1;
     if(typeof(varargin(i))=='string') then
         continue;
     end
     
   [row col]=size(x0);
     if(typeof(varargin(i))==['state-space'] & col==1) then
         [a b c d] =  abcd(varargin(i));
         ///////////checking matrix sizes////////////////
         if(size(b)<>[mtlb_size(a,1),mtlb_size(d,2)]| size(c)<>[mtlb_size(d,1),mtlb_size(a,2)]| size(d)<>[mtlb_size(c,1),mtlb_size(b,2)]|size(varargin(i).a,'r')<>row) then
               error('wrong size of matrices')
         else
            //b1 = zeros(mtlb_size(a,1),mtlb_size(d,2));
             d1 = zeros(mtlb_size(c,1),mtlb_size(x0,2));
             c1 = zeros(mtlb_size(d,1),mtlb_size(a,2));
         end
         
         
             if(((varargin(i).dt)=='c')) then
                sys=syslin('c',a,x0,c1,d1);
        
                 if (k==1) then
                     pq=pole(varargin(i))
                    pq=cell2mat(pq);
                           mm=real(pq);
                         //qq= find(pq>0)
                            t1=0:0.1:100;
                         if or(mm>0) then
                             totalTime=t1;
//                             [aa,x]=csim('impuls',t1,sys)
//                             aa=c*x;
                         else
                            [aa,x]=csim('impuls',t1,sys)
                             aa=c*x;
                             tt=find(aa>0.002);
                            totalTime=0:0.1:t1(tt($))
                         end
                         [aa,x]=csim('impuls',totalTime,sys)
                         aa=c*x;
                    

                  else
        //sys=syslin('c',a,x0,c1,d1);
                      [aa,x]=csim('impuls',totalTime,sys);
                      aa=c*x;
       
               
       
       
                end

         
         else  /////////////////this else is of string and continuous therefore for discrete//////////////////////
  
                if (k==1) then
       
                               pq=pole(varargin(i));
                               pq=cell2mat(pq);
                                t1=0:0.1:100;
                               mm=real(pq);
                               if(varargin(i).dt=='d') then
                                   dtime=1;
                               else
                                   dtime=varargin(i).dt;
                               end
                        if or(mm>0) then
                           // disp('paro')
                           totalTime=0:dtime:100;
                           //disp('paro')
//                           [aa,x]=ltitr(a,b,zeros(1,length(totalTime)),x0);
//                           aa=c*x;
                        elseif and(mm<0) then 
                           [aa,x]=ltitr(a,b,zeros((size(b,'c')),length(t1)),x0);
                           aa=c*x;
                           tt=find(or(aa>=0.002))       
                          totalTime=0:dtime:t1(tt($));
       
                       end
                      [aa,x]=ltitr(a,b,zeros(size(b,'c'),length(totalTime)),x0);
                      aa=c*x;
 
                 
             else
                [aa,x]=ltitr(a,b,zeros(size(b,'c'),length(totalTime)),x0);
                aa=c*x;
            end
        
      end

   elseif(typeof(varargin(i))<>['string']) then
       error('wrong type of initial condition entered');
   end
   ////////////////////////////plotting//////////////
   if (lhs==1) then
       if (size(varargin(i))==[1 1]) then
  if(typeof(varargin(i+1))=='string') then
      plot(totalTime,aa,varargin(i+1));
      hh=gce();
    if(varargin(i).dt<>'c') then
       hh.children.polyline_style=2;
    end
  else
    plot(totalTime,aa)
    hh=gce();
    hh.children.foreground=index;
  if(varargin(i).dt<>'c') then
   hh.children.polyline_style=2;
  end
 end
else
    //xx=size(aa,'r');
    xx=size(varargin(i),'r');
    Q=1;
    if(typeof(varargin(i+1))=='string') then
        
            for oo=1:xx
      //          for pp=1:yy
        subplot(max(xx1),1,Q)
          Q=Q+1;
          plot(totalTime,aa(oo,:),varargin(i+1))
           if(varargin(i).dt<>'c') then
               hh=gce();
               hh.children.polyline_style=2;
            end
        end
        
  // end

      else
         for oo=1:xx
      //for pp=1:yy
       subplot(max(xx1),1,Q)
      Q=Q+1;
     plot(totalTime,aa(oo,:))
     hh=gce();
     hh.children.foreground=index;
if(varargin(i).dt<>'c') then
               
               hh.children.polyline_style=2;
            end 
//end
 
end
end
end
//end    
    
h=gcf();
h.figure_name="INITIAL_RESPONSE";
varargout(1)="INITIAL_RESPONSE"
//end
/////////fot ends here////////////////


/////////////more than one output argument///////////////////////
elseif(lhs>1) then
     if((lhs==2)|(lhs==3)) then
         varargout(1)=aa;
         varargout(2)=totalTime;
         if(lhs==3) then
 
            varargout(3)=x;
          end
    end
end
end

endfunction
//if (size(varargin(i))==[1 1]) then
//  if(typeof(varargin(i+1))=='string') then
//      plot(totalTime,aa,varargin(i+1));
//      hh=gce();
//    if(varargin(i).dt<>'c') then
//       hh.children.polyline_style=2;
//    end
//  else
//    plot(totalTime,aa)
//    hh=gce();
//  if(varargin(i).dt<>'c') then
//   hh.children.polyline_style=2;
//  end
// end
//else
//    xx=size(varargin(i),'r');
//    yy=size(varargin(i),'c');
//    Q=1;
//    if(typeof(varargin(i+1))=='string') then
//        if(varargin(i).dt=='c') then
//            for oo=1:xx
//                for pp=1:yy
//        subplot(xx,yy,Q)
//          Q=Q+1;
//     plot(totalTime,aa,varargin(i+1))
//            end
//            end
//   else
//       for oo=1:xx
//                for pp=1:yy
//                    subplot(xx,yy,Q)
//                Q=Q+1;
//        plot(totalTime,aa,varargin(i+1));
//        hh=gce();
//        hh.children.polyline_style=2;
//end
//end
//end
//else
//    if(varargin(i).dt=='c')  then
//    
//    for oo=1:xx
//      //for pp=1:yy
//    
//    subplot(xx,yy,Q)
//    Q=Q+1;
//    plot(totalTime,aa) 
//end
// else
//    for oo=1:xx
//       //for pp=1:yy
//         subplot(max(oo),max(pp),Q)
//        Q=Q+1;
//        plot(totalTime,aa);
//        hh=gce();
//        hh.children.polyline_style=2;
//end
//end
//end
//end    
//    
//h=gcf();
//h.figure_name="INITIAL_RESPONSE";
//varargout(1)="INITIAL_RESPONSE"
////end
///////////fot ends here////////////////
//
///////////////more than one output argument///////////////////////
//else
//     if((lhs==2)|(lhs==3)) then
//         varargout(1)=aa;
//         varargout(2)=totalTime;
//         if(lhs==3) then
// 
//            varargout(3)=x;
//          end
//    end
//end
//end
//


// 
/////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
//  if(typeof(varargin(1))<>'state-space') then
//     error('wrong input argument .Enter state space model only' )
// end
// k=0;
// for j=1:rhs
//     if(typeof(varargin(j))=='constant')
//         k=k+1;
//     end
// end
// if (k==2) then
//     x0=varargin($-1);
//     if size(varargin($)) == [1 1] then
//                if varargin($) <= 0 then
//                    error("the final time value must be positive no.")
//                end
//                tFinal = varargin($)
//                totalTime =0:0.01:tFinal; 
//            elseif isequal(size(varargin($)),[1 1]) == %f then
//                // finding that the time vector has positive time value
//                tempTimeIndex = find(varargin($) >= 0)
//                if isequal(size(varargin($)),size(tempTimeIndex)) == %t then
//                    t= varargin($);
//                    for u=1:n
//                    if(((varargin(u)).dt<>'c')& typeof(varargin(u))<>'constant') then
//                        dt=varargin(u).dt;
//                        if(t(tempTimeIndex(2))-t(tempTimeIndex(1))<>dt) then                          
//                              error('time interval must be equal to sampling time')
//                          else
//                              t=0:dt:t(tempTimeIndex($));
//                          end
//                          end
//                          
//                      end
//                       
//                else
//                    tempTime = varargin($)
//                    tempTime = tempTime(tempTimeIndex(1):tempTimeIndex($))
//                    totalTime = tempTime
//                end
//            end
//        end
//     
// end
//else if(k==1) then
//        xo=varargin($) ;
//        totalTime=0:0.1:3;   
// end
// for i=1:n-k
//      [a b c d] =  abcd(varargin(i));
//         if(size(b)<>[mtlb_size(a,1),mtlb_size(d,2)]| size(c)<>[mtlb_size(d,1),mtlb_size(a,2)]| size(d)<>[mtlb_size(c,1),mtlb_size(b,2)]) then
//    error('wrong size of matrices')
//else
//            //b1 = zeros(mtlb_size(a,1),mtlb_size(d,2));
//    d1 = zeros(mtlb_size(c,1),mtlb_size(x0,2));
//    c1 = zeros(mtlb_size(d,1),mtlb_size(a,2));
//   
//   end
//    if(varargin(i).dt=='c') then
//       varargin(i)=syslin('c',a,x0,c1,d1);
//       [aa,x]=csim('impuls',totalTime,varargin(i));
//       aa=c*x;
//         
//   else
//       [aa,x]=ltitr(a,b,zeros(1,length(totalTime)),x0);
//       aa=c*x;
////       for i=1:length(totalTime);
////
////end
//end
//   else
//       error('wrong type of initial condition entered');
//end
//plot(totalTime,aa)
//end
//if(lhs>1) then
//     if((lhs==2)|(lhs==3)) then
//         varargout(1)=aa;
//         varargout(2)=totalTime;
//     if(lhs==3) then
//         varargout(3)=x;
//     end
// end
//end
// 
//endfunction
/////////////////////////
//for op=1:rhs-1
// sys1=ss2tf(varargin(i))
// pq=pole(sys)
//if (pq>0) then
//t=100;
//totalTime=0:1:100;
//else
//aa=0.02*aa;
//0.2*aa=csim('impuls',)
//
/// for i=1:2
//[a b c d]=abcd(op(i));
// t=0:op(i).dt:5
//[aa,x]=ltitr(a,b,zeros(1,length(t)),[4;5])
//end