Matlab里 y=[y,C*x];这句代码是什么意思?r=20.4;numo=[1.875e6 1.562e6];deno=[1 54 204 213 64];numh=0.002;denh=1;[num,den]=feedback(numo,deno,numh,denh);[A,B,C,D]=tf2ss(num,den);Tf=5; h=0.02;x=[zeros(length(A),1)];y=0;t=0;for i=1:Tf/hK1=A*x+B*r
来源：学生作业帮助网 编辑：六六作业网 时间：2024/05/01 03:05:48
Matlab里 y=[y,C*x];这句代码是什么意思?r=20.4;numo=[1.875e6 1.562e6];deno=[1 54 204 213 64];numh=0.002;denh=1;[num,den]=feedback(numo,deno,numh,denh);[A,B,C,D]=tf2ss(num,den);Tf=5; h=0.02;x=[zeros(length(A),1)];y=0;t=0;for i=1:Tf/hK1=A*x+B*r
Matlab里 y=[y,C*x];这句代码是什么意思?
r=20.4;
numo=[1.875e6 1.562e6];deno=[1 54 204 213 64];
numh=0.002;denh=1;
[num,den]=feedback(numo,deno,numh,denh);
[A,B,C,D]=tf2ss(num,den);
Tf=5; h=0.02;
x=[zeros(length(A),1)];y=0;t=0;
for i=1:Tf/h
K1=A*x+B*r;
K2=A*(x+h*K1/2)+B*r;
K3=A*(x+h*K2/2)+B*r;
K4=A*(x+h*K3)+B*r;
x=x+h*(K1+2*K2+2*K3+K4)/6;
y=[y,C*x];t=[t;t(i)+h];
end
plot(t,y)
Matlab里 y=[y,C*x];这句代码是什么意思?r=20.4;numo=[1.875e6 1.562e6];deno=[1 54 204 213 64];numh=0.002;denh=1;[num,den]=feedback(numo,deno,numh,denh);[A,B,C,D]=tf2ss(num,den);Tf=5; h=0.02;x=[zeros(length(A),1)];y=0;t=0;for i=1:Tf/hK1=A*x+B*r
在原矩阵y右边加x（两者行数相同）,增广为一个矩阵,并记录在y中