Matlab - Error using trapz

%---------------------------------------------------------------------------------------------------

tengo este codigo de un convertidor, me presenta un error usando la funcion trapz y no he podido corregir el

error, cualquier ayuda les agradeceria.

%---------------------------------------------------------------------------------------------------

V= 220;

R= 10;

L= 0.02;

Ls= 0.002;

Vc= 50;

f= 60;

%------------------------------------------------------------------------

% Constantes

L=L+Ls;

m=Vc/( sqrt (2)*V);

Z= sqrt (R ^2+(2* pi*f*L) ^2);

fi= atan ((2* pi*f*L)/R);

alfa1 = asin (m); % Alfa mínimo

alfa2 =pi - alfa1 ; % Alfa máximo

disp ('Limite de Controlabilidad en grados ');

[ alfa1 *180/ pi , alfa2 *180/ pi]

alfa = 30; % ángulo de disparo en grados

alfa = alfa *pi /180;

% Limite de controlabilidad

if (alfa1 > alfa )|| (alfa > alfa2 )

   disp ('El ángulo de disparo está fuera del límite de contrabilidad se fijo en alfa minimo ');

   alfa = alfa1 ;

end

%------------------------------------------------------------------------

 t0 = pi;

 funx = @(t) sin(t-fi)-m/cos(fi)+(m/ cos (fi)-sin(alfa -fi))* exp (-(t-alfa )/tan(fi));

 b = fzero(funx,t0) ;

 g = b- alfa; % Angulo de Conducción

 a= alfa ;

 %------------------------------------------------------------------------

t= linspace (a,b ,1000) ;

i= sqrt (2)*V/Z*( sin(t-fi)-m/ cos (fi)+(m/cos(fi)-sin(a-fi)*exp (-(t-a)/tan(fi))));

v= sqrt (2)*V*sin(t);

tdes = linspace (b,pi+a ,300) ;

t=[t, tdes ];

v=[v,Vc* ones ( size ( tdes ))];

i=[i, zeros ( size ( tdes ))];

%------------------------------------------------------------------------

I_ave =1/pi* trapz (t,i);

Vo_ave =1/pi* (trapz (t,v));

%------------------------------------------------------------------------

%---------------------------------------------------------------------------------------------------

% tengo este codigo de un convertidor, me presenta un error usando la funcion trapz y no he podido corregir el

% error, cualquier ayuda les agradeceria.

%---------------------------------------------------------------------------------------------------

clear all

V= 220;

R= 10;

L= 0.02;

Ls= 0.002;

Vc= 50;

f= 60;

%------------------------------------------------------------------------

% Constantes

L=L+Ls;

m=Vc/( sqrt (2)*V);

Z= sqrt (R ^2+(2* pi*f*L) ^2);

fi= atan ((2* pi*f*L)/R);

alfa1 = asin (m); % Alfa mínimo

alfa2 =pi - alfa1 ; % Alfa máximo

disp ('Limite de Controlabilidad en grados ');

[ alfa1 *180/ pi , alfa2 *180/ pi]

alfa = 30; % ángulo de disparo en grados

alfa = alfa *pi /180;

% Limite de controlabilidad

if (alfa1 > alfa )|| (alfa > alfa2 )

   disp ('El ángulo de disparo está fuera del límite de contrabilidad se fijo en alfa minimo ');

   alfa = alfa1 ;

end

%------------------------------------------------------------------------

 t0 = pi;

 funx = @(t) sin(t-fi)-m/cos(fi)+(m/ cos (fi)-sin(alfa -fi))* exp (-(t-alfa )/tan(fi));

 b = fzero(funx,t0) ;

 g = b- alfa; % Angulo de Conducción

 a= alfa ;

 %------------------------------------------------------------------------

t= linspace (a,b ,1000) ;

i= sqrt (2)*V/Z*( sin(t-fi)-m/ cos (fi)+(m/cos(fi)-sin(a-fi)*exp (-(t-a)/tan(fi))));

v= sqrt (2)*V*sin(t);

tdes = linspace (b,pi+a ,300) ;

t=[t, tdes ];

v=[v,Vc* ones(1,length( tdes))];

i=[i, zeros(1,length(tdes))];

%------------------------------------------------------------------------

I_ave =1/pi* trapz (t,i)

Vo_ave =1/pi* (trapz (t,v))

%------------------------------------------------------------------------

>> convertidor

Limite de Controlabilidad en grados

ans =

    9.2479  170.7521

I_ave =

   16.0206

Vo_ave =

173.6121