longer_2_change_kp.m

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

kp=[0.05,0.06,0.07,0.08,0.09,0.1];
length=199;
for j=1:6
ki=0.02;
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(j)*e(i)+ki*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