How to plot values from a for loop?

Question

Ljix 2016-2-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/266423-how-to-plot-values-from-a-for-loop

编辑： Star Strider 2016-2-3

在 MATLAB Online 中打开

Maybe my question isn't so specific, but I'll try to explain what's giving me a hard time.

I have to calculate

(C/(s*I-A))*B

where A,B and C are matrices and s is a eigenvalue of some given matrix. I have to calculate it for every eigenvalue and plot it. This is what I did

svA=eig(A);
for P=1:50
fr=(CN/(i*svA(P)*eye(n)-A))*B;
fr_org{P}=fr;
end
fr_original=fr_org{1:50};
plot(fr_original)

but I didn't get anything. Any advice?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Walter Roberson 2016-2-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/266423-how-to-plot-values-from-a-for-loop#answer_208362

在 MATLAB Online 中打开

Try

fr_original = cell2mat(fr_org);

3 个评论
显示 1更早的评论隐藏 1更早的评论

Walter Roberson 2016-2-3

plot(fr_original, '*')

Ljix 2016-2-3

Now it's working. Thank you.

请先登录，再进行评论。

Answer 2

Suresh Garimella 2016-2-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/266423-how-to-plot-values-from-a-for-loop#answer_208367

编辑：Walter Roberson 2016-2-3

在 MATLAB Online 中打开

Try this, Hope it works

svA=eig(A); 
for P=1:50 fr=(CN/(i*svA(P)*eye(n)-A))*B; 
fr_org{P}=fr; 
end 
%fr_original=fr_org{1:50}; comment this line 
fr_original=fr_org; % directly assign 
%plot(fr_original) comment this line 
plot(cell2mat(fr_original)) % convert class using cell2mat to plot

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Ljix 2016-2-3

在 MATLAB Online 中打开

Almost what I want. Using

plot(cell2mat(fr_original),'+') gives what I want.

请先登录，再进行评论。

Answer 3

Star Strider 2016-2-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/266423-how-to-plot-values-from-a-for-loop#answer_208378

编辑：Star Strider 2016-2-3

在 MATLAB Online 中打开

Actually, this is incorrect:

Y = (C/(s*I-A))*B

The correct expression is:

Y = C*((s*I-A)\B)

that you then invert to create:

y = C*expm(A*t)*B

You cannot do what you want to do in the Laplace domain, especially since s = σ + j*ω. You have to invert it to the time domain to do anything with it.

If you want to do a pole-zero plot, there are functions for that such as pzplot in the Control Systems Toolbox.

2 个评论
显示无隐藏无

Ljix 2016-2-3

编辑：Ljix 2016-2-3

Thank you. Isn't frequency response given by Y = C*((s*I-A)\B), where s=i*ω and ω is eigenvalue of A? How to plot y = C*expm(A*t)*B? I tried defining t=linspace(0,100,25) and then calculate y as function (as defined) and I got message about dimension error.

Star Strider 2016-2-3

编辑：Star Strider 2016-2-3

My pleasure.

In a pole-zero plot, σ is the real part, and jω is the imaginary part. This becomes equivalent to the Fourier transform if you set σ=0.

If you want the Bode plot (frequency response only, along the σ=0 line), use the transfer function representation of your system created from the core MATLAB ss2tf function, then plot it using the Signal Processing Toolbox freqs function.

EDIT — I Answered your question about y = C*expm(A*t)*B*u in How to plot y = C*expm(A*t)*B?

请先登录，再进行评论。