matlab的问题，为什么画出iac与t的图形得不到正确的答案，啊啊啊啊崩溃了，速救

代码如下：
clc;
clear; 
w=2*pi*50;
a=0.1;
b= 0.46;
V(1)=sqrt(2)*220;
phi(1)=0;
V(3)=sqrt(2)*20;
phi(3)=1.047;
V(5)=sqrt(2)*7;
phi(5)=1.57;
V(7)=sqrt(2)*3;
phi(7)=-1.047;
c=0.001;
R=15;
H=7;
K=49;
Y=zeros(49,49); 
YY=zeros(49,49); 
U=zeros(49);
U1=zeros(49);
for k=1:2:K
    for h=1:2:H 
        Y(k,h)=sqrt(1+(h*w*R*c)^2)*exp(-1i*(h+k)*(a+b)/2-1i*atan(h*w*R*c))*sin((h+k)*(b-a))/(w*pi*R*(h+k));   
        if h~=k
           YY(k,h)=sqrt(1+(h*w*R*c)^2)*exp(1i*(h-k)*(a+b)/2+1i*atan(h*w*R*c))*sin((h-k)*(b-a))/(w*pi*R*(h-k));    
           else 
            YY(k,h)=sqrt(1+(h*w*R*c)^2)*exp(1i*atan(h*w*R*c))*(b-a)/(2*pi*R);    
        end
    U(h)=abs(V(h))*exp(1i*phi(h));
    U1(h)=abs(V(h))*exp(-1i*phi(h));
    end
end
 I=YY*U+Y*U1; 
 display(I);
 i=zeros(49);
 A=zeros(49);
 B=zeros(49);
 C=zeros(49,49);
 for t=0:0.01:0.4
    iac(t)=0;
   for k=1:2:K  
   A(k)=abs(I(k));    %I的模值
   B(k)=angle(I(k));   %I的相角
   C(k,t)=cos(k*w*t+B(k));
   i(k)=A(k)*C(k);
   iac(t)=iac(t)+i(k);
   end
   display(iac(t));
    plot(t,iac(t));
    hold ;
    grid 
end
大神们最好让代码在自己的电脑上跑跑哈！！！我明天就要交了，拜托了matlab图形plot