longer_3_change_ki.m

来自「通过图形实现OFMD的图像仿真」· M 代码 · 共 32 行

ki=[0.01,0.02,0.03,0.04,0.05,0.06];
length=199;
for j=1:6
kp=0.06;
kd=0.01;
r=150;
c0=20;
current=c0;
u=zeros(1,length);
c=zeros(1,length);
e=zeros(1,length+1);
sum=0;
d=0;
e(1)=r-c0;
for i=1:length
    u(i)=kp*e(i)+ki(j)*sum+kd*d;
    c(i)=current+u(i);
    current=c(i);
    if (u(i)<=0 && u(i-1)>0 )
        c(i)-150
    end

    e(i+1)=r-current;
    d=e(i+1)-e(i);
    sum=sum+e(i+1);
end
t=1:length;
subplot(2,3,j);plot(t,c);
end