syms t
[a,t]=gensig('sin',2*pi);aa=tf([1 1],[1 1 1]);
x1=a; x2=diff(a);
x1=x1(1:length(x1)-1);
x=[x1,x2];
t=t(1:length(t)-1);

sigma11=5;sigma21=5;sigma31=5;sigma12=5;sigma22=5;sigma32=5;
m11=5;m21=5;m31=5;m12=5;m22=5;m32=5;
theta1=3;theta2=3;theta3=3;
omega1=1;omega2=1;omega3=1;
sigma=[sigma11,sigma21,sigma31,sigma12,sigma22,sigma32];
m=[m11,m21,m31,m12,m22,m32];
theta=[theta1,theta2,theta3];
omega=[omega1,omega2,omega3];



ua11=exp(-(x1-m(1)).^2/(sigma(1)).^2);
ua21=exp(-(x1-m(2)).^2/(sigma(2)).^2);
ua31=exp(-(x1-m(3)).^2/(sigma(3)).^2);
ua12=exp(-(x1-m(4)).^2/(sigma(4)).^2);
ua22=exp(-(x1-m(5)).^2/(sigma(5)).^2);
ua32=exp(-(x1-m(6)).^2/(sigma(6)).^2);
ua=[ua11,ua21,ua31,ua12,ua22,ua32];
  
u1=ua(1).*ua(1+j); u2=ua(2).*ua(2+j);u3=ua(3).*ua(3+j);
u=[u1,u2,u3];
h1=u(1).*theta(1);h2=u(2).*theta(2);h3=u(3).*theta(3);
h=[h1,h2,h3];
f1=1/(1+exp(-h(1)));f2=1/(1+exp(-h(2)));f3=1/(1+exp(-h(3)));
f=[f1,f2,f3];
u(1)=f(1).*u(1);u(2)=f(2).*u(2);u(3)=f(3).*u(3);

  

y=u(1).*omega(1)+u(2).*omega(2)+u(3).*omega(3)
plot(t,y)
y=double(y);
lsim(aa,y,t,0)

运行的结果出现ERROR:
Error using ==> rfinputs
Input data U and time vector T must have the same number of samples.

错误的语句在于:lsim(aa,y,t,0) plot(t,y)命令可以运行

请高手指教！
[em17]